| 查看: 401 | 回复: 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
» 猜你喜欢
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有3人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
有考过系统分析师的虫子吗？交流一下你们的话拿来干什么
已经有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