flagrange2.m

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

% flagrange2([2 5],[3 4],[2.5 3 3.5 4 4.5])
function y=flagrange2(x0,y0,z);
syms x
if length(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);
   subs(y,x,z)
else
   disp('input error');
end