In[]:=
(*Findzforfixedvaluesofxandy;forexample(x=0,y=0)*)
fun1=fun/.{x0,y0};
Solve[fun10,z]