Question: Why I receive this peculiar different result?

I have 2 same integrals, but there is difference between  them at the end.WHy?


 

``

restart

N := 8:

I__1 := add(int(zt*ro[i], zt = Z[i] .. Z[i+1]), i = 1 .. N)

(1/2)*ro[1]*(-Z[1]^2+Z[2]^2)+(1/2)*ro[2]*(-Z[2]^2+Z[3]^2)+(1/2)*ro[3]*(-Z[3]^2+Z[4]^2)+(1/2)*ro[4]*(-Z[4]^2+Z[5]^2)+(1/2)*ro[5]*(-Z[5]^2+Z[6]^2)+(1/2)*ro[6]*(-Z[6]^2+Z[7]^2)+(1/2)*ro[7]*(-Z[7]^2+Z[8]^2)+(1/2)*ro[8]*(-Z[8]^2+Z[9]^2)

(1)

ro := Vector(8, {(1) = 2051.904, (2) = 2051.904, (3) = 2051.904, (4) = 2051.904, (5) = 2051.904, (6) = 2051.904, (7) = 2051.904, (8) = 2051.904})

Z := Vector(9, {(1) = -0.4998720000e-6, (2) = -0.3749040000e-6, (3) = -0.2499360000e-6, (4) = -0.1249680000e-6, (5) = 0, (6) = 0.1249680000e-6, (7) = 0.2499360000e-6, (8) = 0.3749040000e-6, (9) = 0.4998720000e-6})

I__2 := add(int(zt*ro[i], zt = Z[i] .. Z[i+1]), i = 1 .. N)

0.1e-19

(2)

I__1-I__2

-0.1e-19

(3)

``

``

``

``

``

``

``


 

Download Ajib.mw

Please Wait...