digui.m

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

function yy=digui(n,u)
% see also  digui_r(revised)
if n==1
   %y(1)=0.5;
   yy=0.5;
else
   %y(n)=digui(n-1,u)-sign(u)*0.1
   yy=digui(n-1,u)-sign(u)*0.1;
end
%yy=y(n);

%*mat_vib亮*<[e-mail]> 2008-05-03 13:37:13
%if index=1
%m(1)=0.5;
%else
%m(n)=m(n-1)-sign(u)*0.1
%end