begin library A0, A6; real sum; integer i, j; array A [1:10000]; for i := 1 step 1 until 10000 do A[i] := 1.0/i; sum := 0.0; for i := 1 step 1 until 10000 do sum := sum + A[i]^2; output(30, sqrt(sum*6.0)); comment the scalar product approximates pi^2 / 6; end |