*harmonics[0]=0.0;
for(int i=1;i<numharmonics;i++){
*harmonics[i]=1.0/sqrt(i);
} *harmonics[0]=0.0;
for(int i=1;i<numharmonics;i++){
*harmonics[i]=1.0/sqrt(i);
}