how do i use fsolve to solve a stream function ?

I have a stream function and am asked to solve it numerically to
calculate velocity at the x and y co-ordinates and then plot the pressure
co-efficient on a rankin oval. But am unsure of how to use the fsolve function to 
solve the stream function . 

the code i have so far:

syms U0 x y real

psi = U0*y + 4/2/pi * atan2(y, x+1) - 4/2/pi * atan2 (y, x-1);

v=-diff(psi, x)

u=diff(psi, y)

U0=1

x=-3:0.01:3;

y=-3:0.01:3;

u=eval(u)

v=eval(v)

any help is appreciated.

[Ameer Hamza]

Can you shows the mathematical form of the equation?

[steven hong]

[Thiago Henrique Gomes Lobato]

Why you need to "solve" the stream function since it is already there? There's nothing to solve. Most likely what was meant was discretize/numeric calculate it, and your code seems fine so far.

[steven hong]

Ive added the full question to make it clear what im trying to achieve.

ive added the following to get the velocity magnitude after calculating u and v.

Z=sqrt(u.^2 + v.^2)

then after this i have to calculate pressure and the pressure co-efficient on the surface of the rankin oval using bernoullis equation.