Динамічні властивості нелінійних локалізованих мод у лінійних молекулярних ланцюжках

$fu = $u[$m][0];//Тотожність, зручно записувати у файл

//просто змінну не масива

$fm = $f[$m][0]*$f[$m][0]+$fi[$m][0]*$fi[$m][0];

//Запис $m $t $fm $fu

fputs($xwail, "$m 0 $fm $fu\n");

fputs($nulli, "$m $fm $fu\n");

}

fclose ($nulli);

print "У початковий момент часу норма = $norm
";

//Накладаємо періодичні умови, "нульова" молекула тотожня до N-ї

$f[0][0] = $f[$N][0];

$fi[0][0] = $fi[$N][0];

$p[0][0] = $p[$N][0];

$u[0][0] = $u[$N][0];

//Так само вчинимо з молекулою N+1

$f[$N+1][0] = $f[1][0];

$fi[$N+1][0] = $fi[1][0];

$p[$N+1][0] = $p[1][0];

$u[$N+1][0] = $u[1][0];

for ($t=$dtau; $t < $tau+$dtau; $t=$t+$dtau) //Зсув за часом

{

$norm = 0;

$pot = fopen ("$t.txt", w);

for ($m = 1; $m <= $N; $m++) //Пішов обрахунок N молекул

{