funtest3.m

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

function f=funtest3(x,y)
% f = -y + cos(2*x) - 2*sin(2*x) + 2*x*exp(-x);
% [x,y]=adams_wjl('funtest3',0,2,1,0.2);
% s= dsolve('Dy=-y +cos(2*x) - 2*sin(2*x) + 2*x*exp(-x)','y(0)=1','x')
% s =
% cos(2*x)+exp(-x)*x^2

%hold on
% plot(x,cos(2*x)+exp(-x).*x.^2,'*r')
f = -y + cos(2*x) - 2*sin(2*x) + 2*x*exp(-x);

%error is y-(cos(2*x)+exp(-x).*x.^2) ,not  y-cos(2*x)+exp(-x).*x.^2
% but y-(cos(2*x)-exp(-x).*x.^2)
%remember the  sign of the expression  , or the  ()