Скачиваний:
44
Добавлен:
10.05.2014
Размер:
1.75 Mб
Скачать

Рассмотрим область D={(x,y): y1(x) y y2(x), x [a,b]}, где y1(x), y2(x) – непрерывные функции на [a,b]. Области такого вида будем называть областями типа A. Области вида D={(x,y): x1(y) x x2(y),y [c,d]}, где x1(y), x2(y) – непрерывные функции на [c,d] называются областями типа B .

Теорема. Если для области типа A существуют ∫∫ f ( x, y)dxdy и для x [a,b]

 

 

D

y2 ( x )

b

y 2 ( x )

существует f ( x, y)dy , то существует dx

f ( x, y)dy и

y1 ( x)

a

y1 ( x )

by2 ( x )

dx

f ( x, y)dy = ∫∫ f ( x, y)dxdy .

a

y1 ( x)

D

Доказательство. Пусть D={(x,y): y1(x) y y2(x), x [a,b]}, где y1(x), y2(x) – непрерывные функции на

[a,b]. Рассмотрим функцию

f ( x, y), ( x, y) D

f *(x,y) =

,

 

0, ( x, y) R \ D

где R=[A,B] [C,E] прямоугольник, содержащий область D. Для функции f * выполнены условия предыдущей теоремы, поэтому

BE

dxf * ( x, y)dy = ∫∫ f * ( x, y)dxdy .

 

 

 

A

C

 

R

B

E

b

E

 

b

y 2 ( x )

Далее dxf * ( x, y)dy = dxf * ( x, y)dy = dx f ( x, y)dy . По теореме 3 из параграфа 4 выполнено

A C a C a y1 ( x )

равенство ∫∫ f * ( x, y)dxdy = ∫∫ f ( x, y)dxdy откуда и следует требуемое равенство. Аналогично

R D

доказывается

Теорема. Если для области типа B существуют ∫∫ f ( x, y)dxdy и y [c,d]

 

 

D

x2 ( y )

d x2 ( y )

существует f ( x, y)dx , то существует dy

f ( x, y)dx и

x1 ( y )

c

x1 ( y )

d

x2 ( y )

 

dy

f ( x, y)dx = ∫∫ f ( x, y)dxdy .

c

x1 ( y )

D

11

Примеры: Расставить пределы интегрирования в интеграле ∫∫ f ( x, y)dxdy в том и другом порядке.

 

D

1. D={(x,y):0 ≤ x ≤ 1, x2≤ y≤ 1+(x-1)2}

2.D={(x,y):0 ≤ x ≤ 1, x2-1≤ y≤ cos( π x ).

 

2

§7. Замена переменных в двойном интеграле

1. Отображение плоских областей. Криволинейные координаты.

Рассмотрим два экземпляра плоскости, плоскость переменных x, y и область D в этой плоскости, плоскость переменных ξ, η и область Σ в этой плоскости

Пусть имеется взаимно однозначное отображение области D на Σ

ξ = ξ ( x, y) ( x, y) D η = η( x, y)

x = x(ξ ,η)

= ξ η (ξ ,η) Σ y y( , )

(1),

(2).

Будем предполагать, что отображения (1), (2) непрерывно дифференцируемы и якобианы этих отображений

12

D(ξ ,η) 0, D( x, y) 0.

D( x, y) D(ξ ,η)

Отметим, что

D(ξ ,η) D( x, y) =1.

D( x, y) D(ξ ,η)

Вобласти Σ рассмотрим некоторую кусочно-гладкую кривую

ξ= ξ (t)

t [α,β].

η = η(t)

Ее образ в D имеет параметризацию

x = x(ξ (t),η(t))

y = y(ξ (t),η(t)) t [α,β]

и будет также кусочно-гладкой кривой. Действительно,

x' = ξx ξ '+ ηx η'

∂ ∂ (3). y' = ξy ξ '+ ηy η'

Если (ξ′,η′)≠(0,0), то и (x′,y′)≠(0,0). Если предположить противное, то система (3) с не вырожденной

матрицей коэффициентов должна будет иметь только тривиальное решение, что противоречит условию (ξ′,η′)≠(0,0).

Определение. Кривая, составленная из точек области D вида

x = x(ξ ,η0 )

x = x(ξ0 ,η)

y = y(ξ ,η0 ) или

y = y(ξ0 ,η)

называется координатной линией

Неявное задание этой линии имеет вид η(x,y)=η0 (соответственно ξ(x,y)=ξ0).

Определение. Числа ξ0 , η0 из области Σ плоскости (ξ , η) определяющие положение точки (x0 ,y0) из области D плоскости (x ,y) называются криволинейными координатами точки (x0 ,y0). Наоборот, на (x0 ,y0) можно смотреть, как на криволинейные координаты точки 0 , η0).

13

Фиксируя значения ξ или η на плоскости (ξ ,η) можно получить два семейства координатных линий. В области D появляется криволинейная координатная сетка. При сделанных предположениях две линии одного семейства не пересекаются между собой и, через любую точку области D проходит по одной линии из каждого семейства

2. Изменение площади при отображениях.

 

Рассмотрим отображение

 

ξ = ξ ( x, y) ( x, y) D и его обратное

x = x(ξ ,η) (ξ ,η) Σ

η = η( x, y)

y = y(ξ ,η)

удовлетворяющее условиям предыдущего пункта и разбиение области D линиями, порожденными линиями ξ=const, η=const плоскости ξ, η

Рассмотрим прямоугольник ξ, ξ +Δξ, , η, η +Δη в плоскости ξ, η и его образ в плоскости x, y.

Обозначим для краткости x=x(ξ,η), y=y(ξ,η), тогда

x(ξ +Δξ ,η)= x +

 

x

 

Δξ + o(ρ),

 

y(ξ +Δξ ,η)= y +

 

y

 

Δξ + o(ρ),

 

 

 

ξ

 

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(ξ ,η +Δη)= x +

 

x

Δη + o(ρ),

 

y(ξ ,η +Δη)= y +

 

y

Δη + o(ρ),

 

 

 

 

η

 

 

η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(ξ +Δξ , η +Δη)= x +

x

Δξ +

x

Δη + o(ρ), y(ξ +Δξ , η +Δη)= y +

y

Δξ +

y

Δη + o(ρ).

 

η

ξ

η

 

 

 

 

 

ξ

 

 

 

 

 

 

 

Для вычисления площади фигуры с вершинами

14

A(x,y), B(x(ξ +Δξ ,η), y(ξ +Δξ ,η)), C(x(ξ +Δξ , η +Δη), y(ξ +Δξ , η +Δη)), E( x(ξ ,η +Δη), y(ξ ,η +Δη))

рассмотрим параллелограмм A=A′, B′, C′, E′ с координатами вершин

A′=A=(x,y), B′=( x + ξx Δξ, y + ξy Δξ ), C′=( x + ξx Δξ + ηx Δη, y + ξy Δξ + ηy Δη ),

E′=( x + ηx Δη, y + ηy Δη ) .

Этот параллелограмм построен на векторах A′B′, A′E′,

a=A′B′ = (

x

Δξ,

y

Δξ), b=A′E′ = (

x

Δη,

y

 

Δη). Поэтому его площадь равна

ξ

ξ

 

 

 

η

 

 

 

 

η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

ξ

 

 

y

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[a,b] =

 

ξ

 

ξ

=ΔξΔη

D( x, y)

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

η

 

y

η

D(ξ ,η)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η

 

η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вершины A,A′, B,B′, C,C′, E,E′ отличаются на o(ρ). Можно показать, что в этом случае площади будут отличаться на o(ρ2)

µ(A′,B′,C′,E′)= D( x, y) ΔξΔη + o(ρ2).

D(ξ ,η)

Отсюда, в свою очередь, следует, что площадь области D будет равна

 

 

D( x, y)

 

 

 

µD= ∫∫dxdy = ∫∫

 

 

dξdη

(4).

D(ξ ,η)

D

Σ

 

 

 

 

 

 

 

 

Докажем последнее равенство для случая, когда область Σ представляет собой квадрат [α,β]×[α,β]

15

Разобьем Σ на равные части линиями ξ=ξi , η=ηj .

В этом случае Δξii+1 - ξi = (β - α)/n , Δηjj+1 - ηj = (β - α)/n , ρ= ξi2 + ηi2 =(β - α)/n,

n −1

 

D( x, y)

 

 

 

n −1

1

 

µD=

 

 

 

ξi

η j + o(

) .

D(ξ ,η)

 

 

2

i , j =0

 

(ξ

,η

)

i, j =0

n

 

 

 

i

j

 

 

 

 

Можно показать, что последнее слагаемое является бесконечно малой при n→∞ , откуда и следует равенство (4).

Замечание. Выражение dxdy иногда называют элементом площади в плоскости x,y, а выражение dξdη - элементом площади в плоскости ξ, η. Равенство (4) позволяет говорить, что модуль якобиана является коэффициентом искажения площади при данном отображении

 

 

 

 

 

dxdy =

D( x, y)

dξdη.

 

 

 

 

 

D(ξ ,η)

 

 

 

 

 

 

 

 

 

 

 

D( x, y)

 

µΣ следует, что в любой точке области M0=(ξ0 0 0 )

Из равенства µD=

 

D(ξ ,η)

 

 

(ξ ',η ')

 

 

 

 

 

 

 

 

 

 

 

 

 

D( x, y)

 

 

 

=

lim V .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D(ξ ,η)

 

 

 

 

 

 

 

 

(ξ

,η

Σ→M 0

µΣ

 

 

 

 

 

)

 

 

 

 

 

 

 

0

0

 

 

 

3.Примеры отображений.

 

 

 

 

 

 

 

 

 

 

 

Экспонента

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u = e

x

cos y

 

 

 

 

 

 

 

 

 

( x, y) Σ , Σ=[-3,1]×[0,π]

 

 

 

v = ex sin y

 

 

 

 

16

Функция Жуковского

 

 

 

 

 

 

 

 

u =

1

 

+

 

 

 

x

 

 

x

 

 

 

 

 

 

 

 

 

2

 

2

 

2

 

 

 

x

+ y

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

y

 

 

( x, y) Σ ,Σ=(

 

 

 

 

 

 

 

 

v =

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

x2 + y2

 

Дробно линейное отображение

 

 

 

 

u = 1 +

 

 

x

 

 

 

 

 

 

 

 

2

+ y

2

 

 

x

 

 

( x,

v = 1 −

 

 

 

y

 

 

 

 

+ y

 

 

 

x2

2

17

в полярных кординатах) [0,π]×[0.25,0.9]

y) Σ ,Σ= [0.25,1]×[0,1]

4. Замена переменных в двойном интеграле.

 

Рассмотрим отображение

 

ξ = ξ ( x, y) ( x, y) D и его обратное

x = x(ξ ,η) (ξ ,η) Σ ,

η = η( x, y)

y = y(ξ ,η)

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

Лемма. Если функция f(x, y) интегрируема на D, то функция F(ξ, η)=f(x(ξ, η),y(ξ, η)) интегрируема на Σ .

Доказательство. Любое разбиение области D порождает разбиение области Σ и наоборот. Таким образом связанные между собой разбиения, будем обозначать D={Dk} , Σ={Σk}. Здесь Dk и Σk – подобласти, переходящие друг в друга при заданном отображении .

Для разбиений D={Dk} , Σ={Σk} можно выписать соотношение между колебаниями функций

ωk(F)= sup

 

f ( x(M ), y(M )) − f ( x( N ), y( N ))

 

=

sup

 

f (M ') − f ( N ')

 

k(f).

 

 

 

 

M , N Σ k

 

∫∫

 

∫∫

 

D(ξ ,η)

 

 

 

 

M ', N ' Dk

 

 

 

 

Далее µΣ

k

=

 

dξdη =

 

 

dxdy CµD .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D( x, y)

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поэтому

 

 

 

Σk

 

Dk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S(F, Σ)-s(F,

 

Σ)= ωk (F )µΣk = ωk ( f )µΣk Сωk ( f )µDk =C(S(f, D)-s(f, D)).

 

 

 

 

 

k

 

 

 

k

k

Откуда и следует интегрируемость функции F(ξ, η) на Σ.

Теорема. Пусть функция f интегрируема в D, тогда

D( x, y)

∫∫D f ( x, y)dxdy = ∫∫Σ f [ x(ξ ,η), y(ξ ,η)] D(ξ ,η) dξdη .

Доказательство. Согласно доказанной лемме интеграл справа существует. Выберем некоторое разбиение области Σ на подобласти Σi и соответствующее ему разбиение области D на множества Di. Тогда по теореме о среднем для каждого i будет существовать точка i , ηi ), для которой

µDi =

∫∫

dxdy =

∫∫

D( x, y)

dξdη =

D( x, y)

 

 

 

µΣi .

D(ξ ,η)

D(ξ ,η )

 

 

 

 

 

 

 

 

Di

 

Σi

 

 

 

 

(ξ

,η

)

 

 

 

 

 

 

 

 

i

i

 

18

Для этих точек jj ) и соответствующих им точек (xj,yj ) можно выписать интегральные суммы

f ( x j , y j )µD j = f [ x(ξ j ,η j

), y(ξ j ,η j

)]

D( x, y)

 

 

µΣ j .

D(ξ ,η)

 

 

j

j

 

 

(ξ

,η

)

 

 

 

 

 

j

j

 

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

∫∫ f ( x, y)dxdy ,

∫∫ f [ x(ξ ,η), y(ξ ,η)]

D( x, y)

 

dξdη ,

D(ξ ,η)

D

Σ

 

 

 

 

соответственно.

Пример 1. Рассмотреть область D={ϕ [α,β],r [r1,r2]} и сделать замену в интеграле ∫∫ f ( x, y)dxdy ,

D

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

Пример 2. Сделать замену переменных u=x+y, v=y – x в интеграле ∫∫ f ( x, y)dxdy для области

D

D={|x|+|y|≤ 1}.

Пример 3 (3959). Сделать замену переменных x=u cos4v, y= u sin4v в интеграле ∫∫ f ( x, y)dxdy для

D

области D, ограниченной кривыми x=0, y=0, x + y = a .

19

Пример 4. Является ли конечной площадь области, заключенной между биссектрисой 2-4 –го координатных углов и кривой (x3+y3)2=x2+y2.

Перейдем к полярным координатам

r6(cos3ϕ+sin3ϕ)2=r2, r =

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos3

ϕ + sin3 ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos3 ϕ +sin 3 ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

1

4

 

 

 

 

 

 

 

 

 

 

 

 

1

4

 

 

 

 

 

µD= ∫∫dxdy =

 

 

rdr

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

=

 

 

 

 

 

 

 

 

 

3

ϕ + sin

3

ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

π

 

 

 

 

0

 

 

 

 

 

2

π

cos

 

 

2

π

(cosϕ + sin ϕ)(1 − cosϕ sin ϕ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

π

sin(ϕ + π )(2 − sin 2ϕ )

 

 

 

 

 

 

2

 

π

sin(ϕ + π )(2 + cos(2(ϕ + π )))

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

π

 

 

d cos(ϕ +

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

π

 

 

 

 

 

 

d cos( ϕ +

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

4

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

4

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

(1 − cos 2 (ϕ + π ))(1 + 2 cos

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

π sin 2 (ϕ + π )(2

+ cos(2(ϕ + π )))

 

 

2

 

π

 

(ϕ + π ))

 

 

 

 

 

 

 

 

 

 

4

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

4

 

4

 

 

 

 

 

1

 

 

−1

 

 

du

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. Последний интеграл расходится.

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 − u

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

1

)(1 +

2u

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Глава 2.

 

 

 

Кратные интегралы. Продолжение

 

 

 

 

§1. Тройные и n-кратные интегралы

1.Определение тройного и n-кратного интеграла

Пусть D кубируема, ее площадь будем обозначать µD , функция f(M)=f(x,y,z) определена и ограничена в области D. Предположим, что D разбита поверхностями на кубируемые подобласти Dk (совокупность {Dk} называется разбиением области D). В каждой из подобластей Dk выберем точку Mk=(ξkkk) Dk. Полученный набор точек обозначим Ξ ={Mk}. Интегральной суммой для набора f, разбиения , набора промежуточных точек Ξ называется выражение

σ ( f , , Ξ) = f (M k )µDk

(1)

k

 

20