Integral anizkoitz

testwikitik
Nabigaziora joan Bilaketara joan

 

Integrala bi kurben arteko azalera gisa.
Integral bikoitza, z=10x2y28 gainazalaren azpiko bolumen gisa. Gorputzaren behealdeko eremu angeluzuzena integralaren definizio-eremua da. Gainazala, integratu beharreko bi aldagaiko funtzioaren grafikoa da.

Matematikan (zehazki, aldagai anitzeko kalkuluan), integral anizkoitza aldagai erreal anitzeko funtzioen integrala da. Bi aldagaiko f(x,y) funtzioaren 2 planoko eremu bateko integralari integral bikoitz deritzo, eta hiru aldagaiko f(x,y,z) funtzioaren 3 espazioko eremu bateko integralari integral hirukoitz. Aldagai bakarreko integral anizkoitzerako, ikus "Integrazio errepikaturako Cauchy-ren formula".

Sarrera

Aldagai bateko funtzio positibo baten integralak funtzioaren grafikoaren eta OX ardatzaren arteko eskualdearen azalera adierazten du. Modu berean, bi aldagaiko funtzio positibo baten integral bikoitzak funtzioak definitutako gainazalaren eta funtzioaren definizio-eremua duen planoaren arteko eskualdearen bolumena adierazten du hiru dimentsioko plano kartesiarrean, non z=f(x,y) den. Funtzioak aldagai gehiago baditu, integral anizkoitzak dimentsio anitzeko funtzioaren hiperbolumena adierazten du.

n aldagaiko f(x1,x2,,xn) funtzioaren integral anizkoitza D eremuan adierazteko, integralaren ikurra hainbat aldiz jartzen da bata bestearen atzetik. Ordena bati jarraituz kalkulatzen dira integralak, ezkerreko integrala kalkulatzen azkena izanik. Ondoren, integratuko den funtzioa idazten da eta amaieran diferentzialen ikurrak, ordena egokian idatzita (eskuineko diferentzialari dagokion integrala kalkulatzen azkena izango da). Integrazio-eremua integrazio-ikur bakoitzean modu sinbolikoan adieraz daiteke, edo bestela modu laburtuan, eskuinean dagoen integral-ikurrean soilik:

Df(x1,x2,,xn)dx1dxn.

Jatorrizko funtzioaren kontzeptua aldagai erreal bakarreko funtzioetarako bakarrik definitzen denez, integral mugagabearen ohiko definizioa ez da zuzenean orokortzen integral anizkoitzera.

Definizio matematikoa

Izan bedi honela definitutako n-dimentsioko T tarte erdiirekia, n>1 izanik:

T=[a1,b1)×[a2,b2)××[an,bn)n.

[aj,bj) tarte bakoitzaren Pj partizio bat kalkula daiteke, hau da, haien artean disjuntuak diren eta guztien artean tarte osoa osatzen duten Ij azpitarteen familia finitu bat, non azpitarteak ezkerreko muturrean itxiak eta eskuinekoan irekiak diren.

Hortaz, T tartearen partizio bat

C=P1×P2××Pn

biderkadura kartesiarra da. C partizioa Ck izendatuko ditugun m azpitarteren familia bat bada, non Ck azpitarteak ez diren gainezartzen (disjuntuak dira) eta guztien bildura T den.

Har dezagun T tartearen C partizio bat,

T=C1C2Cm.

Izan bedi T tartean definitutako f:T funtzioa. Behetik n-dimentsioko T tarteak eta goitik f funtzioaren n-dimentsioko grafikoak bornatzen duten (n+1)-dimentsioko bolumenaren kalkulua Riemann-en batura honen bidez hurbil daiteke[1]:

k=1mf(pk)m(Ck),

non pk hori Ck azpitarteko puntu bat den eta m(Ck) hori Ck azpitartearen neurria den (ingelesetik, measure), hau da, azpitarteen luzeren biderkadura.

Ck azpitartearen diametroa Ck sortzen duten Ij tarteen luzerarik handiena da. T tartearen C partizio baten diametroa partizioko Ckazpitarteen diametrorik handiena da. Intuizioz, C partizioaren diaTxantiloi:Mvaretroa txikiagoa den heinean, azpitarteen m kopurua handitu egiten da, eta azpitarteen m(Ck) neurria txikitu egiten da. f funtzioa Riemann-en zentzuan integragarria dela esaten da, baldin

S=limδ0k=1mf(pk)m(Ck)

limitea existitzen bada, limitea diametroa gehienez δ duten T tartearen partizio posible guztietarako egiten delarik.[2]

f funtzioa Riemann-en zentzuan integragarria bada, T tartearen gaineko f funtzioaren Riemann-en integrala Txantiloi:Mvar dela esaten da, eta honela idazten da:

Tf(x1,x2,,xn)dx1dxn

Askotan, notazio hori laburtu egiten da, honela:

Tf(𝐱)dn𝐱.

non 𝐱=(x1,x2,,xn) n-kotea den eta dn𝐱 n-dimentsioko bolumen-diferentziala den.


Hemendik aurrera, Riemann-en n dimentsioko integralari integral anizkoitz deituko zaio[3].

Propietateak

Aldagai bakarreko funtzioen integraletan betetzen diren propietate asko integral anizkoitzetan ere betetzen dira (linealtasuna, trukakortasuna, monotonia-legea, etab.). Integral anizkoitzen propietate garrantzitsu bat honakoa da: integralaren balioa integrakizunen ordenaren independentea da, baldintza batzuk betetzen badira. Propietate horri Fubini-ren teorema deitzen zaio.

Kasu partikularrak

Txantiloi:Nowrapdenean,

l=Tf(x,y)dxdy

f funtzioaren integral bikoitza da T tartean, eta T3bada,

l=Tf(x,y,z)dxdydz

f funtzioaren integral hirukoitza da T tartean.

Integral bikoitza bi integral-ikurren bidez adierazten da eta integral hirukoitza hiru integral-ikurren bidez. Notaziorako akordio hori egokia gertatzen da integral anizkoitz bat integral berritu gisa kalkulatzen denean, artikulu honetan aurrerago azaltzen den bezala.

Integrazio-metodoak

Integral anizkoitzak ebazteko, gehienetan, integral berritu batera laburtzen da, hau da, aldagai bakarreko hainbat integraletan deskonposatzen da, haietako bakoitza zuzenean ebazteko modukoa delarik. Funtzio jarraituetan, hori Fubini-ren teoremaren bidez justifikatzen da. Batzuetan, zuzenean lor daiteke integral anizkoitzaren emaitza, inolako kalkulurik egin gabe.

Hauek dira integrazio-metodo sinple batzuk:[4][5]

Funtzio konstanteak integratzea

Integratu behar dena c funtzio konstantea denean, integrala c konstantearen eta integrazio-eremuaren neurriaren arteko biderkadura da. Baldin c=1 bada eta integrazio-eremua Txantiloi:Nowrap planoko eremu bat bada, integralak eremuaren azalera ematen du, eta integrazio-eremua Txantiloi:Nowrap espazioko eremua bada, integralak eremuaren bolumena ematen du.

Adibidea. Izan bitez

f(x,y)=2

funtzioa eta integrazio-eremu hau:

D={(x,y)2 : 2x4 ; 3y6}.

Hortaz,

3624 2 dxdy=23624 1 dxdy=2azalera(D)=2(23)=12,

definizioz

3624 1 dxdy=azalera(D)

baita.

Aldagai-aldaketa

Askotan, integrazio-mugak ez dira trukatzen errazak izaten. Aldagai-aldaketa egitean, integrala eremu "erosoago" batean idazten da eta, horrela, formula sinpleago baten bidez adieraztea lortzen da. Horretarako, funtzioa koordenatu berrietara egokitu behar da.

1a adibidea. Izan bedi

f(x,y)=(x1)2+y

funtzioa eta izan bitez

u=x1

eta

v=y

, hau da,

x=u+1

eta

y=v

. Aldagai-aldaketa eginez, honako funtzioa lortzen da:

f2(u,v)=u2+v

.

  • Integrazio-eremuarekin antzera egiten da, aldatuak izan diren jatorrizko aldagaiek mugatzen dutelako (adibidean, x eta y).
  • dx eta dy diferentzialak matrize jacobitarraren determinantearen balio absolutuaren bidez aldatzen dira. Bertan aldaketaren aldagai berriekiko deribatu partzialak daude (kontuan hartu, adibide gisa, transformazio diferentziala koordenatu polarretan).

Aldagai-aldaketen hiru "mota" nagusi daude (bat Txantiloi:Nowrap planoan, bi Txantiloi:Nowrap espazioan); hala ere, aldaketa orokorragoak egin daitezke printzipio bera erabiliz.

Koordenatu polarrak

Koordenatu kartesiarretatik koordenatu polarretara aldatzea.

Txantiloi:Nowrap planoan, integrazio-eremuak simetria zirkularra badu eta funtzioak ezaugarri berezi batzuk baditu, koordenatu kartesiarrak polar bihurtzeko aldaketa egin daiteke (ikus adibidea irudian); horrek esan nahi du, koordenatu kartesiarretako P(x,y) puntu orokorrak koordenatu polarretan dagozkien puntu bihurtuko direla. Aldaketa horrek integrazio-eremuaren forma aldatzeko eta eragiketak sinplifikatzeko aukera ematen du.

Aldaketa hori egiteko funtsezko erlazioa honako hau da:

x=ρcosφ,
y=ρsinφ.

2a adibidea. Izan bedi

f(x,y)=x+y

funtzioa. Koordenatu kartesiarrak polar bihurtzeko aldaketa aplikatuz, funtzioa honela geratzen da:

f(x,y)=f(ρcosφ,ρsinφ)=ρcosφ+ρsinφ=ρ(cosφ+sinφ).

2b adibidea. Izan bedi

f(x,y)=x2+y2

funtzioa. Kasu honetan, honela geratzen da:

f(x,y)=ρ2(cos2φ+sin2φ)=ρ2,

Pitagoras-en identitate trigonometrikoa erabiliz, eragiketa horiek asko sinplifikatzen dira.

Integrazio-eremuaren aldaketa egiteko,

x

eta

y

aldagaietatik abiatuz, erradioaren koroaren luzera eta deskribatutako angeluaren anplitudea definitzen dira

ρ

eta

φ

aldagaien tarteak definitzeko.

Integrazio-eremu baten aldaketaren adibidea, kartesiarretatik polarretara.

2c adibidea. Izan bedi honako integrazio-eremua (2ko erradioa duen zirkulua):

D={x2+y24}

Estaltzen duen angelua zirkuluaren angelua denez,

φ

aldagaia 0tik 2

π

-ra doa eta koroaren erradioa 0tik 2ra.

2d adibidea. Izan bedi honako integrazio-eremua:

D={x2+y29, x2+y24, y0},

hau da, koroa zirkularra y aldagaiaren planoerdi positiboan (ikus irudia 2c adibidean). φ aldagaiaren bidez planoaren angelua deskribatzen da eta ρ aldagaia 2tik 3ra doa. Hortaz, integrazio-eremua aldatu ondoren, honako laukizuzena geratzen da:

T={2ρ3, 0φπ}.

Aldaketaren determinante jacobitarra honakoa da:

(x,y)(ρ,φ)=|cosφρsinφsinφρcosφ|=ρ.

Determinante hori lortzeko, hau egin da: x=ρcosφ eta y=ρsinφ berdintzen ρ-rekiko deribatu partzialak lehen zutabean kokatu dira eta φ-rekiko deribatu partzialak bigarrenean. Horrela, aldaketarekin dxdy diferentziala ρdρdφ bihurtu da.

Funtzioa aldatu eta domeinua ebaluatu ondoren, koordenatu kartesiarretatik koordenatu polarretara aldatzeko formula defini daiteke:

Df(x,y)dxdy=Tf(ρcosφ,ρsinφ)ρdρdφ.

φ

aldagaia

[0,2π]

tartean dago eta

ρ

aldagaiak, aldiz, balio positiboak hartuko ditu, luzeraren neurri bat delako.

2e adibidea. Izan bedi

f(x,y)=x

funtzioa eta 2d adibideko integrazio-eremua.

D

eremuaren aurreko analisitik, tarteak ezagutzen ditugu:

ρ

(2tik 3ra) eta

φ

(0tik

π

-ra). Funtzioaren aldaketa honela egingo dugu:

f(x,y)=xf(ρ,φ)=ρcosφ.

Integrazioa aplikatuz, honela geratzen da:

Dxdxdy=Tρcosφρdρdφ.

Tarteak ezagunak izanik, honela geratzen da:

0π23ρ2cosφdρdφ=0πcosφ dφ[ρ33]23=[sinφ]0π (983)=0.


Koordenatu zilindrikoak

Koordenatu zilindrikoak.

Txantiloi:Nowrap espazioan, oinarri zirkularra duten integrazio-eremuetan integrala kalkulatzeko, koordenatu zilindrikoetara aldatzen da; funtzioaren aldaketa erlazio honen bidez egiten da:

x=ρcosφ,
y=ρsinφ,
z=z.

Eremuaren aldaketa grafikoki lortzen da, oinarriaren forma baino ez baita aldatzen, eta altuerak, berriz, hasierako eremuaren formari jarraitzen dio.

3a adibidea. Izan bedi honako eremua:

D={x2+y29, x2+y24, 0z5}

(hau da, 2d adibideko irudiko koroa oinarri duen eta 5eko altuera duen "hodi" bat); aldaketa aplikatzen bada, eremu hau lortzen da:

T={2ρ3, 0φ2π, 0z5}.

Aldaketarekin z aldagaia ez denez aldatzen, dxdydz diferentziala koordenatu polarretarako aldaketan bezala aldatzen dira; beraz, ρdρdφdz bihurtzen da.

Azkenik, koordenatu zilindrikoetarako formula aplika daiteke:

Df(x,y,z)dxdydz=Tf(ρcosφ,ρsinφ,z)ρdρdφdz.

Metodo hau egokia da integrazio-eremua zilindrikoa edo konikoa denean.

3b adibidea. Izan bedi

f(x,y)=x2+y2+z

funtzioa eta integrazio-eremua honako zilindroa:

D={x2+y29,  5z5}.

D eremua koordenatu zilindrikoetara aldatuta, honela geratzen da:

T={0ρ3, 0φ2π, 5z5}.

Funtzioa honela geratzen da:

f(ρcosφ,ρsinφ,z)=ρ2+z.

Hortaz, integrala honela geratzen da:

D(x2+y2+z)dxdydz=T(ρ2+z)ρdρdφdz;

Formula garatuz,

55dz02πdφ03(ρ3+ρz)dρ=2π55[ρ44+ρ2z2]03dz=2π55(814+92z)dz==405π.

Koordenatu esferikoak

Koordenatu esferikoak.

Txantiloi:Nowrap espazioan, zenbait integrazio-eremuk simetria esferikoa dute; beraz, integrazio-eremuko puntuen koordenatuak bi angeluren eta distantzia baten bidez adieraz daitezke. Orduan, koordenatu esferikoetara pasatzea egokia gertatzen da; funtzioa erlazio honen bidez aldatzen da:

x=ρcosθsinφ,
y=ρsinθsinφ,
z=ρcosφ.

OZ ardatzeko puntuek ez dute karakterizazio zehatzik koordenatu esferikoetan; beraz, θ aldagaia 0 eta 2π artean alda daiteke.

Aldaketa honetarako integrazio-eremurik egokiena esfera da.

4a adibidea. Izan bedi

D={x2+y2+z216}

(zentroa jatorrian duen 4 erradioko esfera). Aldaketa egin ondoren, honako eremua lortzen da:

T={0ρ4, 0φπ, 0θ2π}.

Aldaketa horren determinante jacobitarra honakoa da:

(x,y,z)(ρ,θ,φ)=|cosθsinφρsinθsinφρcosθcosφsinθsinφρcosθsinφρsinθcosφcosφ0ρsinφ|=ρ2sinφ.

Beraz, dxdydz diferentziala honela geratzen da: ρ2sinφdρdθdφ.

Integrazio-formula honela geratzen da:

Df(x,y,z)dxdydz=Tf(ρcosθsinφ,ρsinθsinφ,ρcosφ)ρ2sinφdρdθdφ.

Hobe da metodo hori erabiltzea integrazio-eremua esferikoa denean eta funtzioa erraz sinplifikatzen denean Txantiloi:Nowrap espaziora hedatutako trigonometriaren lehen oinarrizko erlazioen bidez (ikus 4b adibidea). Beste kasu batzuetan, koordenatu zilindrikoak erabiltzea egokiagoa izan daiteke (ikus 4c adibidea).

Tf(a,b,c)ρ2sinφdρdθdφ.

Formulan agertu diren ρ2 eta sinφ jacobitarretik datoz.

Hurrengo adibideetan

φ

eta

θ

aldagaiek jokatzen duten paperak trukatuta daude.

4b adibidea. Izan bedi 4a adibideko

D

eremua eta demagun

f(x,y)=x2+y2+z2

funtzioaren integrala kalkulatu behar dela. Aldaketa oso erraza da:

f(ρcosθsinφ,ρsinθsinφ,ρcosφ)=ρ2.

D eremua aldatuz lortutako T eremuaren tarteak ezagutzen ditugu:

T={0ρ4, 0φπ, 0θ2π}.

Beraz, integrazio-formula aplikatzen dugu:

D(x2+y2+z2)dxdydz=Tρ2ρ2sinθdρdθdφ,

eta, garatuz, honakoa lortzen da:

Tρ4sinθdρdθdφ=0πsinφdφ04ρ4dρ02πdθ=2π0πsinφ[ρ55]04dφ=2π[ρ55]04[cosφ]0π=4096π5.

4c adibidea. Izan bedi

D

eremua erdigunea jatorrian duen

3a

erradioko pilota,

D={x2+y2+z29a2},

eta izan bedi f(x,y)=x2+y2 integratu beharreko funtzioa.

Integrazio-eremua ikusita, koordenatu esferikoetara pasatzea egokia dela dirudi; T eremu berria mugatzen duten aldagaien tarteak hauek dira:

T={0ρ3a, 0φ2π, 0θπ}.

Hala ere, aldaketa eginda, funtzioa honela geratzen da:

f(x,y,z)=x2+y2ρ2sin2θcos2φ+ρ2sin2θsin2φ=ρ2sin2θ.

Integral hau lortzen da:

Tρ2sin2θρ2sinθdρdθdφ=Tρ4sin3θdρdθdφ.

Integral berrituen bidez kalkulatzen da.Tρ4sin3θdρdθdφ=03aρ4dρI0πsin3θdθII02πdφIII.

I=03aρ4dρ=ρ55|03a=2435a5;

II=0πsin3θdθ=0πsin2θd(cosθ)=0π(cos2θ1)d(cosθ)=cos3θ3|0πcosθ|0π=43;

III=02πdφ=2π.

Zati guztiak bilduz, honela geratzen da:Tρ4sin3θdρdθdφ=IIIIII=2435a5432π=6485πa5.

Beste modu batera egin daiteke, koordenatu zilindrikoetara pasatuz. Hauek dira T eremuaren tarte berriak:

T={0ρ3a, 0φ2π, 9a2ρ2z9a2ρ2};

z aldagaiaren tartea lortzeko, pilota bi hemisferiotan zatitu da, D eremuaren formulatik aldagaia askatuz (eta, orduan, x2+y2 zuzenean aldatuz ρ2 lortzeko). Funtzio berria ρ2 besterik ez da. Integrala honela geratzen da:

Tρ2ρdρdφdz.

Hau lortzen da:

02πdφ03aρ3dρ9a2ρ29a2ρ2dz=2π03a2ρ39a2ρ2dρ=2π9a20(9a2t)tdtt=9a2ρ2=2π09a2(9a2ttt)dt=2π(09a29a2tdt09a2ttdt)=2π[9a223t3225t52]09a2=227πa5(6185)=648π5a5.

Integral hirukoitza koordenatu zilindrikoetara aldatuz, oso erraz ebazten den aldagai bakarreko integral bat lortu da.

Adibideak

Integral bikoitza laukizuzen batean

Demagun aldagai anitzeko f funtzioaren integrala kalkulatu nahi dugula A eremu batean:

A={(x,y)𝐑2 : 11x14 ; 7y10} eta f(x,y)=x2+4y.

Honako integrala kalkulatuko dugu.

7101114(x2+4y)dxdy

Lehenik, barneko integrala egiten da, x-rekiko integratuz eta y konstante moduan hartuz, ez baita integrazio-aldagaia. Integral horren emaitza soilik y-ren mendeko funtzio bat da.

1114(x2+4y)dx=[13x3+4yx]x=11x=14=13(14)3+4y(14)13(11)34y(11)=471+12y.

Emaitza hori, ondoren, y-rekiko integratzen da:

710(471+12y) dy=[471y+6y2]y=7y=10=471(10)+6(10)2471(7)6(7)2=1719.

Funtzioaren balio absolutuaren integral bikoitza finitua den kasuetan, integrazio-ordena trukagarria da, hau da, lehenengo x-rekiko integratzeak edo lehenengo y-rekiko egiteak emaitza bera sortzen du (Fubiniren teorema). Adibidez, aurreko kalkulua alderantzizko ordenan eginez gero, emaitza bera lortzen da:

1114710(x2+4y)dydx=1114[x2y+2y2]y=7y=10dx=1114(3x2+102)dx=[x3+102x]x=11x=14=1719.

Integral bikoitza eremu erregular batean

Adibidea: integral bikoitza D eremu erregularrean.

Izan bedi honako eremua (ikus grafikoa):

D={(x,y)2 : x0,y1,yx2}.

Kalkula dezagun honako integrala:

D(x+y)dxdy.

D integrazio-eremua erregularra da OX eta OY ardatzekiko. Integrala kalkulatzeko, D eremua zehazten duten funtzioak eta funtzio horiek definituta dauden tarteak aurkitu behar dira. Kasu honetan, bi funtzioak honakoak dira:

α(x)=x2 eta β(x)=1.

Tarteak kalkulatzeko, funtzioen eta x=0 ardatzaren arteko ebaki-puntuak kalkulatu behar dira. Tartea [a,b]=[0,1] da (OX ardatzarekiko erregulartasuna aukeratu da, hobeto ulertzeko).

Integrala honela geratzen da:

D(x+y)dxdy=01dxx21(x+y)dy=01dx [xy+y22]x21

(hasteko, bigarren integrala kalkulatzen da, x konstante moduan hartuz). Gainerako eragiketak egiteko, integraziorako oinarrizko teknikak aplikatzen dira:

01[xy+y22]x21dx=01(x+12x3x42)dx==1320.

Erregulartasuna OY ardatzarekiko aukeratuz gero, integrala honela geratzen da:

01dy0y(x+y)dx

eta balio bera lortzen da.

OXY planoarekiko erregularra den Txantiloi:Nowrap espazioko eremu bat.

Bolumenen kalkulua

Aipatutako metodoak erabiliz, hainbat solido arrunten bolumenak kalkula daitezke.

  • Zilindroa: R erradioko oinarri zirkularra duen h altuerako zilindro baten bolumena kalkulatzeko, h funtzio konstantea integra daiteke integrazio-eremu gisa oinarri zirkularra hartuz eta koordenatu polarrak erabiliz.
Bolumena=02πdφ0Rhρdρ=2πh[ρ22]0R=πR2h.

Hori bat dator prisma baten bolumena kalkulatzeko formularekin:

Bolumena=oinarriaren azalera×altuera.
  • Esfera: R erradioa duen esfera baten bolumena kalkulatzeko, 1 funtzio konstantearen integrala kalkula daiteke integrazio-eremu gisa esfera hartuz eta koordenatu esferikoak erabiliz.
Bolumena=Df(x,y,z)dxdydz=D1dV=Sρ2sinφdρdθdφ=02πdθ0πsinφdφ0Rρ2dρ=2π0πsinφdφ0Rρ2dρ=2π0πsinφR33dφ=23πR3[cosφ]0π=43πR3.
  • Tetraedroa (piramide triangeluarra edo 3-simplex-a): goi-erpina jatorrian eta l luzerako albo-ertzak OX, OY eta OZ ardatzetan dituen tetraedroaren bolumena kalkula daiteke, 1 funtzio konstantea integratuz eta integrazio-eremu gisa tetraedroa hartuz.
Bolumena=0ldx0lxdy0lxydz=0ldx0lx(lxy)dy=0l(l22lx+x2(lx)22)dx=l3ll2+l33[l2x2lx22+x36]0l=l33l36=l36.
Hori bat dator piramide baten bolumena kalkulatzeko formularekin:
Bolumena=13×oinarriaren azalera×altuera=13×22×=36.
Integrazio-eremu bornegabe baten adibidea.

Integral inpropio anizkoitzak

Integrazio-eremua bornegabea denean edo integratu behar den funtzioa eremuaren mugan bornegabea denean, integral inpropio bikoitza edo hirukoitza erabili behar da.

Ikus, gainera

Erreferentziak

Txantiloi:Erreferentzia zerrenda

Kanpo estekak

Txantiloi:Autoritate kontrola