File: diff.ml

package info (click to toggle)
ocamlgsl 0.3.5-3
  • links: PTS
  • area: main
  • in suites: sarge
  • size: 3,444 kB
  • ctags: 2,901
  • sloc: ml: 7,956; ansic: 6,796; makefile: 303; sh: 87; awk: 13
file content (30 lines) | stat: -rw-r--r-- 658 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
 
open Gsl_fun

let f x = 
  (* raise Exit ;*)
  x ** 1.5

let test () = 
  let gslfun = f in
  Printf.printf "f(x) = x^(3/2)\n" ;
  flush stdout ;

  begin
    let { res=result; err=abserr } = Gsl_diff.central gslfun 2.0 in
    Printf.printf "x = 2.0\n" ;
    Printf.printf "f'(x) = %.10f +/- %.5f\n" result abserr ;
    Printf.printf "exact = %.10f\n\n" (1.5 *. sqrt 2.0)
  end ;

  flush stdout ; 

  begin
    let { res=result; err=abserr } = Gsl_diff.forward gslfun 0.0 in
    Printf.printf "x = 0.0\n" ;
    Printf.printf "f'(x) = %.10f +/- %.5f\n" result abserr ;
    Printf.printf "exact = %.10f\n\n" 0.0
  end

let _ = 
  Gsl_error.init ();
  test ()