solve {y'+2xy=0, y(1)=2} Dormand‐Prince method x=1..3