Per calcolare con questo script l'integrale tra -1 e 6/7 di √(1-x³) calcolo prima il valore numerico di 6/7 (arrotondato a 17 cifre dopo ".") con la calcolatrice (0.8571428571428571).  Ecco il valore dell'integrale che si può ottenere con WolframAlpha e poi quello ottenibile con lo script.

Fermandomi qui posso prendere l'arrotondamento 1.89305726590.

Una alternativa è ricorrere ad un programmino in Basic, facilmente modificabile per calcolare qualunque integrale:

10 a = -1 : b = 6/7 : d = b-a
20 n=54321 : 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
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 = sqr(1-x^3) : return
n = 54321 integrale = 1.8930572659387934
n = 108642  integrale = 1.8930572659113554  variaz.= -2.743805183058612e-11
n = 217284  integrale = 1.8930572659045026  variaz.= -6.852740597196316e-12
n = 434568  integrale = 1.8930572659028158  variaz.= -1.6868728636154628e-12
n = 869136  integrale = 1.8930572659023528  variaz.= -4.629630012686903e-13
n = 1738272  integrale = 1.8930572659022316  variaz.= -1.212363542890671e-13
n = 3476544  integrale = 1.8930572659023248  variaz.= 9.325873406851315e-14

Fino alla penultima uscita le variazioni (negative) si dividono circa per 4, poi intervengono errori di arrotondamento. Assumo l'approssimazione 1.89305726590223.

Per calcolare con questo script l'integrale di √(1-x³) tra -1 e ³√(6/7) posso calcolare la radice cubica usando la funzione cbrt (cube root): vedi l'help della "calcolatrice".

[0.94991425159300 è l'arrotondamento a 17 cifre]

Posso prendere 1.94000795397.  WolframAlpha con integrate sqrt(1-x^3) form -1 to (6/7)^(1/3) avrebbe fornito 1.9400079539688392728741876759981085...

In Basic:

10 a = -1 : b = (6/7)^(1/3)  : d = b-a
...
n = 54321 integrale = 1.94000795410416
n = 108642  integrale = 1.9400079540026602  variaz.= -1.0149969753570076e-10
n = 217284  integrale = 1.9400079539773027  variaz.= -2.5357493882438575e-11
n = 434568  integrale = 1.9400079539709083  variaz.= -6.394440532631052e-12
n = 869136  integrale = 1.9400079539693733  variaz.= -1.5349943538467414e-12
n = 1738272  integrale = 1.940007953968869  variaz.= -5.042632977847461e-13
n = 3476544  integrale = 1.9400079539687958  variaz.= -7.327471962526033e-14

Fino alla penultima uscita le variazioni (negative) si dividono circa per 4, poi intervengono errori di arrotondamento. Assumo l'approssimazione 1.9400079539688.