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

Конспект лекций Высшая математика (Басканова)

.pdf
Скачиваний:
151
Добавлен:
17.01.2018
Размер:
5.08 Mб
Скачать

469

dy1 = f

 

(x, y

, y

2

,..., y

n

),

 

dx

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy2

 

 

 

 

 

 

 

 

 

 

 

 

= f

 

(x, y

 

, y

 

 

,..., y

 

),

 

dx

2

 

2

 

n

 

 

1

 

 

 

 

 

 

 

.....................................

 

dyn

 

 

 

 

 

 

 

 

 

 

 

 

 

= f

 

(x, y

 

, y

 

 

,..., y

 

),

 

 

n

 

2

n

 

dx

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и поставлены начальные условия (4.8.7),

y1= y1(0), y2= y2(0),…, yn= yn(0) при x=x0.

Предположим, что функции, стоящие в правых частях системы (4.8.2), определены в некоторой замкнутой ограниченной области R:

x x a ,

y

k

y(0)

b,

(k=1,2,…,n)

0

 

k

 

 

с точкой (x0, y1(0), y2(0),..., yn(0)) внутри (a и b – заданные положительные числа) и удовлетворяют в этой области следующим двум условиям:

1. Функции fk(x, y1, y2,…,yn) (k=1,2,…,n) непрерывны по всем своим аргументам и, следовательно, ограничены, т. е.:

fk (x, y1, y2 ,..., yn ) M (k=1,2,…,n)

(4.8.8)

где M – постоянное положительное число, а (x, y1, y2,…,yn) – любая

точка, области R.

 

2. Функции fk(x, y1, y2,…,yn) имеют ограниченные частные

производные по аргументам y1, y2,…,yn, т.е.

 

fk (x, y1, y2 ,...yn ) K , (k,j = 1,2,…,n),

(4.8.9)

y j

 

где К – постоянное положительное число, а (x, y1, y2,…,yn) – любая точка области R.

При этих предположениях система (4.8.2) имеет единственное решение

(4.8.5).

y1=y1(x), y2=y2(x),…, yn=yn(x),

удовлетворяющее начальным условиям (4.8.7). Это решение заведомо

определено и непрерывно дифференцируемо в интервале

x x0 h

где

b

 

h = min a,

 

.

 

 

M

Определение 7. Совокупность n функций

470

y С= ϕС(x,

C,

2

,...,

n

),

 

1

1

1

 

 

 

),

 

y

С= ϕС(x,

C,

2

,...,

n

 

2

2

1

 

 

(4.8.10)

 

 

 

 

 

 

 

 

 

......................................

 

y

С= ϕС(x,

C,

2

,...,

n

),

 

n

n

1

 

 

 

определенных в некоторой области изменения переменных x, C1, C2,...,Cn имеющих непрерывные частные производные по x, будем называть общим решением системы (4.8.2) в области D, если система (4.8.10) разрешима относительно произвольных постоянных C1, C2,...,Cn в области D, так что при любых значениях x, y1, y2,..., yn, принадлежащих области D, системой (4.8.10) определяются значения C1, C2,...,Cn, и если совокупность n функций (4.8.10) является решением системы (4.8.2) при всех значениях произвольных постоянных C1, C2,...,Cn, когда точка (x, y1, y2,..., yn) пробегает область D.

Определение 8. Если решение системы (4.8.2) состоит только из точек единственности решения задачи Коши для этой системы, то такое решение называется частным решением.

Решение, получающееся из формулы общего решения при частных числовых значениях произвольных постоянных C1, C2,...,Cn, будет, очевидно, частным решением.

Определение 9. Решение системы (4.8.2), в каждой точке которого нарушается единственность решения задачи Коши для этой системы,

называется особым решением.

4.8.2. Линейные системы дифференциальных уравнений

4.8.2.1. Однородные линейные системы

Рассмотрим один специальный класс нормальных систем дифференциальных уравнений – линейные системы дифференциальных уравнений. Это системы вида:

dy1 = p

(x)y + p

(x)y

2

+ . .+.p

(x)y

n

+ f

(x),

 

 

dx

11

1

12

 

1n

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy2

 

 

 

 

 

 

 

 

 

 

 

 

 

= p21

(x)y1 + p22

(x)y2

+ . .+.p2n (x)yn + f2 (x),

 

 

dx

 

(4.8.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

..............................................................................

 

 

dyn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= p

(x)y + p

 

(x)y

 

+ . .+.p

(x)y

 

+ f

 

(x),

 

 

 

n2

2

n

n

 

 

dx

n1

1

 

nn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Будем предполагать, что в системе (4.8.11) коэффициенты pkj(x) (k=1,2,…n) и функции fk(x) (k=1,2,…n) непрерывны в интервале (a, b). Тогда,

471

согласно теореме Пикара система (4.8.11) имеет единственное решение задачи Коши для этой системы, и это решение будет определено во всем интервале непрерывности функций pkj(x) и fk(x).

Особых решений линейная система (4.8.11) не имеет. Всякое решение этой системы является частным решением.

Если все функции fk(x)0 в (a, b), то система (4.8.11) называется однородной. В этом случае система принимает вид:

dyk

n

 

= pkj (x)y j , (k=1,2,…n).

(4.8.12)

dx

j=1

 

Если же в системе (4.8.11) не все функции fk(x)0 в (a, b), то такая система называется неоднородной.

Общая теория линейных систем во многом аналогична общей теории линейного уравнения n-го порядка. Наша окончательная задача состоит в нахождении всех вещественных решений системы (4.8.12). Однако для реш е- ния этой задачи так же, как и в случае однородного линейного уравнения n- го порядка введем некоторые понятие.

Решения однородной системы (4.8.12) обладают следующими характерными свойствами, аналогичными свойствам решений однородного линейного уравнения n-гo порядка.

1. Если

y1=ϕ1(x), y2=ϕ2(x),…, yn=ϕn(x)

(4.8.13)

есть решение однородной системы (4.8.11), то

 

y1=Cϕ1(x), y2=Cϕ2(x),…, yn=Cϕn(x)

(4.8.14)

где C – произвольная постоянная, тоже будет решением этой системы,

т. е. умножая все функции, составляющие решение однородной системы, на одну и туже постоянную мы снова получаем решение.

2. Пусть дано m решений системы (4.8.12), записанных в виде таблицы:

y

,

y

,

L

y

,

 

 

11

 

 

 

12

 

 

 

 

1n

 

 

y21,

y22,

L y2n ,

(4.8.15)

L

 

 

L

 

 

L

L

 

 

 

 

 

 

 

y

m1

,

y

m2

,

L

y

 

.

 

 

 

 

 

 

 

mn

 

Здесь первый индекс обозначает номер решения, а второй – номер функции. Например, y12 – вторая функция первого решения.

Из линейности системы (4.8.12) следует, что линейная комбинация решений (4.8.15) c любыми постоянными коэффициентами C1, C2,…,Cm, т. е. совокупность функций

n

 

yk = Ci yik (k =1,2,...n) ,

(4.8.16)

i=1

тоже будет решением системы (4.8.12).

472

Предположим теперь, что нам известно n частных решений однородной системы (4.8.12). Поставим основной вопрос: при каком условии линейная комбинация этих решений с произвольными постоянными коэффициентами C1,C2,...,Cn даст общее решение однородной системы. Чтобы ответить на поставленный вопрос, введем понятие о линейной независимости систем функций.

Определение 10. Рассмотрим m систем функций:

y

,

y

,

L

y

,

 

 

11

 

 

 

12

 

 

 

 

1n

 

 

y21,

y22,

L y2n ,

(4.8.17)

L

 

 

L

 

 

L

L

 

 

 

 

 

 

 

y

m1

,

y

m2

,

L

y

 

.

 

 

 

 

 

 

 

mn

 

Эти системы называются линейно независимыми в интервале, (a, b),

если не существует чисел α1, α2,…, αm, не равных одновременно нулю, при которых для всего интервала (a, b) выполнялись бы соотношения:

m

 

αi yik 0 (k=1,2,…n),

(4.8.18)

i=1

впротивном случае системы (4.8.17) называются линейно зависимыми

в(a, b).

Пусть мы имеем n систем функций:

y ,

y , ...

y

,

 

 

 

11

 

12

1n

 

 

 

 

y21,

y22 , ...

y2n ,

 

 

(4.8.19)

...

 

... ...

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yn1, yn2 , ... ynn .

 

 

 

Введем в рассмотрение определитель

 

 

 

 

 

 

y11,

y12 , ... y1n

 

 

 

 

 

 

 

W (x) =

 

y21,

y22 , ...

y2n

.

(4.8.20)

 

...

... ... ...

 

 

 

 

 

 

 

 

yn1,

yn2 ,

...

ynn

 

 

Этот определитель называется определителем Вронского для систем функций (4.8.19) или вронскианом этих систем функций.

Теорема. Если n систем функций (4.8.19) линейно зависимы в интерва-

ле (a, b), то W(x) = 0 в (a, b).

Пусть теперь каждая из систем функций (4.8.19) является решением системы (4.8.12).

Теорема. Если n решений (4.8.19) линейно независимы в интервале (a, b), в котором определены и непрерывны Pkj(x) (k=1,2,…n) то их вронскиан не обращается в нуль ни в одной точке этого интервала.

Из этой теоремы следует, что для того, чтобы n решений системы

473

(4.8.12) были линейно независимы в интервале (a, b), необходимо и достаточно, чтобы их вронскиан не обращался и нуль ни в одной точке этого интервала.

Однако для установления линейной независимости n решений системы (4.8.12) достаточно убедиться, что W(x) отличен от нуля хоть в одной точке интервала (a, b). Это вытекает из следующих двух замечательных свойств вронскиана n решений системы (4.8.12).

1.Если W(x) обращается в нуль хоть в одной точке интервала (a, b), т. е. интервала непрерывности коэффициентов системы (4.8.12), то W(x) равен нулю во всех точках этого интервала.

2.Если W(x) не равен нулю хоть в одной точке интервала, (a, b); то он не обращается в нуль ни в одной точке интервала (a, b).

Таким образом, для линейной независимости n решений системы (4.8.12) в интервале (a, b) необходимо и достаточно, чтобы их вронскиан был отличен от нуля хоть в одной точке этого интервала.

Определение 11. Совокупность n решений однородной системы (4.8.12), определенных и линейно независимых в интервале (a, b), называется

фундаментальной системой решений в этом интервале.

Так же, как и в случае однородного линейного уравнения n-го порядка, знание фундаментальной системы решений дает возможность построить общее решение системы (4.8.12).

Основная теорема. Если

y11,

y12 , ... y1n ,

 

y21,

y22 , ...

y2n ,

(4.8.21)

 

 

 

 

...

...

...

...

 

yn1,

yn2 , ...

 

 

ynn .

 

есть фундаментальная система решений однородной линейной системы (4.8.12) в интервале (a, b), то формулы

n

 

yk = Ci yik (k =1,2,...n) ,

(4.8.22)

i=1

где C1, C2,...,Cn – произвольные постоянные, дают общее решение системы (4.8.12) в области

a<x<b, |yk|< (k=1, 2,…,n),

(4.8.23)

т. е. во всей области задания системы (4.8.11).

4.8.2.2. Неоднородные линейные системы

Рассмотрим теперь неоднородную систему

474

 

 

 

 

dyk

n

 

 

 

 

 

 

 

= pkj (x)y j + fk (x) (k=1,2,…n),

(4.8.25)

 

 

 

 

dx

j=1

 

 

Предположим, что нам известно некоторое частное решение этой сис-

темы:

 

 

y1= y1(1), y2= y2(1),…, yn= yn(1)

(4.8.26)

 

 

 

 

так что мы имеем тождества

 

 

 

 

 

 

(1)

n

 

 

 

 

 

 

dyk

pkj (x)y(1)j

+ fk (x) (k=1,2,…n).

(4.8.27)

 

 

 

 

dx

j=1

 

 

Введем новые неизвестные функции z1, z2,...,zn по формулам

 

 

 

 

 

 

yk = yk(1) + zk

(k =1,2,...n)

(4.8.28)

Подставляя функции (4.8.28) в неоднородную систему (4.8.25), полу-

чим:

 

 

 

 

 

 

 

(1)

 

dzk

n

 

n

 

 

 

dyk

+

= pkj (x)y(1)j + pkj (x)z j + fk (x) (k=1,2,…n), (4.8.29)

 

dx

dx

 

 

j=1

 

j=1

 

 

Отсюда, в силу тождеств (4.8.27), получаем для функций z1, z2,...,zn

следующую однородную систему, дифференциальных уравнений:

 

 

 

 

 

dzk

n

 

 

 

 

 

 

 

= pkj (x)z j

(k=1,2,…n),

(4.8.30)

 

 

 

 

dx

j=1

 

 

Эта система называется однородной системой, соответствующей

неоднородной системе (4.8.25).

 

 

Общее решение однородной системы (4.8.30) дается формулой

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

zk

= Ci zik

(k=1,2,…n)

(4.8.31)

i=1

где {zik} – некоторая фундаментальная система решений этой однородной системы.

Подставляя (4.8.31) в (4.8.28), получаем:

 

n

 

 

yk = yk(1)

+ Ci zik

(k =1,2,...n)

(4.8.32)

i=1

Можно показать, что все решения системы (4.8.25) содержатся в формуле (4.8.32). Эта формула представляет собою общее решение системы (4.8.25) во всей области задания системы (4.8.25).

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

(4.8.30).

475

4.8.2.3. Метод вариации произвольных постоянных (метод Лагранжа)

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

Теорема. Если известна фундаментальная система решений однородной системы (4.8.30), то общее решение неоднородной системы (4.8.25) может быть найдено при помощи квадратур.

Доказательство. Для доказательства этой теоремы применим метод вариации произвольных постоянных. Будем искать решение неоднородной системы (4.8.25) в виде

n

 

yk = Ci (x)zik (k =1,2,...n)

(4.8.33)

i=1

где {zik} – фундаментальная система решений однородной системы (4.8.25), а Ci(x) (i=l,2,...,n) – некоторые непрерывно дифференцируемые функции от x. Выберем эти функции так, чтобы формула (4.8.33) давала решение системы (4.8.25).

Подставляя (4.8.33) в (4.8.25), получаем:

n

n

n

n

Ci

(x)zik + Ci (x)zik = pkj (x)Ci (x)zij + fk (x)

i=1

 

i=1

 

j=1

i=1

(k =1,2,...n)

или

n

n

n

n

Ci

(x)zik + Ci (x)zik = Ci (x)pkj (x)zij + fk (x)

i=1

 

i=1

 

i=1

j=1

(k =1,2,...n)

Переписав эти равенства в виде

(4.8.34)

(4.8.35)

n

n

n

Ci

(x)zik + Ci (x)[zik pkj (x)zij ] = fk (x) (k =1,2,...n) (4.8.36)

i=1

 

i=1

 

j=1

и приняв во внимание, что {zik} – фундаментальная система решений однородной системы (4.8.30), мы приходим к следующей системе n уравнений для определения Ci (i = 1, 2,...,n):

n

 

 

Ci(x)zik = fk (x)

(k =1,2,...n)

(4.8.37)

i=1

Так как определитель этой системы, будучи равным W(x), отличен от нуля во всем интервале (a, b), то разрешая ее относительно Ci(x), находим:

476

n

Wki (x)

 

Ci(x) = fk (x)

(i =1,2,...n)

W (x)

k=1

 

где Wki(x) есть алгебраическое дополнение элемента zki W(x). Интегрируя (4.8.38), находим:

(4.8.38)

вронскиана

n

x

 

 

Ci (x) =

fk (x)Wki (x) dx + Ci

(i =1,2,...n) .

k=1 x

W (x)

 

 

0

 

 

Подставляя эти значения Ci(x) в формулу (4.8.33), получаем:

n

n x

yk = zik ∑∫ fs

i=1

s=1 x

 

0

Полагая здесь C1 = C2

n yk(1) = zik

i=1

 

n

(x)Wsi (x) dx + Ci zik (k =1,2,...n).

W (x)

i=1

=...= Cn = 0, получаем частное решение:

n x

 

fs (x)Wsi (x) dx

(k =1,2,...n),

W (x)

 

s=1 x0

 

(4.8.39)

(4.8.40)

(4.8.41)

так что (4.8.40) можно записать в виде (4.8.32) и, следовательно, решение, определяемое формулой (4.8.40), является общим решением неоднородной системы (4.8.25) в области (4.8.23). Теорема доказана.

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

4.8.3. Линейные системы дифференциальных уравнении с постоянными коэффициентами

Рассмотрим линейную систему уравнений:

dyk

n

 

= akj y j + fk (x) (k=1,2,…n),

(4.8.42)

dx

j=1

 

где коэффициенты akj (k, j=1,2,...,n) – постоянные вещественные числа, a fk(x) (k=l, 2,...,n) – функции от x, непрерывные в интервале (a, b).

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

dyk

n

 

 

= akj y j

(k=1,2,…n),

(4.8.43)

dx

j=1

 

 

Для построения общего решения системы (4.8.43) достаточно построить хоть одну фундаментальную систему решений.

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

эффициентами будем искать частное решение системы (4.8.43) в виде

 

y1=γ1eλx, y2=γ2eλx,…,yn=γneλx,

(4.8.44)

477

где γ1, γ2,…, γn и λ – некоторые постоянные числа, причем числа γ1,

γ2,…, γn не равны нулю одновременно, ибо в противном случае мы получили бы очевидное нулевое решение, которое не может входить в состав фундаментальной системы и, следовательно, не может быть использовано для построения общего решения.

Обратим особое внимание на то, что число λ мы берем одно и то же для всех функций, составляющих решение.

Подставляя функции (4.8.44) в систему (4.8.43), сокращая на eλx и пере-

нося все члены направо, получим для определения чисел γ1, γ2,…, γn следующую систему:

(a11 − λ)γ1 + a12γ2 +... + a1nγn = 0,

 

a21γ1 + (a22 − λ)γ2 +... + a2nγn

= 0,

(4.8.45)

 

 

 

 

 

 

 

 

 

 

....................................................

 

a

γ + a

n2

γ

2

+... + (a

nn

− λ)γ

n

= 0.

 

n1

1

 

 

 

 

 

Нас интересует ненулевое решение этой системы. Такое решение существует лишь при условии, что определитель системы равен нулю, т. е. при условии

 

a11 − λ

a12

L

a1n

 

 

 

 

 

(λ)

a21

a22 − λ

L

a2n

= 0

(4.8.46)

 

L

L

L

L

 

 

 

an1

an2

L ann − λ

 

 

Уравнение (4.8.46) называется характеристическим уравнением сис-

темы (4.8.43), его корни – характеристическими числами, а определитель (λ) – характеристическим определителем.

Рассмотрим сначала случай, когда все характеристические числа λ1,

λ2,...,λn различны и вещественные. В этом случае имеем: (λi) = 0, но ∆′(λi)

≠ 0 (i=1,2,…n). Вследствие этого ранг матрица составленной из коэффициентов системы

(a11 − λi )γ1 + a12γ2

+... + a1nγn

= 0,

 

a21γ1 + (a22 − λi )γ2

 

 

 

 

 

 

 

+... + a2nγn = 0,

(4.8.47)

....................................................

 

 

 

a γ + a

n2

γ

2

+... +

(a

− λ

)γ

n

= 0,

 

 

n1 1

 

 

nn

i

 

 

 

 

которая получается из системы (4.8.45) после замены в ней λ на λi равен n 1.

Поэтому одно из уравнений системы (4.8.47) есть следствие остальных и эта система имеет ненулевое решение, определенное с точностью до произ-

вольного множителя пропорциональности Ai:

 

 

γi1=Aimi1, γi2=Aimi2,…, γin=Aimin

(i=1,2,…n).

(4.8.48)

478

Например, в качестве γik можно взять алгебраические дополнения эле-

ментов любой строки определителя (λi), если не все они равны нулю. Фи к- сируя в формулах (4.8.48) множитель Ai мы получим определенное решение системы (4.8.47).

Подставляв теперь в (4.8.44) вместо λ последовательно характеристи-

ческие числа λi, а вместо γ1, γ2,…, γn – соответствующие им решения системы (4.8.47), определяемые формулами (4.8.48) при фиксированных множителях Ai, получим n решений:

y = γ eλ1x ,

y = γ eλ1x ,

L

y

= γ

 

eλ1x;

 

 

11

 

11

 

12

 

12

 

 

 

1n

1n

 

 

 

y

 

= γ

eλ2x ,

y

 

= γ

 

eλ2x , L

y

 

= γ

 

 

 

21

22

22

2n

2n

eλ2x;

 

 

 

21

 

 

 

 

 

 

 

(4.8.49)

 

 

L

 

 

 

L

 

 

L

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

n1

= γ

eλnx ,

y

n2

= γ

n2

eλnx ,

L

y

nn

= γ

nn

eλnx;

 

 

 

n1

 

 

 

 

 

 

 

 

 

Эти решения линейно независимы в интервале (−∞, + ).

Если при этом все корни λ1, λ2,…,λn вещественны, то все решения (4.8.49) тоже будут вещественными.

Таким образом, в случае различных вещественных корней характеристического уравнения система (4.8.43) имеет n вещественных линейно независимых частных решений вида (4.8.49), так что последние образуют фундаментальную систему решений.

Поэтому, формулы

y

= C

γ eλ1x

+ C

γ

 

eλ2x

+... + C

γ

 

eλnx ,

 

 

 

1

1

11

 

2

 

 

21

 

 

 

 

n

 

 

n1

 

 

 

 

 

 

y

 

 

= C

γ

 

λ x

+ C

γ

 

 

e

λ

x

+ + C

γ

 

 

e

λ

x

 

2

e 1

22

2

 

n2

n

 

,

 

 

1

12

 

2

 

 

 

 

 

 

n

 

 

 

 

 

 

 

(4.8.50)

..................................................................

 

y

n

= C

γ

eλ1x + C

 

γ

2n

eλ2x +... + C

 

γ

nn

eλnx

 

 

 

1

1n

 

2

 

 

 

 

 

n

 

 

 

 

 

 

 

дают общее решение системы (4.8.43) в области

x

 

< +∞,

 

y1

 

< +∞,

 

y2

 

< +∞,...,

 

yn

 

< +∞.

(4.8.51)

 

 

 

 

 

 

 

Если характеристические числа различные, но среди них есть комплексные, то последние входят сопряженными парами. Пусть a + ib и aib – простые корни характеристического уравнения. Корню a + ib соответствует решение

y

 

= eax (γ

cosbx − γ

21

sinbx),

 

y

 

 

= eax (γ

cosbx − γ

22

sinbx),

11

11

 

 

 

 

 

12

 

12

 

 

 

 

 

...,y

= eax (γ

cosbx − γ

 

 

 

sinbx);

 

 

 

 

 

 

2n

 

 

 

 

 

 

 

 

1n

 

1n

 

 

 

 

 

 

 

 

 

 

y

 

= eax (γ

sinbx + γ

 

cosbx),

y

 

 

 

= eax (

γ

sinbx + γ

 

 

 

21

21

22

22

cosbx),

 

11

 

 

 

 

 

 

12

 

 

 

 

 

...,y

 

= eax (γ

sinbx + γ

 

 

 

cosbx).

 

 

 

 

 

 

 

2n

2n

 

 

 

 

 

 

 

 

 

 

1n

 

 

 

 

 

 

 

 

 

(4.8.52)