f=inline('1-6*x.^2*y')
i=dblquad(f,0,2,-1,1)
i