0) Инициализировать графический режим

1. задать начальные условия

t=0; C_old=C; C_oldi=0; i=1,…n

3) найти новые концентрации

4) Вывести на экран зависимость концентрации на данном шаге по времени от координаты (начало 20; 400; масштаб. Множители по х , по у, подпрограммаMoving-curve

5) Переприсвоить новые значения концентраций старым

t=t+dt

да

Тестовый пример: D=1e-10 м²/с; с=1е-3; l=lе-10м; dx=1е-12м.

dt- определить из условия сходимости метода конечных разностей.

Задание : модифицировать алгоритм для ненулевого потока по правой границе и исследовать зависимость c(x;t) от соотношения потоков на левой и правой границе.

Рекомендована література

1. А.М. Гусак, В.Г. Гриценко, Т.В. Запорожець. Статистична фізика – основні положення та моделі.-Черкаси, ЧДУ, 1998. –288с.

2. Задачи по физике твердого тела. Ред. Г. Дж. Голдсмита. –М.; Наука, 1976. – 96с.

3. Тихонов А.Н., Самарский А.А. Уравнения математической физики.-М.; Гостехиздат, 1953.-679с.

4. Новожилов Б.В.Метод Монте-Карло.- М.:Наука,1960.-43с.

5. Хеерман Д.В. Методы комп’ютерного эксперимента в теоретической физике/под ред. С.А.Ахманова.-М.:Наука,1990.-176с. (файл «2017 Д.В.Хеереман Методы компьютерного эксперимента в теоретической физике.pdf в локальній мережі кафедри»)

6. Райченко А.И. Математическая теория диффузии в приложениях.-К.; Наукова думка, 1981. –396с.

7. Гулд Х., Тобочник Я.. Компьютерное моделирование в физике: в 2-х частях. Часть 1; Пер. с англ. М,. Мир, 1990. – 349 с., ил. (файл Х. Гулд Я.Тобочник Компьютерное моделирование в физике Часть 1.pdf в локальній мережі кафедри)