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,2*pi,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