求解二階微分方程，初始條件還需要給出y1"(0)和y2"(0)。這裡暫時按照0處理。

function zd530003514

a=0.1;

b=0.1;

Y0 = [b-1; 0; b; 0];

% 解方程

[t,Y]= ode45(@ode,[0 10],Y0);

y1=Y(:,1);

y2=Y(:,3);

% 繪圖

subplot 211

plot(t,y1);

subplot 212

plot(t,y2);

% 微分方程定義

function dY = ode(t, Y)

L1=5;

L2=0.01;

a0=2;

b0=2;

c0=2;

y1=Y(1);y2=Y(3);

dY = [

Y(2);

-(a0*y2+b0*y2^2+c0*y2^3) - L1^2*L2*y1 - L1^2*y1;

Y(4);

-(a0*y2+b0*y2^2+c0*y2^3) - L1^2*L2*y1;

];