cout<<"Namagnichennost vsey sistemi="<div>
sreden(n, arr);
double T, Efull=0.0, M=0.0, C=0.0, Esum1=0.0, Esum2=0.0;
double Estep[10], Mstep[10], Cstep[10];
double Msum=0.0, Csum=0.0;
double Emist=0.0, Mmist=0.0, Cmist=0.0;
Ef.open("Efull.txt");
Mf.open("M.txt");
Cf.open("C.txt");
for(T = 0.01; T<4; T+=0.1)
{
for(int step=0; step<10; step++)
{
for(long int MK=0; MK<100000; MK++)
{
double E1=0.0, E2=0.0, p=0.0, slch=0.0;
int i=rand()%n;
int j=rand()%n;
E1=e1sp(i, j, n, arr);
arr[i][j]=(-1)*arr[i][j];
E2=e1sp(i, j, n, arr);
if(E2>=E1)
{
p=exp(((-1)*(E2-E1))/T);
p*=100;
slch=rand()%101;
if(p>=slch)
{
arr[i][j]=(-1)*arr[i][j];
}
}
}
Efull=poles(n, arr);
Estep[step]=Efull;
Esum1+=Efull;
Esum2+=(Efull*Efull);
Mstep[step]=namag(n, arr);
Msum+=Mstep[step];
}
cout<<"T="<div>
Efull=poles(n, arr);
Emist=er(Estep, Esum1);
cout<<"E(T)="<div>
Ef<div>
M=namag(n, arr);
Mmist=er(Mstep, Msum);
cout<<"M(T)="<div>
Mf<div>
C=tepl(Esum1, Esum2, T);
cout<<"C(T)="<div>
Cf<div>
}
for(int i=0; idiv>
{
delete[] arr[i];
}
delete[] arr;
Ef.close();
Mf.close();
Cf.close();
}
