import math , random def int(f,a,b,N): sum=0.0 start=float(a) step=float(b-a) for i in range(N): sum+= f(start+step*random.random())*step/N return sum f=lambda x: math.sqrt(1-x**2) a=0 b=1 N=input("Number of steps: ") print int(f,a,b,N)-math.pi/4