ndd x = [0 0.25 0.5 0.75]; f=[1 1.64872 2.71828 4.4816]; p=0.43; n=4; for i=1:n F(i,1)=f(i) end for i=2:n for j=i:n F(j,i)=(F(j,i-1)-F(j-1,i-1)) / (x(j) - x(j-i+1)); end end for i=1:n product(i)=1; for j=1:i-1 product(i)=product(i)(p-x(j)); end end sum = 0; for i=1:n sum=sum+(F(i,i) product(i)); end sum

laint n=8;x=linspace(0,2pi,8);y= @(x) (sin(x))^2;p=input('Enter p of approx: '); for i=1:n l(i)=1; for j=1:n if(i~=j) l(i)=((p-x(j))/(x(i)-x(j)))l(i); end end end sum=0; for i=1:n sum=sum+l(i)*y(x(i)); end sum