SUBROUTINE FCN

In order to be able to utilize the NLEODE program code, the operator must, as a minimum, write the program code for the subroutine FCN, in which the functions for the right-hand side of Equations A3.1 and A3.2 are defined. If the Jacobian matrix of the right-hand side in Equations A3.1 and A3.2 is calculated analytically, an additional subroutine, FCNJ, must be written by the operator. In case the elements in the Jacobian matrix are... 

Parameters ki and k2 obtain constant values. The value of the independent variable x (residence time) should be between 0 and 5. Subroutines FCN and FCNJ are listed below, as well as the numerical results of program code execution. 

The differential equations for the remaining nodes (3 through 9) are given by Eq. 5.113. These equations are written in subroutine FCN. The values of 0 at the boundaries, 0i and 011, are calculated after each time step (or call to DIVPAG) using Eq. 5.115 for 0i and the following equation for 011 (seeEq. 5.117) ... 

---