3.solvepde函数