DSolve[y[x] - (x + 3*y[x])*Derivative[1][y][x] + (5 + 3*x)*Derivative[1][y][x]^2 == 0,y[x],x]