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