# simeq.py   solve simultaneous equations AX=Y 
# solve real linear equations for X where Y = A * X
# usage:  X = simeq(A,Y)

from numpy import array
from numpy.linalg import solve

def simeq(A,Y):
  n = len(Y)
  a = array([[A[i][j] for j in range(n)]for i in range(n)])
  y = array([Y[j] for j in range(n)])
  x=solve(a,y)
  X = [x[j] for j in range(n)]
  return X;