exm030102.m

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

syms t x
ft=exp(-sin(t))
sx=int(ft,t,0,2*pi)

dt=0.05;				
t=0:dt:4;			   
Ft=exp(-sin(t));		
Sx=dt*cumtrapz(Ft);	
Sx(end)				
plot(t,Ft,'*r','MarkerSize',4)
hold on
plot(t,Sx,'.k','MarkerSize',15)
hold off
xlabel('x')
legend('Ft','Sx')