added a normal() to compare the output from lsolve with what is expected