vdp_solve.m

来自「解偏微分方程」· M 代码 · 共 12 行

function vdp_solve

mu  = 1;
[t,y] = ode45(@vdp,[0 20],[2; 0],[],mu);
plot(t,y(:,1),t,y(:,2))

% [t,y] = ode15s(@vdp,[0 3000],[2; 0],[],1000);
% plot(t,y(:,1),'-');

function y_dot = VDP(t,y,mu)
y_dot   = [y(2);
            mu*(1-y(1)^2)*y(2)-y(1)];