jerome/rust-numerical-solvers
1
0
Fork 0
Code Issues 2 Pull requests Projects Releases Packages Wiki Activity

Added more mathematica reference solutions.

Browse source
This commit is contained in:
Jérôme 2023-03-26 00:19:00 +01:00
parent 35ea6fe20b
commit 31c341e808
1 changed files with 5 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
src/main.rs
View file

@ -35,7 +35,9 @@ fn main() {
let tol : f64 = 1e-10; let tol : f64 = 1e-10;
let max_iter : u32 = 100; let max_iter : u32 = 100;
let dx_num : f64 = 1e-6; let dx_num : f64 = 1e-6;
let x_mathematica = -3.26650043678562449167148755288; let x_mathematica: f64 = -3.26650043678562449167148755288;// 30 digits of precision
let x_mathematica_2: f64 = -6.27133405258685307845641527902;// 30 digits of precision
let x_mathematica_3: f64 = -9.42553801930504142668603949182;// 30 digits of precision
let x_newton = univariate_solvers::newton_solve(&(fct as fn(f64) -> f64), &(dfct as fn(f64) -> f64), x0, tol, max_iter); let x_newton = univariate_solvers::newton_solve(&(fct as fn(f64) -> f64), &(dfct as fn(f64) -> f64), x0, tol, max_iter);
let x_newton_num: f64 = univariate_solvers::newton_solve_num(&(fct as fn(f64) -> f64), x0, tol, dx_num, max_iter); let x_newton_num: f64 = univariate_solvers::newton_solve_num(&(fct as fn(f64) -> f64), x0, tol, dx_num, max_iter);
let x_halley: f64 = univariate_solvers::halley_solve(&(fct as fn(f64) -> f64), &(dfct as fn(f64) -> f64), &(ddfct as fn(f64) -> f64), x0, tol, max_iter, false).unwrap(); let x_halley: f64 = univariate_solvers::halley_solve(&(fct as fn(f64) -> f64), &(dfct as fn(f64) -> f64), &(ddfct as fn(f64) -> f64), x0, tol, max_iter, false).unwrap();
@ -44,6 +46,8 @@ fn main() {
let x_secant : f64 = univariate_solvers::secant_solve(&(fct as fn(f64) -> f64), -1.0, 1.0, tol, max_iter); let x_secant : f64 = univariate_solvers::secant_solve(&(fct as fn(f64) -> f64), -1.0, 1.0, tol, max_iter);
let x_ridder : f64 = univariate_solvers::ridder_solve(&(fct as fn(f64) -> f64), -5.0, 1.0, tol, max_iter).unwrap(); let x_ridder : f64 = univariate_solvers::ridder_solve(&(fct as fn(f64) -> f64), -5.0, 1.0, tol, max_iter).unwrap();
println!("Mathematica : x = {}\tf(x) = {}", x_mathematica, fct(x_mathematica)); println!("Mathematica : x = {}\tf(x) = {}", x_mathematica, fct(x_mathematica));
println!("Mathematica 2 : x = {}\tf(x) = {}", x_mathematica_2, fct(x_mathematica_2));
println!("Mathematica 3 : x = {}\tf(x) = {}", x_mathematica_3, fct(x_mathematica_3));
println!("Newton's method : x = {}\tf(x) = {}", x_newton, fct(x_newton)); println!("Newton's method : x = {}\tf(x) = {}", x_newton, fct(x_newton));
println!("Newton's method (num) : x = {}\tf(x) = {}", x_newton_num, fct(x_newton_num)); println!("Newton's method (num) : x = {}\tf(x) = {}", x_newton_num, fct(x_newton_num));
println!("Halley's method : x = {}\tf(x) = {}", x_halley, fct(x_halley)); println!("Halley's method : x = {}\tf(x) = {}", x_halley, fct(x_halley));