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
19-08-2020-JacobiSistemaNewtonaprpx
BoxPlot_Anomaly
Simple Plot
Leap Year or Not
Variance Gamme process 2
Gab1Bdicussão
LDA
此次次
Generating prime numbers in R
print name and age