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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 信息科学 » 程序算法 » RK4（龙格库塔）程序，大家交流
21/1返回列表
查看: 414  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
    • 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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-07-30 11:00:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weivguo

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
Funval程序不足以实现功能
赞一下 回复此楼
2楼2015-07-26 11:08:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 curton 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭