Hi,
I would like to solve a system of 2 PDE with this code, but it does not work.
lin_eqs := {diff(h(x, t), t) + H*diff(u(x, t), x) = 0, diff(u(x, t), t) + g*diff(h(x, t), x) = 0};
sol := pdsolve(lin_eqs);
print('sol', sol);

Can you help me?Thanks
Best regards