Unless I’m missing something, everything you need except time discretisation is covered in the poisson demo.
If you really need to prescribe both a Neumann and Dirichlet boundary condition on the same boundary, then things get tricky. See for example here.