orthpolval.m

来自「用matlab写的一些数值算法」· M 代码 · 共 14 行

function  f = orthpolval(b,g,c,x)
% Value of orthogonal polynomial expansion

% Version 11.12.2003.  INCBOX

n = length(c)-1; 
m = length(x);            % number of simultaneous arguments
P = ones(m,2);   % for values of two consecutive polynomials
f = c(1) * P(:,2);                              % initialize
for  j = 1 : n
  P = [P(:,2)  (x(:) - b(j)).*P(:,2) - g(j)*P(:,1)];
  f = f + c(j+1)*P(:,2);
end