from visad.python.JPythonMethods import *
# construct a 2 x 2 matrix in a VisAD Field
matrix = field([[1, 2], [1, 3]])
# construct a 2 vector in a VisAD Field
vector = field([2, 1])

# solve the linear system
solution = solve(matrix, vector)

# print the solution
print solution[0], solution[1]

# prints 4.0 -1.0
#
# note
# 1  2       4       2
#        *       =
# 1  3      -1       1