#include <S.h>
void my_sum(x, n, sum)
double *x ; long *n ; double *sum ;
{
        long i ;
        *sum = 0 ;
        for (i = 0 ; i < *n ; i++)
                if (is_na(&x[i], DOUBLE)) {
                        na_set(sum, DOUBLE) ; break ;
                } else {
                        *sum += x[i] ;
                }
}