The determination of profiles and bifurcation diagrams of tubular reactors and CSTR cascade is described. This method makes use of the fictitious dynamics and the continuation numerical process. It can be applied for different reactor configurations and chemical reactions. In this method the optimization procedure to calculate the boundary conditions is not necessary.