>> syms x y0>> y=dsolve('Dy=y+1/y','y(0)=y0','x')y =(-1+exp(2*x)*(1+y0^2))^(1/2)-(-1+exp(2*x)*(1+y0^2))^(1/2)>> help dsolveDSOLVE Symbolic solution of ordinary differential equations.DSOLVE('eqn1','eq...