% inline_fun.m
% 求解超越函數的所有零點
% f(t)=sin(t)^2*exp(-a*t)-b*abs(t)
y=inline('sin(t)^2*exp(-a*t)-b*abs(t)','t','a','b');
a=0.2;b=0.6;t=-10:0.01:10;
% 向量化內置函數處理
y_vector=vectorize(y);
yy=feval(y_vector,t,a,b);
clf
plot(t,yy,'r','Linewidth',2);
hold on
plot(t,zeros(size(t)),'k');
xlabel('t');
ylabel('y(t)');
grid on
hold off