// math_64.h  computer dot and cross product of two vectors
double dot(long int n, double x[], double y[]);
void cross(long int n, double x[], double y[], double z[]);