版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

curton

主管区长 (文坛精英)

Hello Kitty

专家经验: +151
应助: 195 (高中生)
贵宾: 49.033
金币: 81392.2
散金: 25181
红花: 225
沙发: 109
帖子: 20300
在线: 15504.3小时
虫号: 950839
注册: 2010-01-30
性别: GG
专业: 电力电子学
管辖: 科研生活区

[交流] 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

2013-07-30 10:59:37, 420 bytes

附件 2 : Funval.m

2013-07-30 10:59:52, 246 bytes

附件 3 : main.m

2013-07-30 11:00:02, 79 bytes

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

有考过系统分析师的虫子吗？交流一下你们的话拿来干什么已经有12人回复
有关fortran的一次而问题，希望大家能帮帮忙，谢谢已经有4人回复
【求助/交流】哪位好人告诉我小麦反转录的程序已经有6人回复
【交流】怎样用mathematica 做积分嵌套积分的运算已经有5人回复
【分享】C/C++常用算法设计方法,fortran其他语言编程也可以参考下已经有8人回复
【讨论】子程序Umat讨论帖已经有12人回复
【交流】VC++, C#, VB答疑专贴已经有145人回复
【交流】Fortran语言答疑专帖已经有157人回复

1楼 2013-07-30 11:00:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

weivguo

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 1821.8
帖子: 160
在线: 28.7小时
虫号: 1300089
注册: 2011-05-19
性别: GG
专业: 偏微分方程

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

Funval程序不足以实现功能

赞一下