精选优质文档-倾情为你奉上新建一个函数脚本文件t0=0; tN=15; y0=5;0; h=0.01;t = t0 : h : tN; N = length (t); j = 1;for i = 1 : N t1 = t0 + h; K1 = R_K(t0, y0); K2 = R_K(t0 + h/2, y0 + h*K1/2); K3 = R_K(t0 + h/2, y0 + h*K2/2); K4 = R_K(t0 + h, y0 + h*K3); y1 = y0 + (h/6)*(K1 + 2*K2 + 2*K3 + K4); Y1(j)=y1(1); Y2(j)=y1(2); t0=t1; y0=y1; j=j+1;end plot (t, Y1, t, Y2); grid; 新建一个函数M-K.m文件function ydot = f (t, y)yd
Copyright © 2018-2021 Wenke99.com All rights reserved
工信部备案号:浙ICP备20026746号-2
公安局备案号:浙公网安备33038302330469号
本站为C2C交文档易平台,即用户上传的文档直接卖给下载用户,本站只是网络服务中间平台,所有原创文档下载所得归上传人所有,若您发现上传作品侵犯了您的权利,请立刻联系网站客服并提供证据,平台将在3个工作日内予以改正。