function y = DEYCJZ_ml(f,h,a,b,y0,varvec,type,s)
format long;
N = (b-a)/h;
y = zeros(N+1,1);
x = a:h:b;
y(1) = y0;
y(2) = y0+h*Funval(f,varvec,[x(1) y(1)]);
y(3) = y(2)+h*Funval(f,varvec,[x(2) y(2)]);
y(4) = y(3)+h*Funval(f,varvec,[x(3) y(3)]);
if type == 1
for i=5:N+1
v1 = Funval(f,varvec,[x(i-3) y(i-3)]);
v2 = Funval(f,varvec,[x(i-2) y(i-2)]);
v3 = Funval(f,varvec,[x(i-1) y(i-1)]);
t = y(i-4) + 4*h*(2*v3 - v2 + 2*v1)/3;
for m=1:s
ft = Funval(f,varvec,[x(i) t]);
y(i) = y(i-2)+h*(4*v3 + v2 + ft)/3;
t = y(i);
end
end
else
p0 = 0;
c = 0;
for i = 5:N+1
v1 = Funval(f,varvec,[x(i-3) y(i-3)]);
v2 = Funval(f,varvec,[x(i-2) y(i-2)]);
v3 = Funval(f,varvec,[x(i-1) y(i-1)]);
t = y(i-4) + 4*h*(2*v3 - v2 + 2*v1)/3;
ft = Funval(f,varvec,[x(i) t]);
p = y(i-4)+4*h*(2*v3 - v2 + 2*v1)/3;
M = p - 28*(p0 - c)/29;
F = Funval(f , varvec, [x(i) ,M]);
c = y(i-2)+h*(4*v3 + v2 + F)/3;
y(i) = c + ( p - c)/29;
p0 = p;
end
end
format short;
没有合适的资源?快使用搜索试试~ 我知道了~
matlab 解方程算法程序(部分)
共26个文件
m:26个
需积分: 9 20 下载量 123 浏览量
2011-04-20
12:37:30
上传
评论
收藏 10KB ZIP 举报
温馨提示
常微分方程的初值问题 里面有很多的解方程的算法程序 可以自己学用
资源推荐
资源详情
资源评论
收起资源包目录
常微分方程的初值问题.zip (26个子文件)
DELGKT4_jer.m 479B
DEYCJZ_myds.m 702B
DEYDS.m 628B
DELGKT3_suen.m 360B
DELGKT4_lungkuta.m 419B
DEYCJZ_mid.m 1KB
DEMS.m 572B
DELSBRK.m 586B
DEYCJZ_adms2.m 607B
DELGKT3_kuta.m 358B
DEWT.m 612B
DEWT_glg.m 991B
DELGKT2_suen.m 313B
DELGKT2_mid.m 303B
DEimpEuler1.m 325B
DEYCJZ_adms.m 406B
DEYCJZ_hm.m 1KB
DELGKT4_qt.m 423B
DEYCJZ_ml.m 1KB
Funval.m 253B
NewtonRoot.m 810B
DEYCJZ_yds.m 604B
DEMiren.m 501B
DEimpEuler.m 272B
DEEuler.m 200B
DEModifEuler.m 321B
共 26 条
- 1
资源评论
haoshilin
- 粉丝: 2
- 资源: 7
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功