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])
