inline_fun.m

来自「基本的源程序」· M 代码 · 共 16 行

% 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