exm05622_1.m

来自「这是本人收集的一些程序源代码」· M 代码 · 共 15 行

%exm05622_1.m
y=inline('sin(t)^2*exp(-a*t)-b*abs(t)','t','a','b');				%<1>  

a=0.1;b=0.5;
t=-10:0.01:10;
y_char=vectorize(y);														%<4>
Y=feval(y_char,t,a,b);
clf,plot(t,Y,'r');hold on,plot(t,zeros(size(t)),'k');
xlabel('t');ylabel('y(t)'),hold off  

zoom on,disp('先把零点区域局部放大；然后按任意键，继续。'),pause
[tt,yy]=ginput(5);tt
zoom off

[t4,y4,exitflag]=fzero(y,tt(4),[],a,b)								%<11>