Mathematica script implementing the Riccati System method for solving Linear PDEs