Facebook
From cubic Interpolation, 6 Months ago, written in Plain Text.
Embed
Download Paste or View Raw
Hits: 295
  1. A=[1,10,100,1000;1,15,225,3375;1,20,400,8000;1,22.5,506.25,11390.625];
  2. b=[227.04;362.78;517.35;602.97];
  3. n=length(b);
  4. for k=1:n-1
  5.   for i=k+1:n
  6.     factor=A(i,k)/A(k,k);
  7.     A(i,k:n)=A(i,k:n)-factor*A(k,k:n);
  8.     b(i)=b(i)-factor*b(k);
  9.   endfor
  10. endfor
  11. A
  12. x=zeros(n,1);
  13. for i=n:-1:1
  14.   x(i)=(b(i)-A(i,i+1:n)*x(i+1:n))/A(i,i);
  15. endfor
  16. y=x(1)+x(2)*16+x(3)*16*16+x(4)*16*16*16;
  17. disp(y)