daily pastebin goal
55%
SHARE
TWEET

Untitled

a guest Mar 26th, 2019 74 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. clear all, clc,
  2. % punkt początkowy
  3. x=2;
  4. h0=0.2;
  5. % funkcja badana
  6. y0=@(x) exp(-x).*cos(x);
  7. % pierwsza pochodna
  8. y1=@(x) -exp(-x).*cos(x)-exp(-x).*sin(x);
  9. % druga pochodna
  10. y2=@(x) 2.*(exp(-x)).*sin(x);
  11. % trzecia pochodna
  12. y3=@(x) 2.*(-exp(-x)).*sin(x)+cos(x).*exp(-x);
  13. % krok
  14. i=1:40
  15. h =h0*2.^(1-i);
  16. % różnica progresywna
  17. rp=(y0(x+h)-y0(x))./h;
  18. % błąd względny // I BEZWZGLEDNY
  19. bw=abs((rp-y1(x))/y1(x));
  20. bgp=abs(h./2*y2(x)./y1(x));
  21. % różnica centralna
  22. rc=(y0(x+h)-y0(x-h))./(2*h);
  23. % błąd wzgldny // I BEZWZGLEDNY
  24. ec=abs((rc-y1(x))/y1(x));
  25. bgc=abs(h.^2/6*y3(x)/y1(x));
  26.  
  27. % wykresy
  28. loglog(h,bw,'r.');
  29. hold on
  30. grid on
  31. loglog(h,ec,'b.');
  32. loglog(h,bgp);
  33. loglog(h,bgc);
  34.  
  35. legend('r. progresywna', 'r. centralna', 'blad wzgledny odciecia r. progr.', 'blad wzgledny odciecia r.centr.')
  36.  xlabel('log(h)');  
  37.  ylabel('log(blad wzgledny)');
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top