thrbsvx

xp=0.0:0.2:1.0;
yp=[0.00000, 0.20271, 0.422793, 0.684137, 1.029639, 1.557408];
N=5;

for j=0:N

a(j+1,1)=yp(j+1);

end

for j=1:N

for i=j:N

a(i+1,j+1)=(a(i+1,j)-a(i,j))/(xp(i+1)-xp(i-j+1));

end

end

for j=0:N

c(j+1)=a(j+1,j+1);

end

x=input('Podaj wartośc x w przedziale [0.0 ,1.0] =');

cn=0;

pom=1;

for j=0:N

cn=cn+c(j+1)*pom;

pom=pom*(x-xp(j+1));

end

disp('Wartość obliczona metodą Newtona:;')

disp(cn)

disp('Wartość wyznaczona z funkcji interp:')

d=splin(xp,yp);

disp(d)

y=interp(x, xp, yp, d);

disp(y)

p=tan(x);

disp('Wartość obliczona z rówania funkcji:')

disp(p)

disp('Różnica między interpolacją a funkcją interp:'+string(abs(y-p)));