Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- >> format long
- >> %zadanie 8.1
- >> function y=f(x);
- y = 1/(1+x);
- endfunction
- >> a = 0;
- >> b = 0.6;
- >> %Wynik Rzeczywisty obliczony recznie
- >> wynik_rz = 0.470004;
- >> %dk to dokladnosc jest takze epsilonem
- >> dk1 = (10^(-4));
- >> dk2 = (10^(-6));
- >> %Obliczanie kwadratury i bledu_teo za pomoca funkcji
- >> [wynik_rbg(1) , blad_teo(1),numer_przyb(1)] = oblicz_Rbg(@f,a,b,dk1);
- >> [wynik_rbg(2),blad_teo(2),numer_przyb(2)] = oblicz_Rbg(@f,a,b,dk2);
- >> %Wyznaczenie kwadratury za pomoca funkcji quad
- >> wynik_q = quad(@f,a,b)
- wynik_q = 0.470003629245736
- >> %Numer przyblizenia (10^(-4))
- >> numer_przyb(1)
- ans = 4
- >> %Wynik dla (10^(-4))
- >> wynik_rbg(1)
- ans = 0.470003638311235
- >> %Blad Kwadratury
- >> blad_teo(1)
- ans = 1.839718980656269e-06
- >> %Blad rzeczywisty
- >> blad_rz1 = abs(wynik_rbg(1) - wynik_q)
- blad_rz1 = 9.065499462224125e-09
- >> %Sprawzenie czy |Qr - I(f) <= Epsilon
- >> blad_rz1<= dk1
- ans = 1
- >> %Numer przyblizenia (10^(-6))
- >> numer_przyb(2)
- ans = 5
- >> %Wynik dla (10^(-6))
- >> wynik_rbg(2)
- ans = 0.470003629264101
- >> %Blad Kwadratury
- >> blad_teo(2)
- ans = 9.047133875395019e-09
- >> %Blad rzeczywisty
- >> blad_rz2 = abs(wynik_rbg(2) - wynik_q)
- blad_rz2 = 1.836558682910550e-11
- >> %Sprawzenie czy |Qr - I(f) <= Epsilon
- >> blad_rz2 <= dk2
- ans = 1
Add Comment
Please, Sign In to add comment
