solving PDE