exm030101.m

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

syms t x				
ft=t^2*cos(t)			
sx=int(ft,t,0,x)	

dt=0.05;				
t=0:dt:5;				
Ft=t.^2.*cos(t);		
Sx=dt*cumtrapz(Ft);	
t(end-4:end)			
Sx(end-4:end)			
plot(t,Sx,'.k','MarkerSize',12)
xlabel('x'),ylabel('Sx'),grid on