A Função Erro e Outras Funções Relacionadas

A função erro $erf(x)$ tem uma longa história que começa com os artigos de De Moivre (1718-1733) e Laplace (1774), onde foi expressa através da seguinte integral:

\begin{equation*} \int e^{-t^2}dt \end{equation*}

Mais tarde, Kramp (1799) utilizou esta integral para a definição da função de erro complementar $erfc(x)$.

A integral de probabilidade foi assim chamada porque é amplamente utilizada na teoria da probabilidade e estatística, aparecendo também em equações diferenciais parciais.

A função erro imaginário $erfi(z)$ também chamada de Função Erro de Gauss, foi desenvolvida para calcular a integral da distribuição normal. Gauss mostrou como a probabilidade pode ser representada por uma curva em forma de sino ou normal, onde o erro se distribui simetricamente com picos na média e caindo rapidamente para mais e menos infinito. Em 1809, Gauss usou esta função para analisar dados astronômicos.

A função erro é definida por:

\begin{equation} erf(x)=\frac{2}{\sqrt{\pi}}\int_0^x e^{-t^2}dt \end{equation}

A função erro complementar é definida por:

\begin{equation} erfc(x)=1-erf(x)=\frac{2}{\sqrt{\pi}}\int_x^\infty e^{-t^2}dt \end{equation}

A função erro imaginário introduzida por Gauss é definida por:

\begin{equation} erfi(z)=-ierf(iz)=-\frac{2i}{\sqrt{\pi}}\int_0^{iz}e^{-t^2}dt \end{equation}

Fazemos uma mudança de variável:

$$t=iy \Rightarrow dt=i\ dy$$

Vejam que na integral dada em $(3)$, $t$ assume valores de $0$ a $iz$ e com a mudança de variável, teremos que:

\begin{cases} \text{se } & t=0 & \Rightarrow & y=0\\ \text{se } & t=iz & \Rightarrow & iz=iy & \Rightarrow z=y \end{cases}

Obtemos então:

\begin{equation} erfi(z)=-\frac{2i}{\sqrt{\pi}}\int_0^z e^{-(iy)^2}i\ dy \end{equation}

Como $i^2=-1$, temos:

\begin{equation} erfi(z)=\frac{2}{\sqrt{\pi}}\int_0^z e^{y^2}dy \end{equation}

Podemos expressar $(5)$ em termos de $t$:

\begin{equation} erfi(z)=\frac{2}{\sqrt{\pi}}\int _0^z e^{t^2}dt \end{equation}

Proposição:

Vamos provar que:

\begin{equation*} \lim_{x \rightarrow +\infty}erf(x)=1 \end{equation*}

Demonstração:

Seja $I$ o limite acima. Assim:

$$ I=\lim_{x \rightarrow +\infty} erf(x) $$

$$ I=\lim_{x \rightarrow +\infty}\frac{2}{\sqrt{\pi}}\int _0^xe^{-u^2}du $$

\begin{equation} I=\frac{2}{\sqrt{\pi}}\int _0^{+\infty}e^{-u^2}du \end{equation}

e

$$ I=\lim_{x \rightarrow +\infty} erf(x) $$

$$ I=\lim_{x \rightarrow +\infty}\frac{2}{\sqrt{\pi}}\int _0^xe^{-v^2}dv $$

\begin{equation} I=\frac{2}{\sqrt{\pi}}\int _0^{+\infty}e^{-v^2}dv \end{equation}

Multiplicando as expressões $(7)$ e $(8)$, obtemos:

$$ I^2=\frac{4}{\pi}\cdot \int_0^{+\infty}e^{-u^2}du \cdot \int_0^{+\infty} e^{-v^2}dv $$

\begin{equation} I^2=\frac{4}{\pi}\cdot \int_0^{+\infty}\int_0^{+\infty}e^{-u^2-v^2}du\ dv \end{equation}

Observe que a região de integração é o $1^\circ$ quadrante. Para resolver a integral dupla $(9)$, De Moivre teve a ideia de usar coordenadas polares, fazendo $u=r\ \cos(\theta)$ e $v=r\ \text{sen}(\theta)$. O elemento de área $du\ dv$ está relacionado com o elemento $dr\ d\theta$ através da expressão:

\begin{equation*} du\ dv=\left | J \right |dr\ d \theta \end{equation*}

onde $J$ é o Jacobiano da transformação dado por:

\begin{equation*} J=\frac{\partial (u,v)}{\partial (r, \theta)} \\ \ \\ J=\begin{vmatrix} \displaystyle \frac{\partial u}{\partial r} & \displaystyle \frac{\partial u}{\partial \theta}\\ \\ \displaystyle \frac{\partial v}{\partial r} & \displaystyle \frac{\partial v}{\partial \theta} \end{vmatrix} = \begin{vmatrix} \cos(\theta) & -r\ \text{sen}(\theta)\\ \\ \text{sen}(\theta) & r\ \cos(\theta) \end{vmatrix} \\ \ \\ J = r\cos^2(\theta)+r\ \text{sen}^2(\theta)=r \end{equation*}

Assim, $du\ dv=r\ dr\ d\theta$, de modo que:

\begin{equation*} I^2= \frac{4}{\pi}\int_0^{\pi/2}\int_0^{+\infty}e^{-r^2}\cdot r\ dr\ d\theta \\ \ \\ I^2= \frac{4}{\pi}\int_0^{\pi/2}d\theta \cdot \lim_{p \rightarrow +\infty}\int_0^pe^{-r^2}\cdot r\ dr \\ \ \\ I^2=\frac{4}{\pi}\cdot \frac{\pi}{2}\cdot \lim_{p \rightarrow +\infty}\left[-\frac{1}{2}e^{-r^2}\right]_0^p \\ \ \\ I^2=2\cdot \lim_{p \rightarrow +\infty} \left[-\frac{1}{2}e^{-p^2}+\frac{1}{2}\right]\Rightarrow \\ \ \\ I^2=1 \\ \ \\ I=1 \end{equation*}

Observação: Os limites de integração foram $0\leq\theta\leq\pi/2$ e $0\leq r <+\infty$, pois desta forma, preenchemos o $1^\circ$ quadrante, que é a região de integração das variáveis antigas.

Corolário:

\begin{equation*} \lim_{x \rightarrow -\infty}erf(x)=1 \end{equation*}

Demonstração:

\begin{equation*} \lim_{x \rightarrow -\infty}erf(x)=\lim_{y \rightarrow +\infty}erf(-y)=\lim_{y \rightarrow +\infty}erf(y)=1 \end{equation*}

onde $y=-x$. Se $x \rightarrow -\infty$, então $y\rightarrow +\infty$.

No penúltimo passo, usamos o fato de que a função $erf(x)$ é par.

Exemplo 1:

A função gama denotada por $\Gamma(x)$ é definida por:

\begin{equation*} \Gamma(x)=\int_0^{+\infty}t^{x-1}e^{-t}dt\text{ , }x>0 \end{equation*}

Mostre que: $\Gamma (1/2)=\sqrt{\pi}$.

Façamos $t=u^2$ na integral acima. Note que $dt=2udu$ e que os limites de integração são os mesmos. Assim:

\begin{equation*} \Gamma(1/2)=\int_0^{+\infty}(u^2)^{\frac{1}{2}-1}e^{-u^2}\cdot 2u\ du \\ \ \\ \Gamma(1/2)=2\cdot \int_0^{+\infty}u^{-1}\cdot e^{-u^2}\cdot u\ du \\ \ \\ \Gamma(1/2)=2\cdot \int_0^{+\infty}e^{-u^2}du \\ \ \\ \Gamma(1/2)=\sqrt{\pi}\cdot \frac{2}{\sqrt{\pi}}\int_0^{+\infty}e^{-u^2}du \\ \ \\ \Gamma(1/2)=\sqrt{\pi} \cdot \lim_{p \rightarrow +\infty}\frac{2}{\sqrt{\pi}}\int_0^pe^{-u^2}du \\ \ \\ \Gamma(1/2)=\sqrt{\pi} \cdot \lim_{p \rightarrow +\infty}erf(p)=\sqrt{\pi}\cdot 1=\sqrt{\pi} \end{equation*}

Exemplo 2:

Vamos calcular a integral $\displaystyle \int x^2 e^{x^2}dx$. Temos que:

\begin{equation*} \int x^2e^{x^2}dx=\int x \cdot xe^{x^2}dx \end{equation*}

Aplicamos o método de integração por partes, de modo que:

\begin{equation*} u=x \Rightarrow \frac{du}{dx}=1 \Rightarrow du=dx \end{equation*}

e

\begin{equation*} dv=xe^{x^2}dx \Rightarrow v=\frac{e^{x^2}}{2} \end{equation*}

Observação: Para integrarmos $\displaystyle dv=xe^{x^2}dx$, aplicamos o método de integração por substituição. Seja $u=x^2$ e $du=2x\ dx$. Assim, se $dv=xe^{x^2}$ então:

\begin{equation*} v=\int xe^{x^2}dx \\ \ \\ v=\int \frac{xe^u}{2x}du \\ \ \\ v=\frac{1}{2}\int e^udu \end{equation*}

Como a integral de $e^u=e^u$, temos:

\begin{equation*} v=\frac{1}{2}e^u+C \\ \ \\ v=\frac{e^u}{2}+C \\ \ \\ v=\frac{e^{x^2}}{2}+C \end{equation*}

Retomando as ideias, temos pela integração por partes que:

\begin{equation*} \int udv=uv-\int vdu \end{equation*}

Assim:

$$ \int x^2e^{x^2}dx=x\cdot \frac{e^{x^2}}{2}-\int \frac{e^{x^2}}{2}\cdot dx+C $$

$$ \int x^2e^{x^2}dx=\frac{xe^{x^2}}{2}-\frac{1}{2}\int e^{x^2} dx+C $$

$$ \int x^2e^{x^2}dx=\frac{xe^{x^2}}{2}-\frac{1}{2}\cdot \frac{\sqrt{\pi}}{2}\cdot \frac{2}{\sqrt{\pi}}\int e^{x^2} dx+C $$

$$ \int x^2e^{x^2}dx=\frac{xe^{x^2}}{2}-\frac{\sqrt{\pi}}{4} erfi(x)+C $$

Autores:

• Kleber Kilhian
• Paulo Sérgio Costa Lino

Veja mais:

• Método de Integração por Partes
• Método de Integração por Substituição
• Sobre as Funções Gama e Beta