flagrange2l.m

来自「自己编写的matlab程序 数值分析里的程序」· M 代码 · 共 17 行

% flagrange2l([2 5])
function y=flagrange2l(x0);
syms x
if length(x0)==2
   y0(1)=exp(x0(1));
   y0(2)=exp(x0(2));
   l1=(x-x0(2))/(x0(1)-x0(2));
   l2=(x-x0(1))/(x0(2)-x0(1));
   y=y0(1)*l1+y0(2)*l2; 
   ezplot(y,[0 5]);
   hold on
   fplot('exp(k)',[0 5])  
   hold off
else
   disp('input error');
end