Q = quad(FUN,A,B) tries to approximate the integral of scalar-valued
function FUN from A to B to within an error of 1.e-6 using recursive
adaptive Simpson quadrature.
Example: Q = quad(@(x)myfun2(x,5),0,2);
Where are your A and B Parameters?
Are -500 and 500 supposed to be A and B? You have them trapped inside of exp(ixw-abs(w)^2,-500,500) so that means you don't have enough input arguments