在matlab中,怎么编写隐式RK方法求解线性微分方程的程序,急需解决啊,谢谢!
微分方程为:
dy=f(y,t)=y+2(t+1)-(t+1)^2; t属于[0,1.5]
y(0)=1;
隐式龙格库塔方法为:A=[0.222 -0.3003 -0.2778;-0.0360 -0.01389 0.0098;0.4804 -0.2680 0.1389];b=[0.4444 -0.2778 0.2778];c=[0.5000,0.1127,0.8873];
取h=1.5,算一步就可以,谢谢啊