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:17] – [Integração Numérica no R] adalardo | exercicios:calc1 [2024/01/09 18:18] (atual) – edição externa 127.0.0.1 | ||
---|---|---|---|
Linha 218: | 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[, | ||
- | |||
- | </ |