This study shows how the optimal control theory for distributed parameter systems can be implemented for a problem of tubular reactor with axial dispersion described by partial differential equations. Two methods are implemented. One is based on differential equation approach and the other is based on integral equation approach. It was found that the approach with partial differential equations is preferable to the one with integral equations for the type of problems treated in this study. Computation algorithms and programs for both cases are developed.