

Una alternativa è ricorrere ad un programmino in Basic, facilmente modificabile per calcolare qualunque integrale:
10 a = 0 : b = 100 : d = b-a 20 n=4321 : gosub 400 : print "n = "; n ; " integrale = "; I 30 I0=I: n=n*2: gosub 400: print "n = "; n ; " integrale = "; I; " variaz.= ";I-I0 40 I0=I: n=n*2: gosub 400: print "n = "; n ; " integrale = "; I; " variaz.= ";I-I0 50 I0=I: n=n*2: gosub 400: print "n = "; n ; " integrale = "; I; " variaz.= ";I-I0 60 I0=I: n=n*2: gosub 400: print "n = "; n ; " integrale = "; I; " variaz.= ";I-I0 70 I0=I: n=n*2: gosub 400: print "n = "; n ; " integrale = "; I; " variaz.= ";I-I0 80 I0=I: n=n*2: gosub 400: print "n = "; n ; " integrale = "; I; " variaz.= ";I-I0 90 I0=I: n=n*2: gosub 400: print "n = "; n ; " integrale = "; I; " variaz.= ";I-I0 150 end 400 h=d/n: s=0: for j=1 to n: x=a+(j-1/2)*h: gosub 500: s=s+y: next: I=s*h: return 500 y = 34.7-0.611*exp(-0.05*x)-0.0244*x*exp(-0.05*x) : return n = 4321 integrale = 3448.496912006572 n = 8642 integrale = 3448.4969119181005 variaz.= -8.847155186231248e-8 n = 17284 integrale = 3448.4969118959652 variaz.= -2.2135282051749527e-8 n = 34568 integrale = 3448.4969118904346 variaz.= -5.530637281481177e-9 n = 69136 integrale = 3448.4969118890626 variaz.= -1.3719727576244622e-9 n = 138272 integrale = 3448.4969118886875 variaz.= -3.751665644813329e-10 n = 276544 integrale = 3448.496911888631 variaz.= -5.638867150992155e-11 n = 553088 integrale = 3448.4969118886597 variaz.= 2.864908310584724e-11
Le variazioni (negative) si dividono via via circa per 4, fino alla penultima uscita. Posso prendere l'approssimazione 3448.4969118886.


