《龙格库塔方法matlab实现》由会员分享,可在线阅读,更多相关《龙格库塔方法matlab实现(1页珍藏版)》请在金锄头文库上搜索。
龙格库塔方法matlab实现function ff=rk(yy,x0,y0,h,a,b)%yy为y的导函数,x0,y0,为初值,h为步长,a,b为区间c=(b-a)/h+1;i1=1; %c为迭代步数;i1为迭代步数累加值y=y0;z=zeros(c,6); %z生成c行,5列的零矩阵存放结果; %每行存放c次迭代结果,每列分别存放k1k4及y的结果for x=a:h:b if i1 rk(yy,0,1,0.2,0,1)%将得到结果%结果矩阵,第一列为x(n),第二列第五列为k1k4第六列为y(n+1)的结果%z = % 0 1.0000 1.2000 1.2200 1.4440 1.2428 % 0.2000 1.4428 1.6871 1.7115 1.9851 1.5836 % 0.4000 1.9836 2.2820 2.3118 2.6460 2.0442 % 0.6000 2.6442 3.0086 3.0451 3.4532 2.6510 % 0.8000 3.4510 3.8961 3.9407 4.4392 3.4365 % 1.0000 4.4365 4.9802 5.0345 5.6434 4.4401