"Solución de la ecuación diferencial exacta: y dx+x dy = 0, bajo la condición de que y(1)=1,5 y 9";
"La solución es";
sol=DSolve[y[x]+x *y'[x]==0,y[x],x];
sol[[1,1,2]]

"El valor de la función en y(1) es";
sol[[1,1,2]]/.x->1

"Las trayectorias solución de la anterior ecuación para valores iniciales y(1)=1,3,5 son las que siguen";
trajectorias=Plot[Evaluate[sol[[1,1,2]]/.C[1]->{1,5,9},{x,0,10}],PlotRange->{0,10},AxesOrigin->{0,0},AxesLabel->{x,y},AspectRatio->1,PlotStyle->Thick]