octave:2> quad(@integranda1,-20,20) warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/integranda1.m' is in the future ans = 0.88623 octave:3> clear all octave:4> quad(@integranda1,-20,20) warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/integranda1.m' is in the future ans = 0.88623 octave:5> quad(@integranda1,-30,30) ans = 0.88623 octave:6> quad(@integranda1,-40,40) ans = 1.6654e-13 octave:7> quad(@integranda1,-40,41) ans = 0.88623 octave:8> quad(@integranda1,-40,40,1e-12) ans = 1.6654e-13 octave:9> trapezi(@integranda1,-4,4,1e-12) warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/trapezi.m' is in the future warning: Impossibile raggiungere la tolleranza richiesta warning: con il numero di nodi di quadratura consentito. ans = 7.2023e-06 octave:10> trapezi(@integranda1,-4,4,1e-2) warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/trapezi.m' is in the future ans = 7.2023e-06 octave:11> trapezi(@integranda1,-4,4,1e-2,10) warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/trapezi.m' is in the future ans = 0.88623 octave:12> trapezi(@integranda1,-4,4,1e-2,4) warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/trapezi.m' is in the future warning: Impossibile raggiungere la tolleranza richiesta warning: con il numero di nodi di quadratura consentito. ans = 0.81672 octave:13> trapezi(@integranda1,-4,4,1e-2,2) warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/trapezi.m' is in the future ans = 7.2023e-06 octave:14> [I,stimaerrore]=trapezi(@integranda1,-4,4,1e-2,2) warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/trapezi.m' is in the future I = 7.2023e-06 stimaerrore = 3.6011e-06 octave:15> [I,stimaerrore]=trapezi(@integranda1,-4,4,1e-2,4) warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/trapezi.m' is in the future warning: Impossibile raggiungere la tolleranza richiesta warning: con il numero di nodi di quadratura consentito. I = 0.81672 stimaerrore = 0.39289 octave:16> [I,stimaerrore,x]=trapeziadatt1(@integranda3,-pi,pi,1e-3); warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/trapeziadatt1.m' is in the future warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/integranda3.m' is in the future octave:17> I I = -0.088253 octave:18> stimaerrore stimaerrore = 3.8456e-04 octave:19> length(x) ans = 1025 octave:20> plot(x,integranda3(x)) warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/integranda3.m' is in the future octave:21> axis([0,pi,-pi,pi]) octave:22> plot(x,integranda3(x),'x') warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/integranda3.m' is in the future octave:23> axis([0,pi,-pi,pi]) octave:24> [I,stimaerrore,x]=trapeziadatt2(@integranda3,-pi,pi,1e-3); warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/trapeziadatt2.m' is in the future octave:25> I I = -0.088384 octave:26> stimaerrore stimaerrore = 0.0032899 octave:27> length(x) ans = 849 octave:28> [I,stimaerrore,x]=trapeziadatt1(@integranda3,-pi,pi,1e-3); warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/trapeziadatt1.m' is in the future octave:29> stimaerrore stimaerrore = 3.8456e-04 octave:30> length(x) ans = 1025 octave:31> [I,stimaerrore,x]=trapeziadatt2(@integranda3,-pi,pi,1e-3); warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/trapeziadatt2.m' is in the future octave:32> plot(x,integranda3(x),'x') warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/integranda3.m' is in the future octave:33> axis([0,pi,-pi,pi]) octave:34> plot(x,integranda3(x)) warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/integranda3.m' is in the future octave:35> plot(x,zeros(size(x)),'*',x,integranda3(x)) warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/integranda3.m' is in the future octave:36> axis([0,pi,-pi,pi]) octave:37> octave:37> octave:37> testgauss warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/integranda1.m' is in the future warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/r_jacobi.m' is in the future warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/quadratura/gauss.m' is in the future error = 2.7756e-15 octave:38> help quad -- Loadable Function: [V, IER, NFUN, ERR] = quad (F, A, B, TOL, SING) Integrate a nonlinear function of one variable using Quadpack. The first argument is the name of the function, the function handle or the inline function to call to compute the value of the integrand. It must have the form y = f (x) where Y and X are scalars. The second and third arguments are limits of integration. Either or both may be infinite. The optional argument TOL is a vector that specifies the desired accuracy of the result. The first element of the vector is the desired absolute tolerance, and the second element is the desired relative tolerance. To choose a relative test only, set the absolute tolerance to zero. To choose an absolute test only, set the relative tolerance to zero. The optional argument SING is a vector of values at which the integrand is known to be singular. The result of the integration is returned in V and IER contains an integer error code (0 indicates a successful integration). The value of NFUN indicates how many function evaluations were required, and ERR contains an estimate of the error in the solution. You can use the function `quad_options' to set optional parameters for `quad'. It should be noted that since `quad' is written in Fortran it cannot be called recursively. /usr/local/octave-3.0.2/libexec/octave/3.0.2/oct/i686-pc-linux-gnu/quad.oct Additional help for built-in functions and operators is available in the on-line version of the manual. Use the command `doc ' to search the manual index. Help and information about Octave is also available on the WWW at http://www.octave.org and via the help@octave.org mailing list. octave:39> [I,ier,nfun]=quad(@integranda1,-1,1); octave:40> nfun nfun = 21 octave:41> diary off