Implementation code
syms y(x),D2y=diff(y,2); Dy=diff(y, 1);
"Analytical solutions of ordinary differential equations"
y=dsolve(( 1+x^2)*d2y==2*x*dy,y(0)== 1,dy(0)==3)