bab4

R Markdown

This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.

When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:

summary(cars)

##      speed           dist       
##  Min.   : 4.0   Min.   :  2.00  
##  1st Qu.:12.0   1st Qu.: 26.00  
##  Median :15.0   Median : 36.00  
##  Mean   :15.4   Mean   : 42.98  
##  3rd Qu.:19.0   3rd Qu.: 56.00  
##  Max.   :25.0   Max.   :120.00

Including Plots

You can also embed plots, for example:

Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot. HORIZONTAL

library(mosaicCalc)

## Loading required package: mosaic

## Registered S3 method overwritten by 'mosaic':
##   method                           from   
##   fortify.SpatialPolygonsDataFrame ggplot2

## 
## 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.

## 
## Attaching package: 'mosaic'

## The following objects are masked from 'package:dplyr':
## 
##     count, do, tally

## 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: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

## Loading required package: mosaicCore

## 
## 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

g <- makeFun(sin(x^5)*cos(sin(x^1 - 2 )-x^2) - x + 5 ~ x)
slice_plot(g(x) ~ x, domain(x = -6 : 5)) %>%
  gf_hline(yintercept  = 0, color = "red")%>%
  gf_vline(xintercept = 0, color = "blue")

VERTIKAL

slice_plot(g(x) ~ x, domain(x = 1:4)) %>%
  gf_hline(yintercept = 0, color = "red")%>%
  gf_vline(xintercept = 0, color = "blue")

findZeros( sin(x) - 0.35 ~ x, xlim=range(-10,10) )

##          x
## 1 -12.2088
## 2  -9.7824
## 3  -5.9256
## 4  -3.4992
## 5   0.3576
## 6   2.7840
## 7   6.6408
## 8   9.0672

g <- makeFun(4 + exp(k*t) - 2^(b*t) ~ b, k=0.00050, t=1)
findZeros( g(b) ~ b , xlim=range(-1000, 1000) )

##        b
## 1 2.3221

findZeros( 3*exp(-t/5)*sin(pi*t) ~ t, xlim=range(7 ,10))

##   t
## 1 6
## 2 7
## 3 8
## 4 9

findZeros( sin(cos(x^2) - x) -x - 0.52~ x, xlim=range(-1,4))

##        x
## 1 0.1981

g <- makeFun(4 + exp(k*t) - 2^(b*t) ~ b, k=0.00035, t=1)
findZeros( g(b) ~ b , xlim=range(-1000, 1000) )

##       b
## 1 2.322

findZeros( 7*x^4 -2*x^3 - 4*x^2 - 3*x - 9 ~ x, xlim=c(-9,10))

##         x
## 1 -1.0328
## 2  1.3928

findZeros( 7*x^4 -2*x^3 - 4*x^2 - 3*x - 10 ~ x, xlim=c(-10,10))

##         x
## 1 -1.0628
## 2  1.4123

findZeros( g(b, t=2) ~ b, xlim=range(-1000,1000) )

##        b
## 1 1.1611