看了看,运行了一下确实出问题 原因是在用ode数值求解时,x并不是1:0.01:3均匀分散的 解决方法：1.在画解析解和欧拉解时横轴用x的转置；在画数值解时横轴用x,已运行成功2.在用ode数值求解时设定分散点的个数,与x相同即可