leehcaabral

T2 - CalcNum

Jun 23rd, 2019
559
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. x = [-2,-1,0,1,2,3,4];  
  2. y = [-14,0.5,3.1,0,-3,0,16];  
  3. n=length(x);
  4. #inicializando a tabela  
  5. T = zeros(n,n);  
  6. #primeira coluna  
  7. T(:,1)=y;  
  8. #segunda coluna  
  9. for j = 2:n
  10.   T(j,2)=(T(j,1)-T(j-1,1))/(x(j)-x(j-1));  
  11.   for k = 3:n
  12.   #terceira coluna  
  13.     T(k,3)=(T(k,2)-T(k-1,2))/(x(k)-x(k-2));  
  14.     for m = 4:n
  15.       #quarta coluna  
  16.       T(m,4)=(T(m,3)-T(m-1,3))/(x(m)-x(m-3));
  17.       for q = 5:n
  18.         #quinta coluna
  19.         T(q,5)=(T(q,4)-T(q-1,4))/(x(q)-x(q-4));
  20.         for r = 6:n
  21.           #sexta coluna
  22.           T(r,6)=(T(r,5)-T(r-1,5))/(x(r)-x(r-5));
  23.           for s = n
  24.             #setima coluna
  25.             T(s,7)=(T(s,6)-T(s-1,6))/(x(s)-x(1));
  26.           endfor
  27.         endfor
  28.       endfor
  29.     endfor
  30.   endfor
  31. endfor
  32. T
  33.  
  34. 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))
  35. p
RAW Paste Data

Adblocker detected! Please consider disabling it...

We've detected AdBlock Plus or some other adblocking software preventing Pastebin.com from fully loading.

We don't have any obnoxious sound, or popup ads, we actively block these annoying types of ads!

Please add Pastebin.com to your ad blocker whitelist or disable your adblocking software.

×