Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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)));
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement