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

$fi[$m][$t] = $fi[$m][$tm1] + ($dtau/6)*($w1[$m][$tm1] + 2*$w2[$m][$tm1] + 2*$w3[$m][$tm1] + $w4[$m][$tm1]);

$p[$m][$t] = $p[$m][$tm1] + ($dtau/6)*($y1[$m][$tm1] + 2*$y2[$m][$tm1] + 2*$y3[$m][$tm1] + $y4[$m][$tm1]);

$u[$m][$t] = $u[$m][$tm1] + ($dtau/6)*($z1[$m][$tm1] + 2*$z2[$m][$tm1] + 2*$z3[$m][$tm1] + $z4[$m][$tm1]);

if ($m == $N)

{

$v1[0][$tm1] = $v1[$N][$tm1];

$v2[0][$tm1] = $v2[$N][$tm1];

$v3[0][$tm1] = $v3[$N][$tm1];

$v4[0][$tm1] = $v4[$N][$tm1];

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

$w1[0][$tm1] = $w1[$N][$tm1];

$w2[0][$tm1] = $w2[$N][$tm1];

$w3[0][$tm1] = $w3[$N][$tm1];

$w4[0][$tm1] = $w4[$N][$tm1];

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

$y1[0][$tm1] = $y1[$N][$tm1];

$y2[0][$tm1] = $y2[$N][$tm1];

$y3[0][$tm1] = $y3[$N][$tm1];

$y4[0][$tm1] = $y4[$N][$tm1];

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

$z1[0][$tm1] = $z1[$N][$tm1];

$z2[0][$tm1] = $z2[$N][$tm1];

$z3[0][$tm1] = $z3[$N][$tm1];

$z4[0][$tm1] = $z4[$N][$tm1];

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