%exm045_1.m
%（1）使用字符串表示被處理函數(shù)
P1=0.1;P2=0.5;		%按泛函指令要求，這里參數(shù)必須用P1,P2,...表示
y_C='sin(x).^2.*exp(-P1*x)-P2*abs(x)';	%這里自變量必須用x表示	<2>  

%（2）作圖法觀察函數(shù)零點(diǎn)分布
x=-10:0.01:10;	    %對自變量采樣，采樣步長不宜太大。
Y=eval(y_C);		%在采樣點(diǎn)上計(jì)算函數(shù)值。
clf,shg
plot(x,Y,'r');hold on,plot(x,zeros(size(x)),'k');	%畫坐標(biāo)橫軸
xlabel('t');ylabel('y(t)'),hold off  

disp('按任意鍵繼續(xù)！')
disp('下面將通過鼠標(biāo)操作從圖上獲取 5 個(gè)零點(diǎn)的近似值。') 
pause

zoom on						            %在MATLAB指令窗中運(yùn)行，獲局部放大圖
[tt,yy]=ginput(5);zoom off	     %在MATLAB指令窗中運(yùn)行，用鼠標(biāo)獲5個(gè)零點(diǎn)猜測值。
tt			                                 %顯示所得零點(diǎn)初始猜測值（該指令可在Notebook中運(yùn)行）。  

disp('按任意鍵繼續(xù)！')
disp('下面求靠近tt(4) 的精確零點(diǎn)。') 
pause

[t4,y4,exitflag]=fzero(y_C,tt(4),[],P1,P2)