|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我想调用四阶龙格库塔算法:
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;
以上程序是书上的,我的调用程序是:
z=1+log(x+1);
yy=DELGKT4_lungkuta(z,0.1,0,1,1,[x y])
为什么老出错啊?应该怎样修改?我用的是matlabR2008a 我是新手。谢谢。
[ 本帖最后由 sigma665 于 2008-12-26 09:43 编辑 ] |
|