odefun5.m

来自「(有源代码)数值分析作业,本文主要包括两个部分,第一部分是常微分方程(ODE)的」· M 代码 · 共 14 行

function dy = ODEfun5(t,y)
% dy -> dy(1), and dy(2), is output
% y and t is input , 
% y -> y(1), y(2)
% is vector to  form the system of eqution

% y(1) is x, dy(1) is dx/dt, and y(2) , is y ; dy(2) is dy/dt
mu = 0.1;
dy =zeros(2,1);% initialize 
dy(1) = y(2) - mu * ( y(1) ^ 3 / 3 - y(1) );
dy(2) = - y(1) ;

return