def METHOD_NAME():
    # Initialize
    problem      = MyProblem()
    initGuess    = problem.getSolution()
    initGuess[:] = [.5,.5]
    jacobian     = problem.getJacobian()
    solver       = NOX.Epetra.defaultSolver(initGuess,
                                            problem,
                                            problem,
                                            jacobian,
                                            maxIters=40)
    print(solver.statusTest)
    # Solve
    solveStatus = solver.solve()
    print()
    if solveStatus != NOX.StatusTest.Converged:
        print("Nonlinear solver failed to converge")
        print("End Result: TEST FAILED")
    else:
        group = solver.getSolutionGroup()
        soln  = group.getX()
        print("Solution =", soln)
        print("End Result: TEST PASSED")