Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Гуревич, А. Г. Магнитный резонанс в ферритах и антиферромагнетиках

.pdf
Скачиваний:
138
Добавлен:
23.10.2023
Размер:
31.33 Mб
Скачать

400

СПИНОВЫЕ в о л н ы

[ГЛ . 8

(7.1.21),

позволит определить эти частоты с учетом

магнитного

взанмо действия.

Ограничимся рассмотрением первого — коллииеарного основ-

лого

состоящія феррпмагнетика, которое,

как

было показано в

§ 4.4,

имеет место при Я 0 ^ Л ( М х 0 — М

2 0)

(здесь Л — кон­

станта однородного обменного взаимодействия между подрешетка­ ми, М х о п М 2 о — статические намагниченности подрешеток) . Уравнения движения в этом состоянии с учетом эффективных полей одноосной кристаллографической анизотропии и эффек­

тивных полей (8.2.36) в циркулярных переменных mli2 і- =

r +

+ imU2y, как

нетрудно

убедиться, примут вид1)

 

[+ ю — у2 (Яі0 -f- Н е-і +

77.u + ffm)l mi+

 

[+ СО— у2( Н і о

— Ti (HEI 4" Hkl) ra, -|; =

XiMi Л ъ

(8.2.37)

I I E l I I А і —■I I кг)]т і + +

 

 

 

+ Т г Е -i + I I кг) т і +

— Т г 372 a h+ .

 

Здесь Н і0 — внутреннее постоянное поле; 77а і,г — эффективные поля одноосной кристаллографической анизотропии (в предпо­ ложении аддитивности энергий анизотропии подрешеток и с уче­ том только первых констант они имеют вид (4.4.38)); величины

77к і,2 п Як 1,2, аналогично (8.2.6), обозначают

зз

Я киг = 3 7 1.2о 2 2 Чіаps kpkSi

Р = 1 s = l

(8.2.38)

3 3

77KI,2= 371>202 2 4pshpks.

P = T s = l

Из уравнений (8.2.37) видно, что учет неоднородного обмен­ ного взаимодействия свелся к замене, аналогичной (8.2.7),

Яи 1.2 —> Яд ij2 + 77к 1.2,

(8.2.39)

Ял1.2—>77аі,2 77к],2•—IIк іл-

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

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

!) Аналогичные уравнения (4.4.37) были записаны для собственных од­ нородных колебаний. Теперь, поскольку мы рассматриваем колебаппя под действием заданного внутреннего переменного поля (ищем внутренний тепзор восприимчивости), из уравнений следует исключить члены Ti,2Nj_/l'/i,2 0

(но оставить члены Т і ЛИг{М1 0 — Мг „), учитывающие постоянные размаг­ ничивающие поля).

§ 8.2] ВОЛНЫ В АНТИФЕРРОМ АГНЕТИКАХ И ФЕРРИМ АГНЕТИКАХ 40]

породного обменного взаимодействия из соответствующих выра­ жений для однородных колебаний (§4.4).

Произведя замену (8.2.39) в выражении (4.4.39) (без члена N ± (М± о — М 2 о)). мы получим частоту первой (низкочастотной) ветви спектра спиновых воли ферримагнетика без учета магнит­ ного взаимодействия

®х = Тэфф (Н і0 + Н а эфф + Н кэфф).

( 8 .2 . 4 0 )

Здесь т'аффи Н а офф определяются выражениями (4.4.16) и (4.4.40),

а

1(1)(Ь

М10(Нк1 - Н 1 г) + М2 0 (Н L-2 _и'

(8.2.41)

Н н о

 

 

Если учесть вид тензоров qx, q2 и q', определяемый одноосной сим­ метрией задачи, и пренебречь, как и раньше, различием продоль­ ных и поперечных компонент этих тензоров, то

М \0*71-

■2Л/і оЛ/о oq' -)- M ^ 0 q 2

£2 = Д ф Ф/с2. (8.2.42)

Пи офф —

Міо — М 2 о

И, наконец, если принять (хотя трудно привести какие-либо до­

воды в пользу этого) = q2 = q' =

q, то

 

Ни офф = (М г о -

М 8 о) д/<2.

(8.2.43)

Сравним выражения (8.2.40) с выражением (8.1.42) для часто­ ты спиновых волн в одноосном ферромагнетике при постоянном поле, направленном, как и для ферримагнетика, вдоль оси1). Если в (8.1.42) положить Ѳ^- = 0 (что равносильно неучету маг­ нитного взаимодействия), то указанные выражения совпадают при

о 7\Іо о ■= M Q I I аі — Н а эфф, Т Тэфф и D = ПЭфф.

Таким образом, в рассматриваемом частном случае (низкочастот­ ная ветвь спиновых воли, одноосная анизотропия п другие огово­ ренные выше условия) ферримагнетик эквивалентен ферромаг­ нетику с эффективными параметрами. Легко убедиться, что такая

эквивалентность имеет место и без предположения о равенстве

*-У <-*■

продольных и поперечных компонент тензоров qb 2 и q'. Имеет она место, конечно, и для других ориентаций постоянного поля,

идля кристаллов с другой (не одноосной) симметрией.

1)R выражении (8.1.42) (как п в главах 2 и 7) индекс 1 указывает на связь /7.И с первой константой анизотропии. В этом параграфе (как и в главе 4)

учитывается только первая константа анизотропии, а индексы 1 и 2 обозна­ чают номера подрешеток,

402 СПИНОВЫЕ ВОЛНЫ [ГЛ. 8

Производя замену (8.2.39) в выражениях компонент тензора

Ч-»

%для низкочастотных однородных колебаний J) или же непосред­

ственно

решая систему

(8.2.37) при И л і,2, Ни1Л,Нк і, 2 <С

< С е і

— НЕі) и Я 0< ^ ( Я в і

— Я е 2), получим для низкочастотной

области

(со << töjj == т2 НЕі — ГіПе2)

 

 

(8.2.«)

где % определяется выражением (8.2.40). Как видно из (8.2.44), тензор восприимчивости ферримагнетика с учетом неоднородного обменного взаимодействия в низкочастотной области совпадает с тензором восприимчивости ферромагнетика с эффективными пара­ метрами. Отсюда ясно, что для низкочастотной («ферромагнитной») ветви спектра спиновых волн при достаточно малых внешних по­ лях (Н 0 HEX Нео) закон дисперсии п с учетом магнитного взаимодействия будет таким же, как для ферромагнетика.

Исследуем теперь вторую (высокочастотную) ветвь спектра спиновых волн. При этом рассмотрим для простоты изотропный ферримагнетик (HAI = Я л2 == 0). Тогда, приравнивая нулю оп­ ределитель системы (8.2.37) ло-прежпему в предположении, что

Но, Я п . 2 , н 'кг, г *< Н ы НЕг, но не для низкочастотной обла­

сти, а, наоборот, для высоких частот (to

Т0фф#о), получим

coo = сс>Е ТэффЯіо “Ь ТлФфЯ^офф.

(8.2.45)

Здесь сод и Трфф определяются

соответственно формулами (4.4.18)

и (4.4.19), а величина

 

 

 

 

И к яфф —

 

»kа + " £Л'2

Ны - Ны

. (8.2.46)

Л і Ң

_

ң

 

 

 

T i

Я і

T2

11 Ег

 

Если принять ?і =

Я.’ = Q

т0

 

 

 

И к эфф — (Ті — Т2 ) 2 M l аМч о . а

(8.2.47)

 

Ч\Міо — і\Мго

 

Выражение (8.2.45) получено без учета магнитного взаимодей­ ствия. Для учета его следовало бы подставить восприимчивость для второй ветви спиновых воли, полученную в результате реше­ ния системы (8.2.37), в магнитостатическое уравнение (7.1.2). Но восприимчивость эта, как отмечалось в § 4.4, имеет малый числи­ тель — пропорциональный (ух — Тг)2- Частота же второй ветви (мы рассматриваем область небольших Н 0) велика. Поэтому по­

!) Эти выражения с учетом кристаллографической анизотропии нс были приведены в § 4.4, но их нетрудно получить, решая уравнении (4.4.37) без членов с N ± .

§ 8.3] ВОЛНЫ В ОГРАНИЧЕННЫ Х ТЕЛАХ И НЕОДНОРОДНЫХ СРЕДАХ 403

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

В заключение заметим, что кристаллические и магнитные струк­ туры антиферромагнетиков и ферримагнетиков отличаются большим разнообразием. Столь же разнообразны и спектры спиновых волн в них. Небольшое число рассмотренных выше простых примеров мог­ ло лишь дать представление о некоторых особенностях этих спектров и познакомить с методом их расчета, основанным на решении урав­ нений овижепия намагниченностей1). Результаты расчета спект­ ров для различных антиферромагнетиков приведены в [21] и [244].

§ 8.3. Спиновые волны в ограниченных телах

инеоднородных средах

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

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

Задачи о распространении электромагнитных волн в ограничен­ ных телах рассматривались в главах С и 7. Но в главе 6 это рассмот­ рение велось без учета зависимости параметров среды от частоты и волнового вектора. В главе 7 зависимость магнитной проницае­ мости от частоты принимались во внимание, но зависимость от волнового вектора (пространственная дисперсия), обусловленная неоднородным обменным взаимодействием, не учитывалась. Как уже отмечалось в главе 7, совместный учет неоднородного обмен­ ного взаимодействия и граничных условий необходим только в тех случаях, когда размер тела (в каком-либо одном измерении) становится сравнимым пли меньшим, чем величина

‘' - т г - Ѵ

' Т

<8 -З Л )

(&! — такое волновое число, при

котором вклад неоднородного

обменного взаимодействия в частоту спиновых

волн становится

существенным). Оценка для иттрий-железиого

граната (г| = 0,1)

и сантиметрового диапазона волн дает

 

dx— ІО-5

ІО-6 см.

 

J) Других методов расчета спектров спиновых волн в ферро- и антпфер- ро-магпетиках мы коснемся в § 8.5.

404

СПИНОВЫЕ ВОЛНЫ

ИГЛ. 8

Таким образом, улет граничных условий (в слабо проводящих средах) становится необходимым лишь в случае мелких частиц, а также тонких стержней (нитей) и топких пленок, поперечные размеры которых измеряются долями микрона 1). Следует подчерк­ нуть, что это относится лишь к влиянию граничных условий на спектр спиновых волн. Что же касается таких задач, как, напри­ мер, задача о прохождении и отражении спиновых волн через границы раздела, то для них, конечно, учет граничных условий на поверхности раздела необходим независимо от размеров тел.

Аналогичное положение плюет место и для сред с плавно из­ меняющимися параметрами, в частности, для случая неоднород­ ного постоянного магнитного поля. Если изменение параметров — достаточно медленное, т. е. расстояния, на которых оин заметно изменяются, гораздо больше, чем 1/7с (где к — волновое число рассматриваемых волн), то спектр волн в каждой точке определя­ ется значениями параметров в этой точке. Задачи о распростране­ нии спиновых воли в такой среде с медленно изменяющимися па­ раметрами аналогичны задачам геометрической оптики. Если же расстояния, на которых происходят заметные изменения парамет­ ров, сравнимы или меньше, чем 1/к, то магнитиые колебания не являются плоскими волнами, и зависимость m (г, t) может быть найдена лишь путем совместного решения уравнений движения намагниченности и уравнений ноля (обычно — в магнитостати­ ческом приближении) при заданной зависимости параметров от координат.

Возбуждение спиновых волн. Исследование распространения спиновых волн в неоднородных средах тесно связано с проблемой возбуждения спиновых воли. Некогерентные (тепловые) спиновые волны возбуждаются под воздействием тепловых колебаний кри­ сталлической решетки. При этом возникает некоторое равновесное распределение таких спиновых волн с различными а>и к. Когерент­ ные спиновые волны, которые нас больше всего интересуют, по­ могут возбуждаться в результате нелинейных процессов под воз­ действием однородных колебаний намагниченности [522] или непосредственно переменного магнитного поля [528] достаточно большой амплитуды (см. также [537, 15]). Рассмотрение такого возбуждения выходит за пределы этой книги, посвященной ли­ нейным проблемам 2). Не будем мы рассматривать и возбужде­ ния спиновых волн упругими колебаниями, которое происходит благодаря связи между магнитной системой и решеткой в магнитоупорядоченных кристаллах [347, 3]. Остановимся лишь на линей­

х) В металлах, в которых глубина проникновения поля (толщина скинслоя) сравнима с ІД і, учет граничных условий в задачах о спиновых вол­ нах необходим для образцов любых размеров.

2) Вопрос о нелинейном возбуждении спиновых волн (спннопых неста­ бильностях) будет лишь слегка затронут в § 9.2.

§ 8.3І ВОЛНЫ В ОГРА Н И ЧЕН Н Ы Х ТЕЛАХ И Н ЕО Д Н О РО ДН Ы Х СРЕДА Х 405

ном (имеющем место при сколь угодно малых амплитудах) воз­ буждении когерентных спиновых волн переменным магнитным полем.

Не проводя строгого исследования задачи о возбуждении спиновых волн, отметим лишь, что ее решение в случае дискрет­ ного спектра спиновых волн *) можно искать, так же как и для магнитостатических (безобменных) волн, в виде ряда (7.3.30). Величины т ѵ(г) в (7.3.30) будут теперь комплексными амплиту­ дами спиновых колебаний с собственными частотами со ѵ, а для ко­ эффициентов Сѵ при соответствующей нормировке собственных функций піу (г) будет справедлива формула (7.3.31).

Если же спектр спиновых волн — непрерывный, как например, в неограниченной среде или волноводе 2),* то ряд ѵ.7.3.30) должен быть заменен интегралом. И если собственными функциями задачи являются плоские волны іщ (г), то интегрирование может вестись по волновым векторам к этих волн:

(8.3.2)

Величина ск (к) характеризует интенсивность возбуждения нор­ мальной волны с волновым вектором к. Для нее справедлива фор­ мула, аналогичная (7.3.31),

(8.3.3)

где — заданная частота возбуждающего поля, а а к — параметр диссипации; интегрирование в (8.3.3) производится по всему про­ странству, где h (г) ф 0 и тк (г) ф 0.

Рассмотрим в качестве самого простого примера возбуждение бегущей спиновой волны в однородной неограниченной среде. Тогда

пц (г) = me~lkr

(8.3.4)

а он- собственная частота, связанная с к дисперсионными соотно­ шениями, которые подробно исследовались выше. В этом случае возбуждение происходило бы наиболее интенсивно, если бы воз­ буждающее поле представляло собой бегущую волну с тем же вол­ новым вектором, что и спиновая волна, для которой со я- = <»• Наибольший практический интерес мог бы представить случай, когда возбуждающее переменное поле однородно. Однако тогда интеграл в (8.3.3) обращается в нуль. Действительно, в случае

*) Дискретный спектр спиновых волн, как и в случае безобменных маг­ нитостатических волн (§ 7.3), электромагнитных волн (§ 6.3) п других волпопых процессов, имеет место для ограниченных тел.

2) См. примечание 1 на стр. 340.

406 СПИНОВЫЕ ВОЛНЫ [ГЛ. 8

однородного поля h (ѵ) может быть вынесено на знак интеграла и ии_

теграл ^ in* (г) dV, который представляет собой величину,

V'

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

переменным полем (если не касаться нелинейных процессов) не­ возможно.

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

Граничные условия для намагниченности. Учет неоднородно­ го обменного взаимодействия приводит, как мы видели, к появле­ нию в уравнениях движения производных от намагниченности по координатам. Порядок дифференциальных уравнений, которые получаются при совместном решении уравнений движения и урав­ нений электромагнитного поля, повышается, и обычных электро­ динамических граничных условий становится недостаточно для однозначного решения граничных задач. Электродинамические граничные условия (непрерывность касательных составляющих магнитного и электрического полей и нормальных составляющих магнитной и электрической индукции, см. § 5.1) должны быть дополнены теперь некоторыми условиями, накладываемыми на намагниченностьнаграницераздела. Они должны быть следствиями уравне ний движения, так как уравнения поля уже использованы при вы воде электродинамических граничных условий. Поскольку в уравнения движения входят теперь (см., например, выражения

(8.1.3)

и (8.2.3)) вторые производные переменной намагниченности

по

координатам, то можно

полагать, что искомые грапичпые

условия

будут содержать

составляющие намагниченности и

их

первые производные.

 

 

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

рассмотрим наиболее важный случай, когда магнитная среда пред­ ставляет собой ферромагнетик *) и граничит с немагнитной (ц0 = = 1) средой. Следуя Радо и Уиртмену [256], мы получим допол­ нительные граничные условия путем интегрирования уравнения

Ц Очевидно, что граничіше условия, которые мы получим, будут справед­ ливы и для низкочастотпой («ферромагнитной») ветви магнитных колебаний в ферримагнетике.

§ 8.3] ВОЛНЫ В ОГРА Н И ЧЕН Н Ы Х ТЕЛАХ И НЕОД Н О РО ДН Ы Х СРЕДАХ 407

движения аналогично тому, как электродинамические граничные условия получаются интегрированием уравнений Максвелла (см., например, [29]).

Считая для простоты неоднородное обменное взаимодействие

изотропным, т. е. тензор q в (2.1.7) скалярной величиной (qvs = = q), запишем исходное уравнение движения в виде

- і ® + М Х Н е„ + g M X Ѵ 2М + М X H s = 0. (8 .3 .5)

Здесь Herr включает в себя внешнее поле, размагничивакщее по­ ле и все эффективные поля (поле объемной анизотропии, эффектив­ ное поле, учитывающее диссипацию), не зависящие или слабо зави­ сящие от координат; поле Hs есть эффективное поле поверхност­ ной анизотропии, которое очень сильно зависит от координаты в направлении, перпендикулярном к граничной поверхности.

Поле Hs существует только в поверхностном слое толщиной d. Оно учитывает особые условия, в которых находятся магнитные моменты вблизи границы ферромагнетика. Не обсуждая подробно природу поверхностной анизотропии, заметим лишь, что источни­ ком ее, согласно Неелхо (см., например, [244]), может явиться то, что ионы у поверхности находятся в другом окружении, чем внутри ферромагнетика. Поверхностная анизотропия может воз­ никать также из-за наличия на поверхности тонкой пленки дру­ гого химического состава, чем внутри, например, антиферромаг­ нитной окисной пленки на границе ферромагнитного металла. Существенно, что интеграл от поля Hs по объему поверхностного слоя не является, вообще говоря, малым, несмотря на малую толщину слоя.

Поверхностная анизотропия может являться либо одноосной с осью — легкой или трудной, совпадающей с нормалью к поверх­ ности, либо однонаправленной (когда, в отличие от одноосной анизотропии, направления к поверхности и от нее не эквивалент­ ны). Рассмотрим для определенпости первый случай. Тогда, со­ гласно (2.2.5), ограничиваясь первой константой анизотропии, можно записать

H s = n 0 ^ j ( M n 0),

( 8 .3 . 6 )

где п0 — единичный вектор внешней нормали к граничной поверх­ ности S (рис. 8.3.1). Величина Hs зависит — причем, как уже от­ мечалось, сильно — от координаты z' в направлении, перпенди­ кулярном поверхности, и обращается в нуль вне поверхностного слоя.

Проинтегрируем уравнение (8.3.5) по объему V (рис. 8.3.1)— тонкому диску с площадью а и толщиной 2d1 (где d^, хотя и ма­

Рис. 8.3.1. К выводу граппчпых усло­ вий для намагниченности. S — граница между ферромагнетиком (z' < 0) и немагннтпой средой.

408

С П И Н О В Ы Е в о л н ы

[ Г Л . 8

ло, но больше, чем толщина поверхностного слоя d), «разрезанно­ му пополам» граничной поверхностью S. При достаточно малом dx первые два члена (8.3.5) дадут пренебрежимо малые вклады и результат интегрирования можно будет записать в виде

Іі + і2 = о,

(8.3.7)

І! = ЛVМ X Ѵ2М dV,

(8.3.8)

V

 

Ія = ( м х HsdV.

(8.3.9)

Применяя вторую формулу Грина [35] к проекциям (8.3.8) на оси декартовой системы коорди­ нат, нетрудно преобразовать Іх в пптеграл по поверхности 2 , ограничивающей объем диска У:

Іі = ? J М X (noaV)MdS.

(8.3.10)

Е

 

Здесь п0е — единичный

вектор

внешней нормали к 2 , а (поЕу)М = = дШ/дп-z — производная М по направлению этой нормали. Ин­ тегралом по боковой поверхности диска в (8.3.10) можно пренебречь. На внешнем (лежащем вне ферро­ магнетика) основании диска М —0, а на внутреннем noS = —п0. Тог­ да, считая площадь основания о

малой (хотя и большей, чем боковая поверхность диска), получим окончательно

І1 = — одМ X (п0Ѵ)М,

(8.3.11)

где (п0у)М — производная М по направлению нормали к гранич­ ной поверхности 5. Подчеркнем, что в (8.3.11) входит величина М внутри ферромагнетика — на расстоянии от границы, превышаю­ щем толщину поверхностного слоя.

Займемся теперь интегралом (8.3.9). Учтем, что Hs существует лишь в поверхностном слое и изменяется быстро только в направ­ лении z', а площадь основания диска а является малой. Тогда, принимая для Hs выражение (8.3.6), получим

§ 8.3] В О Л Н Ы В О Г Р А Н И Ч Е Н Н Ы Х Т Е Л А Х И Н Е О Д Н О Р О Д Н Ы Х С Р Е Д А Х

409

В (8.3.12) от z' может зависеть не только Hs, но и вектор М — его направление и длина. Если изменением М в поверхностном слое можно пренебречь, то (8.3.12) примет вид

І2 = ^ М X п0(Мп0),

(8.3.13)

где

 

и.

 

Ks = ^ H s c l z '

(8.3.14)

о

 

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

Можно рассмотреть и более общий случай, когда М изменя­ ется в поверхностном слое, однако таким «симметричным» обра­ зом, что направления вектора І2 и вектора М х н0 (вне поверхно­ стного слоя) совпадают. В этом случае выражение (8.3.13) можно по-прежнему считать справедливым, принимая в нем, как и в (8.3.11), величину М вне поверхностного слоя. Но теперь опреде­ лением поверхностной энергии K s будет уже не формула (8.3.14), а равенство выражений (8.3.12) и (8.3.13).

Подставляя (8.3.11) и (8.3.13) в (8.3.7), получим граничное

условие

 

 

gM X (п„Ѵ) М +

n0 X М(п0М) = 0.

(8.3.15)

Это условие накладывается на вектор М вблизи границы, но вне поверхностного слоя, а свойства поверхности учитываются един­ ственным параметром Ks-

Рассмотрим теперь случай малых колебаний намагниченно­

сти:

 

М = М0 + т е '“' , т < ^М 0

(8.3.16)

и получим граничное условие для комплексной амплитуды пере­ менной намагниченности т . Подставляя (8.3.16) в (8.3.15) и при­ равнивая нулю линейные по щчлены, найдем в случае М0 = const, т. е. отсутствия доменной структуры,

£z0 X (п0Ѵ) m + n0 X z0 (n0m) + (n0z0) n0 X m = 0, (8.3.17)

где z0 — единичный вектор в направлении М0, а

(8.3.18)

Выражение (8.3.17) эквивалентно двум скалярным условиям для составляющих вектора т . Направив оси координат, как по­ казано на рис. 8.3.1 (ось х касательиа к граничной поверхности),

Соседние файлы в папке книги из ГПНТБ