% script du corrigé du TP3 % ATTENTION, calculs symbolique longs !!! Nmin=1; Nmax=5; N=1000; syms x; A=sym(0); B=sym(pi)/2; f='sin'; disp(int_fcn(4,3,A,B,f)); pause; Iexa=double(int(feval(f,x),A,B)); logerreur(double(A),double(B),f,Nmin,Nmax,N,Iexa); pause; A=0; B=sym(1); f='sqrt'; Iexa=double(int(feval(f,x),A,B)); logerreur(double(A),double(B),f,Nmin,Nmax,N,Iexa);