lab7

clc
clear
% a=2;
% b=3;
m=12;
syms x
% % tutaj wyrazenie
% % wzor naszej funkcji
% fun1=3*x^2 -x*log(x);

a=2;
b=3;
%fun1 = cos(9*x)*exp(x);

dokladny=2.399044155917648846433327704453978787367128532409215707476694552237890741233994118308931425819765129;
%tu \/ kopiujemy na 1 pozycje
%int( 'mojafunkcja(x)',x,a,b)
%dokladny=int( 'mojafunkcja',a,b)
%\/ tu kopiujemy jw.
%  cos(9*x).*exp(x)
r1=rect('mojafunkcja(x)',a,b,12)
1-r1/dokladny
%\/ jw.
t1=trap('mojafunkcja(x)',a,b,12)

1-t1/dokladny
r2=rect('mojafunkcja(x)',a,b,120)
1-r2/dokladny
%\/ jw.
t2=trap('mojafunkcja(x)',a,b,120)
1-t2/dokladny
s2 = simp('mojafunkcja(x)',a,b,6)
1-s2/dokladny
s21 = simp('mojafunkcja(x)',a,b,60)
1-s21/dokladny
s3 = simp38('mojafunkcja(x)',a,b,4)
1-s3/dokladny
s31 = simp38('mojafunkcja(x)',a,b,40)
1-s31/dokladny


pause
clear
'zad 2'
%zad2 nowa funkcja - podstawiamy otrzymana

% 'mojafunkcja2(x)'
%fun1 = 7*x^4 - 5*x^3*log(x) + 38.234*x*sin(x/2)
%%trzeba wstawic do pliku...
% podstawiamy c i d
a = -10
b = 10
niedokladnoscQ = 0.001
%dokladny=eval(int( 'mojafunkcja2(x)',x,a,b))
dokladny = 1.253314137315500251207882642405522626503493370304969158314961788171146827303920987473297850165907899;
%%myfun - dobra funkcja...
[ q ] = quad('mojafunkcja2(x)',a,b,niedokladnoscQ,1)
% % %poprawka na niedokladnosc quad.m - dodatkowy arg funkcji
% % if q>dokladny
% %     q=q+niedokladnoscQ;
% % else
% %     q = q- niedokladnoscQ;
% % end

bq =q % blad wzgledny
1- bq/dokladny
n=8;
s2 = simp ('mojafunkcja2(x)',a,b,n/2)
1-s2/dokladny
pause
clear
'zad 3'
%trzecie i znowu nowa funkcja...
% tak jak a,b

%fun1 =2*x^2
% podstawiamy c i d
a = 0
b = 1
epsilon = 0.0000001
a=a+epsilon
b=b-epsilon
%dokladny=eval(int( fun1,x,a,b))
dokladny=0.9460830703671830149413533138231796578123379547381117904714547735666870365407979180887021330817407112;
s10=simp('mojafunkcja3(x)',a,b,5)
1-s10/dokladny
s100=simp('mojafunkcja3(x)',a,b,50)
1-s100/dokladny