Solves laplace in 1 dimension: -d2u/dx2 =  (lambda)u using petsc4py and slepc4py, a matrix-free method.