Recent content by Ally1

Thank for your respond. To clear my question, let me take an example of the Exponential Time Differencing-Runge-Kutta Method [F. de la Hoz and F. Vadillo Computer Physics Communications 179 (2008) 449–456 ], I tried to code the stability regiong of this method as you can see below but the code...

Ok thanks but I am more concern about point 1).

In many cases one finds code such as x = linspace(-3,1.5,200); y = linspace(-3.5,3.5,200); [X,Y] = meshgrid(x,y); Z = X +Y*i; %Euler's Method M = abs(1+Z); [c,h] = contour(X,Y,M,[1,1]); set(h,'linewidth',2,'edgecolor','g') to plot the stability region of the Euler's Method, where in fact the...