Aqui você vê as diferenças entre duas revisões dessa página.
Ambos lados da revisão anteriorRevisão anteriorPróxima revisão | Revisão anterior | ||
exercicios:calc1 [2012/05/11 17:14] – [Método de Euler] adalardo | exercicios:calc1 [2024/01/09 18:18] (atual) – edição externa 127.0.0.1 | ||
---|---|---|---|
Linha 138: | Linha 138: | ||
A aproximação foi boa? Tente repetir o mesmo código com dt = 0.01 e 0.001 e compare. | A aproximação foi boa? Tente repetir o mesmo código com dt = 0.01 e 0.001 e compare. | ||
- | ====== Integração Numérica no R ====== | + | ====== Integração Numérica no R ====== |
- | Vamos integrar numericamente algumas equações usando o pacote deSolve e a função //ode//. | + | Não precisamos fazer todo o procedimento anterior para fazer a integração numérica no {{: |
+ | Vamos integrar numericamente algumas equações usando o pacote | ||
Antes de tudo precisa instalar e carregar o pacote. Para instalar você pode usar o menu do RGui ou pela linha de comando, digite: | Antes de tudo precisa instalar e carregar o pacote. Para instalar você pode usar o menu do RGui ou pela linha de comando, digite: | ||
Linha 217: | Linha 218: | ||
*1. $ \frac{dy}{dt} = y-y^2*f(t)$ | *1. $ \frac{dy}{dt} = y-y^2*f(t)$ | ||
- | * sendo: $f(t)= | + | * sendo: $f(t)= |
* M=15 | * M=15 | ||
*2. $\frac{dn}{dt} = r(t)*n $ | *2. $\frac{dn}{dt} = r(t)*n $ | ||
- | * sendo: $r(t)= | + | * sendo: $r(t)= 0.1 - 100 sin(2 \pi t)$ |
- | Uma solução: | ||
- | < | ||
- | ### primeiro caso | ||
- | fy3 <- function(time, | ||
- | { | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | } | ||
- | y0 = 10 | ||
- | prmt=c(M=1) | ||
- | st=seq(0.1, | ||
- | res.fy3= ode(y=y0, | ||
- | plot(res.fy3[, | ||
- | |||
- | ### segundo caso | ||
- | |||
- | fy4 <- function(time, | ||
- | { | ||
- | | ||
- | rt= 0.1-100*(sin(2*pi*time)) | ||
- | | ||
- | | ||
- | } | ||
- | y0 = c(.1) | ||
- | # | ||
- | st=seq(0.1, | ||
- | res.fy4= ode(y=y0, | ||
- | plot(res.fy4[, | ||
- | |||
- | </ |