Numerical solver for a specific type of PDEs