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

(вводяться з клавіатури у форму soliton.php і передаються цьому файлу)

*/

// Початкові умови

// Обрахунок параметрів рівнянь

$tau = $T;

$dtau = $dt;

//вивід параметрів на екран

print "Параметри:
N=$N, tau=$tau, dtau=$dtau, j=$j, g0=$g0, c0=$c0, hamma=$hamma
";

$xwail=fopen($filo, w);

fputs($xwail, "№ t f^2+fi^2 u(x,t)\n");

$norm = 0;

$nulli = fopen ("0.txt", w);

for ($m=1; $m <= $N; $m++)

//Задаємо початкове збудження для всіх молекул

{

if ($m == $N/2)

$f[$m][0] = 0.8;

if (($m == $N/2+1)||($m == $N/2-1))

$f[$m][0] = 0.4;

if (($m == $N/2+2)||($m == $N/2-2))

$f[$m][0] = sqrt(0.02);

$fi[$m][0] = 0;

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

$p[$m][0] = 0;

$u[$m][0] = 0;