Home > K25104 > ODEs > RK_stab_fun_script_4.m

RK_stab_fun_script_4

PURPOSE ^

Usage example for RK_stab_fun

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 Usage example for RK_stab_fun
 Boundary of stability domain of Nystrom's method

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 % Usage example for RK_stab_fun
0002 % Boundary of stability domain of Nystrom's method
0003 
0004 b = [1/4; 3/8; 3/8];
0005 c = [0; 2/3; 2/3];
0006 A = [0 0 0;
0007      2/3 0 0;
0008      0 2/3 0];
0009 [ R ] = RK_stab_fun( b,c,A )
0010 [X,Y] = meshgrid(-3:0.1:0.5,-3:0.1:3);
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(['Nystrom',char(39),'s method']);

Generated on Mon 18-Jan-2016 10:25:49 by m2html © 2005