Como você resolveu a EDO usando (ax+b) ao invés de (ax-b) (e analogamente para o termo envolvendo c e d), para efeitos de comparação com o Maple, deve-se considerar a solução a seguir.
# declara variáveis simbólicas: variável independente e parâmetrosx,a,b,c,d=var('x,a,b,c,d')# impõe que todos os parâmetros são estritamente positivosassume(a>0)assume(b>0)assume(c>0)assume(d>0)# declara a função incógnitap=function('p')(x)# define a equação diferencialed=((a*x+b)^2)*diff(p,x,2)+(c*x+d)*diff(p,x,1)==1# obtém a solução da EDO e exibesol=desolve(ed,p,ivar=x)show(sol)# simplifica a solução (elimina os logaritmos, essencialmente) e exibesol1=sol.canonicalize_radical()show(sol1)forget()