得到一个错误使用θ
5视图(30天)
显示旧的评论
运行这段代码:
函数[tvals, wvals] = RK4SYS (rhsf、init_t init_y, stop_t, Nsteps)
h = (stop_t-init_t) / Nsteps;
tvals (1) = init_t;
wvals (: 1) = init_y;
为指数= 1:Nsteps
tvals(指数+ 1)= tvals(指数)+ h;
k1 = h * rhsf (tvals(指数),wvals(:,指数));
k2 = h * rhsf (tvals(指数)+ h / 2, wvals(:,指数)+ k1/2);
k3 = h * rhsf (tvals(指数)+ h / 2, wvals(:,指数)+ k2/2);
k4 = h * rhsf (tvals(指数)+ h, wvals(:,指数)+ k3);
wvals(:,指数+ 1)= wvals(:,指数)+ (k1/6 + k2/3 + k3/3 + k4/6);
结束
结束
要做的:
Rhsprob = @ (t, U) [U(2), 32.17 / 2 *罪(θ)]
initU1 =[π/ 6;0]
电视,西弗吉尼亚州]= RK4SYS (Rhsprob 0 initU1 2 5);
但是我得到一个错误当θ说做这条线需要在15我的代码和错误的东西