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

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

disp('按任意鍵繼續！')
disp('下面將通過鼠標操作從圖上獲取 5 個零點的近似值。') 
pause

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

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

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