Revision history [back]

Thanks for comments!

x=var('x')
f = lambda x : x^2
def defintegral_viaMathFree(f,x,a,b):
F(x)=integrate(f(x),x,algorithm='mathematica_free')
return F(b)-F(a)

And

defintegral_viaMathFree(f,x,0,1)

output: 1/3