Приложение. Комментарии к коду программы

Прилагаемая программа на языке С снабжена подробными комментариями. Программа реализует моделирование движения броуновской частицы и сбор статистики.

Особое внимание следует обратить на функцию langevene_solution, реализующую расчет координат и скоростей для одного шага по времени. Здесь параметр dt задает шаг по времени. В структуре P хранятся координаты и проекции скоростей частицы.

Далее следуют функции, которые собирают данные для построения графиков и осуществляют запись и чтение данных о координатах и скоростях в файл. Эти функции можно не изучать подробно.

Основной является функция main, где собственно и происходит выполнение моделирования и накопление статистики. Переменная count определяет количество проходов. Переменная eps задает точность расчетов. При запуске программы следует ввести число шагов по времени N, шаг по времени t, максимальное значения случайной силы F. Начнется процесс моделирования. При нажатии комбинации клавиш Ctrl+C, вычисления приостановятся. Далее следуйте инструкциям на экране.

Список литературы

[1]Хеерман Д.В. Методы компьютерного эксперимента в теоретической физике. М.: Наука, 1990.

[2]Гулд Х., Тобочник Я. Компьютерное моделирование в физике. М.: Мир, 1990. т.1,2.

[3]Метод молекулярной динамики в физической химии. Под ред. Ю.К. Товбина. М.: Наука, 1996.

[4]Frenkel D., Smit B. Understanding Molecular Simulation. From Algorithms to Applications. San Diego, Academic Press, 2002.

[5]Федоренко Р.П. Введение в вычислительную физику. М.: Изд-во Моск. физ.-техн. ин-та, 1994.

49

[6]Allen M.P., Tildesley D.J. Computer Simulation of Liquids, Oxford, Oxford University Press, 1990.

[7]Ландау Л.Д.. Лифшиц Е.М. Статистическая физика. Часть 1. М.: Наука, 1995.

[8]Metropolis N., Ulam S. Journal of the American Statistical Association, Vol. 44, No. 247, pp. 335-34 (1949).

[9]Metropolis N., Rosenbluth A.W., Rosenbluth M.N., Teller A.H., Teller E., Equation of state calculations for fast computing machines, J. Chem. Phys.

6, 1087 (1953).

50