Module Gsl.Poly

Polynomials

type poly = float array

Polynomial Evaluation

external eval : poly ‑> float ‑> float = "ml_gsl_poly_eval"

eval p x returns p.(0) +. p.(1) *. x +. p.(2) *. x**2 +. ... +. p.(n) *. x**n where n = Array.length p.

Quadratic Equations

type quad_sol =
| Quad_0
| Quad_2 of float * float
external solve_quadratic : a:float ‑> b:float ‑> c:float ‑> quad_sol = "ml_gsl_poly_solve_quadratic"
external complex_solve_quadratic : a:float ‑> b:float ‑> c:float ‑> Gsl_complex.complex * Gsl_complex.complex = "ml_gsl_poly_complex_solve_quadratic"

Cubic Equations

type cubic_sol =
| Cubic_0
| Cubic_1 of float
| Cubic_3 of float * float * float
external solve_cubic : a:float ‑> b:float ‑> c:float ‑> cubic_sol = "ml_gsl_poly_solve_cubic"
external complex_solve_cubic : a:float ‑> b:float ‑> c:float ‑> Gsl_complex.complex * Gsl_complex.complex * Gsl_complex.complex = "ml_gsl_poly_complex_solve_cubic"

General Polynomial Equations

val solve : poly ‑> Gsl_complex.complex_array