应用Runge_kutta法的程序怎么改正?
请大家帮我看一下错误在哪里,我找不出来:
function dydt=vdp1(t,y)
dydt=[y(2);(1-y(1)^2)*y(2)-y(1)];
[t,y]=ode23(@vdp1,[0 20],[0 -0.5]);
% plot(t,y(:,1));
??? function dydt=vdp1(t,y)
|
Error: Function definitions are not permitted at the prompt or in scripts.
|