Part 1

Solució per diferents valors de a (demanat a la pràctica) i de c (no demanat).

library(mosaic)
## Warning: package 'mosaic' was built under R version 3.5.1
## Loading required package: dplyr
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
## Loading required package: lattice
## Loading required package: ggformula
## Warning: package 'ggformula' was built under R version 3.5.1
## Loading required package: ggplot2
## Warning: package 'ggplot2' was built under R version 3.5.1
## Loading required package: ggstance
## Warning: package 'ggstance' was built under R version 3.5.1
## 
## Attaching package: 'ggstance'
## The following objects are masked from 'package:ggplot2':
## 
##     geom_errorbarh, GeomErrorbarh
## 
## New to ggformula?  Try the tutorials: 
##  learnr::run_tutorial("introduction", package = "ggformula")
##  learnr::run_tutorial("refining", package = "ggformula")
## Loading required package: mosaicData
## Warning: package 'mosaicData' was built under R version 3.5.1
## Loading required package: Matrix
## 
## The 'mosaic' package masks several functions from core packages in order to add 
## additional features.  The original behavior of these functions should not be affected by this.
## 
## Note: If you use the Matrix package, be sure to load it BEFORE loading mosaic.
## 
## Attaching package: 'mosaic'
## The following object is masked from 'package:Matrix':
## 
##     mean
## The following object is masked from 'package:ggplot2':
## 
##     stat
## The following objects are masked from 'package:dplyr':
## 
##     count, do, tally
## The following objects are masked from 'package:stats':
## 
##     binom.test, cor, cor.test, cov, fivenum, IQR, median,
##     prop.test, quantile, sd, t.test, var
## The following objects are masked from 'package:base':
## 
##     max, mean, min, prod, range, sample, sum
library(mosaicCalc)
## Warning: package 'mosaicCalc' was built under R version 3.5.1
## Loading required package: mosaicCore
## Warning: package 'mosaicCore' was built under R version 3.5.1
## 
## Attaching package: 'mosaicCore'
## The following objects are masked from 'package:dplyr':
## 
##     count, tally
## 
## Attaching package: 'mosaicCalc'
## The following object is masked from 'package:stats':
## 
##     D
soln <- integrateODE(dy~a*y*(c-y), a=.001, c=300, y=200, tdur=50)
plotFun(soln$y(x)~x, xlim=c(0,10))

soln <- integrateODE(dy~a*y*(c-y), a=.1, c=300, y=200, tdur=50)
plotFun(soln$y(x)~x, xlim=c(0,10))

soln <- integrateODE(dy~a*y*(c-y), a=.5, c=300, y=200, tdur=50)
plotFun(soln$y(x)~x, xlim=c(0,10))

soln <- integrateODE(dy~a*y*(c-y), a=1, c=300, y=200, tdur=50)
plotFun(soln$y(x)~x, xlim=c(0,10))

soln <- integrateODE(dy~a*y*(c-y), a=.001, c=100, y=200, tdur=50)
plotFun(soln$y(x)~x, xlim=c(0,10))

soln <- integrateODE(dy~a*y*(c-y), a=.1, c=100, y=200, tdur=50)
plotFun(soln$y(x)~x, xlim=c(0,10))

soln <- integrateODE(dy~a*y*(c-y), a=.5, c=100, y=200, tdur=50)
plotFun(soln$y(x)~x, xlim=c(0,10))

soln <- integrateODE(dy~a*y*(c-y), a=1, c=100, y=200, tdur=50)
plotFun(soln$y(x)~x, xlim=c(0,10))

Part 2

Per un sol valor de c (no especificat a la pràctica).

soln <- integrateODE(dy~y*(c-y/40)-10, c=6, y=200, tdur=50)
plotFun(soln$y(x)~x, xlim=c(0,5))