odefun5.m

(有源代码)数值分析作业,本文主要包括两个部分,第一部分是常微分方程(ODE)的三个实验题,第二部分是有关的拓展讨论,包括高阶常微分的求解和边值问题的求解(BVP).文中的算法和算例都是基于Matla

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