Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function wynik = rich()
- h = 0.1;
- x=1;
- A = zeros(3);
- for j=1:2;
- #metoda progresywna
- if j == 1
- for i=1:3;
- A(i,1)=(sin((x+h).^2)- sin(x.^2))./h;
- h=h/2;
- end;
- end;
- #metoda centralna
- if j == 2
- for i=1:3;
- A(i,1)=(sin((x+h).^2) - sin((x-h).^2))./2*h;
- h=h/2;
- end;
- end;
- A(2,2) = A(2,1) + (A(2,1) - A(1,1))/(2.^2 -1);
- A(3,2) = A(3,1) + (A(3,1) - A(2,1))/(2.^2 -1);
- A(3,3) = A(3,2) + (A(3,2) - A(2,2))/(2.^4 -1);
- W(j)=A(3,3);
- end;
- display(W)
- end;
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement