See the following example:
from dolfin import *
# Set up some simple demonstration nonlinear problem:
mesh = UnitIntervalMesh(10)
V = FunctionSpace(mesh,"CG",1)
u = Function(V)
v = TestFunction(V)
x = SpatialCoordinate(mesh)
F = (u+1)**2*dot(grad(u),grad(v))*dx + u*v*dx - sin(x[0])*v*dx
problem = NonlinearVariationalProblem(F,u,J=derivative(F,u),bcs=[])
# Create solver:
solver = NonlinearVariationalSolver(problem)
stype = 'newton'
solver.parameters['nonlinear_solver']=stype
sprms = solver.parameters[stype+'_solver']
# Set maximum iterations:
sprms['maximum_iterations'] = 100
# Force solver to reach maximum iterations for demonstration:
sprms['relative_tolerance'] = 0.0
sprms['absolute_tolerance'] = 0.0
# Print out other available options:
info(solver.parameters,True)
# Run solver:
solver.solve()