solve {y'+2xy=0, y(1)=2} implicit midpoint method x=1..3