It comes out OK.
Write y as a poynomial
[tex]y = a_0 + a_1.x + a_2.x^2 + a_3.x^3 + a_4.x^4 + ...[/tex]
Differentiate
[tex]y' = a_1 + 2a_2x + 3a_3x^2 + 4a_4x^3 + ...[/tex]
Substitute
[tex]xy' - 3y = a_1x + 2a_2x^2 + 3a_3x^3 +4a_4x^4 + ... - 3a_0 - 3a_1x - 3a_2x^2 - 3a_3x^3 - 3a_4x^4 + ...[/tex]
[tex]xy' - 3y = (-3a_0) + (a_1 -3a_1)x + (2a_2 - 3a_2)x^2 + (3a_3 - 3a_3)x^3 + (4a_4 - 3a_4)x^4 + ...[/tex]
[tex]\mbox{substituting for } xy' - 3y = k,[/tex]
[tex]k = (-3a_0) + (a_1 -3a_1)x + (2a_2 - 3a_2)x^2 + (0)x^3 + (4a_4 - 3a_4)x^4 + ... ----------------------(1)[/tex]
Analysis
[tex]\mbox{For the lhs to equal the rhs, } k = -3a_0 \mbox{ (a constant) and all the other coefficents must be zero: }a_n = 0, n \in N, n \neq 0,3.[/tex]
ergo,
[tex]y = a_0 + a_3x^3[/tex]
or
[tex]y = -k/3 + Cx^3[/tex]
============
As long as the original DE is made up of powers of x only (I think it would be difficult for trig functions and exponentials) then you're just manipulating the coefficents in the rhs of (1) to give your answer.