dampmass_test.m

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

% function dampmass_test
plot(x,sindy)
o = odeset('MaxStep',1e-4) ;
[x,y]=ode45('dampmass',[0 1],0.001,o);
dydt=feval('dampmass1',y);
sindy=dydt.*(sin(2*pi*75*x));