static Standard_Integer findplanarsurface(Draw_Interpretor&, Standard_Integer n, char** a)
{
  if (n<3) {
    cout << "bad number of arguments" <<endl;
    return 1;
  } 
  
  // try to read a shape:
  TopoDS_Shape inputShape=DBRep::Get(a[2]);
  if (inputShape.IsNull() || inputShape.ShapeType() != TopAbs_WIRE) {
    cout << "Invalid input shape"<< endl;
    return 1;
  }
  Standard_Real toler = atof(a[3]);
  TopoDS_Wire aWire = TopoDS::Wire(inputShape);
  BRepLib_FindSurface FS(aWire, toler, Standard_True);
  if(FS.Found()) {
    cout<<"Planar surface is found"<<endl;
    Handle(Geom_Surface) aSurf = FS.Surface();
    BRepBuilderAPI_MakeFace aMakeFace(aSurf,aWire,Standard_True);
    if(aMakeFace.IsDone()) {
      TopoDS_Face aFace = aMakeFace.Face();
      DBRep::Set(a[1],aFace);
    }
  }
  else cout<<"Planar surface is not found with toler = "<<toler <<endl;
  return 0;
  
}
theCommands.Add("findplanarsurface","findplanarsurface Result wire Tol",findplanarsurface,g);