Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- x = [-2,-1,0,1,2,3,4];
- y = [-14,0.5,3.1,0,-3,0,16];
- n=length(x);
- #inicializando a tabela
- T = zeros(n,n);
- #primeira coluna
- T(:,1)=y;
- #segunda coluna
- for j = 2:n
- T(j,2)=(T(j,1)-T(j-1,1))/(x(j)-x(j-1));
- for k = 3:n
- #terceira coluna
- T(k,3)=(T(k,2)-T(k-1,2))/(x(k)-x(k-2));
- for m = 4:n
- #quarta coluna
- T(m,4)=(T(m,3)-T(m-1,3))/(x(m)-x(m-3));
- for q = 5:n
- #quinta coluna
- T(q,5)=(T(q,4)-T(q-1,4))/(x(q)-x(q-4));
- for r = 6:n
- #sexta coluna
- T(r,6)=(T(r,5)-T(r-1,5))/(x(r)-x(r-5));
- for s = n
- #setima coluna
- T(s,7)=(T(s,6)-T(s-1,6))/(x(s)-x(1));
- endfor
- endfor
- endfor
- endfor
- endfor
- endfor
- T
- p = T(1,1) + T(2,2)*(x-x(1)) + T(3,3)*(x-x(1))*(x-x(2)) + T(4,4)*(x-x(1))*(x-x(2))*(x-x(3)) + T(5,5)*(x-x(1))*(x-x(2))*(x-x(3))*(x-x(4)) + T(6,6)*(x-x(1))*(x-x(2))*(x-x(3))*(x-x(4))*(x-x(5)) + T(7,7)*(x-x(1))*(x-x(2))*(x-x(3))*(x-x(4))*(x-x(5))*(x-x(6))
- p
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement