odr := diff(y(x),x)=x+y(x);
dsolve(odr,y(x));
dsolve({odr,y(0)=1},y(x));
sol := dsolve({odr,y(0)=2},y(x));
plot(rhs(sol),x=-1..3);
sol := dsolve({odr,y(0)=c},y(x));
Y := unapply(rhs(sol),x,c);
plot({Y(x,c)$c=-2..4},x=0..3);
odr1 := diff(y(x),x$2) = - y(x);
sol1 := dsolve({diff(y(x),x$2) = - y(x),y(0)=1,D(y)(0)=4},y(x));
plot(rhs(sol1),x=0..40);
odr2 := diff(y(x),x) = x*y(x)/(1+x*y(x)^2);
dsolve(odr2,y(x));
odr3 := diff(y(x),x) = y(x)/(1+x*y(x)^2);
dsolve(odr3,y(x));
odr4 := diff(y(x),x) = x*y(x)/(1+y(x)^2);
dsolve(odr4,y(x));
odr5 := diff(y(x),x) = x*y(x)/(1+x*y(x));
ae5 := dsolve(odr5,y(x));
ae50 := subs(_C1=0,ae5);
solve(ae05,y(x));