let error a b =
  let e = ref 0. in
    for i = 0 to (Array.length a)-1 do
      e := !e +. (a.(i)-.b.(i))*.(a.(i)-.b.(i))
    done ;
    sqrt !e