Constant a=1,b=0.5,c=2;
Variable x=[0,5],y1=25,y3=0,y2=[,800];
Plot x[x],y1,y2;
ODEFunction y1'=a*(y2-y1);
            y2'=y3;
            y3'=b*(y2-y1)+c*(y2-30);

2楼: Originally posted by shikang999 at 2014-12-01 16:50:43
如果是作业的话，建议将y3约去，得到一个二阶的y2与一阶的y1，然后把导数项展开，然后组装到一个方程组里面，根据边界条件设置相应的值，最后求解一个线性方程组即可。当然为了保证精度，可以采用不同步长求解，前后 ...