Rambler's Top100Astronet    
  по текстам   по ключевым словам   в глоссарии   по сайтам   перевод   по каталогу
 

На первую страницу
Теория относительности для астрономов

<< 4. Эксперименты по ОТО | Оглавление | 6. Анализ в неэвклидовой >>

Разделы


5. Неэвклидова геометрия

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

В произвольной геометрии рассматриваются произвольные преобразования координат:

\begin{displaymath}
x^{\mu} = f^{\mu}(\tilde x^{\nu}),
\end{displaymath}

например,

\begin{displaymath}
\begin{array}{l}
x^0 = \tilde x^0, \\
\qquad \\
x^1 = \til...
...2, \\
\qquad \\
x^3 = \tilde x^1 \cos \tilde x^2,
\end{array}\end{displaymath}

как видно из приведенных преобразований, они описывают просто переход от сферических к декартовым координатам.

Дифференциал в нетильдованной системе связан с дифференциалом в системе координат с тильдой уравнениями вида:

\begin{displaymath}
d x^{\mu} = {\displaystyle\partial f^{\mu}\over\displaystyle...
...mu}\over\displaystyle\partial \tilde x^{\nu}} d\tilde x^{\nu}
\end{displaymath}

В геометрии вводится понятие геометрического объекта. Не перечисляя всех объектов геометрии назовем только скаляры, векторы и тензоры. Строго говоря, все перечисленные величины являются тензорами, так скаляр является тензором нулевого ранга, а вектор тензором первого ранга. Однако исторически принято выделять скаляры и векторы.

Скалярная величина (или скаляр) при преобразованиях системы координат не преобразуется,

\begin{displaymath}
\varphi(x^{\mu}) = \varphi(\tilde x^{\mu})
\end{displaymath}

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

\begin{displaymath}
A^{\mu} = {\displaystyle\partial x^{\mu}\over\displaystyle\partial \tilde x^{\nu}} \tilde A^{\nu}
\end{displaymath} (5.1)

для контравариантных компонент и

\begin{displaymath}
A_{\mu} = {\displaystyle\partial \tilde x^{\nu}\over\displaystyle\partial x^{\mu}} \tilde
A_{\nu}
\end{displaymath} (5.2)

для ковариантных компонент.

Как можно заметить, из закона преобразования (5.1) преобразование контравариантных компонент какого - либо вектора подобно преобразованию дифференциалов координат, представляющих разность положений двух точек. Следовательно контравариантные компоненты вектора можно представить как систему из двух точек, с указанием какая из точек является первой, а какая является второй. Примером контравариантных компонент является четырехмерная скорость какого - либо тела.

Примером ковариантных компонет вектора является градиент некоторого поля:

\begin{displaymath}
A_{\mu} = {\displaystyle\partial \varphi(x^{\alpha})\over\displaystyle\partial x^{\mu}}
\end{displaymath}

Образом градиента от поля являются наборы двумерных поверхностей определяющих постоянство поля $\varphi(x^{\alpha})$. Следовательно геометрическим образом ковариантных компонент тензора являются двумерные поверхности "раскрашенные" в различные "цвета" для того, чтобы определить какая из поверхностей является первой, а какая - второй.

Выпишем также правило по которому преобразуются компоненты тензора второго ранга:

\begin{displaymath}
\begin{array}{l}
A_{\mu \nu} = {\displaystyle\partial \tilde...
...aystyle\partial x^{\nu}}
\tilde A^{\alpha}_{\beta}
\end{array}\end{displaymath}

Геометрическая величина $A^{\mu_1 \mu_2 ... \mu_n}_{\nu_1 \nu_2 ...
\nu_m}$, которая имеет $n$ верхних ипндексов и $m$ нижних индексов и преобразуется по закону:

\begin{displaymath}
\begin{array}{l}
A^{\mu_1 \mu_2 ... \mu_n}_{\nu_1 \nu_2 ... ...
...ha_2 ... \alpha_n}_{\beta_1 \beta_2 ...
\beta_m}
\end{array}
\end{displaymath}

называется тензором ранга $n+m$ с $n$ контравариантными $m$ ковариантными индексами.

В алгебре тензоров определяется суммирование тензоров одинакового ранга и с одинаково расположенными индексами:

\begin{displaymath}
C^{\mu}_{\nu} = A^{\mu}_{\nu} + B^{\mu}_{\nu},
\end{displaymath}

умножение определено для тензоров любого ранга с произвольно расположенными индексами:

\begin{displaymath}
C^{\mu \alpha \beta}_{\nu} =A^{\mu}_{\nu} \cdot B^{\alpha \beta}
\end{displaymath}

В тензорной алгебре определяют также свертку двух тензоров по правилу:

\begin{displaymath}
C^{\mu \nu} = \sum_{\alpha} A^{\mu \alpha} \cdot B^{\nu}_{\alpha}
\end{displaymath}

для произвольного расположения индексов важно только, чтобы индекс по которому проводится свертка в первом тензоре был верхним, во втором - нижним (или наоборот, в первом - нижним, а во втором -верхним). Положение и количество остальных индексов переносится в свертку по их месту. Свертки могут быть образованы различными способами, однако каждое суммирование должно выполняться по верхнему и нижнему индексу. Например, сумма

\begin{displaymath}
\sum_{\alpha} A^{\mu}_{\alpha} \cdot B^{\nu}_{\alpha}
\end{displaymath}

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

В современной тензорной алгебре, анализе и физике пользуются правилом суммирования Эйнштейна, который для удобства ввел следущее обозначение:

\begin{displaymath}
\sum_{\alpha} A^{\mu \alpha} \cdot B^{\nu}_{\alpha} \equiv A^{\mu \alpha}
\cdot B^{\nu}_{\alpha}
\end{displaymath}

при этом индекс $\alpha$ называют мертвым индексом и по повторяющимся индексам производится суммирование. Так запись $C_{\mu \nu}A^{\mu}B^{\nu}$ обозначает:

\begin{displaymath}
C_{\mu \nu}A^{\mu}B^{\nu} \equiv \sum_{\mu=0, \nu=0}^{\mu=3, \nu=3}
C_{\mu \nu}A^{\mu}B^{\nu}
\end{displaymath}

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

Определим также обратные тензоры согласно уравнению:

\begin{displaymath}
A_{\mu \alpha} \cdot B^{\nu \alpha} = \delta^{\mu}_{\nu}
\end{displaymath}

Здесь $\delta^{\mu}_{\nu}$ - символ Кронекера, который определяется следущими условиями $\delta^{\mu}_{\nu} =1$ если $\mu = \nu$ и $\delta^{\mu}_{\nu}=0$ если $\mu \ne \nu$. Можнл также выписать матрицу символа Кронекера

\begin{displaymath}
\delta^{\mu}_{\nu} =\left(
\begin{array}{rrrr}
1 & 0 & 0 & 0...
...& 0 \\
0 & 0 & 1 & 0 \\
0 & 0 & 0 & 1 \\
\end{array}\right)
\end{displaymath} (5.3)

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

\begin{displaymath}
A^{\mu} B_{\mu} = {\displaystyle\partial x^{\mu}\over\displa...
...laystyle\partial x^{\mu}} \tilde A^{\alpha}
\tilde B_{\beta}.
\end{displaymath} (5.4)

Различные авторы [8], [9], [12] вводят также операцию симметрирования, албтернирования, выключения индексов и т.п. Для нас важной остается только последняя операция над тензорами - образование следа, образование поперечной части тензора и образование поперечно - бесследовой части тензора, которые будут применяться во второй части книги. Однако для их определения нам придется использовать метрику, поэтому отложим введения этих понятий после обсуждения метрики.

5.1 Метрика искривленного пространства - времени

Геометрия четырехмерного пространства - времени полностью определяется десятью функциями, которые являются компонентами симметричного тензора второго ранга. Метрика четырехмерного интервала есть:


\begin{displaymath}
d s^2 = g_{\mu \nu} d x^{\mu} d x^{\nu}
\end{displaymath} (5.5)

Здесь $g_{\mu \nu}$ - ковариантные компоненты метрического тензора или, как обычно говорят, метрический тензор второго ранга.

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

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

\begin{displaymath}
\eta_{\mu \nu} =\left(
\begin{array}{rrrr}
1 & 0 & 0 & 0  ...
...0 \\
0 & 0 & -1 & 0 \\
0 & 0 & 0 & -1 \\
\end{array}\right)
\end{displaymath} (5.6)

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

\begin{displaymath}
A_{\mu} = g_{\mu \alpha} A^{\alpha}
\end{displaymath}

Теперь можно определить норму вектора, аналогично тому, как это делается в эвклидовой геометрии. Норомой вектора по определению полагают величину [8], [9], [12]:

\begin{displaymath}
\vert\vert A \vert\vert^2 = A_{\mu} A^{\mu} = g_{\mu \nu} A^{\mu} A^{\nu}
\end{displaymath} (5.7)

В зависимости от того, чему равна норма ветора его называют:

Норма называется также длиной вектора. Как видно из определения изотропного вектора он обладает нулевой длиной.

Рассмотрим свертку двух векторов. Два вектора, которые удовлетворяют условию:

\begin{displaymath}
A_{\mu} B^{\mu} = g_{\mu \nu} A^{\mu} B^{\nu} = 0
\end{displaymath}

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

Кроме этого важного свойства приведем также уравнение описывающее угол $\varphi$ между двумя векторами $A^{\mu}$ и $B^{\nu}$ в неэвклидовой геометрии:

\begin{displaymath}
\cos \varphi = {\displaystyle g_{\mu \nu} A^{\mu} B^{\nu}\ov...
...isplaystyle\vert\vert A\vert\vert\cdot \vert\vert B\vert\vert}
\end{displaymath} (5.8)

Здесь следует заметить, что оба вектора $A$ и $B$ являются либо пространственноподобными, либо изотропными.

5.1.1 Примеры метрики искривленного пространства

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

Приведем в качестве примера метрику на поверхности сферы. Для этого выпишем вначале метрику эвклидова пространства в сферической системе координат:

\begin{displaymath}
dl^2=dr^2 + r^2 \left(d\theta^2 + \sin^2 \theta d\varphi^2\right)
\end{displaymath}

Будем считать, что мы вычисляем расстояние только между теми точками, которые лежат на поверхности сферы радиуса $R$. Тогда можно полагать, что $d R =0$, а значит метрика получает вид:

\begin{displaymath}
dl^2= R^2 \left(d\theta^2 + \sin^2 \theta d\varphi^2\right)
\end{displaymath} (5.9)

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

5.2 Геодезические линии в неэвклидовой геометрии

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

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

Длина линии в геометрии общей теории относительности есть интеграл по интервалу между двумя событиями разделенными по времени и пространству:

\begin{displaymath}
S = \int d s = \int_{(L)} \sqrt{g_{\mu \nu} {\displaystyle d...
...le d s}
{\displaystyle d x^{\nu}\over\displaystyle d s}} d s
\end{displaymath}

Линии экстремального пути удовлетворяет условию:

\begin{displaymath}
\delta S = \int \delta d s= 0
\end{displaymath}

Здесь символ $\delta$ означает первую вариацию. Из этого условия получается уравнение геодезической линии.

Для вычисления вариации величины $S$ применим следущий прием. Проварьируем не дифференциал интервала $d s$, а квадрат дифференциала $ds^2$. Вариации самого дифференциала по вариациям квадрата можно получить согласно уравнению $\delta ds^2 = 2ds \delta ds$. С другой стороны, вариации квадрата $ds^2$ есть:

\begin{displaymath}
\delta ds^2 = \delta \left( g_{\mu \nu} d x^{\mu} d x^{\nu}\right)
\end{displaymath}

правую часть этого равенства можно представить в виде суммы:

\begin{displaymath}
\delta ds^2 = \delta g_{\mu \nu} d x^{\mu} d x^{\nu} + g_{\m...
...d
x^{\mu} d x^{\nu} + g_{\mu \nu} d x^{\mu} \delta d x^{\nu}
\end{displaymath}

Воспользуемся равенствами $\delta d x^{\mu} =d \delta x^{\mu}$ и $\delta
g_{\mu \nu} = {\displaystyle\partial g_{\mu \nu}\over\displaystyle\partial x^{\alpha}} \delta
x^{\alpha}$ и получим уравнение для вариации дифференциала интервала:

\begin{displaymath}
\delta(ds) = {\displaystyle\delta ds^2\over\displaystyle 2ds...
...}
{\displaystyle d x^{\nu}\over\displaystyle d s}
\right) ds
\end{displaymath}

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

\begin{displaymath}
\int ds g_{\mu \nu} {\displaystyle d(\delta x^{\mu})\over\di...
...tyle d x^{\nu}\over\displaystyle d s} \right)
\delta x^{\mu}
\end{displaymath}

Воспользуемся теперь этими равенствами и вычислим первую вариацию полного пути:

\begin{displaymath}
\delta S = \int \delta ds =
g_{\mu \nu} \delta x^{\mu} {\di...
...{\nu}\over\displaystyle d s}
\right)
\rbrace \delta x^{\alpha}
\end{displaymath}

В этом уравнении первый член после второго знака равенства представляет из себя вариации в конечных точках пути, по определению эти вариации равны нулю [13], [14]. Член под знаком интеграла в фигурных скобках должен обращаться в нуль на геодезической линии. Это и дает нам уравнение геодезической линии:

\begin{displaymath}
{\displaystyle d u_{\alpha}\over\displaystyle d s} = \frac{1...
...mu
\nu}\over\displaystyle\partial x^{\alpha}} u^{\mu} u^{\nu}
\end{displaymath} (5.10)

Здесь $u^{\mu}$ - вектор касательный к геодезической линии.

Уравнение, которое мы получили является уравнением второго порядка. В случае эвклидовой геометрии, когда все компоненты метрического тензора являются диагональными и равны единице, правая часть равенства (5.10) равна нулю. Тогда решение этого уравнения есть $u_{\alpha} =
const$. Это решение описывает прямые линии, которые являются линиями кратчайшего пути в эвклидовой геометрии.

Приведем это уравнение к стандартной форме, принятой в неэвклидовой геометрии. Для этого воспользуемся тождеством:

\begin{displaymath}
{\displaystyle d\over\displaystyle d s}\left( g_{\alpha \nu}...
...\alpha \nu}\over\displaystyle\partial x^{\mu}} u^{\mu} u^{\nu}
\end{displaymath}

Подставив его в левую часть уравнения (5.10) получаем уравнение вида:

\begin{displaymath}
{\displaystyle d u^{\alpha}\over\displaystyle d s} + \frac{1...
...ver\displaystyle\partial x^{\beta}}\right) u^{\mu} u^{\nu} = 0
\end{displaymath}

Величина

\begin{displaymath}
\Gamma^{\alpha}_{\mu \nu} = \frac{1}{2} g^{\alpha \beta}\lef...
...artial g_{\mu \nu}\over\displaystyle\partial x^{\beta}}\right)
\end{displaymath} (5.11)

называется символом Кристоффеля. Символ Кристоффеля имеет три индекса и не является тензором. В алгебре такие величины носят название символы первого рода.

Стандартный вид уравнения геодезической линии в неэвклидовой геометрии записывается с помощью символа Кристоффеля:

\begin{displaymath}
{\displaystyle d u^{\alpha}\over\displaystyle d s} + \Gamma^{\alpha}_{\mu \nu} u^{\mu} u^{\nu} = 0
\end{displaymath} (5.12)

Кроме этого, стандартного вида уравнения геодезической линии, можно также записать как меняется дифференциал касательного вектора $u^{\mu}$ при переносе вдоль геодезической линии:

\begin{displaymath}
d u^{\alpha} =- \Gamma^{\alpha}_{\mu \nu} u^{\mu} d x^{\nu}
\end{displaymath} (5.13)

5.2.1 Геодезические линии на поверхности сферы

Рассмотрим в качестве примера геодезические линии на сфере - двумерной неэвклидовой поверхности. Метрический тензор на поверхности сферы имеет вид:

\begin{displaymath}
g_{\mu \nu} =\left(
\begin{array}{cc}
1 & 0 \\
0 & \sin^2 \theta \\
\end{array}\right)
\end{displaymath} (5.14)

Контравариантные компоненты метрического тензора на поверхности сферы есть:

\begin{displaymath}
g^{\mu \nu} =\left(
\begin{array}{cc}
1 & 0 \\
0 & {\displaystyle 1\over\displaystyle\sin^2 \theta} \\
\end{array}\right)
\end{displaymath} (5.15)

Вычисляя символ Кристоффеля получаем, что все компоненты равны нулю за исключением двух:

\begin{displaymath}
\begin{array}{l}
\Gamma^1_{22} = - \sin \theta \cos \theta \\
\qquad \\
\Gamma^2_{12} = \ctg \theta
\end{array}\end{displaymath}

Уравнения геодезических линий имеют вид:

\begin{displaymath}
{\displaystyle d^2 \theta\over\displaystyle d\lambda^2} - \s...
...playstyle d
\varphi\over\displaystyle d \lambda} \right)^2 =0
\end{displaymath} (5.16)

и


\begin{displaymath}
{\displaystyle d^2 \varphi\over\displaystyle d\lambda^2} + 2...
...bda} {\displaystyle d \varphi\over\displaystyle d \lambda} =0
\end{displaymath} (5.17)

Здесь $\theta$, $\varphi$ сферические координаты наклонения и азимута, а $\lambda$ - афинный параметр вдоль геодезической линии.

При приведении к первому интегралу уравнение (5.17) получает вид:

\begin{displaymath}
{\displaystyle d\over\displaystyle d \lambda}\left( \ln {\di...
...aystyle d \lambda}\left( \ln \vert\sin \theta\vert \right) = 0
\end{displaymath}

или


\begin{displaymath}
\sin^2 \theta {\displaystyle d \varphi\over\displaystyle d \...
...ta_0 {\displaystyle d
\varphi\over\displaystyle d \lambda}_0
\end{displaymath} (5.18)

Здесь $\sin^2 \theta_0 {\displaystyle d \varphi\over\displaystyle d \lambda}_0$ обозначает константу интегрирования. Первый интеграл уравнения (5.16) можно легко привести к удобной форме, если выбрать в качестве афинного параметра длину вдоль геодезической $d s = d \lambda$. Тогда интеграл первого уравнения становится:

\begin{displaymath}
\left( {\displaystyle d \theta\over\displaystyle d \lambda} ...
...splaystyle d \varphi\over\displaystyle d \lambda} \right)^2 =1
\end{displaymath} (5.19)

Рассмотрим два частных случая геодезических линий.

Движение по экватору

Пусть наша линия начинается на экваторе $\theta_0 = {\displaystyle\pi\over\displaystyle 2}$, причем направление этой линии совпадает с направлением экватора, т.е. ${\displaystyle d
\theta\over\displaystyle d \lambda} _0=0$ и ${\displaystyle d \varphi\over\displaystyle d \lambda} _0=1$.

Решение для геодезической линии получается из уравнений (5.18, 5.19). Действительно, правая часть равенства (5.18) есть 1. Подставляем выражение для ${\displaystyle d \varphi\over\displaystyle d \lambda}$ в уравнение (5.19) и получаем уравнение:

\begin{displaymath}
({\displaystyle d \theta\over\displaystyle d \lambda})^2 = - \ctg^2 \theta
\end{displaymath}

Удовлетворяет этому уравнению только значения $\theta= {\displaystyle\pi\over\displaystyle 2}$. Таким образом доказано, что экватор является геодезической линией на сфере.

Движение по меридиану

Рассмотрим теперь движение вдоль меридиана и покажем, что меридиан также является геодезической линией. Теперь вектор, касательный к линии имеет компоненты:

\begin{displaymath}
{\displaystyle d \varphi\over\displaystyle d \lambda} _0=0 \...
...m} {\displaystyle d \theta\over\displaystyle d
\lambda} _0=1
\end{displaymath}

Начальную точку также выберем на экваторе $\theta_0 = {\displaystyle\pi\over\displaystyle 2}$. Из уравнения (5.18) получаем, что

\begin{displaymath}
{\displaystyle d \varphi\over\displaystyle d \lambda} = 0
\end{displaymath}

откуда следует, что $\varphi(\lambda) = const$. Из второго интеграла системы (5.19) получаем, что

\begin{displaymath}
{\displaystyle d \theta\over\displaystyle d \lambda} = \pm 1
\end{displaymath}

откуда непосредственно следует, что

\begin{displaymath}
\theta (\lambda) =\pm \lambda + {\displaystyle\pi\over\displaystyle 2}
\end{displaymath}

Не рассматривая общего случая решения системы (5.18, 5.19) упомянем только, что экватор и меридиан на поверхности сферы являются частными случаями большого круга - сечения поверхности сферы плоскостью, которая проходит через центр сферы. Рассматривая трехмерные вращения, можно показать, что только отрезок большого круга на сфере представляет кусок геодезической линии. Полная геодезическая линия представляет из себя большой круг.

Движение по малому кругу

Продемонстрируем отличие произвольной линии от геодезической линии на примере малого круга на сфере. Малый круг на сфере - это сечение сферы плоскостью, которая не проходит через центр сферы.

Выберем малый круг так, чтобы он был сечением с $\theta = const$. Тогда уравнения, которые описывают малый круг есть:

\begin{displaymath}
\begin{array}{l}
\theta(\lambda) = \theta_0 \ne {\pi \over 2...
...qquad \\
\varphi( \lambda) = \lambda + \varphi_0 ,
\end{array}\end{displaymath}

отсюда легко найти первые и вторые производные от координат по афинному параметру:

\begin{displaymath}
\begin{array}{l}
{\displaystyle d \varphi\over\displaystyle ...
...yle d^2
\theta\over\displaystyle d \lambda^2} = 0
\end{array}\end{displaymath}

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

\begin{displaymath}
{\displaystyle d^2 \theta\over\displaystyle d \lambda^2} -\s...
...aystyle d \varphi\over\displaystyle d \lambda}\right)^2 \ne 0
\end{displaymath}

Таким образом малый круг на сфере не является геодезической линией.



<< 4. Эксперименты по ОТО | Оглавление | 6. Анализ в неэвклидовой >>

Публикации с ключевыми словами: Общая теория относительности - специальная теория относительности - системы отсчета - измерение времени
Публикации со словами: Общая теория относительности - специальная теория относительности - системы отсчета - измерение времени
См. также:
Все публикации на ту же тему >>

Мнения читателей [21]
Оценка: 3.1 [голосов: 127]
 
О рейтинге
Версия для печати Распечатать

Астрометрия - Астрономические инструменты - Астрономическое образование - Астрофизика - История астрономии - Космонавтика, исследование космоса - Любительская астрономия - Планеты и Солнечная система - Солнце


Астронет | Научная сеть | ГАИШ МГУ | Поиск по МГУ | О проекте | Авторам

Комментарии, вопросы? Пишите: info@astronet.ru или сюда

Rambler's Top100 Яндекс цитирования