W[j]=U[j]+V[j];
int j;
double U[DIMENSION], V[DIMENSION];
#define DIMENSION 3
double W[DIMENSION];
for (j=0 ; j<DIMENSION ; ++j)