Update solvers.f90

error blows up if phi(nx+1) is not computed