@mmcdara

Thank you very much for your work, that you were not too lazy to code the program.

I expected that the solution of the system of PDEs would be the graph z(x, b_n) in space R3.

As I understand it, having studied your code, you simply depict sections of this figure on a plane.

As I understand it, the Maple solver cannot solve differential equations of this kind numerically, so I decided to use Matlab in my work.

Do you have experience in solving similar problems in Matlab?

Thx u very much.