% Section 3.1, problem 6(a) format long s=[ 0.3 0.4 0.5 0.6]; x= exp(-2*s) - s; y= s; xstar = 0.0; n = length(x)-1; Q = zeros(n+1); Q(:,1) = y'; R = Q; for i = 1:n for j = 1:i Q(i+1,j+1) = ( (xstar-x(i-j+1))*Q(i+1,j) - (xstar-x(i+1))*Q(i,j) ) ... /( x(i+1)-x(i-j+1) ); end end Q(n+1,n+1) for j = 1:n for i = n:-1:j R(i+1,j+1) = ( (xstar-x(i-j+1))*R(i+1,j) - (xstar-x(i+1))*R(i,j) ) ... /( x(i+1)-x(i-j+1) ); end end R(n+1,n+1)