IVANA MARIA FERNANDES MACHADO

MATEMÁTICA APLICADA: O USO DAS EQUAÇÕES

DIFERENCIAIS ORDINÁRIAS EM MODELOS
MATEMÁTICOSDESISTEMAS FÍSICOS E BIO-QUÍMICOS.

ANÁPOLIS, FEVEREIRO
2016
 IVANA MARIA FERNANDES MACHADO

MATEMÁTICA APLICADA: O USO DAS EQUAÇÕES

DIFERENCIAIS ORDINÁRIAS NA MODELAGEM
MATEMÁTICADE SISTEMAS FÍSICOS E BIO-QUÍMICOS.

Trabalho de Conclusão do Curso de Licenciatura em

Química apresentado à Coordenação de

Química do Instituto Federal de Educação,

Ciência e Tecnologia de Goiás.

Orientador: Prof. Msc. Thársis Silva Souza.

Co-Orientador: Prof. Dr. Júlio César Saavedra Vasquez.

ANÁPOLIS, FEVEREIRO
2016
 FICHA CATALOGRÁFICA

Machado, Ivana Maria Fernandes

M149 Matemática aplicada: o uso das equações diferenciais ordinárias em
modelos matemáticos de sistemas físicos e bio-químicos. / Ivana Maria
Fernandes Machado. – – Anápolis: IFG, 2016.
58 p. : il.
Inclui CD- Rom.

Orientador: Prof. Me. Thársis Silva Souza

Trabalho de Conclusão do Curso de Licenciatura Plena em Química,

Instituto Federal de Goiás, Campus Anápolis, 2016.

1. Matemática aplicada. 2. Sistemas químicos. 3. Equações

diferenciais ordinárias - EDO.
I. Título

CDD 540.7

Código 003.2016
Ficha catalográfica elaborada pelo bibliotecário Helio Lino Delfino,
CRB-1/3031.
Biblioteca Clarice Lispector, Campus Anápolis
Instituto Federal de Educação, Ciência e Tecnologia de Goiás
 DEDICATÓRIA

Ao meu querido Pai que me apoiou, amou e deu coragem em todos os momentos.
A quem amo sem que palavras possam descrever.

Á minha melhor amiga, minha Mãe, que me inspira todos os dias a ser e fazer o
meu melhor.

Á meus queridos irmãos, de sangue ou coração, que foram meus companheiros,

conselheiros e amigos. Que me deram motivos para me esforçar cada dia mais.

Ao Pedro pela amizade que transcende os laços de sangue, e que também ama a
Química. Ao meu pequeno Josias por me fazer uma pessoa melhor com seu belo sorriso e à
Aline minha querida companheira de quarto.

Aos professores que me acompanharam em cada etapa da vida, e àqueles que me

ensinaram mesmo sem o dever formal.

Ao meu querido Lucas, que com carinho e paciência me acompanha. E a quem

amo de coração.
 AGRADECIMENTOS

Agradeço a meu Deus e pai Jeová pela saúde, cuidado, oportunidade, sobretudo
pelo amor.

Agradeço a minha família que se esforçou em fazer dessa uma caminhada feliz e
de sucesso.

Agradeço ao Instituto Federal de Goiás pela oportunidade, a todo seu corpo

docente que acompanhou com dedicação essa minha caminhada e aos técnicos
administrativos que com presteza e agilidade tanto colaboraram.

Agradeço o auxilio financeiro concedido pela Coordenação de Aperfeiçoamento

de Nível Superior- CAPES, pelo projeto PIBID, que apoiou essa conclusão do curso. E aos
meus queridos orientadores, supervisores e colegas.

Agradeço em especial ao meu orientador pela paciência, dedicação, atenção e

confiança desde as matérias que com ele cursei até a conclusão desse trabalho.

Agradeço a Giovanna S. Spagnol, que me fez desacelerar e aproveitar cada

pedacinho do curso de uma maneira muito mais feliz, e a quem tenho um imenso carinho.

Agradeço ao companheirismo e apoio de queridos colegas.

 - Se enxerguei mais longe [...] é porque me
apoiei em ombros de gigantes. Sir. Isaac Newton
(1642- 1727) (BOYER, 1974).
 RESUMO

A modelagem de muitos problemas do mundo real, utilizando as equações

diferenciais ordinárias (EDO), é uma ferramenta importante que tem grande potencial e pode
descrever inúmeros fenômenos. Seguindo este contexto, o presente trabalho objeta abordar
aplicações de modelos matemáticos que usam as EDO em sistemas químicos em algumas
condições especificas. A metodologia usada envolveu estudo de literatura padrão na área,
tanto de EDO como de sistemas químicos, com o objetivo de tornar clara a relação
matemática-Química do sistema, fomentando material útil para um ensino menos superficial
de cálculo em cursos de Química. Os resultados indicam que o meio matemático utilizado é
de grande emprego em sistemas químicos na determinação de fatores complexos.

Palavras-chave: Matemática Aplicada, Sistemas Químicos, Equações diferenciais Ordinárias

(E.D.O.).
 Sumário
1. INTRODUÇÃO. ......................................................................................................... 09

2. REVISÃO BIBLIOGRÁFICA. . ............................................................................... 11

2.1. Aplicação Matemática nos Cursos de Licenciatura Em Química................................. 11

2.2. Introdução as Equações Diferenciais Ordinárias. ........................................................ 12
2.3. Equações Diferenciais de Primeira Ordem. .................................................................. 15
2.4. Equações Diferenciais de Segunda Ordem. . ................................................................ 18
2.5. Equações Diferenciais de Ordens mais Altas. ............................................................. 24
2.6. Introdução a Análise Qualitativa das Equações Diferencias. ...................................... 26

3. APLICAÇÕES DE EDO EM FENÔMENOS QUÍMICOS, FÍSICOS E

BIOLOGICOS. ....................................................................................................................... 30

3.1. APLICAÇÕES EDO DE PRIMEIRA ORDEM. ......................................................... 30

3.1.1. Exemplo-1 Concentrações em uma mistura. ................................................................ 30

3.1.2. Exemplo-2 Concentração de resíduos químicos em um açude..................................... 32

3.1.3. Exemplo-3 Cinética Química: Lei de velocidade. ....................................................... 34

3.1.4. Exemplo-4 Decaimento radioativo: Carbono 14. ......................................................... 37

3.2. APLICAÇÃOPARA EDO DE SEGUNDA ORDEM. ................................................ 39

3.2.1. Exemplo- Pendulo......................................................................................................... 39

3.3. APLICAÇÃOPARA EDO „QUALITATIVA‟. ........................................................... 44

3.3.1. Exemplo- Catalise enzimática. ..................................................................................... 44

3.4. CONTEXTUALIZAÇÃO E FINALIZAÇÃO ............................................................. 51

3.4.1. Aplicação de sistemas harmônicos simples no ensino da mecânica quântica pela

Equação de Schördinger. . ........................................................................................................ 51

3.4.2. A equação de Michaelis-Menten. ................................................................................ 52

CONCLUSÕES. ...................................................................................................................... 54

REFERÊNCIAS...................................................................................................................... 55
 9

1. INTRODUÇÃO

O que é mais importante ensinar a um químico, decorar uma equação pronta ou

dominar o aparato matemático que a originou? Independente das muitas discussões que traz
esse tema no âmbito educacional algo é uníssono entre os pesquisadores do ensino de
Química, a aprendizagem significativa e motivadora é um desafio. (CARVALHO, 2011;
FERNANDES, 2014; FERNANDES e SALDANHA, 2014).

Entre as muitas dificuldades no processo de ensino-aprendizagem pesquisastem

mencionado que a abstração da disciplina de cálculo se mostra um grande desafio, e maior
ainda é a sua aplicação em outras disciplinas. FERNANDES, 2014; FERNANDES e
SALDANHA, 2014).

Conforme Bassanezi (2002), as ciências naturais, como a Química, tem o

desenvolvimento de muitos conceitos respaldados pelo formalismo matemático. Para esse
mesmo autor a modelagem matemática é o processo de atacar um fenômeno real por teorias
matemáticas resultando em um modelo que o descreve. A modelagem envolvida neste
processo é caracterizada pela multidisciplinaridade, tornando tênue a distância entre a
matemática e outras ciências como a Química, a Física e a Biologia.

Em cursos superiores de Química a matemática aplicada a sistemas químicos

permite métodos mais contextualizados ao curso e por sua vez a diminuição da distância entre
a ciência central do curso, a Química, e a matemática envolvida na validação de seus
sistemas. (BASSANEZI, 2002; FERNANDES, 2014; FERNANDES e SALDANHA, 2014)

O estudo da interação entre a matéria e a energia, a ciência Química, depende de

conhecer as propriedades físico-Químicas de muitos componentes. Com base nas variações
sofridas por esses padrões físico-químicos se torna possível determinar por que tipo de
fenômenos, Físico ou Químico, e em que proporção ocorrem. (ATKINS e DE PAULA, 2012).

Essas propriedades podem ser relacionadas como equações segundo as variações,

ou taxas, em que ocorrem durante um determinado tempo. O formalismo matemático que as
descrevem são ferramentas do cálculo, essas taxas de variação são derivadas. (BOYCE e
DIPRIMA, 2002).

Quando as tendências de um sistema são regidas por relações em que os

coeficientes envolvem taxas como varáveis o aparato matemático envolvido é um pouco mais
 10

rebuscado. Em muitos casos compreender as teorias envolvendo as equações diferenciais

ordinárias, EDO, se torna imprescindível. (LUZ e CORRÊA, 2001).

Historicamente o desenvolvimento matemático esteve extremamente interessado

em desbravar essa área. Entre os nomes de peso relacionados ao início desse desenvolvimento
estão Isaac Newton (1642- 1727) e Gottfried Wilhelm Leibniz (1646-1716) durante o século
XVII. Os problemas por eles desbravados foram além da álgebra, às aplicações de Newton
generalizando as ideias de Galileu Galilei (1564-1642) vieram a ser chamadas de “princípios
do movimento de Newton” e Leibniz por sua vez formulou junto a Christian Huygens (1629-
1695) a noção de energia cinética. (BOYER, 1974).

Uma parte da matemática aplicada se interessa em compreender o comportamento

de uma equação diferencial. Segundo Boyer (1974), com tal objetivo durante os séculos XIX
e XX muitos matemáticos detiveram suas pesquisas a investigações teóricas de existência e
unicidade e estudos de equações diferenciais mais complexas. Algumas equações resistiram
ainda, e graças aos séculos da computação (sec. XX e XXI), muitos métodos complexos de
aproximação e integração numérica foram desenvolvidos. Muitos problemas ainda perduram e
fazem das equações diferenciais um tema fértil. (LUZ e CORRÊA, 2001).

O desenvolvimento da área de modelagem de sistemas com uso de equações

diferenciais ordinárias é importantíssimo para compreensão de fenômenos físicos, biológicos
e químicos. E ainda se mostra um campo produtivo para pesquisa. (MURRAY, 1993)

Esta pesquisa tem por objetivo exemplificar aplicações do uso de EDO em

sistemas químicos, que em geral matematicamente são apenas superficialmente trabalhados
no ensino superior. Os sistemas aqui propostos advêm de um estudo da literatura padrão na
área, com o objetivo de tornar clara a relação matemática-Química do sistema, fomentando
material útil ao ensino de cálculo em cursos de Química. (ABREU e LOPES-CASTILHO,
2012; PEIXOTO, 1978).

O texto abrange uma revisão teórica do tema de EDO, e aplicações fundamentadas

nas teorias Químicas pertinentes aos diferentes sistemas químicos analisados.
 11

2. REVISÃO BIBLIOGRÁFICA

Apresentam-se aqui alguns elementos básicos imprescindíveis para compreender

as demais partes deste trabalho. Contém alguns detalhamentos teóricos relevantes quanto às
características individuais das EDO com os principais teoremas a elas relacionados e às leis
naturais aplicadas às ciências Químicas.

2.1. APLICAÇÃO MATEMÁTICA NOS CURSOS DE

LICENCIATURA EM QUÍMICA

Em pesquisas recentes se observa que os cursos de Química enfrentam uma

grande barreira no progresso de seus alunos em várias disciplinas que envolvem aplicações
matemáticas. Por exemplo, a pesquisa de Fernandes e Saldanha (2014) detectou que 67,8%
dos alunos entrevistados têm dificuldade com aplicações matemáticas no curso sendo ainda
que 70,97% afirmam que Cálculo é a disciplina mais dificultosa.
Para Santos e Borges (1993, p. 4-5) isso ocorre pela defasagem no ensino de
matemática básica, o que é refletido nos resultados no ensino superior. E para Nery, Liegel e
Fernandez (2007, p.588) a forma como são abordadas as matérias de cálculo trazem ainda
mais dificuldades. O autor justifica que isso se dá devido ao caráter inter-disciplinar da
Química, que é composta com uma gama de conteúdos inter-relacionados a outras disciplinas.
Entretanto, infelizmente dificilmente os alunos compreendem de maneira significativa essa
interdependência.
O estudo da Química exige que os discentes possuam habilidades de domínio das
ferramentas matemáticas. Segundo Morais e Teixeira Junior (2014) muitas vezes os discentes
são alvo de um ensino mecânico e sobrecarregado de conteúdos informativos que eles não
conseguem relacionar de maneira significativa ao contexto de disciplinas específicas do curso.
Nesse ponto os autores concluem que esse desafio exige uma evidente
necessidade de repensar a postura do docente, as ementas das disciplinas e até mesmo os
livros didáticos adotados, além da função desempenhada pelos conteúdos. (MORAIS E
TEIXEIRA JUNIOR, p. 205, 2014).
Quanto a essa necessária modificação muitos docentes acreditam o fato de
docentes matemáticos e químicos do mesmo curso de graduação trabalham em separado
contribui em muito para essas dificuldades. Na pesquisa de Morais e Teixeira Junior (2014)
 12

fica clara a necessidade de aumentar a contextualização, e nisto a matemática aplicada é uma

grande aliada. Entretanto os docentes matemáticos do curso alegam não dominarem os
conteúdos Químicos para vencer essa barreira em contra partida os Químicos recebem alunos
que não conseguem relacionar as duas disciplinas.
Uma alternativa mencionada por um dos entrevistados seriam trabalhar em
matérias optativas e cursos de extensão de maneira contextualizada com o apoio dos docentes
das duas áreas.
Devido à grande aplicabilidade das E.D.O. nas disciplinas químicas este trabalho
traz uma série de aplicações possíveis para os conteúdos previstos em ementas de cursos de
Licenciatura em Química. Apesar de E.D.O. nem sempre ocuparem lugar nessas ementas o
contato do aluno ocorrerá, mesmo que sem a devida percepção deste, nas disciplinas Físico-
Químicas e Bioquímicas. Os exemplos mais complexos poderiam fazer parte de matérias
optativas e cursos de extensão.

2.2. INTRODUÇÃO AS EQUAÇÕES DIFERENCIAIS ORDINÁRIAS

Muitos fenômenos químicos, físicos e biológicos possuem resultados modelados

por equações e funções, por exemplo, na Física as funções de primeiro e segundo grau
modelando o espaço, velocidade e aceleração em função do tempo no Movimento
Uniformemente Variado (MVU); na Química as equações polinomiais modelando a
concentração de reagentes de uma reação. (BOYCE e DIPRIMA, 2002).

O objeto de estudo nas ciências naturais são fenômenos, que não raro, apresentam
tendências regidas por relações de variáveis de estado com parâmetros estabelecidos. Muitos
fenômenos são modelados com equações envolvendo coeficientes que correspondem
principalmente a taxas de variação e acabam por exigir um pouco mais das funções
matemáticas que os exemplos anteriores. (BOYCE e DIPRIMA, 2002).

Do mesmo modo para Bassanezi (2004), durante o processo de modelagem as

quantidades, ou variáveis de estado, podem ser obtidas como equações das variações. Quando
estas variações são instantâneas, e as oscilações do fenômeno se dão continuamente as
equações matemáticas são chamadas por Equações Diferenciais e quando são constituídas de
uma gama de pontos das variações, utilizando cálculo de médias, são definidas como
Equações de Diferenças. A compreensão do comportamento de ambas possibilita seu uso em
modelagens prevendo eventos do fenômeno modelado.
 13

Para Boyce e Diprima (2002) os chamados princípios ou leis são muitas vezes
modelados como proporções, ou relações, que envolvem a taxa segundo a ocorrência. Essas
relações, no sentido de igualdade, são equações. E a partir do conceito de derivada em sua
principal utilização, como à taxa de variação, há a possibilidade de modelar diversas situações
físicas e Químicas com um tipo especial de equações as Equações Diferenciais Ordinárias
(EDO).

As equações diferenciais podem ainda ser classificadas segundo suas variáveis.

Segundo a descrição de Boyce e Diprima (2002), essa classificação diz respeito à quantidade
de variáveis envolvidas. Quando uma única variável está relacionada à derivação se designa
como Equação Diferencial Ordinária, ao passo que o envolvimento de duas ou mais variáveis
com a derivação resulta nas chamadas Equações Diferenciais Parciais.

O mesmo autor ainda acresce a essas designações o ordenamento segundo o

número de funções desconhecidas. Na ocorrência de envolvimento da determinação de uma
única função o uso de uma única equação é suficiente, enquanto o envolvimento da
determinação de mais de uma função faz-se necessário um sistema de equações. E possível
ainda à descrição segundo o termo ordem dependendo da maior ordem de diferenciação
apresentada.

E mencionada ainda nesta mesma obra que há referência a Equações como

Lineares ou Não-Lineares designação que exprime a condição de uma função linear de
variáveis do tipo: 𝑦, 𝑦’, . . . 𝑦(𝑛) e não linear para funções de variáveis como 𝑦𝑦’.

Segundo a definição de Beltrão (2009) para o termo “cálculo”, na visão

matemática, são as ferramentas de análise das variações ocorridas em fenômenos com alguma
natureza física, ou palpável. Sendo que o desenvolvimento do cálculo esteve intimamente
relacionado ao avanço de pesquisas das chamadas ciências naturais, especialmente a Física.
Os nomes relacionados ao início do desenvolvimento do cálculo são Isaac Newton (1642-
1727) e Gottfried Wilhelm Leibniz (1646-1716). [Boyer (1974)].

Muitas vezes no estudo de equações, de forma mais básica, são propostos

problemas onde se faz necessário definir um ou mais valores reais que satisfaça a relação
entre as incógnitas da equação. Não muito longe disso no estudo de uma equação diferencial o
objetivo é definir uma função tal que ela, e suas derivadas, satisfaçam as relações da equação
original, tal como afirma Machado (1988, p. 153),
 14

De maneira geral, uma equação diferencial é uma pergunta do tipo: “Qual a função
cuja derivada satisfaz a seguinte relação?” Ou seja, uma equação diferencial é uma
equação (no sentido de igualdade envolvendo uma incógnita) onde a incógnita é
uma função, sendo que as informações disponíveis para a determinação da função
desconhecida envolvem sua derivada.

Para possibilitar a análise de uma equação diferencial é importante compreender

suas tendências. As tendências das E.D. são analises qualitativas por comparações de
possibilidades perante a proximidade ou não proximidade há um ponto de equilíbrio, ou ponto
crítico, do sistema. Em busca da exatidão se constroem sistemas mais complexos e robustos
que são mais exigentes do aparato matemático. A análise de tendências é muito útil, pois há
uma dificuldade em analisar com maior exatidão explicita sistemas de construção mais
simples, devido às aproximações necessárias. Boyce e Diprina (2002) destacam que há
“sempre há uma troca entre precisão e simplicidade” (Boyce e Diprima, 2002, p. 10). Ainda é
possível dizer que a coexistência entre estas é algo muito raro e difícil, portanto é necessário
contrabalançar as percas e a utilidade qualitativa do modelo.

Nesse sentido os campos de direções são muito úteis, pois possibilitam analisar as
tendências de uma equação diferencial de maneira bem clara. Campos de direção são meios
de analisar uma tendência sem achar uma solução determinada. Usando da definição de
derivada, encontrada em Boyce e DIPRIMA (2002):

Definição 1- Se𝑦 = 𝑦(𝑥) é uma função derivável em um intervalo aberto 𝐼 = (𝑎, 𝑏) e seja
𝑥0 ∈ 𝐼, então 𝑦’(𝑥0 )que apresenta a inclinação da reta tangente t ao gráfico de 𝑦 = 𝑦(𝑥),
no ponto ( 𝑥0 , 𝑦’(𝑥0 )).

Sendo assim quando não é obtida a solução da diferencial há como recurso usar
um diagrama auxiliador. O campo de direções da equação diferencial é um gráfico da função
g de forma que para cada ponto 𝑄 = (𝑥, 𝑦) do domínio, g apresenta uma direção de qualquer
solução y que passe pelo ponto 𝑄 = (𝑥, 𝑦). Direção essa representada por um segmento de
reta cujo coeficiente angular é o valor da função g naquele ponto 𝑄 = (𝑥, 𝑦) e cuja origem
do segmento é o ponto 𝑄 = (𝑥, 𝑦).[7]

DAS EQUAÇÕES DIFERENCIAIS ORDINÁRIAS

Como mencionado às equações diferenciais ordinárias podem ser classificadas

segundo a ordem da sua mais alta diferencial. Elas têm a forma geral:

y(n)= f(t, y, y’, y”,...... y( n-1)). (01)

 15

Reside nessa formulação o questionamento quanto às garantias de haver solução

ou mesmo a quantidade de possíveis soluções. A fórmula geral não propicia tal garantia, mas
há teoremas que fomentam tal necessidade. Em especial para a modelagem estes teoremas
tem valor imprescindível, pois permite a generalização dos princípios, ou leis, modelados.
Entretanto obter explicitamente essas soluções é extremamente raro no âmbito das E. D. O., e
depende exclusivamente das características de cada equação. (BOYCE e DIPRIMA, 2002).

O estudo de E.D.O. pode ser organizado segundo a ordem da equação desde as

mais baixas em direção às de ordens superior ou mesmo segundo as dimensões da função
solução. Ou seja, levando em conta a ordem da derivada y(n) envolvida na E.D.O. e a
dimensão da variável y, que pode estar em R ou no Rn.(BOYCE e DIPRIMA, 2002).

Porquanto segue um estudo mais aprofundado das E.D.O. segundo seus graus.

2.3. Equações Diferenciais Lineares de Primeira Ordem

A forma geral de uma E.D.O. de primeira ordem é dada por,

f x, y, y'0 . (02)

Onde a solução é qualquer função y= ϕ(t) que satisfaça essa equação em todo t
num dado intervalo. Encontrar um método geral para resolvê-las, entretanto não é tão simples.
Essas mesmas podem ser dispostas em subclasses que exigem métodos próprios.

Se a função f, respeitando a equação (02), depender linearmente da variável

dependente y pode ser chamada de equação linear de primeira ordem. Que pode ser descrita
pela equação geral;

𝑑𝑦
+ 𝑝 𝑡 𝑦 = 𝑔(𝑡) (03)
𝑑𝑡

Sendo que p e g são funções da variável independente t, quando g(t) é nula

nomeia-se como equação linear homogênea. Quando p e g são constantes pode-se utilizar
um método chamado de integração direta desde que 𝑎 ≠ 0 e 𝑦 ≠ 𝑏/𝑎(BOYCE e
DIPRIMA, 2002). Reescrevendo assim a equação (03) como;

𝑑𝑦
= −𝑎𝑦 + 𝑏 (04)
𝑑𝑡

𝑑𝑦 /𝑑𝑡
𝑏 = −𝑎 (05)
𝑦−
𝑎
 16

Após a integração;

𝑏
ln 𝑦 − = −𝑎𝑡 + 𝑐 (06)
𝑎

Obtendo a solução geral para a equação (04)

𝑦 = 𝑏/𝑎 +𝑐𝑒 −𝑎𝑡 (07)

Esse é um método muito útil. Por exemplo, a função que modela o movimento de
um objeto em queda segue uma estrutura semelhante. Mas esse método não é geral para uma
equação como (03), para essa há um método derivado de Leibniz. Esse método se baseia em
multiplicar uma função μ(t) que será o fator integrante. Ao descobri-lo se torna possível
responder a função cuja derivada satisfaz a relação da E.D.O. modelada.

O teorema fundamental a seguir, encontrado em BOYCE e DIPRIMA (2002)

denota algo importante para as E.D.O., à existência e unicidade, e se aplica especificamente as
equações lineares.

Teorema 1- (Existência e Unicidade) Sejam p e g funções contínuas definida

num intervalo aberto I a< t < b contendo o ponto t= t0, então existe uma única função
y=ϕ(t) que satisfaz a equação diferencial;
𝑑𝑦
+𝑝 𝑡 𝑦 =𝑔 𝑡 . (03)
𝑑𝑡

Para cada y(t0)= y0 onde y0 problema de valor inicial (P.V.I.).

Esse teorema diz que uma E.D.O. dada pela Eq. (03) e uma condição inicial
𝑑𝜙 (𝑡)
y(t0)=y0, existe uma única função ϕ(t) tal que + 𝑝 𝑡 𝜙(𝑡) = 𝑔(𝑡) e inicial ϕ(t0)=y0 será a
𝑑𝑡

única solução.
Um método de encontrar soluções de E.D.O. de primeira ordem, já mencionado,
consiste na multiplicação de um fator μ(t), chamado de fator integrante na equação (03).
Neste caso,
𝑑𝑦
𝑑𝑡
𝜇 t + 𝑝 𝑡 𝑦 𝜇(t) = 𝑔 𝑡 𝜇(t) (08)

Com isso busca-se uma função μ(t) tal que a primeira parte da equação (08) possa
ser [𝜇 t y]′ , ou seja,

𝜇 t + 𝑑𝜇𝑑𝑡 t 𝑦 .
𝑑𝑦
𝑑𝑡
𝜇 t + 𝑝 𝑡 𝑦 𝜇 t = [𝜇 t y]′ = 𝑑𝑦
𝑑𝑡
(09)
 17

Essa equação propiciará a equação;

𝑑𝜇 t
=𝑝 𝑡 𝜇 t , (10)
𝑑𝑡

Que por integração direta gera,

𝑝 𝑡 𝑑𝑡
𝜇 t =c∙e (11)

Onde a constante c da Eq. (10) pode ser desconsiderada por ser a mesma que
aparece em ambos os lados da igualdade da Eq (08).

Nestas condições a Eq (08) fica da forma,

[𝜇 t ∙ y′] = 𝑔 𝑡 ∙ 𝜇(t) (12)

Que pode ser facilmente integrada, e em seguida isolar a função y.A seguir um
exemplo que demonstra a utiliza este método.

Considerando o P.V.I. y′ + 2𝑦 = t; 𝑦0 = 1;

i. Multiplicação por um fator integrante μ(t),

y′ 𝜇 t + 2𝑦𝜇 t = t𝜇 t (13)

ii. Conforme a Eq. 02,

𝑑𝜇 t
=2𝜇 t (14)
𝑑𝑡

𝑑
ln 𝜇 t = 2𝑑𝑡 (15)
𝑑𝑡

iii. Por esse processo isola-se o fator integrante,

𝜇 t = 𝑒2𝑡+c (16)

iv. Substituindo-o conforme a Eq. 12, que permite isolar y,

𝑒 2𝑡 ∙ y′ = 𝑡𝑒 2𝑡 (17)

′
y = 𝑡𝑒 −2𝑡 𝑑𝑡 (18)
 18

1
y = − 4 e−2t (2t + 1) + c (19)

v. Substituindo a condição inicial, y0= 1, se obtém a constante c,

1
1 = − e−2∙0 (2 ∙ 0 + 1) + c,
4

1=1+𝑐

c = 0, (20)

Chegando assim a solução da E.D.O. dada,

1
y = − 4 e−2t (2t + 1) . (21)

Esse tópico trouxe uma revisão dos conceitos básicos e mais necessários de
E.D.O. de primeira ordem. O método do fator integrante mencionado e o teorema de
existência e unicidade são pontos importantes que permitiram compreender os aspectos mais
complexos, e gerais, das E.D.O. de ordens mais altas.

2.4. Equações Diferencias Lineares de Segunda Ordem

A forma geral de uma E.D.O. de segunda ordem é dada,

𝑑2 𝑦 𝑑𝑦
= 𝑓 𝑡, 𝑦, 𝑑𝑡 , (22)
𝑑𝑡 2

Quando f é uma função similar a equação (03) está é designada linear. Sendo este
o caso de f ser linear em y e y’ e ter que g, p e q como funções continuas da variável t e
independentes de y. Deste modo temos a equação (23);

𝑑𝑦 𝑑𝑦
𝑓 𝑡, 𝑦, 𝑑𝑡 = 𝑔 𝑡 − 𝑝 t − 𝑞 t 𝑦, (23)
𝑑𝑡

Muitas vezes também encontrada na forma,

𝑃 𝑡 𝑦" + 𝑄 t 𝑦′ + 𝑅 t 𝑦 = 𝐺 t , (22)

com um par de condições iniciais tais como:

 19

y(t0)=y0;y’(t0)=y’0, (24)

onde são dados valores numéricos a y0e y’0 , tem-se um problema de valor inicial.

Similar ao Teorema 1 de unicidade e existência pra E.D.O. de primeira ordem o

teorema a seguir relaciona o resultado teórico fundamental para P.V.I. das E.D.O. lineares de
segunda ordem.
O problema de valor inicial para uma E.D.O. de segunda ordem é tal que:
𝑑2 𝑑1
 Resolva: 𝑎2 (𝑥) 𝑑𝑥 2 + 𝑎1 (𝑥) 𝑑𝑥 1 + 𝑎0 (𝑥)𝑦 = 𝑓(𝑥)

 Satisfazendo: y(t0)=y0e y’(t0)=y’0

Sendo que o resultado é igual para equações não-homogêneas e homogêneas e
pode ser encontrado em BOYCE e DIPRIMA (2002) e a prova pode ser observada em
CODDINTON (1961).
Teorema 2- (Existência e Unicidade de soluções)Sejam p, q e g funções
contínuas definida num aberto I. Então existe uma única função y=ϕ(t) que satisfaz esse
problema, e a solução existe em todo o intervalo I;
𝑦′′ + 𝑝 𝑡 𝑦′ + 𝑞 𝑡 𝑦 = 𝑔(𝑡) (25)

para o problema de valor inicial (P.V.I.) y(t0)= y0 e y’(t0)= y’0.

Esse teorema traz três pontos essenciais:

1] O problema de valor inicial tem uma solução; ou seja, existe uma solução;
2] Esta solução é única;
3] A solução ϕ está definida em todo o intervalo I, neste todos os coeficientes são
contínuos e, ao menos, duas vezes diferenciáveis.

 Soluções de Equações Diferenciais Lineares de Segunda Ordem.

Para que a E.D.O. de segunda ordem seja dita homogênea função g(t)tem de valer
zero para todo t. Sendo porquanto esta equação muito conhecida como o termo não-
homogêneo. Para as equações lineares vale o princípio da superposição, sendo assim no caso
de se ter várias soluções a combinação linear entre elas também e uma solução. Em relação à
solução das homogêneas têm-se as seguintes definições e teoremas.
Sejam as funções A, B e C constantes tal que,
𝑎𝑦′′ + 𝑏𝑦′ + 𝑐𝑦 = 0 (26)
lhe é atribuída a termologia autônoma, uma vez que independe da variável t.
 20

Serão agora analisadas as soluções, e métodos, para essas características das

E.D.O. de segunda ordem;

 Solução Geral para as E. D. O. Lineares Homogêneas.

Teorema 3- (Princípio da Superposição):Se y1 e y2 são soluções da equação diferencial,
𝒅𝟐 𝒅𝟏
𝒂𝟐 (𝒙) 𝒅𝒙𝟐 + 𝒂𝟏 (𝒙) 𝒅𝒙𝟏 + 𝒂𝟎 (𝒙)𝒚 = 𝟎, (27)

então a combinação linear c1y1+ c2y2 é solução para quaisquer valores das constantes c1e
c2.(BOYCE E DIPRIMA, 2002).

Quando duas funções são LD e LI e o Wronskiano da seguinte maneira;

Definição 2:
(i) Duas funções φ1, φ2: (α,β) →R são linearmente dependentes (L.D) se existe
uma constante k tal que φ2(x)= k φ1(x), ∀ x∈ (α, β).
(ii) Duas funções φ1, φ2: (α,β) →R são linearmente independentes (L.I) se a
condição α1φ1(x)+ α2φ2(x)= 0∀ x∈ (α, β) a implicar que α1= α2= 0.
Definição 3:
Dadas duas funções diferenciáveis φ1, φ2: (α,β) →R, o determinante

𝝋𝟏, 𝝋𝟐
𝑾[φ𝟏, φ𝟐] = (28)
𝝋′𝟏, 𝝋′𝟐

é chamado o Wronskiano das funções φ1e φ2.

O Teorema 3 demonstra portanto que por base de apenas um par de soluções de

y’’(x)+ py’(x) +q(x)y= 0, é possível elaborar uma família de soluções duplamente infinitas.

Teorema 4: (Solução Geral)(encontrado em LUZ e CORRÊA, 2001)

Se y1 e y2 são duas soluções particulares da equação linear homogênea
y''+ p(x)y'+ q(x)y=0, (29)
num intervalo (α, β), e se num ponto x o ∈ (α, β), o Wronskiano das duas soluções é diferente
de zero, então o Wronskiano será diferente de zero em qualquer outro ponto no intervalo(α,
β) e as soluções serão linearmente independentes no intervalo.

Teorema 5.(encontrado em LUZ e CORREA, 2001) Sejam φ1e φ2:(a,b) duas soluções L.I de
y''+ p(x)y'+ q(x)y=0. Então qualquer solução φ para a equação é da forma
φ=α1 φ1+ α2 φ2 (30)
 21

com α1 e α2 constantes escolhidas convenientemente.

Observação: Pelos teoremas da existência e unicidade e o da superposição

conclui-se que toda E.D.O.. de segunda ordem homogêneas tem duas soluções linearmente
independentes (L.I.).

 Soluções para as E. D. O. Lineares Homogêneas com Coeficientes

Constantes.
Nas E.D.O. lineares homogêneas as mais simples são as com coeficientes
constantes (a,b,c) tal qual:
𝑎𝑦′′ + 𝑏𝑦′ + 𝑐𝑦 = 0 ; (26)
Uma base de soluções é um par de funções soluções y1 e y2, que são linearmente
independentes em I. Quanto às soluções há o teorema a seguir encontrado em Luz e Correia,
2001.

Teorema 6. Se y1 e y2 são soluções de 𝑎𝑦′′ + 𝑏𝑦′ + 𝑐𝑦 = 0, linearmente independentes em I

então, toda solução y = ϕ(x) é uma combinação linear de y1 e y2, isto é
y = ϕ(x) = C1y1 + C2y2. (31)
Prova: O objetivo é mostrar que a equação (31)e a forma geral de qualquer solução de y’’(x)+
py’(x) +q(x)y= 0.Seja y = ϕ(x) uma solução de y’’(x)+ py’(x) +q(x)y= 0e seja x = x0 um ponto
do intervalo I para o qual W(y1, y2)(x0)≠ 0; Suponha-se ainda que
ψ(x0) = ψ0 ; ψ (x0) = ψ0. (32)
Procure-se provar então que existem constantes C1 e C2 tais que a solução y = ϕ(x)
C1y1 + C2y2 satisfaz as mesmas condições iniciais ψ(x0) = ψ0 ; ψ (x0) = ψ0de ψ e assim,ψ ≡ ϕ
pelo Teorema 2 (T. E. U). Para que isto aconteça deve-se ter:
C1y1(x0) + C2y2(x0) = ψ0 (33)
C1y‟1(x0) + C2y‟2(x0)= ψ0 (34)
Esse sistema nas incógnitas C1 e C2 tem solução desde que
𝑦1 𝑦2
𝑑𝑦 1 𝑑𝑦 2 (x0) = W(y1, y2)(x0) ≠ 0 (35)
𝑑𝑥 𝑑𝑥

Saindo de que a E.D.O. de primeira ordem y’ (x) – y(x)= 0 tem solução 𝑦 = 𝑒 𝑎𝑡 e

supondo que uma solução similar tal que 𝑦 = 𝑒 𝑟𝑡 com r ∈R, segundo encontrado em BOYCE
e DIPRIMA (2002), sendo r o padrão a ser definido. Tem-se que 𝑦′ = 𝑟𝑒 𝑟𝑡 e 𝑦" = 𝑟 2 𝑒 𝑟𝑡 .
Substituindo em (26) obtém-se;

(𝑎𝑟 2 + 𝑏𝑟 + 𝑐)𝑒 𝑟𝑡 = 0 (36)

 22

Sendo que 𝑒 𝑟𝑡 ≠ 0, com 𝑎𝑟 2 + 𝑏𝑟 + 𝑐 = 0chamada de equação característica.

Que por ser uma equação polinomial de segunda ordem tem duas raízes que podem
dependendo do valor do discriminante ∆ = b2 − 4ac, as raízes do polinômio
característico, e o modelo de combinação linear de equações para solução, podem ser:
 Reais distintas se ∆= b2 – 4ac> 0; com solução:

y(t) = C1𝑒 𝑟1 𝑡 + C2𝑒 𝑟2 𝑡 ; (37)

 Reais iguais se ∆= b2 − 4 ac = 0; com solução:

y(t) = ϕ(x) = 𝑒 𝑟𝑡 ; (38)

 Complexas conjugadas se ∆= b2 − 4 ac < 0:

y(t) = ϕ(x) =𝑒 𝛼𝑡 𝑐1 𝑐𝑜𝑠 𝛽𝑡 + 𝑐2 𝑠𝑖𝑛 𝛽𝑡 ,onder =𝛼 ± 𝑖𝛽. (39)

Observação: Existem ainda outros casos particulares não autônomos em que é

possível obter soluções explicitas, entretanto para as aplicações à frente a parte homogênea é
sempre autônoma.

 Solução Geral da E. D. O Linear Não-Homogênea com Coeficientes

Constantes.

Partindo da definição de uma E.D.O. homogênea é possível observar que as não-

homogêneas tem um comportamento definido segundo suas partes. Um por sua parte
homogênea e outro pela parte a torna não homogênea, do coeficiente de não homogeneidade.

Quanto a isso é possível encontrar o seguinte teorema encontrado em Boyce e

DIPRIMA (2002),

Teorema 7: Se Y1 e Y2 são duas soluções da equação não-homogênea

L 𝑦 = 𝑦" + 𝑝(𝑡)𝑦′ + 𝑞(𝑡)𝑦 = 𝑔(𝑡) (40)

então sua diferença entre Y1 e Y2 é uma solução da equação homogênea associada, com g(t)=
0. Se além disso, y1 e y2 formam um conjunto de soluções para a equação associada, então
c1y1 (t) + c2y2 (t) = Y1(t) – Y2(t), (41)
Onde C1 e C2 são constantes determinadas.
 23

Então se tem que a solução geral da equação não homogênea segue o teorema a
seguir; disponível na mesma obra anteriormente mencionada;

L 𝑦 = 𝑦" + 𝑝(𝑡)𝑦′ + 𝑞(𝑡)𝑦 = 𝑔(𝑡) (40)

Teorema 8:A solução geral da referida equação pode ser escrita tal que,

y= ϕ(t)=c1y1 (t) + c2y2 (t)+ Y(t), (42)

Onde Y1 e Y2 formam um conjunto fundamental de soluções da equação

homogênea associada, c1 e c2 são constantes arbitrárias e Y é alguma solução específica da
equação não homogênea. (BOYCE e DIPRIMA, 2002)

Das proposições destes teoremas entende-se que para solucionar a equação não
homogênea,

L 𝑦 = 𝑦" + 𝑝(𝑡)𝑦′ + 𝑞(𝑡)𝑦 = 𝑔(𝑡) (40)

faz-se necessário encontrar a solução de c1y1 (t) + c2y2que é chamada de solução

complementar e denotada por yc(t). Também se necessita encontrar uma única solução Y(t) da
equação não-homogênea, que é por sua vez chamada de solução particular. E por ultimo se
deve somar estes resultados.

As proposições para encontrar a solução complementar são as mesmas para uma

equação homogênea, como as do tópico anterior. Dois métodos são necessários compreender
para encontrar a solução particular, são estes o método dos coeficientes indeterminados e o
método da variação dos parâmetros. (BOYCE e DIPRIMA, 2002).

O método dos coeficientes indeterminados consiste na proposição de uma

hipótese inicial sobre a forma da solução particular, mas com coeficientes não especificados.
A hipótese é substituída na equação,

L 𝑦 = 𝑦" + 𝑝(𝑡)𝑦′ + 𝑞(𝑡)𝑦 = 𝑔(𝑡) (40)

e se busca determinar os coeficientes. Se houver sucesso tem-se uma solução da equação (40),
e se denota que se pode usar essa solução particular. Do contrario se tem de formular uma
nova hipótese, pois o insucesso não demonstra a inexistência de solução para e equação.
(BOYCE e DIPRIMA, 2002).
 24

Este método é extremamente útil devido a sua simplicidade, entretanto é limitado.

Depende de a equação homogênea associada ter termos constantes e o termo não-homogêneo
vai pertencer a uma classe restrita de funções, geralmente polinômios, funções exponenciais,
senos e co-senos. (BOYCE e DIPRIMA, 2002).

Já o método de variação dos parâmetros, que é derivado dos estudos de calculo

de variações de Lagrange (1736-1813), complementa bem o método anterior. Sua maior
vantagem é a generalidade, entretanto muitas das integrais necessárias envolvendo o
coeficiente não-homogêneo trazem grandes dificuldades. Quanto a esse método se tem
enunciado o seguinte teorema, encontrado em Boyce e Diprima, 2002.

Teorema 9- Se as funções p, q e g são contínuas em um intervalo aberto I e suas funções y 1e

y2 são soluções linearmente independentes da equação homogênea,

𝑦" + 𝑝(𝑡)𝑦′ + 𝑞(𝑡)𝑦 = 0 (29)

associada à equação não- homogênea

𝑦" + 𝑝(𝑡)𝑦′ + 𝑞(𝑡)𝑦 = 𝑔(𝑡) (43)

então uma solução particular da equação não- homogênea é

𝑦2 𝑔(𝑠) 𝑦1 𝑔(𝑠)
𝑌 𝑡 = −𝑦1 𝑑(𝑠) + 𝑞𝑦2 𝑑(𝑠) (44)
𝑊(𝑦1 ,𝑦2 )(𝑠) 𝑊(𝑦1 ,𝑦2 )(𝑠)

onde t0 é qualquer ponto em I escolhido convenientemente,

A solução geral é

y=c1y1 (t) + c2y2 (t)+ Y(t), (45)

conforme o teorema 8.

Esse subtema trouxe elementos básicos imprescindíveis para a compreensão do

tema de E.D.O.. como o princípio da superposição, os métodos de solução mencionados para
encontrar as soluções da parte homogênea e os de coeficientes indeterminados e método de
variação do parâmetros para as não-homogêneas.

2.4 Equações Diferencias Lineares de Ordem Mais Altas

 25

As equações diferenciais lineares de ordem mais altas (ordem n) têm, segundo

Boyce e Diprima (2002), a metodologia e a estrutura teórica compatível diretamente com a de
segunda ordem. Todos os elementos de desenvolvimento de equações diferenciais de segunda
ordem como: princípio da superposição; métodos para encontrar as soluções da parte
homogênea; método dos coeficientes indeterminados; método de variação dos parâmetros.
Tudo vale para as equações diferenciais de ordem mais alta, variando o número de soluções
independentes de acordo com a ordem da equação.

A forma geral de uma equação desse tipo é;

𝑑𝑛 𝑦 𝑑 𝑛 −1 𝑦 𝑑𝑦
𝑝0 𝑡 + 𝑝1 𝑡 + ⋯ + 𝑝𝑛 −1 𝑡 + 𝑝𝑛 (𝑡)𝑦 = 𝑔(𝑡) (46)
𝑑𝑡 𝑑𝑡 𝑛 −1 𝑑𝑡

É mais difícil perceber diretamente eventos naturais que podem ser expressos em
ordens superiores a dois devido, principalmente a dificuldade geométrica de “enxergar” as
derivadas de ordem maior que dois. Há eventos naturais que são expressos em ordens
superiores a dois e os estudar permite uma melhor compreensão das tendências, pois quando
expressos em E.D.O. mais simples, aqueles com soluções explicitas, dificilmente
possibilitariam.

Como mencionado anteriormente, em casos que não é obtida a solução da

diferencial é possível recorrer a um diagrama auxiliador. Cada um dos segmento dereta tem
em seu coeficiente angular o valor da função naquele ponto. Um gráfico como esse quando
tem uma malha relativamente fina pode permitir uma boa visão da tendência de
comportamento da função, uma visão qualitativa.

2.5 Introdução a Análise Qualitativa das Equações Diferencias

Equações diferenciais ordinárias, como as anteriores, são funções envolvendo

derivadas com uma ou mais variáveis. Com isso todas as características aplicadas às derivadas
são aplicadas as E.D.O..

Se𝑦 = 𝑦(𝑥) é uma função Y:ℜ→ ℜ, seja t0 ∈ ℜno domínio da

função 𝑦. Diretamente da derivada como taxa de variaçãose tem que:

𝑑𝑦
a. Se (t0 ) > 0, y é crescente em t0.
𝑑t
 26

𝑑𝑦
b. Se (t0 ) < 0, y é decrescente em t0.
𝑑t
𝑑𝑦
c. Se (t0 ) = 0, y não é decrescente nem é crescente em t0.
𝑑t
𝑑𝑦
Geometricamente, pensando em t0 como a taxa de variação, ou seja, como o
𝑑t

coeficiente angular da reta que passa no ponto 𝑑𝑦(t , t0 )e é tangente a y em t0 se tem que;

Figura 1: Gráfico- Da função y(t), e retas que a tangencia nos pontos 𝑑𝑦(t , t 0 )

Com isso dada uma equação diferencial qualquer pode-se buscar as características
da variável onde a função cresce, decresce ou está em equilíbrio. Como no exemplo a seguir:

𝑑𝑦
= 4 − 2y (47)
𝑑t

Analisando o sinal dessa equação se tem;

𝑑𝑦
a. Se y> 2, então < 0 a função y decresce.
𝑑t
𝑑𝑦
b. Se y< 2, então > 0 a função y cresce.
𝑑t

𝑑𝑦
Se y= 2, então = 0 , sendo assim y= 2, e está é uma solução de equilíbrio.
𝑑t

Geometricamente se tem;
 27

Figura 2: Gráfico-Diagrama de direções e soluções da função y(t).

Isso denota que toda solução em que 𝑦(0) < 2 tende a crescer para atingir o valor
2 e toda solução começando em 𝑦(0) > 2tende a decrescer para atingir o valor 2. Se 𝑦 0 =
2 se tem a solução trivial constante ou solução de equilíbrio.

Muitos fenômenos ocorrem com a existência de muitos fatores interagindo entre

si. Esses casos, não raro, são modelados com E.D.O.. com varias variáveis.

Formando um sistema de equações, como mencionado no capítulo 2 e seção 2.1.

Observando um sistema de equações de duas variáveis e duas equações, para

análise das taxas de crescimento e descimento de soluções, tal que;

𝑑𝑥
=𝑥−𝑦 (48)
𝑑𝑡

𝑑𝑦
= 𝑥+𝑦 (49)
𝑑𝑡

Para compreender seu comportamento se deve analisar as variações no ambiente

do plano (x, y) como varáveis.

𝑑𝑥
1º caso: Analisando ;
𝑑𝑡

𝑑𝑥
=0 →𝑥−𝑦 =0→𝑥 = 𝑦 (50)
𝑑𝑡

𝑑𝑥
>0 →𝑥−𝑦 >0→𝑥 >𝑦 (51)
𝑑𝑡

𝑑𝑥
<0 →𝑥−𝑦 <0→𝑥 <𝑦 (52)
𝑑𝑡

e geometricamente;
 28

Figura 3:Diagrama de Direção dx/dt..

𝑑𝑦
2º caso: Analisando 𝑑𝑡
;

𝑑𝑦
= 0 → 𝑥 + 𝑦 = 0 → 𝑦 = −𝑥 (53)
𝑑𝑡

𝑑𝑦
> 0 → 𝑥 + 𝑦 > 0 → 𝑦 > −𝑥 (54)
𝑑𝑡

𝑑𝑦
< 0 → 𝑥 + 𝑦 < 0 → 𝑦 < −𝑥 (55)
𝑑𝑡

e geometricamente;

Figura 4:Diagrama de Direção dy/dt.

Unindo-os se tem geometricamente o “campo de vetor”;

 29

Figura 5 :Diagrama de Direção.

É esperado então que a solução y(t) ou x(t) faça trajetórias circulares, da forma da
figura abaixo.

Figura 6: Gráfico- soluções em trajetórias circulares.

De fato não se atinge a clareza das soluções x(t) e y(t) por esta analise de direções,
mas há uma noção de como deve se comportar cada uma das soluções pelo comportamento
qualitativo.

O campo da análise qualitativa das equações diferenciais é muito vasto, mas aqui
se limita a exemplificação dessas analise que demonstram o caminho para o exemplo de
catálise enzimática.
 30

3. APLICAÇÕES DE E.D.O.. EM FENÔMENOS QUÍMICOS, FÍSICOS E

BIOLOGICOS.

A modelagem, segundo já definido, exige a formulação de um modelo

matemático que descreva um evento da melhor forma possível com o uso das ferramentas
matemáticas disponíveis. Ao trabalho de análise que se origina num modelo matemático em
direção ao evento que este descreve define-se como „Aplicação Matemática‟. Dizendo então
que “uma situação do mundo real que pode ser atacada por meio da matemática é chamada
uma aplicação matemática” (APPLICATIONS, 2002, p. 5).

Muitos são os exemplos de aplicação matemática nas áreas físico-Químicas e

bioquímicas. Compreender satisfatoriamente estes fenômenos e suas variáveis exige dominar
a modelagem envolvida. Os exemplos a seguir são a reunião de alguns modelos matemáticos
já conhecidos, que usam as E.D.O. em sua formulação,e algumas aplicações físico-Químicas e
bioquímicas possíveis.

3.1. APLICAÇÕES DE E.D.O. DE PRIMEIRA ORDEM

Neste tópico incluem exemplos de aplicações das E.D.O.de primeira ordem e

algumas considerações analíticas.

3.2.1 Exemplo 01- Concentração em umaMisturas1

Um tanque contém 100 galões (cerca de 455 litros) de água e 50 onças (cerca de
1,42 kg) de sal. Água contendo uma concentração de sal de ¼ (1 + ½ sen t) oz/gal entra no
tanque a uma taxa de 2 galões por minuto e a mistura no tanque sai á mesma taxa. Denotando
por x a massa de sal (em onças) na vasilha, e o tempo (em minutos) por t.

𝑑𝑥
A razão (instantânea) da variação de x com relação a t, , é igual a razão da
𝑑𝑡

variação da quantidade de sal que entra na vasilha com relação ao tempo, menos a razão da
variação da quantidade de sal que sai.

𝑑𝑥 𝑥
= 2 ¼ 1 + ½ sen 𝑡 − 2 (56)
𝑑𝑡 100

𝑑𝑥 1 1 𝑥
= + 4 sen 𝑡 − (57)
𝑑𝑡 2 50

1
Esse exemplo hipotético segue a linha de exemplos semelhantes aos propostos em Boyce e DIPRIMA, 2002.
 31

Essa é uma equação diferencial ordinária de primeira ordem. Inclui uma variável
diferenciada, x, e uma variável de diferenciação, t.

Pelo método do fator integrante, μ(t), é possível encontrar uma quantidade de sal
para qualquer instante t, encontrando assim o fator integrante:

𝑑𝑥 𝑥
∙𝜇 t + 2 𝜇 t = 2 ¼ 1 + ½ sen 𝑡 𝜇 t (58)
𝑑𝑡 100

onde conforme a Eq 10,

𝑑𝜇 t 𝑥
= 2 𝜇 t , (59)
𝑑𝑡 100

então,

𝜇 t = 𝑒𝑡/50. (60)

Substituindo-o conforme a Eq. 12, que permite isolar x,

𝑑𝑥
𝑒 𝑡/50 ∙ 𝑑𝑡 = 2 ¼ 1 + ½ sen 𝑡 𝑒 𝑡/50 , (61)

𝑑𝑥
𝑒 𝑡/50 ∙ 𝑑𝑡 = 2 ¼ 1 + ½ sen 𝑡 𝑒 𝑡/50 𝑑𝑡 ,

625 25
𝑥 = 25 − 2501 cos 𝑡 + 5002 sin t + 𝑐 ∙ 𝑒 −𝑡/50 . (62)

Substituindo a condição inicial, x0= 50, obtém-se a constante c,

625 25
50 = 25 − 2501 cos 0 + 5002 sin 0 + 𝑐 ∙ 𝑒 −0/50 ,(63)

63150
c = . (64)
2501

Encontrando assim a relação,

63150 625 25
𝑥= 𝑒 −𝑡/50 + 25 − 2501 cos 𝑡 + 5002 sin t. (65)
2501

É possível observar que esta equação reflete características do ponto de vista

físico, que são descritos pelo modelo, em t(0)tem-se que a massa de sal, x, é igual a 50 oz.
 32

A Figura-7 é um gráfico de solução da E.D.O.. que demonstra o

comportamento limite da solução. Observa-se que há oscilação em torno de um nível
25 2501
constante de 25 oz com uma amplitude de oscilação de ≅ 0,24995.
5002

Figura 7:Gráfico- Variação da massa de sal pelo tempo.

As concentrações em misturas, como neste exemplo, que envolvem a

solubilização de um dado composto sólido, um sal, em um solvente, água, foram mensuradas
considerando uma propriedade física constante, no caso o volume. Nesse exemplo, ao
trabalhar como uma E.D.O.. linear se despreza possíveis interações intermoleculares que o
possam variar.

3.2.2 Exemplo 02- Concentração de Resíduos Químicos em um Açude

Do mesmo modo uma hipotética solução de água com resíduos químicos que
contamina um açude de água potável foi trabalhada em condição ideal, ou seja, sem levar em
conta volumes parciais ou outras interações intermoleculares. Esse exemplo hipotético segue
a linha de exemplos semelhantes aos propostos em Boyce e Diprima, 2002.

Tomou-se assim um dado açude contendo inicialmente 10 milhões de galões de

água potável. O açude recebe um fluxo indesejável de resíduos químicos contaminantes de
origem industrial a uma taxa de 5 milhões de galões por ano e a mistura sai a uma mesma
taxa. A concentração γ(t) de resíduos químicos na água que está entrando varia
periodicamente com o tempo de acordo com a fórmula γ(t) = 2+ sen 2t g/ 32al.

Supondo que estes resíduos permanecem inertes no meio e que suas interações
intermoleculares não influenciem no volume total. Sendo assim visto o fluxo de entrada e
 33

saída de água ser iguais o volume total do açude permanece constante em 10 7 galões. Usando
a notação t para o tempo, medido em anos, e a quantidade de resíduos por ϕ(t), medido em
gramas, tem-se,

𝑑𝜙
= 𝑡𝑎𝑥𝑎 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 − 𝑡𝑎𝑥𝑎 𝑑𝑒 𝑠𝑎í𝑑𝑎 (66)
𝑑𝑡

Sendo que “taxa de entrada” e “taxa de saída” referem-se respectivamente as taxas

pelas quais ocorrem o acréscimo e decréscimo de resíduos no açude.

Taxa de entrada = [(5 . 106 )𝑔𝑎𝑙/𝑎𝑛𝑜(2 + 𝑠𝑒𝑛 2𝑡)𝑔/𝑔𝑎𝑙] (67)

Taxa de saída = [(5 . 106 )𝑔𝑎𝑙/𝑎𝑛𝑜[𝜙(𝑡)/107 ] 𝑔/𝑔𝑎𝑙 (68)

= [ 𝜙(𝑡)/2] 𝑔𝑎𝑙/𝑎𝑛𝑜

Tem-se então a E.D.O.. que modela o sistema em g/ano,

𝑑𝜙
= [(5 . 106 )(2 + 𝑠𝑒𝑛 2𝑡)] − [𝜙(𝑡)/2] 𝑔𝑎𝑙/𝑎𝑛𝑜 (69)
𝑑𝑡

Se a variável dependente for trabalhada como:

𝜙 𝑡 = 𝜑. 106 (70)

Então,

𝑑𝜑
+ 1/2 𝜑(𝑡) = 10 + 5 𝑠𝑒𝑛 2𝑡 (71)
𝑑𝑡

Inicialmente não havia contaminantes,sendo assim a condição inicial q(0)=0.

O fator integrante, µ 𝑡 , que será encontrado é𝑒 𝑡/2 . A equação resultante com o

uso deste é:

40 10
𝜑 𝑡 = 20 − 17 𝑐𝑜𝑠2𝑡 + 17 𝑠𝑒𝑛 2𝑡 + 𝑐𝑒 −𝑡/2 . (72)

A condição inicial que fomenta encontrar a constante c= -300/17. A solução do

P.V.I. portanto é:

𝑡
40 10 300
𝜑 𝑡 = 20 − 𝑐𝑜𝑠2𝑡 + 𝑠𝑒𝑛 2𝑡 − 𝑒 −2 . (73)
17 17 17
 34

A Figura 8- traz o gráfico da solução do P.V.I., onde em valores de t

pequenos a parte exponencial é importante o que diminui com o crescimento de t, mostra
ainda que a partir de certa etapa a solução consiste numa oscilação resultante das parcelas
sen2t e cós 2t presentes.

Figura 8: Gráfico- Variação da quantidade de resíduos pelo tempo.

3.2.3 Exemplo 4- Cinética Química.

Um campo em que as E.D.O.. estão muito presentes são os estudos físico-

químicos de cinética Química.Cinética Químicaé o estudo das velocidades da reação.Segundo
Atkins e de Paula, 2012,a velocidade total da reação pode se dada pelo somatório do produto
das concentrações molares pelas velocidades parciais de cada componente.Considerando os
reagentes como taxa de saída e os produtos como a taxa de entrada, e sabendo que estes se
anulam no ponto de equilíbrio,momento em que as reações direta (reagentes → produtos) e
inversa (produtos → reagentes) coexistem com velocidades iguais. Tal que para uma reação,

𝑣1 𝐴1 + 𝑣2 𝐴2 + ⋯ → 𝑣𝑛 𝐴𝑛 + 𝑣𝑛+1 𝐴𝑛+1 + ⋯ (74)

pode ser reescrita tal que é;

𝑣𝑖 𝐴𝑖 = 0 (75)
 35

A cinética pode ser estudada conforme o número de estágios e fases de uma

reação. Um campo importante no estudo da cinética Química é aquele em que se modela uma
equação em termos da concentração das espécies reagentes e das constantes de velocidade
chamada de lei de velocidade da reação.(MOORE, 1976).

Uma Lei de Velocidade ou lei cinéticade reação é uma equação em que há a

velocidade da reação é uma variável dependente das concentrações das espécies envolvidas
pela equação global Química para reação num dado instante. Em muitos casos é possível
obter soluções analíticas, ou experimentais, conhecidas como lei de velocidade integrada
que são muito úteis. (MOORE, 1976).

O acompanhamento analítico exige um termostato para manter a temperatura

constante e um bom cronômetro. A maior dificuldade, entretanto reside em acompanhar a
concentração do reagente ou do produto. Moore, 1976, menciona que segundo as
características de cada reação constantes físicas podem ser analisadas, como volume, pressão
e índice de absorção ou refração. (ATKINS E DE PAULA, 2012).

Os estudos de L. Wilhermy em 1850 demonstraram que a velocidade do

decréscimo da concentração de um reagente [A] com o tempo é proporcional ao reagente não
modificado. (MOORE, 1976).

Para as chamadas reações de „primeira ordem‟, uma reação do tipo A→B +

C,uma reação desse tipo é a da decomposição do pentóxido de dinitrogênio,

[2 𝑁2 𝑂5 (𝑔) → 4 𝑁𝑂2 (𝑔) + 𝑂2 (𝑔) ]. (76)

a forma integrada da lei de velocidade,

𝑑 𝐴
= −𝑘[𝐴] (77)
𝑑𝑡

é,

𝐴 = 𝐴 0 𝑒 −kt . (78)

onde k é a constante de reação que depende apenas da temperatura e na equação representa a

constante de proporcionalidade, e [A] é a concentração molar de uma espécie da reação.
(MOORE, 1976).
 36

A justificativa para as leis integradas baseia-se no fato de que asleis de velocidade

são E.D.O.., e neste caso tem estrutura semelhante a Eq. (04).

A integração assim se dá,

𝑑[𝐴]
= −k[A] (79)
𝑑𝑡

A equação pode ser integrada diretamente, pois k é independente de t. Para o

instante inicial t=0, a concentração de [A] é 𝐴 0 . Serão considerados como limites de
integração,

𝐴 𝑑𝐴 𝑡
𝐴0 𝐴
= −k 0
𝑑𝑡

𝐴
𝑙𝑛
𝑒 𝐴 0 = 𝑒 −kt

𝐴 = 𝐴 0 𝑒 −kt (80)

A Figura-9 mostra um gráfico da tendência de decaimento do reagente numa

reação de primeira ordem. Onde a concentração do reagente diminui exponencialmente com o
tempo numa velocidade determinada pela constante k, sendo que quanto maior k mais rápido
o decaimento.

Figura 9: Gráfico-Decaimento do reagente segundo a variação de k.

 37

3.2.4 Exemplo 4- Decaimento radioativo

A datação por carbono radioativo é uma ferramenta muito importante para a

arqueologia e foi desenvolvida pelo químico americano Willard F. Libby (1908-1980). Essa
datação se baseia em que algumas madeiras plantas e restos animais contêm quantidades
residuais de carbono-14, um isótopo radioativo do carbono. Esse isótopo é acumulado durante
a vida do individuo e passa a decair, ou desintegrar-se, na sua morte. (BELTRÃO, 2009).

A desintegração (variação) desse átomo que ocorre por unidade de tempo é

proporcional à quantidade de átomos radioativos presentes em cada instante. Assim, se R =
R(t) representa o número de átomos radioativos presente em cada instante t, uma equação
matemática, que pode representar o fenômeno, é dado pela Equação diferencial (81)

𝑑𝑅 (𝑡)
= −𝑘𝑅 𝑡 . (81)
𝑑𝑡

onde k, chamado de constante dedesintegração, é uma constante que na equação apresenta-se

como coeficiente de proporcionalidade. (BELTRÃO, 2009).

A meia-vida do carbono é de aproximadamente 5730 anos, por ser tão longa

permite a medição de resíduos datados de milhares de anos. Mesmo pequenas medidas de
carbono identificadas permitem a descoberta proporcional da quantidade original. Ou seja, se
for determinada uma quantidade R de C14 numa amostra, mesmo pequenas, torna-se possível
identificar a quantidade inicial, R0,pela proporcionalidade R/ R0.(JAVARONI, 2009).

Observando que a Eq. (81) apresenta estrutura semelhante a Eq. (04)

𝑑𝑅 (t)
= −kR(t) (82)
𝑑𝑡

A integração pode ser feita diretamente, pois k é independente de t. Para o instante

inicial t =0, a quantidade de isótopos C14 é R0, e segundo a meia-vida depois de 5730 anos
esse número cai pela metade, ou seja, em t=5730 obtém-se R(5730)= 𝑅0 /2. Considerando estas
relações a constante k é identificada, tal que;

𝑑𝑅(𝑡)
= −k𝑑𝑡
𝑅(𝑡)

𝑑
ln 𝑅(𝑡) 𝑑𝑡 = −k 𝑑𝑡
𝑅(𝑡)
 38

𝑅 𝑡 = 𝑐𝑒 −𝑘𝑡 . (83)

Tomando que R(5730)= 𝑅0 /2, então;

𝑅0
= 𝑅0 𝑒 −𝑘5730
2

𝑘 = 𝑅0 ∙ 0,00012097 𝑎𝑛𝑜 −1 . (84)

Sendo assim a E.D.O.. tem como solução a função,

𝑅(t) = 𝑒 −0,00012097∙t . (85)

A Figura 10 demonstra a do tendência de decaimento isótopo segundo o tempo,

sendo que quanto maior o R0mais acentuado é o decaimento.

Figura 10:Gráfico-Decaimento do isótopo C14 segundo a variação de R 0 .

 39

3.3 APLICAÇÕES DE E.D.O.. DE SEGUNDA ORDEM

3.3.1 Exemplo- Pendulo.

Pendulo pode ser caracterizado comoum oscilador harmônico. Um oscilador

harmônico é o modelo matemático que descreve o comportamento de movimento retilíneo de
uma partícula sujeita a uma força que adireciona periodicamente em direção à origem, tal que
𝐹 = −𝑘𝑟 2 . A magnitude é descrita segundo essa força por um múltiplo k (constante positiva)
da distância até origem.(ABREU E LOPEZ-CASTILHO, 2012)

Seja o pendulo conforme a figuraabaixo onde → denota a força peso tal que→=
P P

m massa . a (aceleração); →= α. V velocidade ; a resistência; e → tração de modulo T.

R T

Figura 11: Esquema- Pendulo de massa m.

A partícula, conforme a Figura-11 tem massa m, e o comprimento da corda l, onde

a corda tem massa desprezível.

Tomando-se em coordenadas polares, pelas propriedades do triangulo retângulo se

tem que;

𝑥= 𝑙𝑠𝑖𝑛𝜃
(86)
𝑦= − 𝑙𝑐𝑜𝑠𝜃
 40

A partir da derivação das coordenadas acima em relação ao tempo, t, a velocidade

é explicitada por;

𝑣 = 𝑥 ′ , 𝑦 ′ , 𝑥 ′ = 𝑙𝑐𝑜𝑠𝜃′ 𝑦 ′ = 𝑙𝑠𝑖𝑛𝜃′ (87)

A segunda derivada fomenta a aceleração, tal que;

𝑎 = 𝑥 " , 𝑦 " , 𝑥 " = −𝑙𝑠𝑖𝑛 𝜃 ′ 2

+ 𝑙𝑐𝑜𝑠𝜃 ′ . 𝜃" e 𝑦 " = 𝑙𝑐𝑜𝑠(𝜃 ′ )2 + 𝑙𝑠𝑖𝑛𝜃. 𝜃 (88)

A separação de todas as forças atuantes na partícula nos componentes x e y resulta

em;

→
𝑃𝑥 = 0; 𝑃𝑦 = −𝑚𝑔; 𝑇𝑥 = −𝑇𝑠𝑖𝑛𝜃; 𝑇𝑦 = 𝑇𝑐𝑜𝑠𝜃; 𝑅𝑥 = 𝛼𝑥 ′ ; 𝑅𝑦 = 𝛼𝑦′ (89)

Pela segunda lei de Newton a força resultante (Fr = m.a), é igual à soma das forças
aplicadas num ponto.

Aplicando este conceito em cada direção de x e y.

Em x:

𝑚. 𝑥" = −𝑇𝑠𝑖𝑛𝜃 + 𝛼𝑥 ′ ;

−𝑚.𝑥"+𝛼𝑥 ′
𝑇= (90)
𝑠𝑖𝑛𝜃

Em y:

𝑚. 𝑦" = −𝑇𝑐𝑜𝑠𝜃 + 𝑚𝑔 + 𝛼𝑦′;

𝑚.𝑦"+𝑚𝑔 −𝛼𝑦 ′
𝑇= (91)
𝑐𝑜𝑠𝜃

Igualando os módulos da tração “T” nas equações acima:

−𝑚 .𝑥"+𝛼𝑥 ′ 𝑚.𝑦"+𝑚𝑔 −𝛼𝑦 ′

= (92)
𝑠𝑖𝑛𝜃 𝑐𝑜𝑠𝜃

Substituindo os valores de x’, y’, x”, y”:

−𝑚. −𝑙𝑠𝑖𝑛 𝜃 ′ 2
+ 𝑙𝑐𝑜𝑠𝜃 ′ . 𝜃" + 𝛼 𝑙𝑐𝑜𝑠𝜃 ′ . 𝑐𝑜𝑠𝜃
= 𝑚. 𝑙𝑐𝑜𝑠(𝜃 ′ )2 + 𝑙𝑠𝑖𝑛𝜃. 𝜃" + 𝑚𝑔 − 𝛼 𝑙𝑠𝑖𝑛𝜃′ . 𝑠𝑖𝑛𝜃

𝛼 𝑔
𝜃" = − 𝑚 𝜃′ + 𝑙 𝑠𝑖𝑛𝜃 = 0 (93)
 41

Para θ muito próximo de zero, sinθ se aproximando de θ, trabalha-se a equação

anterior tal que;

𝛼 𝑔
𝜃" − 𝑚 𝜃′ + 𝑙 𝜃 = 0. (94)

𝑔 𝛼
Chamando 𝑤 2 = 𝑙; 𝛾 = 𝑚; a equação (82) se transforma em;

𝜃" − 𝛾 𝜃′ + 𝑤 2 𝜃 = 0. (95)

A equação (95) é a E.D.O.. do pendulo amortecido na vizinhança 𝜃 = 0 onde o

parâmetro 𝛾é muito pequeno.

AE.D.O.. (95) modela o comportamento de oscilação do pendulo. A resolução

depende primariamente de caracterizar o tipo de E.D.O., nesse caso está é homogênea e pode
ser resolvida pelo método explanado comosoluçõespara as E. D. O. lineares homogêneas
com coeficientes constantes.

Resolvendo a equação (95) se tem então;

𝜃 𝑡 = 𝑒 𝑟𝑡 , 𝜃′ 𝑡 = 𝑟𝑒 𝑟𝑡 , 𝜃" 𝑡 = 𝑟 2 𝑒 𝑟𝑡 (96)

𝑟 2 𝑒 𝑟𝑡 − 𝛾 𝑟𝑒 𝑟𝑡 + 𝑤 2 𝑒 𝑟𝑡 = 0. (97)

Tem-se então por solução do polinômio de segundo grau;

2
𝛥 = −𝛾 − 4 1. 𝑤 2 (98)

− −𝛾 ± −𝛾 2 −4 𝑤 2
𝑟= (99)
2

Esse problema pode existir com três comportamentos e foram abertos a seguir;

1º caso:É o caso em que o pendulo não é amortecido,para tanto tem-se queα=0 por quanto
γ=0eΔ< 0. As raízes da equação característica são 𝑟 = ±𝑤𝑖 e a solução geral é:𝜃 𝑡 =
𝑐1 cos 𝑤𝑡 + 𝑐2 sin 𝑤𝑡 ; com as condições iniciais (CI) 𝜃 0 = 1 e 𝜃′ 0 = 0.
Observação: Esse pêndulo não amortecido tem a mesma E.D.O.. do movimento
harmônico simples, usado no movimento da mola e outros movimentos periódicos.
 42

Figura 12: Gráfico- Pendulo α=0,γ=0, w= 1 e CI(𝜃 0 =

1 e 𝜃′ 0 = 0).

2º caso: É o caso em que é considerado o pendulo sofrendo amortecimento, γ< 0,

ou seja, vai dissipando energia do sistema e este vai diminuindo a intensidade em direção a
V= 0. Neste a E.D.O.. é semelhante aos que tem movimentos “periódicos” que vão
diminuindo com movimentos senoidais.
𝛾
As raízes são complexas e conjugadas e são 𝑟 = 2 ± 𝛽𝑖, onde 𝛽 = −𝛾 2 + 4𝑤 2 ,

e a solução geral é;
𝛾
𝜃 𝑡 = 𝑒 2 𝑡 𝑐1 cos 𝛽𝑡 + 𝑐2 sin 𝑤𝛽𝑡 (100)
com as condições iniciais (CI) 𝜃 0 = 1 e 𝜃′ 0 = 0e Pendulo com γ< 0, α<0 , Δ< 0.
 𝟏
43
Figura 13: Gráfico- Pendulo com γ=− 𝟖 , w= 1, Δ< 0;CI: 𝛉 𝟎 =
𝟏 𝐞 𝛉′ 𝟎 = 𝟎.

3º caso: É o caso em que é considerado um impulso, ou excitação, sobre o

𝛾
pendulo, γ> 0. As são raízes complexas e conjugadas e são 𝑟 = 2 ± 𝛽𝑖 e a solução geral é;
𝛾
𝜃 𝑡 = 𝑒 2 𝑡 𝑐1 cos 𝛽𝑡 + 𝑐2 sin 𝑤𝛽𝑡 (101)

com as condições iniciais (CI) 𝜃 0 = 1 e 𝜃′ 0 = 0 e Pendulo com γ> 0, α>0 , Δ< 0.

1
Figura 14: Gráfico- Pendulo com γ=+ 8 e w= 3
CI(θ 0 = 1 e θ′ 0 = 0.
 44

3.4APLICAÇÕES E.D.O.. ‘QUALITATIVAS’

3.4.1Exemplo- Catalise Enzimática.

O potencial de catalisar reações Químicas é uma das bases para possibilitar a vida.
Um catalisador, segundo Atkins e de Paula (2012) p. 259, é uma substância que acelera uma
reação, mas não sofre, modificação de natureza Química, nesse processo. Diminui assim a
energia de ativação, esseparâmetrodetermina a mínima energia cinética que se necessita dos
reagentes para a formação de produtos.

Boa parte das enzimas são proteínas, e não raro são catalisadores muito mais
poderosos que os sintéticos. O reconhecimento e a discrição dessas complexas moléculas se
deu inicialmente no final dos anos 1700 por estudos dirigidos a digestão de carne por
secreções estomacais. Durante os séculos seguintes muitas descobertas nesse campo foram
imprescindíveis para o desenvolvimento da bioQuímica.

As condições biológicas relevantes possibilitam condições de fator de

hidrogênio(pH), meio aquoso e temperatura desfavoráveis à cinética de reações não
catalisadas. As enzimas, segundo Nelson e Cox (2011), fomentam um ambiente especifico e
adequado para determinada reação.

Essas moléculas são estruturadas de forma a conter um sítio ativo,um “bolsão”

em que a reação ocorre isolada,e é ali que ela recebe o substrato, o reagente, fomentando o
ambiente necessário.

O complexo enzima-substrato teve sua existência proposta por Charles-Adolphe

Wutrz em 1880. Sendo este o ponto crucial para a modelagem matemática que defina o
comportamento reacional, pois as reações de velocidade dependem do estudo deste.

A cinética enzimática é um dos métodos mais importantes para compreender um

mecanismo enzimático, e consiste em determinar as variações de velocidade frente a variação
dos parâmetros. Um fator chave é a concentração de substrato, [S], entretanto a constante
mudança deste durante o processo dificulta o estudo.

É necessário então acompanhar a velocidade inicial, V 0. A teoria geral de ação das

enzimas foi proposta por Leonor Michaelis e Maud Menten em 1913. Segundo MURRAY,
1993, eles postularam a seguinte reação total:
 45

(102)

Sendo que a enzima combina-se inicialmente, e de maneira reversível, com o

substrato, e então o complexo é rompido em uma reação mais lenta que libera o produto e a
enzima.

Trabalhando matematicamente se tem uma lei que norteia o processo. Segundo a

lei de conservação da massa a concentração dos produtos e proporcional a dos reagentes.
Então trabalhando as concentrações:

[S]= s [E]= e [SE]= c [P]= p

Os postulados denotam que as velocidades de formação e quebra do complexo são

determinadas pelas constantes k1 (formação), k-1+ k2 (quebra em reagentes e produtos
respectivamente), que remetem a observação da Lei da conservação da massa já mencionada.

Aplicando a lei se tem então as seguintes equações:

𝑑𝑠
= −k1 𝑒𝑠 + k −1 𝑐, (103)
𝑑𝑡

esta diz que a taxa de concentração do substrato é proporcional a quebra do complexo em

reagentes subtraído da formação do complexo,

𝑑𝑒
𝑑𝑡
= −k1 𝑒𝑠 + (k −1 +k2 )𝑐 , (104)

que por sua vez diz que a taxa de concentração da enzima é proporcional a quebra do
complexo em reagentes e produtos subtraído da formação do complexo,

𝑑𝑐
= k1 𝑒𝑠 − (k −1 +k2 )𝑐, (105)
𝑑𝑡

por sua vez diz que a taxa de concentração do complexo é proporcional a quebra do complexo
produtos,

𝑑𝑝
= k2 𝑐 (106)
𝑑𝑡

esta diz que a taxa de concentração do produto é igual a quebra do complexo em reagentes
subtraído da formação do complexo. Tem-se assim um sistema não-linear de equações. Nem
 46

todas são independentes, a formação dos produtos depende da equação de concentração do

complexo, trabalhando-as tem-se;

𝑡
𝑝 𝑡 = 𝑘2 0 𝑐(𝑠)𝑑𝑠, (107)
Quanto à concentração inicial se tem que;

s = s0 e = e0 c=0 p=0

Isso porque a reação inicia-se com substratos e enzimasinteragindo com

concentrações determinadas e objetivo de liberar um produto, entretanto no momento inicial a
concentração de complexo e produto conhecidas e igual a zero.

Analisando as equações (103) e (104) se observa que a concentração do complexo

e enzimas na solução é inversa, equações com módulos iguais, mas sinais opostos. Pois;

𝑑𝑐 𝑑𝑒
+ 𝑑𝑡 = 0, (108)
𝑑𝑡

Essas se anulando implica que;

𝑐(𝑡) + 𝑒(𝑡) = 𝑒0 . (109)

e isolando e(t);

𝑒(𝑡) = 𝑒0 − 𝑐(𝑡) (110)

Por substituir a eq. (109) nas (103) e (105);

𝑑𝑠
= −k1 𝑒0 𝑠 + (k1 𝑠+k−1 )𝑐, (111)
𝑑𝑡

𝑑𝑐
= k1 𝑒0 𝑠 − (k1 𝑠+k −1 +k2 )𝑐. (112)
𝑑𝑡

com condições iniciais s(0) = s0 e e(0) = 0. Como as equações acima [(111) e (112)] são dadas
apenas em termos de “s” e “c” a análise qualitativa inicial é feita apenas com essas duas
equações. As analises aqui feitas só são validas para 𝑡 ≥ 0 e as concentrações das variáveis
acima ou iguais a zero.

Segue-se assim uma analise qualitativa das direções do sistema com as equações
[(111) e (112)]no plano (c, s) para deduzir as formas das soluções, valendo apenas para
𝑐 e 𝑠 > 0.
 47

𝑑𝑠
1º- Análise de 𝑑𝑡 = −k1 𝑒0 𝑠 + k −1 𝑐;

𝑑𝑠 k1 𝑒0 𝑠
>0 ⇒𝑐> k1 𝑠 + k −1 , (113)
𝑑𝑡

𝑑𝑠 k1 𝑒0 𝑠
<0⇒𝑐< k1 𝑠 + k−1 , (114)
𝑑𝑡

𝑑𝑠 k1 𝑒0 𝑠
=0⇒𝑐= k1 𝑠 + k −1 . (115)
𝑑𝑡

k1 𝑒0 𝑠
E denomine-se 𝑐 = k1 𝑠 + k −1 ; como curva “a”,

Figura 15:Diagrama de Direção ds/dt.

𝑑𝑐
2º- Análise de 𝑑𝑡 = k1 𝑒0 𝑠 − (k −1 k1 k2 )𝑐;

𝑑𝑐 k1 𝑒0 𝑠
>0 ⇒𝑐< k1 𝑠 + k−1 , (116)
𝑑𝑡

𝑑𝑐 k1 𝑒0 𝑠
<0⇒𝑐> k1 𝑠 + k −1 , (117)
𝑑𝑡

𝑑𝑐 k1 𝑒0 𝑠
=0⇒𝑐= k1 𝑠 + k −1 . (118)
𝑑𝑡
 48

k1 𝑒0 𝑠
E denomine-se 𝑐 = k1 𝑠 + k2 + k −1 ; como curva “b”

Geometricamente se obtém as seguintes direções;

Figura 16:Diagrama de Direção dc/dt.

Unindo os gráficos anteriores se obtém;

Figura 17:Diagrama de Direções.

 49

O estudo desse gráfico permite observar que na condição inicial s(0)=s 0 e c(0)= 0
𝑑𝑠
se tem que < 0 em todo t>0, portanto s(t) é sempre decrescente e tende a zero quando t
𝑑𝑡

→∞. Quimicamente isso corresponde a dizer que o substrato tende a ser totalmente
𝑑𝑐
consumido atingindo a concentração zero. Em relação a c(t) se observa que > 0 até c(t)
𝑑𝑡

k1 𝑒0 𝑠
atingir a curva “b” no valor 𝑐 = k1 𝑠 + k2 + k −1 ; onde s depende da concentração de

s0. Apos esse valor c decresce e tende a zero quando t →∞. Quimicamente isso diz que o
complexo tende à zero depois atingir a curva b, pois o substrato passa a ser mais escasso.

Assim se formula o seguinte gráfico para a concentração do substrato e do

complexo (SE);

Figura 18: Gráfico- concentração s e c pelo tempo.

Pela equação (110), o gráfico da concentração da enzima é o oposto da do

complexo ES elevando a uma altura “e 0” da concentração inicial da enzima.
𝑑𝑝
Pensando no produto se tem pela equação (105). Em t= 0 se tem = 𝑘2 𝑐 0 =
𝑑𝑡

𝑘2 0. Mas a partir do conceito básico de analises qualitativas, capitulo 2, 𝑐 𝑡 > 0 para todo
t>0, ou seja, a concentração do produto é sempre crescente. Mas como, 𝑡 → ∞, então
𝑑𝑝
= 𝑘2 𝑐, c(t) → 0 quando 𝑡 → ∞, logo a concentração do produto tende a estabilizar em
𝑑𝑡

algum valor p0. Observam-se as concentrações na Figura 19;

 50

Figura 19: Gráfico- concentração s, p, ee cpelo tempo.

As concentrações dos componentes da reação são demonstradas conforme a

Figura 19, onde s e einiciam em sua concentração máxima, e s tende a zerar sendo totalmente
consumido. Pode-se observar também que as concentrações de e e c são inversase que após
um pico de existência de c e consumo de e, o complexo tende a zero e e tende a concentração
inicial e0.
 51

3.5 CONTEXTUALIZAÇÃO E FINALIZAÇÃO

São aqui apresentadas mais possibilidades de contextualizações dos exemplos

anteriormente trabalhados, com foco nas E.D.O., com uma maior profundidade em sistemas
químicos abordados durante a graduação. Os aqui comentados são os exemplos do pendulo e
o de velocidade de reações catalisadas por enzimas. E por fim algumas considerações finais.

Aplicação de sistemas harmônicos simples no ensino da mecânica quântica pela Equação

de Schördinger.

A importância de compreender sistemas de oscilações harmônicas para o estudo

da Química é a aplicação útil no estudo de Química quântica. Uma aplicação admirável é no
ensino do conceito de radiação.

Segundo Abreu e Lopez-Castilho, 2012, a Química quântica é introduzida nos

cursos de graduação em Química, geralmente, de forma abrupta e sem preparar um caminho
de conhecimentos matemáticos e físicos necessários. O que origina muitos problemas para o
ensino-aprendizagem do tema.

Para Peixoto, 1978, muitos problemas de natureza quântica não tem uma solução
analítica explicita. Entretanto modelos hipotéticos de aproximação de fenômenos quânticos
são muito úteis para compreender as tendências dos sistemas naturais.

Esses problemas são aplicados a Equação de Schrödinger, que pode ser

considerada uma análoga na mecânica quântica da segunda lei de Newton, da mecânica
clássica. Esta é uma equação diferencial homogênea de segunda ordem, independente do
tempo.Considerando a oscilação harmônica unidirecional pode ser descrita tal; (PEIXOTO,
1978).

𝑑2 8𝜋𝜇 1
=− [ 𝑘𝑥 2 − 𝐸] 𝜓 𝑥 = 0 (119)
𝑑𝑥 2 𝑕2 2

Sendo k a constante de força do oscilador associada à energia potencial do

sistema, x é o deslocamento apresentado pelo oscilador em relação à posição de equilíbrio e µ
é a massa reduzida do sistema.

Peixoto, 1978, sugere o uso de um modelo a partir do oscilador harmônico

simples, citado no exemplo-1 do subtema de aplicações de E.D.O.. de segunda ordem, como
um ponto de partida que permite conectar o conhecimento da mecânica clássica do aluno e a
 52

teoria de mecânica quântica. Sendo útil por permitir uma analise muito precisa dos espectros
vibracionais causados pelas vibrações moleculares.

3.5.1 A equação de Michaelis-Menten

Tomando a equação de velocidade de transformação do complexo SE (c), após o

estagio transitório de complexação;

𝑑𝑐
= k1 𝑒𝑠 − (k −1 +k2 )𝑐 = 0, (120)
𝑑𝑡

segundo estabelecido pelos estudos de L. Michaelis e M. L. Menten (1913) a constante Km

k−1 +k2
(constante de Michaelis)= k1 .

es
𝑐(𝑡) = 𝐾 , (121)
m

a velocidade pode ser determinada tal que;

es 𝑑𝑠
v = k2 𝑐 = k2 𝐾 = − , (122)
m 𝑑𝑡

Sendo assim a velocidade da reação é dada pela concentração restante de

𝑑𝑠
substrato,𝑣 = − . Entretanto essa concentração depende da concentração de [E] e
𝑑𝑡

experimentalmente só haverá o dado de [E] 0, ou e0, que é a concentração inicial da enzima, na

presença de um excesso de substrato. Pela relação anterior, que determinou c(t)+e(t)=e0,

e s
𝑐 = 𝐾 0+𝑠, (123)
m

A famosa equação de hipérbole retangular de L. Michaelis e M. L. Menten (1913)

é dada tal que;

𝑘2 e 0 s 𝑑𝑠
v = k2 𝑐 = k2 =− (124)
𝐾m +𝑠 𝑑𝑡

e portanto para a velocidade inicial da reação;

(125)
 53

Onde k2[E]0 = Vmax é a velocidade máxima da reação para uma dada concentração
inicial da enzima, [E]0. A constante de velocidade k2 é também chamada de kcat ou número de
renovação. (CARVALHO, 2010).

Espenson, 2002 p. 91,o melhor modo para obtenção dos valores de Km, Vmax e kcat
é o uso de um ajuste não linear por mínimos quadrados. Ajustes menos confiáveis
estatisticamente, mas mais simples, são os ajustes lineares.

Carvalho, 2010, menciona cinco manipulações de ajustes lineares entre eles dentre
as quais se tem as equações de Lineweaver-Burk, de Hanes e de Eadie-Hofstee.

Tabela 1: Equações de manipulações de ajustes lineares da

equação Michaelis-Menten.

Essas equações são mais confiáveis, entre si, conforme a necessidade de maior
precisão de um parâmetro ou outro.2 A sua importância reside em ser didaticamente mais
simples, e a maior rapidez que oferece na percepção desvios de um comportamento no
mecanismo do tipo Michaelis-Menten.

2
Carvalho, 2010, apresenta em seu trabalho uma discussão interessante da capacidade dessas equações em
produzir parâmetros de Michaelis-Menten confiáveis.
 54

CONCLUSÕES

Diante do exposto pode-se concluir que o ensino de cálculo apesar de

imprescindível tem sido trabalhado nos cursos de química com uma carência no âmbito das
aplicações e da contextualização. Muitos desafios contribuem para esse estado tais como as
dificuldades prévias dos discentes em matemática básica, carência de professores matemáticos
que dominem conhecimentos químicos pertencentes à grade dos cursos, bem como de
Químicos que usem o aparato matemático em suas explicações.

Neste trabalho foi possível destacar exemplos de aplicações que podem ser
desbravados de maneira interdisciplinar, bem como individualmente pelos docentes de ambas
as áreas, Matemática e Química. O que possibilitou determinar que as E.D.O., mesmo as mais
simples, têm grande emprego em sistemas químicos na determinação de fatores complexos. E
pode-se concluir que o uso aplicado, mesmo que demande mais tempo de estudo, é possível e
de grande importância para um ensino menos superficial de cálculo em cursos de Química

Da pergunta que iniciou essa discussão, “O que é mais importante ensinar a um

químico, decorar uma equação pronta ou dominar o aparato matemático que a originou?”, se
pode destacar que o “decorar” acaba realmente sendo, infelizmente, usado quando não se
domina o aparato matemático. Mas para um Químico o aparato matemático só pode ter
utilidade direta se for aplicado de maneira clara à ciência central de seu curso. Por isso
aplicações são imprescindíveis. Sendo assim ações práticas que tenham em foco à melhoria
do ensino de Matemática nos cursos de graduação em Química se mostram necessárias e
urgentes.
 55

REFERÊNCIAS

[1] ABREU, L.; LOPES-CASTILHO, A.; Sintonização de estados quânticos: um

estudo numérico do oscilador harmônico quântico. Quim. Nova, Vol. 35, No. 8,
1692-1695, 2012. Disponível em:
<http://quimicanova.sbq.org.br/imagebank/pdf/Vol35No8_1692_32-
ED11942.pdf>, acesso em: 18 jan. 2016.
[2] APPLICATIONS and Modelling in Mathematics Education. International
Reviews on Mathematical Education, v. 34, n. 5 ISSN 1615-679X. Disponível
em: <http://subs.emis.de/journals/ZDM/zdm025i1.pdf>. Acesso em: 18 nov.
2014.
[3] ATKINS, P. W., PAULA, J. Atkins, Físico-Química V. 2, Tradução: Edilson
Clemente da Silva, 8ª ed. – [Reimpr.].- Rio de Janeiro: LTC, 2012.
[4] BASSANEZI R. C e W. C. FERREIRA JR, Equações Diferenciais com
Aplicações. 1ª ed., Edit. Harbra, 1988.
[5] BASSANEZI, R. C. Ensino-Aprendizagem com Modelagem Matemática: uma
nova estratégia. São Paulo: Contexto, 2002.
[6] BELTRÃO, M. E. P. Ensino de Cálculo pela Modelagem e Aplicações- Teoria e
Prática.2009. 320 f. Tese (Doutorado em Educação Matemática)- Pontifícia
Universidade Católica de São Paulo, São Paulo, 2009.
[7] BOYCE, W. E. ; DIPRIMA, R. C. Equações Diferenciais Elementares e
Problemas de Valores de Contorno. Tradução Valéria de Magalhães Iorio. 7. ed.
Rio de Janeiro: LTC, 2002.
[8] BOYER, C. B. História da Matemática. Tradução Elza F. Gomide. São Paulo:
Editora: Edgar Blücher, Editora da Universidade de São Paulo, 1974.
[9] CARVALHO, N. M. F. I; etc e tal, Uso de equações lineares na determinação dos
parâmetros de Michaelis- Menten. Quím. Nova vol.33 No.7 São Paulo, 2010.
Disponível em: <http://dx.doi.org/10.1590/S0100-40422010000700033>; acesso
em: 18/ 01/ 2016.
[10] CARVALHO, R. B. F.; ALMEIDA, A. A. C.; DA SILVA, A. A. C. A.; LEITE,
D. B. C.;ATAIDE, M.A. C. S. Dificuldades apresentadas por alunos do ensino
superior do curso de licenciatura Química diante de situações-problema e questões
retiradas de listas tradicionais, ANNQ, 2011. Disponível
 56

em:<http://www.annq.org/congresso2011/arquivos/1300242646.pdf>. Acesso em:

14 jan. 2016.
[11] CODDINGTON, E. A. Na introduction to ordinary differential equations
(Englewood Cliffs, NJ: Prentice-Hall, 1961; New York: Dover, 1989) Ensino-
aprendizagem com modelagem matemática: uma nova estratégia. 2 ed. São Paulo:
Contexto, 2004.
[12] ESPENSON, J. H.; Chemical Kinetics and Reaction Mechanisms, 2nd ed.;
McGraw-Hill: Nova Iorque, 2002.
[13] FERNANDES, D. M. da S. Utilização de Softwares Educacionais na diminuição
da abstraçãono ensino-aprendizagem de Química em nível superior. 88 f.
Monografia (Graduação em Licenciatura e Química) – Instituto Federal do Ceará,
Quixadá,2014.
[14] FERNANDES, D. M. da S; SALDANHA, G. C. B. Dificuldades de
Aprendizagem no Nível Superior: estudo de caso com graduandosde licenciatura
em Química. 2014. Disponível
em:<http://enalic2014.com.br/anais/anexos/7704.pdf>. Acesso em: 14 jan. 2016.
[15] JAVARONI, S. L. Abordagem geométrica: possibilidades para o ensino e
aprendizagem de Introdução às Equações Diferenciais Ordinárias. 2007. 321 f.
Tese (Doutorado em Educação Matemática)- Universidade Estadual Paulista,
Instituto de Geociências e Ciências Exatas, Rio Claro- SP, 2009.
[16] LUZ, A. M. S.; CORRÊA, F. J. S. A. EQUAÇÕES DIFERENCIAIS
ORDINÁRIAS E APLICAÇÕES. Revista Virtual de Iniciação Acadêmica da
Universidade Federal Pará. Vol. 1, N o 1, março 2001. Disponível em: <
http://www.cultura.ufpa.br/rcientifica/ed_anteriores/pdf/ed_01_amsl.pdf> Acesso
em: 18 nov. 2014.
[17] MACHADO, N.J. Matemática por assunto: noções de cálculo. São Paulo:
Scipione, 1988.
[18] MICHAELIS, L.; MENTEN, M. L.; Biochem. Z. 1913.
[19] MORAES; TEIXEIRA JUNIOR REFLEXÕES SOBRE O ENSINO DE
CÁLCULO DIFERENCIAL E INTEGRAL EM CURSOS DE GRADUAÇÃO
EM QUÍMICA. DiversaPrática, v. 2, n. 1, p. 194-217 – 1° semestre 2014.
[20] MOORE, W. J. Físico-Química. Tradução da 4ª ed. Estadunidense/ Helena Li,
São Paulo: Editora Edgar Blücher, 1ª ed,1976.
[21] MURRAY, J. D. Mathematical Biology, Springer Verlag. 2ª ed., 1993, USA.
 57

[22] PEIXOTO, E. M. A. Modelos simples em mecânica quântica. Quim. Nova, Vol.

1, No.4 08/ 08/ 1978, Disponível em:
<http://quimicanova.sbq.org.br/imagebank/pdf/Vol1No4_21_v01_n4_(6).pdf>,ac
esso em: 18 jan. 2016.
[23] Nelson, D. L. and Cox, M. M. (Editor). Princípios de bioQuímica de Lehninger. 5ª
ed. Artmed, PortoAlegre.2011.
[24] NERY, A. L. P.; LIEGEL, R. M.; FERNANDEZ, C. Um olhar crítico sobre o uso
de algoritmos no Ensino de Química no Ensino Médio: a compreensão das
transformações e representações das equações químicas. Revista Electronica de
Enseñanza de las Ciencias, 6 (3), p. 587-600, 2007.
[25] SANTOS, R. M.; BORGES, H. B. Avaliação do Desempenho no Processo de
Ensino Aprendizagem de Cálculo Diferencial e Integral I (O Caso da UFC). 1993.
Disponível em:
<http://www.multimeios.ufc.br/arquivos/pc/artigos/artigoavaliacao-do-
desempenho-no-processo-de-ensinoaprendizagem.pdf>. Acesso em: Janeiro de
2016.