I am working on differential equation .

The solution with no initial conditions is

dsolve(‘Dy = (x-y)/(x+y)’, ‘x’)

ans =

(-x*C1+(2*x^2*C1^2+1)^(1/2))/C1
(-x*C1-(2*x^2*C1^2+1)^(1/2))/C1

Given the initial condition ‘y(0) = 3′, and setting my solution equal to “sol”,

sol = dsolve(‘Dy = (x-y)/(x+y)’, ‘y(0) = 3′, ‘x’)

sol =

-x+3*(2/9*x^2+1)^(1/2)

Using ezplot i can show the graph with the initial condition.

ezplot(sol, [0 2])