1. 개발/실행 플랫폼 
  (1) OS : win7
  (2) 개발툴 : delphi7
  (3) 사용 Database : 

2. 개발 중인 프로그램 

3. 질문 내용

아주 작은값의 평균값과 표준편차 구합니다.   

var a , b , c : array[0..1000] of double  ;
   i : integer ;
   v_mean , v_Sd : extended ;
  
........
........
  for i := 0 to 1000 do
    begin
      c[i] := ( a[i] - b[i] )/ b[i] ;
    end ; 

// MeanAndStdDev( c :array of double ; var  v_mean , v_Sd : extended )

   MeanAndStdDev( c , v_mean , v_Sd ) ; <-- v_mean , v_Sd 오차가 너무 커서 사용할 수 없습니다.
                                        특히 c 가 0보다 작은값일 때 오차가 더 큽니다.(쓰레기 수준) 

임시로 지금 아래처럼 사용하고 있습니다.

  for i := 0 to 1000 do
    begin
      c[i] := ( ( a[i] - b[i] )/ b[i] ) * 1000000 ;
    end ; 

   MeanAndStdDev( c , v_mean , v_Sd ) ;

  v_mean := v_mean / 1000000 ;
  v_Sd := v_Sd / 1000000 ;

이 값도 만족할 정도는 아닙니다. 
더 좋은 해결책은 없나요?
작은 수를 다뤄본 경험이 있는 고수님 들의 조언을 바랍니다...
읽어주셔서 감사 합니다.