integral

I am getting the following error code when trying to evaluate a definite integral in Python. AttributeError Traceback (most recent call last) <ipython-input-7-2be8718c68ec> in <module>() 7 x, n = symbols('x n') 8 ----> 9 f = (cos(n*x))*(x**2-pi**2)^2 10 integrate(f,(x,-n*pi,n*pi)) 11 AttributeError: 'Mul' object has no attribute 'cos' I have copied my input code below. Thanks for any help. from pylab import * from sympy import * from numpy import * init_printing(use_unicode=False, wrap_line=False, no_global=True) x, n = symbols('x n') f = (cos(n*x))*(x**2-pi**2)^2 integrate(f,(x,-n*pi,n*pi))

I'm looking for a solution for a double integral that is faster than integrate(function(y) { sapply(y, function(y) { integrate(function(x) myfun(x,y), llim, ulim)$value }) }, llim, ulim) with eg myfun <- function(x,y) cos(x+y) llim <- -0.5 ulim <- 0.5 I found an old paper that referred to a FORTRAN program called quad2d , but I couldn't find anything else but some help pages for matlab for the rest. So I'm looking for a C or FORTRAN library that can do double integrals quick (i.e. without the sapply loop), and that can be called from R. All ideas are very much appreciated, as long as they're