exm040109.m

来自「matlab电子书籍」· M 代码 · 共 9 行

tspan=[0,30];				
y0=[1;0];					
[tt,yy]=ode45(@DyDt,tspan,y0);									
plot(tt,yy(:,1))
xlabel('t'),title('x(t)')

plot(yy(:,1),yy(:,2))	
xlabel('位移'),ylabel('速度')