digui.m

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

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