exm040210.m

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

S=solve('sin(t)^2*exp(-0.1*t)-0.5*abs(t)','t')  

y_C=inline('sin(t).^2.*exp(-0.1*t)-0.5*abs(t)','t');  

t=-10:0.01:10;	
Y=y_C(t);			
clf,
plot(t,Y,'r');
hold on
plot(t,zeros(size(t)),'k');	
xlabel('t');ylabel('y(t)')
hold off 

zoom on						
[tt,yy]=ginput(5);zoom off	 

tt

[t4,y4]=fzero(y_C,0.1)