第2个回答 2009-12-05
%令:
%x'=y1
%x=y2
function dx=myfunc(t,x)
dx(1)=cos(t)-x(1)+x(2)-x(2).^3;%x''
dx(2)=x(1);
end
%%%%%%%%%%保存为myfunc.m%%%%%%%%%%%%%%%%
%在命令窗口输入:
clear
t=0:0.1:10; %t范围
x0=[1,0];%零点:x'(0)=1,x(0)=0,自定义
[t,x]=ode45('myfunc',t,y0);
plot(t,x(:,2),'r',t,x(:,1),'b')%t,x; t,x'关系图