% MATLAB function that coordinates the calculations for Example 33.1 of "A
% First Course on Kinetics and Reaction Engineering."
%
function Example_33_1
    % Solve the PFR design equation
    [PFR_length,PFR_conv] = Example_33_1_pfr;
    % Solve the axial dispersion design equation
    [ax_length,ax_conv] = Example_33_1_ad;
    % Make the plot
    figure
    plot(PFR_length,PFR_conv,ax_length,ax_conv)
    xlabel('Length (m)')
    ylabel('Percent Conversion')
    legend('PFR','Axial Dispersion Model')

end % of Example_33_1.m

