xy(dy)/(dx)+4x^2+y^2=0
⚠ 4x is read as the product 4*x: for a power, write 4^x
⚠ xy is read as the product x*y: for a power, write x^y; for the variable xy, write xy_
⚠ xy(dy)/(dx)+4x^2+y^2=0 is read as xy*y'+4x^2+y^2=0: use dif(expression, variable) for a derivative with respect to that variable, use (numerator)/(denominator) for a quotient
Let z = y2: z' = 2y·y'. On intervals where x ≠ 0, the equation becomes z' + 2x·z = -8x
The integrating factor is x2, nonzero for x ≠ 0: the product rule gives (x2z)' = -x3·8
⇥Let's compute the integral: ∫(-x3·8)dx
⇥ -8·x44
⇥ℹ-8·x44 = -8x44‖-8x44 = -2x4
⇥ -2x4
Integrating gives x2z = -2x4 + C: y2 = -2x4 + Cx2, where C is an arbitrary real constant
The real branches have a fixed sign on each interval where -2x4 + Cx2 > 0 and x ≠ 0; these restrictions exclude singular points
▶Solution:   y = { --2x4 + Cx2, -2x4 + Cx2 }