Para analisar a influência da ordem da reação \(n\), os demais parâmetros - do modelo (\(\Phi\)) e dos métodos numéricos (\(\Sigma\) e \(N\)) - são considerados constantes. O modelo de interesse é dado pela seguinte equação:

\[\begin{gathered} \frac{d^2C}{dR^2}+\frac{2}{R}\frac{dC}{dR}=\Phi^2C^n\\ \frac{dC}{dR}=0,\;\;\;\;\;R=0\\ C=1,\;\;\;\;\;R=1\end{gathered}\]

A entrada \(n\) deverá ser um vetor com elementos separados por vírgula com tamnho variando entre 2 e 5.