Run Code
|
API
|
Code Wall
|
Misc
|
Feedback
|
Login
|
Theme
|
Privacy
|
Patreon
root of a real number (bisection)
# Square root of a real number - Bisection root <- function(number,tol = 0.00001){ bis <- function(f,a=0,b=number,tol){ while( (b-a)/2 > tol){ fa = f(a) fb = f(b) pm = (a+b)/2 fpm = f(pm) # Some especial cases if(f(b) == 0){ return(b) break; } if(fpm == 0){ return(pm) break; } if(fpm*fb < 0){ a <- pm }else{ b <- pm pm <- (a+b)/2 } } cat("Root: ", pm) curve(f, pm - 10, pm + 10) abline(h=0,v=0) points(pm,0, pch = 19, col = 'red') } f <- function(x){ x^2 - number } return(bis(f,0,number,tol)) } root(3) sqrt(3)
run
|
edit
|
history
|
help
0
Análise de Variâncias DIC e teste de Tukey 1
19-08-2020-GaussSeidelSistema-menos-calc
CDF- Gaussian (trapezoids method)
untitled
Regressão linear
Zac2
Call option formula for Geometric Laplace distributed stock
123
G
Implied Volatility Surface