How to get the solution obtained by fenics as a vector?

2

Consider the following minimal example:

from dolfin import *

mesh = UnitSquareMesh(2, 2)
V = FunctionSpace(mesh, "CG", 2)
dof_coordinates = V.tabulate_dof_coordinates()

expr = Expression("x[0]", degree=1)
uh = Function(V)
uh.interpolate(expr)

print(dof_coordinates)
print(uh.vector().get_local())