(* eulerstab.wl -- The linear stability domain for Euler's method *)

r=1+z
ra=Abs[r]
ra2=ra/.z->a+I*b
ContourPlot[ra2,{a,-3,1},{b,-2,2},
   Contours->{1},ContourShading->{Red,Blue}]
