Numerical integration

Define a function f(x) to integrate on [0,1] .

> f := x->(16-16*x^2)^(1/2);

[Maple Math]

> readlib(`evalf/int`):

> Digits := 13;

[Maple Math]

> F1 := `evalf/int`(f,0..1,Digits);

[Maple Math]

Try again at higher precision

> Digits := 36;

[Maple Math]

> F2 := `evalf/int`(f,0..1,Digits);

[Maple Math]

Compare with the exact answer .

> pi := evalf(Pi);

[Maple Math]

> zero_3 := pi-F2;

[Maple Math]