Реферат: Спектры элементарных возбуждений в двупериодических одномерных системах

// Перегрузка /

Complex Complex::operator/(Complex &fp1)

{

double k, i, j;

k = fp1.real * fp1.real + fp1.image * fp1.image;

i = (real * fp1.real + image * fp1.image) / k;

j = (fp1.real * image - real * fp1.image) / k;

fp1.real = i;

fp1.image = j;

return fp1;

}

int main()

{

Complex Gkp;

double m;

int i,j;

for(i=0;i<N;i++)

for(j=0;j<N;j++)

{

Gkp.real=0;

Gkp.image=0;

Gkp.real=Gkp.real+1/(double)N*exp(-1/a*sqrt(pow(i*a-j*b,2)+d*d))*cos(6.28*i-6.28*j);

Gkp.image=Gkp.image-1/(double)N*exp(-1/a*sqrt(pow(i*a-j*b,2)+d*d))*sin(6.28*i-6.28*j);

}

Gkp.real=pow(Gkp.absolute(),2);

cout<<"Gkp"<<" "<<Gkp.real<<"\n";

getchar();

}