| 查看: 476 | 回复: 1 | |||||
curton主管区长 (文坛精英)
Hello Kitty
|
[交流]
RK4（龙格库塔）程序，大家交流 已有1人参与
|
|
function fv=Funval(f, varvec, varval) var= findsym(f); if length(var)<4 if var(1)==varvec(1) fv=subs(f,varvec(1),varval(1)); else fv=subs(f,varvec(2),varval(2)); end else fv=subs(f,varvec,varval); end function y = DELGKT4_lungkuta(f, h,a,b,y0,varvec) format long; N = (b-a)/h; y = zeros(N+1,1); y(1) = y0; x = a:h:b; var = findsym(f); for i=2:N+1 K1 = Funval(f,varvec,[x(i-1) y(i-1)]); K2 = Funval(f,varvec,[x(i-1)+h/2 y(i-1)+K1*h/2]); K3 = Funval(f,varvec,[x(i-1)+h/2 y(i-1)+K2*h/2]); K4 = Funval(f,varvec,[x(i-1)+h y(i-1)+h*K3]); y(i) = y(i-1)+h*(K1+2*K2+2*K3+K4)/6; end format short; syms x y; z=1+log(x.*x+1)+sin(x)+y yy=DELGKT4_lungkuta(z,0.1,0,1,1,[x,y]) |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : DELGKT4_lungkuta.m
- 附件 2 : Funval.m
- 附件 3 : main.m
2013-07-30 10:59:37, 420 bytes
2013-07-30 10:59:52, 246 bytes
2013-07-30 11:00:02, 79 bytes
» 猜你喜欢
时间戳又变了
已经有13人回复
-
科研孤儿太难了
已经有4人回复
-
今天系统多次维护，明天很可能放榜！
已经有10人回复
-
放榜前的不淡定
已经有12人回复
-
言之凿凿，8月21号(明天)放榜
已经有16人回复
-
时间戳变了，能看出什么问题？
已经有21人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
只有每年这种时候来逛逛小木虫
已经有21人回复
-
今天基金会出结果吗？20260819
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有考过系统分析师的虫子吗？交流一下你们的话拿来干什么
已经有12人回复
- 有关fortran的一次而问题，希望大家能帮帮忙，谢谢
已经有4人回复
- 【求助/交流】哪位好人告诉我小麦反转录的程序
已经有6人回复
- 【交流】怎样用mathematica 做积分嵌套积分的运算
已经有5人回复
- 【分享】C/C++常用算法设计方法,fortran其他语言编程也可以参考下
已经有8人回复
- 【讨论】子程序Umat讨论帖
已经有12人回复
- 【交流】VC++, C#, VB答疑专贴
已经有145人回复
- 【交流】Fortran语言答疑专帖
已经有157人回复
weivguo
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1821.8
- 帖子: 160
- 在线: 28.7小时
- 虫号: 1300089
- 注册: 2011-05-19
- 性别: GG
- 专业: 偏微分方程
2楼2015-07-26 11:08:57
40