Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- y = @(x) (2./(sqrt(x-2)) -4).^2;
- yzab = @(x) y(x).*(1+2*(randn(1,1)-0.5)*1e-4);
- x0 = 4;
- kroki = logspace(-10, 1, 250);
- y1 =@(x) 8./sqrt((x-2).^3) - 4./(x-2).^2;
- y2 = @(x) 8./(x-2).^3 - 12./sqrt((x-2).^5);
- D = @(x, h) (y(x + h) - y(x))./h;
- Dzab =@(x,h) D(x, h).*(1+2*(randn(1,1)-0.5)*1e-4);
- i = 0;
- for krok = kroki
- i = i+1;
- blad(i) = abs(y1(x0) - D(x0, krok));
- bladzab(i) = abs(y1(x0) - Dzab(x0, krok));
- blad2(i) = 2.*sqrt(abs(y(x0)./y2(x0))*eps);
- end
- loglog(kroki, blad, 'r'); hold on;
- loglog(kroki, bladzab, 'b'); hold on;
- loglog(kroki, blad2, 'g');
- [minimum, ind] = min(blad);
- loglog(kroki(ind), minimum, 'or');
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement