Advertisement
Guest User

Untitled

a guest
Dec 18th, 2018
51
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.65 KB | None | 0 0
  1. y = @(x) (2./(sqrt(x-2)) -4).^2;
  2. yzab = @(x) y(x).*(1+2*(randn(1,1)-0.5)*1e-4);
  3. x0 = 4;
  4. kroki = logspace(-10, 1, 250);
  5.  
  6. y1 =@(x) 8./sqrt((x-2).^3) - 4./(x-2).^2;
  7. y2 = @(x) 8./(x-2).^3 - 12./sqrt((x-2).^5);
  8. D = @(x, h) (y(x + h) - y(x))./h;
  9. Dzab =@(x,h) D(x, h).*(1+2*(randn(1,1)-0.5)*1e-4);
  10. i = 0;
  11.  
  12. for krok = kroki
  13. i = i+1;
  14. blad(i) = abs(y1(x0) - D(x0, krok));
  15. bladzab(i) = abs(y1(x0) - Dzab(x0, krok));
  16.  
  17. blad2(i) = 2.*sqrt(abs(y(x0)./y2(x0))*eps);
  18. end
  19.  
  20. loglog(kroki, blad, 'r'); hold on;
  21. loglog(kroki, bladzab, 'b'); hold on;
  22. loglog(kroki, blad2, 'g');
  23.  
  24. [minimum, ind] = min(blad);
  25. loglog(kroki(ind), minimum, 'or');
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement