Double_t fcn(Double_t *xi, Double_t *par)
{
   Double_t x = xi[0];
   Double_t y = xi[1];
   Double_t z = xi[2];
   Double_t f1 = x**2+y**2+sqrt((par[0]**2+par[1]**2)*(x**2+y**2)+4*par[0]*par[1]*x*y)-z;
   Double_t f2 = x**2+y**2-sqrt((par[2]**2+par[3]**2)*(x**2+y**2)+4*par[2]*par[3]*x*y)-z;
   return f1*f2;
}

void fun3func()
{
   TF3 *fun3 = new TF3("fun3", fcn,-6,6,-6,6,-6,6,4);  
   fun3->SetParameter(0, 2);
   fun3->SetParameter(1, 0);
   fun3->SetParameter(2, 2);
   fun3->SetParameter(3, 0);

   fun3->Draw(); 
}