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

{

$result = 0-$j*($fi_np1 + $fi_nm1)+$g0*$u*$fi;

return $result;

}

function f2 ($j,$f_np1, $f_nm1, $g0, $u, $f)

{

$result = $j*($f_np1 + $f_nm1)-$g0*$u*$f;

return $result;

}

function f3 ($g0, $c0, $fi, $f, $u, $hamma, $p)

{

$result = 0-($g0/$c0)*($fi*$fi+$f*$f)-$u-$hamma*$p;

return $result;

}

function f4 ($p)

{ return $p; }

/* Тут починається тіло скрипта

Параметри, що передаються:

$j, $g0, $c0, - параметри рівнянь

$hamma - коефіцієнт тертя

$dt - крок

$T - кінцевий час, бажано задати цілим

$N - кіль-сть молекул в ланцюжці

$filo – ім’я файла, в який слід записувати результати

Всі ці параметри в даному скрипті є константами