Usage example for RK_stab_fun related to Exercise 6.15
0001 % Usage example for RK_stab_fun 0002 % related to Exercise 6.15 0003 0004 alpha = 1/3; 0005 b = [1/2; 1/2]; 0006 c = [0; 1]; 0007 A = [0 0; 0008 1-alpha alpha]; 0009 [ R ] = RK_stab_fun( b,c,A ) 0010 [X,Y] = meshgrid(-7:0.1:0.5,-5:0.1:5); 0011 Z = X+1i*Y; 0012 Rmod = abs(R); 0013 Rmod = matlabFunction(Rmod); 0014 boundary = feval(Rmod,Z); 0015 contour(X,Y,boundary,[1 1]); 0016 grid on; 0017 title(['Exercise 6.15']);