I have system of differential equations, for which defined boundary condition. I do not know how implement this in maxima. What how I try do this:
F: desolve(['diff(Ep(r), r)=(Er(r)-Ep(r))/r, 'diff(Sr(r), r)=(Sp(r)-Sr(r))/r-r, Sr(1)=0, Sp(1)=3.625], [Sr(r), Sp(r)]);
This return: desolve: can't handle this case.