I found a case where PDEtools:-Solve fail and solve works. I do not understand why.
Is this a bug or what is going on? I thought they should give same result.



> |
eqs:=[-2*arctanh((_C1-1)^(1/2)/_C1^(1/2)) = _C2, 0 = _C1^(3/2)*tanh(1/2*_C2)*sech(1/2*_C2)^2];
unknowns:={_C1,_C2};
|
![[-2*arctanh((_C1-1)^(1/2)/_C1^(1/2)) = _C2, 0 = _C1^(3/2)*tanh((1/2)*_C2)*sech((1/2)*_C2)^2]](/view.aspx?sf=238811_question/170a2f7292cb5210526f44472c65acb8.gif)

> |
sol1:=solve(eqs,unknowns);
|

> |
sol2:=PDEtools:-Solve(eqs,unknowns);
|

|
Download why_PDEtools_Solve_fail_august_14_2024.mw