flagrangen.m

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

% flagrangen([2 3 4 5],[1 10 5 8],[2.5 3.5 4.5 5 5.5 10])
function flagrangen(x0,y0,z);
n=length(x0);
x=sym('x');
s=0.0;
for k=1:n
   p=1.0;
   for j=1:n
     if j~=k
         p=p*(x-x0(j))/(x0(k)-x0(j));
     end
   end
   s=p*y0(k)+s;
   s=simplify(s);
end
s
%画符号函数的图形
ezplot(s,[0 10]);
grid on
subs(s,x,z)