def METHOD_NAME(self):
    prob = om.Problem()
    prob.model.add_subsystem('comp', om.ExecComp('y=3*x+1'), promotes=['*'])
    balance = prob.model.add_subsystem('balance', om.BalanceComp(),
                                       promotes=['*'])
    balance.add_balance('x', lower=-.1, upper=10, rhs_val=0, lhs_name='y')
    nl_solver = prob.model.nonlinear_solver = om.BroydenSolver()
    nl_solver.options['stall_limit'] = 3
    nl_solver.options['stall_tol'] = 1e-8
    nl_solver.options['maxiter'] = 100
    prob.model.linear_solver = om.DirectSolver()
    prob.setup()
    prob.set_solver_print()
    prob.run_model()