Научная статья на тему 'Потенциалы Лиенара-Вихерта для модели расширенного пространства'

Потенциалы Лиенара-Вихерта для модели расширенного пространства Текст научной статьи по специальности «Физика»

CC BY
310
26
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по физике, автор научной работы — Д Ю. Ципенюк, В А. Андреев

В расширенном (1+4)-мерном пространстве G(T] X, S) рассматривается модель, объединяющая электромагнитное и гравитационное поля. Для уравнений, описывающих эти поля, найдены точные решения, отвечающие точечным движущимся источникам. Такие решения выражаются через потенциалы Лиенара-Вихерта в пространстве G{T\X,S). Проанализированы их свойства.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Потенциалы Лиенара-Вихерта для модели расширенного пространства»

УДК 535.32

ПОТЕНЦИАЛЫ ЛИЕНАРА-ВИХЕРТА ДЛЯ МОДЕЛИ РАСШИРЕННОГО ПРОСТРАНСТВА

Д. Ю. Ципенюк, В. А. Андреев

В расширенном (1+4)-мерном пространстве С(Т~, X, Б) рассматривается модель, объединяющая электромагнитное и гравитационное поля. Для уравнений, описывающих эти поля, найдены точные решения, отвечающие точечным движущимся источникам. Такие решения выражаются через потенциалы Лиенара-Вихерта в про-странстве С(Т\Х,8). Проанализированы их свойства.

В работах [1, 2] была предложена модель, являющаяся обобщением специальной

теории относительности Эйнштейна (СТО) на 5-мерное расширенное пространство —*

Сг(Т; X, 5) с метрикой

(+,-,-,-,-)•

В рамках этой модели расширенного пространства (МРП) рассматривалась возможность объединения гравитации и электромагнетизма в одно единое поле. Различные детали, относящиеся к структуре МРП, содержатся в работах [3 - 7]. Был поставлен эксперимент по проверке предсказаний МРП [8, 9]. Для его количественного анализа необходимо знать точные решения уравнений поля, которые были сформулированы в [1]. В данной работе найдены решения, отвечающие точечным источникам.

Главное отличие МРП от СТО состоит в том, что МРП позволяет изучать процессы, при которых меняется масса частицы т. Под массой частицы т, следуя рекомендациям обзора [10], будем понимать ее массу покоя, которая является лоренцевским скаляром. В нашей модели масса частицы - компонента 5-вектора в пространстве С(Т'; X, 5), другими компонентами которого служат ее энергия и импульс. Таким образом, масса частицы

_ —*

меняется при преобразованиях в пространстве Сг(Т; X, Я). Что касается ее заряда, то

предполагается, что он остается постоянным. Мы исходим также из того, что заряд

е всегда связан с некоторой массой т и сопоставляем заряженной частице параметр

—*

те. Поскольку заряд е - скаляр в расширенном пространстве С(Г; X, 5), параметр те - компонента вектора в этом пространстве. В качестве 5-й дополнительной координаты используется та величина, которая уже существует в пространстве Минковского, а именно, интервал 5:

з2 = (сг)2-х2-у2-г2. (1)

Эта величина сохраняется при обычных преобразованиях Лоренца в пространстве Мин-—* —♦

ковского М{Т; X), но меняется при поворотах в расширенном пространстве (7(Т'; X, 5).

—♦

Таким образом, пространство Минковского М(Т; X) - это конус в расширенном про-—* —♦

странстве Сг(Т; X, 5"). 5-вектора в пространстве С(Т; X, 5) будем обозначать черточкой (г), вектора в 4-мерном евклидовом пространстве с координатами (ж, г, в) и вектора в его 3-мерных подпространствах обозначаются стрелочкой сверху и цифрой внизу, указывающей размерность пространства, к которому относится этот вектор: ^(4))-Как правило, из контекста бывает понятно, о каком 3-мерном подпространстве идет речь и его данные в номенклатуре вектора мы не указываем.

В работе [1] было показано, что единое поле, объединяющее гравитацию и электромагнетизм, порождается током, который в случае одной частицы является изотропным 5-вектором

т / • • \ ( ЯС Я™ ^ /г.\

^ = С(2)

где

/З2 = ^(з^/с2; Я — те; дй = тег?(£, х, у, г). (3)

Если имеется несколько частиц, то их общий ток есть сумма токов отдельных частиц. Дополнительные преобразования, которые имеются в расширенном пространстве Сг(Т; X, меняют величину д аналогично тому, как они меняют массу покоя т [1 - 3].

В [1] была построена расширенная система уравнений Максвелла, описывающая это поле. Как и в случае обычной электродинамики, его напряженности выражаются через компоненты 5-вектор-потенциала

(<р, А, А„) = (Аи Ах, Ау, А„, А„). (4)

Потенциал (4) и ток (2) связаны уравнениями

4я*

о Л, =--jt, (5)

с

— 47Г -

оА =--j, (6)

с

47Г .

оЛ, =--Js. 7

с

Здесь

_ д2 д2 д2 ° ~ дз2 + дх2 + ду2 + дг2 с2 дг2' (8)

Для того, чтобы получить из системы (5) - (7) обычную систему уравнений для вектор-потенциала, которая используется в электродинамике, надо учесть, что масса т - скаляр относительно преобразований Лоренца, поделить каждое из уравнений (5), (6) на эту массу т, а уравнение (7) на заряд е, а также учесть то, что теперь все величины не зависят от переменной 5. Тогда оператор (8) превращается в обычный даламбертиан. После этого система (5) - (7) распадется на систему уравнений для потенциалов обычной классической электродинамики и уравнение для гравитационного поля, которое является скаляром по отношению к преобразованиям Лоренца. Эти системы уравнений имеют вид:

oAt = -4тг/9, (9)

-» 47Г

□ А =--р, (10)

с

4я"

аА3 =--р„. (11)

с

В правых частях уравнений (9) - (11) стоят компоненты тока

ее еи \ ,,

mcj . (12)

г Уу/Г^' У/Т

Уравнения (9), (10) описывают электромагнитное поле движущегося точечного заряда, а уравнение (11) - скалярный потенциал гравитационного поля. Таким образом, в случае, когда масса постоянна, т.е. нет зависимости от пятой координаты з, гравитационное и электромагнитное поля можно рассматривать по отдельности, разделяя уравнения (9), (10) и (11). Если же такая зависимость имеет место, то их надо рассматривать как единое поле.

При рассмотрении любой системы уравнений большое значение имеет нахождение ее точных решений. Прежде всего, сам факт существования такого нетривиального решения указывает на совместность такой системы, ее содержательность и внутреннюю непротиворечивость. Изучая структуру точного решения, можно обнаружить такие свой ства и закономерности, которые не видны, если исходить из самой системы. Особенно важны точные решения уравнений, описывающих какое-либо физическое явление. Га кое решение фактически дает нам модель процесса, динамика которого известна во всех деталях. Для системы уравнений Максвелла важнейшую роль играет решение, называемое потенциалами Лиенара-Вихерта. Это решение дает явный вид электромагнитного поля, которое создает заряд, движущийся произвольным, но вполне определенным образом. Закон движения заряда является внешним условием, определяющим структуру решения. В данном разделе построим потенциалы Лиенара-Вихерта для системы (5)

- (7). По своей формальной структуре они имеют вид запаздывающих потенциалов.

Идея построения запаздывающих потенциалов основывается на том, что величина напряженности поля, создаваемого движущимся зарядом в точке (х, у, г) в момент времени t, задается положением заряда в некоторый предыдущий момент времени I' < В обычной электродинамике расстояние между источником поля и точкой наблюдения

- это расстояние в 3-мерном евклидовом пространстве. В нашей модели двум точкам

—»

в расширенном пространстве Сг(Т; X, Б) сопоставляется 4-мерный пространственный интервал

являющийся частью полного интервала (1).

Обсудим вопрос о том, какие параметры следует подставлять в формулу (13) при вычислении потенциалов Лиенара-Вихерта.

Пространственные координаты (х, у, г), (х',у',г') определяют соответственно положение точки, в которой измеряют напряженность поля в момент времени и положение заряда в момент времени

В работах [1 - 3] было показано, что пятую координату й можно связать с показателем преломления п. В соответствии с этим расстояние по 5-й координате (5 — л')2 определим следующим образом. Если частица движется в пустом пространстве и при распространении ее поля ему не приходится пересекать участки с п > 1, то в (13) (5 — з') = 0. Если же такие участки встречаются, т.е. отрезок, соединяющий заряд с точкой наблюдения, лежит в 4-мерном евклидовом пространстве, то пространственный

Д?4) = {х- х'У + (у - уГ + (г- г'У + (* - в')

2

(13)

б

интервал имеет вид (13). При этом величина компоненты (5 — 5') = Rs зависит не от значения показателя преломления в этих точках, а от его значения во всех точках пути, их соединяющего. Отметим, что этот путь может и не быть прямолинейным отрезком, его геометрия зависит от распределения коэффициента преломления n(x,y,z). Пусть точка наблюдения поля и его источник находятся друг от друга на 3-мерном расстоянии

R\3) = (х - х'У + (у - у')2 + (z - г')2 (14)

и свет проходит это расстояние за время Ai, тогда по определению интервала

(S - S'Y = c2(AtУ - Rl3) = R]. (15)

В этом выражении s и 5' не являются значениями параметра s в точке наблюдения поля и в точке - его источнике. Поскольку разность (s — s') характеризует весь отрезок, соединяющий источник с точкой наблюдения, мы будем использовать в дальнейшем для него обозначение Rs. В этих обозначениях скорость заряда по переменной s есть us, а его ускорение ûs. При этом скорость vs распространения света вдоль координаты s определяется с помощью соотношения

Для определения 4-мерного расстояния Л(4) надо величину (15) подставить в формулу (13). 4-мерное расстояние i?(4), скорость света с, 3-мерное расстояние i?(3) и скорость света в среде v связаны друг с другом соотношением

V С

или

7?(3)П = /2(4),

где п - показатель преломления среды.

Связь момента испускания поля t' с моментом его наблюдения t определяется с помощью соотношения

с

которое имеет тот же вид, что и аналогичное соотношение в обычной электродинамике, с той лишь разницей, что вместо 3-мерного расстояния в нем стоит 4-мерное. Выражение (18) является уравнением на время испускания выражающим его через время наблюдения I. В пустоте оно имеет только одно решение [11], однако известно, что в среде это не так и оно может иметь несколько решений [12], далее обсудим этот вопрос подробнее.

Перейдем к непосредственному построению потенциалов Лиенара-Вихерта. Для этого, как и в обычном случае, сначала запишем потенциалы неподвижной заряженной частицы, а затем перепишем их в форме, инвариантной относительно преобразований в расширенном пространстве (7(Т; X, 5).

Итак, пусть имеется точечная частица с зарядом е и массой т. Ей отвечает 5-вектор тока

3 — ет6(г — г') ( —====, —===. с ) .

V ; Уу/Т^Р1 \/1~Р2 )

(19)

Подставляя эти значения компонент тока в правые части уравнений (5) - (7), получим, что при нулевой скорости частицы (и = 0) потенциалы имеют вид:

ет „ , ет .__.

Здесь 72(4) ~ расстояние в 4-мерном пространстве от источника поля до точки наблюдения. Используя (18), выражения для потенциалов можно записать в виде

А1 = ц>= —-Л = 0, А, = —-(21)

Г ф-*') ' сЦ-Г) К '

—*

Поскольку потенциал - это вектор в расширенном пространстве 0(Г; X, 5"), выражение (21) легко обобщается на случай движущихся частиц:

А - етй (22)

л-(й,йу (22)

Здесь - 5-вектор скорости заряженной частицы; И(^) - радиус-вектор, соединяющий источник поля с точкой наблюдения; (Я, й) - их скалярное произведение. Значения обоих этих векторов берутся в момент времени 1'\

Я = {с(г - о, (г-г'Ю), Д.},

й = {«,(<'), «„(О, «.(О}- (23)

Оба эти вектора изотропны. По своей структуре радиус-вектор Й похож на вектор энергии-импульса-массы фотона [1 - 3], а вектор скорости й на вектор энергии-импульса-массы массивной частицы. В пустом пространстве

R = {c(t - О, (r-r'(i')), 0}, й = с} •

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(24)

Расписывая покомпонентно скалярное произведение (Л, и), выражение (22) для потенциалов Лиенара-Вихерта можно записать в виде

Аг = <р(1) = -^-. (25)

сД(4)(*') - (Д(4)(*'), *{4){*))

—* —

Здесь величины Д(4), гГ(4) обозначают пространственные части 5-векторов Я, й

Д(4) = ((« - Л (У - У% (* - Л (» - О)» (26)

«(4) = («а;,^,^,«,) = а их скалярное произведение (Я(4), гГ(4)) имеет вид

(Л(4), гГ(4)) = (Л(4), Д(4)) = (27)

-((х - XV + (у - у')у + (г - + (в -Мы будем пользоваться обозначением

£(4) = Я<4)(0 - "(4)(0)- (28)

С его помощью потенциалы (22) можно записать в виде

л'(<»" 2Ш' <29>

А т - етЦд а т - е77Шу

А - еГПиг А П\ - етиз

~ СС{4)(ГУ М*> ~ СС{4)(ГУ Параметры < и ¿' связаны друг с другом соотношением (18).

Отметим, что найденные таким способом потенциалы Лиенара-Вихерта в рас-

—+

ширенном пространстве (5(Т; X, 5) как частный случай дают обычные потенциалы

Лиенара-Вихерта в сплошной среде [12, 13]. Их связь устанавливается с помощью соотношения (17).

Потенциалы (29) являются решениями системы (5) - (7) вида

/ = /((#( 4), "(4)))-

Это не сферически-симметричное решение, поэтому на него не распространяется утверждение теоремы о том, что п-мерный лапласиан имеет решения вида

/ = г_п+2, при п > 2, и / = 1пгприп = 2.

Тот факт, что для неподвижных частиц решение (29) дает обычный кулоновскин потенциал, непосредственно следует из формулы (21). Действительно, решение (29) мы получили, переходя от неподвижной частицы к движущейся, поэтому останавливая частицу, мы вернемся от общего решения (29) к частному решению (21). Это решение можно переписать в форме (20), которая в случае пустого пространства (и — с, Яз = Я4) дает обычный кулоновский потенциал.

Найдем теперь напряженности полей, соответствующих потенциалам Лиенара Вихерта (29). Для этого воспользуемся формулой

= ¿,¿ = 0,1,2,3,4, (30)

ОХ к ОХ{

связывающей напряженности с потенциалами.

—*

Сначала вычислим напряженность электрического поля Е. Для этого следует воспользоваться формулой

Е = -—-(31) тс оъ т

аналогичной той, по которой вычисляется электромагнитное поле в обычной электродинамике [11]. Отличие состоит в том, что в формуле (31) правая часть поделена на массу покоя источника тока т. В обычной электродинамике такое деление осуществляется еще на стадии построения 4-вектора тока из 4-вектора энергии-импульса, в нашей модели

мы относим эту процедуру на момент вычисления напряженностей электромагнитно-

—♦

го поля. Отметим, что так определенное поле Е является 2-компонентным тензором

—+

в пространстве Минковского М(Т;Х), но это не тензор в расширенном пространстве Сг(Т; X, 51). В этом пространстве тензорами являются величины (30), а напряженность

(31) преобразуется при вращениях в пространстве С(Т; X, 5") по гораздо более сложному закону.

Подставляя сюда выражения для потенциалов (29), получим

С2С\4}

((с2 - af4)) (Я(3) - -Т?(3)Л(4)) + i Д(3) - (Я(4),И(4))-

е / („о i?(4)

- ^г ^ ^(4) - -f (%), ад j «о) j. (32)

Магнитное поле имеет вид:

Я = —rot А. (33)

т

Подставляя сюда значения потенциалов, получаем

"(3), ^Л(З) - Д(4)И(3) + (с(Л(4), «(4)) + Д(4)(Д(4), ¿(4)) - RWU(*)) ■ (34)

—♦ —*

Как и в обычном случае, поля Е vi Н связаны соотношением

Н = -±-[Щ3),Ё]. (35)

(4)

Формулы (32) - (35) внешне похожи на аналогичные формулы для электрического и

магнитного полей в случае обычной электродинамики [11]. Разница состоит в том, что

—*

в формулы (32) - (35) входят как 3-мерные вектора R3, так и модуль 4-мерного вектора i?4, в обычной же электродинамике в соответствующих формулах вместо модуля 4-мерного вектора R4 стоит модуль 3-мерного вектора R3.

Формулы (32), (34) дают выражения для электрических и магнитных полей, но с их помощью можно получить выражения и для всех других напряженностей (30). Для этого

следует учесть, что вектора гГ(3), -К(з) в формулах (32), (34) - это 3-мерные вектора, 3

—♦

компоненты которых выбираются из 4-х компонент 4-мерных векторов гГ(4), Я(4). Такие тройки компонент можно выбирать разным способом. Если мы выбрали компоненты х, у, z, то получим электрическое и магнитное поля, которые отличаются от обычных только тем, что теперь их компоненты зависят еще и от переменной s. Если выбрать другую тройку переменных, то получим другие компоненты тензора напряженностей

(30). Так, например, если выбираются переменные у, г, з, то в формулах (31), (33) будет стоять 3-вектор-потенциал Д(3) с компонентами (Ау, А2, А3). Соответственно и в левой части формулы (31) будут стоять те же самые компоненты 4-вектора Ё(4) : (Еу, Ег, Е, =

Ф = 7*40)• Формула (33) дает в этом случае напряженности НХ,СУ = (?2 = ^43-—* —♦

Новые поля (7 появляются за счет того, что теперь имеется 5 потенциалов (А*, А, Аа) и все они зависят еще и от новой переменной з. Приведем явный вид этих полей:

<2 = с2£3 (V ~ «(4)) (в* ~ + - "^#(4)) (%),«(4))) -

с2£(4)

(Л(4),г?(4)) й5

(36)

Для того, чтобы найти явный вид компонент вектора Сг, введем обозначения

Е> = (Яж, Сг), д[з) = (Я», Я«, Д,), и}*| = К, и„ ив). (37)

Верхний индекс (ж) у 3-мерных векторов указывает на ту координату 4-мерного пространства, которая отсутствует у этого 3-мерного вектора.

—*

Используя эти обозначения, получим выражение для вектора И:

и; '

(3)> I С^(3) - Д(4)И(3) + №(4), "(4)) + #(4)(Я(4)>(4)) - #(4)^(4

Б =

С2£(4)

(з)' и(3)

+

С2Я(4)£(4)

X

(38)

Если в формуле (38) заменить верхний индекс (х) на (у), то она вместо вектора Е) с

—*

компонентами (37) даст вектор О' с компонентами

V = (Ну, (Тг, Сгг).

(39)

При этом

Я§ — (Я2, Я3, Ях), — (и2, и3, их).

Таким образом видно, что в зависимости от того, какие компоненты 4-векторов —*

гГ(4),Д(4) входят в формулы (32), (34), они позволяют вычислять различные наборы компонент тензора напряженностей (30).

Рассмотрим конкретный пример частицы с массой т и зарядом е, у которой постоянны пространственные координаты х, у, г, но может меняться координата 5. Иначе

говоря, ее 4-вектора положения, скорости и ускорения имеют вид:

—*

Я(4) = -г, К*), "(4) = (о, 0,0, и), ¿(4) = (0,0,0, (40)

Здесь мы используем введенное ранее обозначение для расстояния в 4-мерном евклидовом пространстве между источником поля и точкой наблюдения. Отметим, что при движении заряда со скоростью и3 происходит изменение расстояния Я3, но связь между ними достаточно сложна и определяется распределением показателя преломления п(г) во всем пространстве между ними. При ее вычислении следует учитывать соотношения (16) - (18) и свойства среды.

Используя формулу (32), получаем выражение для напряженностей электрического поля

с2 - и] + В.3й, -

Я = еС(сД(<)-Л,и,)зД(з)- (41)

Знаменатель этого выражения всегда больше нуля, поскольку скорость заряда иа всегда меньше скорости света в пустоте с, а расстояние по координате 5 всегда меньше полного расстояния Д3, поскольку источник излучения и точка наблюдения всегда пространственно разделены.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Но числитель выражения (41) может менять знак в зависимости от того, какова величина и знак ускорения й5. Если это ускорение достаточно велико, т.е. велика скорость

изменения оптической плотности среды вокруг неподвижного заряда, то напряжение

—#

электрического поля Е поменяет свой знак, т.е. одноименные заряды начнут притягиваться, а разноименные отталкиваться.

Интересно сравнить выражение (41) с выражением для напряженности электрического поля, которое получается с помощью обычных потенциалов Лиенара-Вихерта [12]. Пусть 3-мерная скорость направлена вдоль оси X

из = (и*, 0,0). (42)

В этом случае компонента Ех электрического поля имеет вид:

Ех —

¿¿и

((с2 - и1) (л, - - (Д2 + В,1)йх) . (43)

Здесь величина £(3) имеет вид аналогичный (28), только вместо 4-мерных векторов стоят 3-мерные:

Аз)=%)(о - «(в)(0)-

Из формулы (43) видно, что компонента напряженности Ех может менять свой знак в зависимости от величины и знака ускорения их. Однако на практике такое изменение знака, т.е. переход притяжения в отталкивание и наоборот, наблюдать довольно затруднительно, поскольку сверхбольшие величины ускорений можно создать только на очень короткие промежутки времени. Причем это ускорение должно быть положи тельным, т.е. скорость движения заряженной частицы должна резко возрасти.

Кроме того, эффект смены знака у силы взаимодействия возникает только в том случае, когда скорость частицы направлена не по линии, соединяющей частицу с точкой наблюдения. Иными словами, величина (Щ + Щ) в формуле (43) не должна быть равна нулю. Отсюда следует, что при лобовом движении заряженных частиц друг к другу смены знака силы взаимодействия не происходит, она имеет место только в том случае, когда прицельный параметр не равен нулю. Такая ситуация может возникнуть, например, при падении заряженной частицы на кулоновский центр. Вблизи этого центра напряженность поля велика, что может вызвать большие ускорения падающей частицы и смену притяжения отталкиванием.

Вернемся к анализу напряженностей (30) при условии (40).

Поле ф имеет вид:

® = ^ ((°2 " ^ (К° ~ ~ {К1 + Е2у + П1)й°) • (44)

Из формулы (44) видно, что поле <5 также, как и поле Е, может менять свой знак в

зависимости от знака и величины ускорения й3.

—* —* —*

Поскольку напряженности Н, С? выражаются через поля Е, ф с помощью формул (35), (38), они также меняют свой знак в тот момент, когда это делают поля

Такое изменение знака напряженностей полей и, как следствие, смену знака силы Лоренца можно связать с реакцией излучения этих полей, которое возникает, когда заряженные частицы движутся с ускорением.

Полученные в данной работе результаты показывают, что при обобщении электродинамики на расширенное пространство удается сохранить все основные элементы

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

ЛИТЕРАТУРА

[1] Ц и п е н ю к Д. Ю., Андреев В. А. Краткие сообщения по физике ФИАН, N 6, 23 (2000).

й ц и п е н ю к д. ю„ Андреев В. А. Препринт ИОФАН N 5, М., 1999.

[3] ц и п е н ю к д. ю„ Андреев В. А. Препринт ИОФАН N 9, М., 1999.

[4] ц и п е н ю к д. ю„ Андреев В. А. Препринт ИОФАН N 2, М., 2000.

[5] ц и п е н ю к д. Ю„ Андреев В. А. Препринт ИОФАН N 1, М., 2001.

[6] ц и п е н ю к д. Ю., Андреев В. А. Препринт ИОФАН N 4, М., 2001.

[7] ц и п е н ю к д. ю. Краткие сообщения по физике ФИАН, N 7, 39 (2001).

[8] О к у н ь Л. Б. УФН, 158, 511 (1989).

[9] Л а н д а у Л. Д., Л и ф ш и ц Е. М. Теория поля, М., Наука, 1967.

[10] Болотовский Б. М., Быков В. П. УФН, 160, 511 (1999).

[11] Г и н з б у р г В. Л. Теоретическая физика и астрофизика, М., Наука, 1981.

Институт общей физики РАН Поступила в редакцию 20 марта 2002 г.

i Надоели баннеры? Вы всегда можете отключить рекламу.