>>>>> "David" == David Billinghurst <billingd at users.sourceforge.net> writes:
David> Update of /cvsroot/maxima/maxima/share/contrib/diffequations/tests
David> In directory sc8-pr-cvs7.sourceforge.net:/tmp/cvs-serv16616
David> Modified Files:
David> rtestode_kamke_2_4.mac
David> Log Message:
David> Improvements to sin.lisp rev 1.18 have resulted in better answers to
David> Kamke 2.337 2.270 2.375 and 2.376.
David> (pn_(337),ans:contrib_ode(eqn:'diff(y,x,2)=-(3*x+2*b+a)*'diff(y,x,1)/(2*(x+a)*(x+b))-(a-b)*y/(4*(x+a)^2*(x+b)),y,x));
David> -[y=%k1*sin(sqrt(a-b)*'integrate(1/sqrt(4*x^3+4*b*x^2+8*a*x^2+8*a*b*x+4*a^2*x+4*a^2*b),x))+%k2*cos(sqrt(a-b)*'integrate(1/sqrt(4*x^3+4*b*x^2+8*a*x^2+8*a*b*x+4*a^2*x+4*a^2*b),x))];
David> +[y=%k1*sqrt(x+b)/(sqrt(a-b)*sqrt((x+b)/(a-b)+1))+%k2/sqrt((x+b)/(a-b)+1)];
Do you have radexpand:all here?
The first integral can be written
integrate(1/sqrt(4*(x+a)^2*(x+b)),x), and this can only be simplified
if (x+a)>0.
I just want to make sure integration is still doing the right thing.
Ray