Switch to numpy and add the start of a simple solver.