Usage example for RK_stab_fun Example of stability domain boundary for an implicit method
0001 % Usage example for RK_stab_fun 0002 % Example of stability domain boundary for an implicit method 0003 0004 b = [1/4; 3/4]; 0005 c = [0; 2/3]; 0006 A = [1/4 -1/4; 0007 1/4 5/12]; 0008 [ R ] = RK_stab_fun( b,c,A ) 0009 [X,Y] = meshgrid(-1:0.1:7,-5:0.1:5); 0010 Z = X+1i*Y; 0011 Rmod = abs(R); 0012 Rmod = matlabFunction(Rmod); 0013 boundary = feval(Rmod,Z); 0014 contour(X,Y,boundary,[1 1]); 0015 grid on; 0016 title(['implicit example']);