dufenfangcheng.m

这是一个求杜芬方程解的小程序。开发环境为matlab。通过仿真曲线可以很容易的看到解的曲线。

TSPAN=[0,50];
x0=[0.3,0.5]';
[t,x]=ode23('xprime',TSPAN,x0);
subplot(2,2,1);
plot(t,x(:,1),':b',t,x(:,2),'-r')
subplot(2,2,2);
%figure(2)
plot(x(:,1),x(:,2));
axis([-10 10 -12 12]);

m1=x(:,1);
m2=x(:,2);




TSPAN=[0,50];
x0=[0.3,0.5]';
[t,x]=ode23('xprime2',TSPAN,x0);
subplot(2,2,3);
plot(t,x(:,1),':b',t,x(:,2),'-r')
subplot(2,2,4);
%figure(2)
plot(x(:,1),x(:,2));
axis([-10 10 -12 12]);

%axis([-10 10 -12 12]);