I am trying to solve the following system of equations in MATLAB:
S = solve([real(integral(@(r)sqrt(qn2(x+iy,r)),0.6,1)/pi) == 1, imag(integral(@(r)sqrt(qn2(x+iy,r)),0.6,1)/pi)== 0], [x, y])
Where qn2 is a function handle. MATLAB has trouble because it is not happy evaluating the integral due to the fact that $x$ and $y$ have to be symbolic.
I am using MATLAB as it is a lot faster at numerically evaluating the integral that Mathematica, Maple etc.
Any advice on how to best solve these would be great!