Kapitola 1 Neurčitý integrál 1.1 Primitivní funkce a neurčitý integrál V kurzu Matematická analýza 1 jsme se seznámili s pojmem derivace, naučili se derivovat některé elementární funkce a řekli si též některé aplikace derivací. Jednou z nich byla i aplikace fyzikální, a to v(t) = s'(t)—slovy okamžitá rychlost je první derivací dráhy podle času. Je-li nám známo jak se mění dráha pohybujícího se tělesa v závislosti na čase, snadno spočítáme, jakou rychlostí se hmotný bod v daném čase pohybuje. Co když však známe závislost okamžité rychlosti na čase a neznáme dráhu? Z fyziky víme, že u pohybu rovnoměrně zrychleného platí pro rychlost vztah v = at, kde a je zrychlení; pro tento typ pohybu se jedná o konstantu. Víme, že dráha pohybu rovnoměrně zrychleného je dána vztahem s = ^at2 + s0, kde s0 je konstanta, jejíž fyzikální význam je dráha v čase t = 0. Snadno se přesvědčíme, že s'(t) = v(t) = at. Vím, že vám aplikace moc pod fousy nelezou, proto se abstrouhám od pohybu a problém postavím čistě matematicky. Otázka zní, zda např. funkce / : y = x není derivací nějaké funkce. Bez velkých problémů zjistíme, že kýžené x obdržíme derivací funkce F : y = \x2. Je zřejmé, že tato nalezená funkce není jediná možná, problém řeší jakákoliv funkce F : y = ^x2 + c, kde c je konstanta. Zeptám-li se, kterou funkci musím zderivovat, abych dostal cosrr, tak mi asi každý řekne, že je to y = sin x + c. Avšak ruku na srdce, kdo je mi schopen během jedné minuty říci, že potřebuji derivovat funkci F : y = |rr — \ sin2x + ^ sin4x + c, abych obdržel poměrně jednoduchou složenou funkce / : y = sin4 x? Tušíme, že tyto výpočty budou složitější než derivování, které v jednodušších případech zvládá i medvěd Brumla mého přítele RK, takto domptéra v jednom nejmenovaném cirkuse. My se však nenecháme odradit a jsouce povzbuzeni slovy básníka Jana Nerudy že malý ten kdo má jen malý cíl shrneme dosavadní poznatky do první definice. Def. 1.1 Nechí F(x) a f(x) jsou funkce definované v otevřeném intervalu J. Funkci F{x) nazýváme funkcí primitivní k funkci f{x), jestliže pro všechna x E I platí F'(x) = f(x). Z úvodních poznámek nám plyne rovněž naše první věta, kterou si jistě sami dokážete coby cvičení. 1 2 KAPITOLA 1. NEURČITÝ INTEGRÁL Věta 1.1 Nechi F(x) je primitivní funkcí k funkci f (x). Pak i libovolná funkce G (x) = F (x) + c, kde c G M je libovolná konstanta, je funkcí primitivní k funkci f (x). K této větě uvedu dvě poznámky. Graf primitivní funkce G(x) vznikne posunutím grafu funkce F(x) o c ve směru osy x. Mezi množinou všech primitivních funkcí k f (x) a množinou IR existuje bijekce. V předdmětu Teorie množin se pak dozvíte, že tyto množiny mají stejnou mohutnost, ale to poněkud předbíhám. Je namístě situaci z předchozí věty nějak pojmenovat. Def. 1.2 Množinu všech funkcí primitivních k funkci f{x) nazýváme neurčitým integrálem funkce f{x) v intervalu I. Píšeme /í(l)iI = ř(l)+r' xeI Funkci f{x) nazýváme integrandem, symbol j integračním znakem. Celý proces budeme nazývat integrováním (integrací) funkce f{x). Záhadný symbol dx tam zatím pište a nepřemýšlejte o něm. Jakmile se probojujeme k integrálu určitému, tak jeho význam ozřejmíme. Otázku, zda uvedený proces má smysl řeší následující věta. Věta 1.2 Ke každé funkci spojité na intervalu I existuje v tomto intervalu funkce primitivní. Pusťme se tedy chutě do hledání funkcí primitivních, tedy do integrování. Někteří národové čtou zprava doleva, vezmeme-li si z nich příklad a přečteme-li vzorce pro derivování tímto směrem, dostaneme sadu základních vzorců pro neurčitý integrál. 1.1. j O.dx = c 1.2. j dx = x + c 1.3. fxndx = ^ + c, ri^-l J n+l ' ' 1.4. f -dx = ln Irrl + c J X i i 1 1.5. j exdx = ex + c 1.6. J cŕdx = if^ + c, o, > 1, a ^ 1 1.7. j sin xdx = — cos x + c 1.8. j cos xdx = sin x + c 1.9. f —\-dx = tsx + c 1.10. f .\ dx = — cotg x + c J sir x ° 1.11. f ,} , dx = arcsin x + c J yl—xÁ 1.1. PRIMITIVNÍ FUNKCE A NEURČITÝ INTEGRÁL 3 1.12. j j^dx = arctgrr + c U všech vzorců předpokládáme, že platí pro všechna přípustná x. Sami cítíme, že podle uvedených vzorců bychom moc integrálů nevyřešili. Zkusíme tedy pátrat, zda bychom nemohli obrátit i jiné vzorce a věty z derivací. V případě lineární kombinace funkcí budeme úspěšní, neb ta se integruje na stejném principu jako se i derivuje, jak nám to říká následující věta. Věta 1.3 Nechť v intervalu I existuji neurčité integrály funkcí f\(x), f2(2), ■ ■ ■, fn{x) a nechť c1; c2,..., cn jsou libovolné konstanty. Pak existuje i neurčitý integrál funkce f(x) = c1f1(x) + c2f2(x) H-----h cnfn(x) a platí cifi(x) + c2f2(x)-\-----h cnfn(x)dx = ci / f1(x)dx+c2 / f2{x)dx-\-----hcn / fn(x)dx Některé integrály můžeme jednoduchou úpravou převést na tzv. integrály tabulkové, jak se o tom přesvědčíme v následujícím příkladu. Příklad 1.1 2 , i sin2 x f 1 — cos2 x i, ± 1 tg xdx = /--—dx = /---dx = /----1 \ dx cos2 x I cos2 x /V cos2 x Nyní s využitím vzorce 1. 9. a vety 1. 3. máme 1 J dx = ľ--—dx — I ldx = tgx — x + c cos2 X J J cos2 X Bohužel výše uvedená je věta je naším posledním úspěchem v tomto směru. Na rozdíl od derivací neexistují vzorce pro derivaci součinu a podílu. Malou útěchou nám budiž skutečnost, že ze vzorce pro derivaci součinu lze odvodit vzorec pro tzv. derivaci per partes, po česku po částech. Toto pro snadnost ponecháváme jako domácí cvičení a budeme pokračovat patřičnou větou. Věta 1.4 Nechť funkce u(x) a v(x) jsou spojité v intervalu (a,b) (může být i a = — 00 či b = 00) a nechť mají v tomto intervalu spojité derivace. Pak zde platí u(x)v'(x)dx = u(x)v(x) — j u'(x)v(x)dx či stručnéji uv'dx = uv — l u'vdx Použití této metody si ukážeme na několika příkladech. 4 KAPITOLA 1. NEURČITÝ INTEGRÁL Příklad 1.2 Vypočtěte integrál J xexdx. Jelikož máme integrovat součin funkcí, můžeme to zkusit metodou per partes. Nadějné je i to, že oba součinitele umíme derivovat a integrovat. Jelikož se ex nemění ani při integraci ani při derivaci, zaměříme svou pozornost na druhého činitele. Zatímco j xdx = ^- + c, tak derivace x' = 1. Volba je tedy jasná— u = x, u' = 1, v' = ex, v = ex. Podle vzorce je xexdx = xex — / exdx = xex — ex + c. Metodu per partes lze použít i opakovaně, musíme však vidět světlo na konci tunelu, nikov však od protijedoucího vlaku. Příklad 1.3 Vypočtěte integrál j x2 sinxdx. Volíme u = x2, u' = 2x, v' = sin x, v = — cos x. Podle vzorce je x2 sin xdx = — x2 cosx + 2 / x cos xdx Na integrál j x cos xdx použijeme opět metodu per partes s volbou u = x, u' = 1, v' = cosx, v = sinx a obdržíme x cos xdx = x sin x — I sin xdx = x sin x + cosx. Je tedy /,^inrf, = -^cos, + 2(,Sin, + coS,) + e Další způsob použití metody per partes ukazuje následující příklad. Příklad 1.4 Vypočtěte integrál jinxdx. Zpočátku nás asi zarazí, že zde není žádný součin. To nás však nesmí odradit, podle Palackého věty si ho prostě vytvoříme—1.x. Pak už je volba jednoznačná, a to u = \nx, u' = ^, v' = 1 a v = x. Aplikací metody per partes máme ľ 1 ln xdx = x ln x — / — .x = x \nx — x + c. J x Obdobným způsobem můžeme integrovat i funkce cyklometrické, jejichž integrály jste v přehledu základních vzorců jistě postrádali. Závěrem ještě jeden způsob, jak lze využit tuto metodu. Příklad 1.5 Vypočtěte integrál j ex sinxdx. Zde součin máme, obě funkce lze snadno derivovat i integrovat, tady je těžké si vybrat. Po delším přemýšlení zvolíme následující možnost: u = sinx, u' = sinx, v' = v = ex. Aplikujeme patřičný vzorec a máme ex sinxdx = ex — ex cos xdx. Z výsledku jsme trochu rozpačití, máme však pevnou vůli a z nastoupené cesty nesejdeme. Volba bude u = cosx, u' = — sinx a v' = v = ex. Výsledkem pak bude 1.1. PRIMITIVNÍ FUNKCE A NEURČITÝ INTEGRÁL 5 Ocitli jsme se ve stejné situaci, v jaké byl jistý Australan, který na dotaz svého kamaráda, proč má na hlavé bouli, odpovédél: Ale koupil jsem si nový bumerang a ten starý jsem odhodil. Nám tato rána do hlavy ale nevadí, naopak se nám v ní rozsvítí, protože si všimneme, že po odstránení závorek mají integrály opačná znaménka. Výše uvedený vztah můžeme tedy chápat tak, že se jedná o rovnici, v níž je neznámou zadaný integrál, či chcete-li zavedeme substituci j ex smxdx = t. Po vyřešení rovnice zjistíme, že je ex(sinx + cos x) e únxdx =--h c. 2 Podobně neexistuje univerzální vzorec pro integraci funkce složené, na základě věty o derivaci složené funkce však můžeme odvodit větu o substituci v integrálu. Věta 1.5 Nechť funkce F(t) je primitivní funkcí k funkci f(t) v intervalu (a, (3). Nechť funkce t = f{x) má derivaci tp'(x) v intervalu (a,b) (intervaly mohou být i nekonečné). Pro každé x G (a, b) nechť je číslo f{x) G («,/?)• Pak v intervalu (a,b) je funkce F(ip(x)) primitivní funkcí k funkci f( O, x E (—a, a). Zde použijeme substituci x = asinŕ, t G (—f, f); dx = acostdt. Dosadíme do integrandu a upravíme. V a2 — x2dx = / Va2 — a2 sin21 cos tdt = a2 / y 1 — sin21 cos tdt = a2 cos2 tdt K výpočtu tohoto interálu můžeme použít více způsobů, nejlepší a nejkratší cesta k výsledku vede přes vzorec pro poloviční úhel. a2 J cos2 tdt=a- j (1 + cos 2ť)dt = y (t + ^ sin 2t | + c Abychom ukázali, že se nám dostalo slušného vychování, vzpomeneme si na vzorec pro dvojnásobný úhel a tzv. goniometricou jedničku. Výsledek upravíme na tvar a2 ( 1 \ a2 -2 — I t + - sin2ŕ ] = — (t + sin t cos t) = —(t + sin ŕ • V 1 - sin2 ŕ) 2 y 2 y 2 2 Ze substituční rovnice x = a sin ŕ máme t = arcsin |; po dosazení a drobné úpravě, kterou jistě odhalíte sami máme a,2 x x Va2 — x2 = — arcsin —|—Va2 — x2 + c 2 a 2 1.2 Integrace racionální lomené funkce Velice často jsme postaveni před úkol integrovat racionální lomenou funkci, proto se této problematice budeme věnovat podrobněji. Nejprve si vezměte vaše poznámky z algebry a zopakujte si vše, co víte o polynomu a racionální lomené funkci. Definice a věty budu zmiňovat, jen někdy je však budu uvádět v jejich klasické podobě. Začnu připomínkou skutečnosti, že každá neryze lomená racionální funkce může být vyjádřena jako součet polynomu a ryze lomené racionální funkce. Jelikož integrace polynomu je brnkačka, budeme se věnovat pouze integraci ryze lomené racionální funkce. Již staří Egypťané používali tzv. kmenové zlomky, tedy zlomky typu -. ostatní zlomky vyjadřovali právě pomocí těchto kmenových, tedy například ^ = 2| + l|. U ryze lomené racionální funkce platí něco obdobného. Zopakujte si pojem kořen polynomu a vězte, že každý polynom lze rozložit na součin polynomů stupně nejvýš dva. Pak lze formulovat následující větu. Věta 1.7 Nechi je dána ryze lomená racionální lomená funkce a nechi platí Q(x) = a(x - ...(x- am)km(x2 + ptx + gx)řl ... (x2 + pnx + qn)ln, přitom cti je kořen polynomu Q(x) násobnosti k,i a kvadratické polynomy odpovídají komplexně sdruženým kořenům násobnosti /j. Pak existují jednoznačně určená reálná čísla označená velkými písmeny tak, že platí P(x) A1 A2 Akl Q(x) x — ct\ ' (x — cii)2 {x — a^kl 1.2. INTEGRACE RACIONÁLNÍ LOMENÉ FUNKCE 9 | Bx B2 ^_____Bkm -řrirr + Li -rř2^ + L2 Kh x + Li x2 + pix + gi (x2 + + gx)2 (x2 + + g^1 Mix + iVi M2x + iV2 Mlnx + iVZn x2 + pnx + gn (x2 + pnx + gn)2 (x2 + pnx + gn)ř Uznávám, že ta věta je strašná a dokáže odpudit většinu lidí od studia matematiky. Leč neházejte flintu do žita, mohl by ji tam někdo najít a vy byste měli problémy. Na následujících příkladech uvidíte, že je to věta velmi jednoduchá a porozumí jí každý z vás. Než se pustíme do vlastního rozkladu, uvedu ještě dvě věty. Věta 1.8 Nechť P{x) a Q{x) jsou mnohočleny stupně n a nechi se shodují pro n + 1 hodnot proměnné x. Pak jsou tyto mnohočleny totožné. Věta 1.9 Nechi P{x) a Q{x) jsou mnohočleny stupně n a nechť se shodují všechny koeficienty u týchž mocnin Pak jsou tyto mnohočleny totožné. Nyní se s optimismem pustíme do řešení příkladů. Budeme hned řešit integrál. Příklad 1.15 Vypočtěte integrál J ^.3-11^+34^-24^x■ ^e ^° ^egra/ ryze lomené racionální funkce, nejdříve rozložíme jmenovatel na součin polynomů stupně nejvýš dva. Po jistém úsilí, například pomocí Homérova schématu zjistíme, že polynom má tři reálné kořeny x = l,x = 4ax = 6. Věřím, že jste začali hledat kořeny jen mezi děliteli čísla 24- Patřičný rozklad bude 5x2 -39^ + 64 A B C x3 — llx2 + 3Ax — 24 x — 1 x — 4 x — 6 Zatím neurčité koeficienty A, B a C určíme následujícím způsobem. Rovnici vynásobíme společným jmenovatelem a obdržíme 5x2 - 39x + 64 = A(x - 4)(x - 6) + B(x - l)(x - 6) + C(x - l)(x - 4) Vzpomeneme si na větu 1. 8. a přinutíme polynomy na levé a na pravé se rovnnaly dosazením tří různých hodnot x . Samozřejmě můžeme dosadit libovolné hodnoty, ale podíváme-li se na pravou stranu, mělo by nás trknout, že bude velmi výhodné dosadit právě kořeny. x = 1 =>. 30 = 15A =>■ A = 2 x = 4 =>■ -12 = -6B =>■ B = 2 a konečně x = 6 =>• 10 = 10C =>• C = 1 Je tedy 5x2 - 39x + 64 ,dx=--1---1--) dx = \n\(x - l)2(x - A)2(x - Q)\+c x3 - llx2 + 34x - 24 J J \x- 1 x-A x-6J IV J K 1 y n Pokud koukáte na výsledek poněkud nedůvěřivě, zopakujte si pravidla pro logaritmování. 10 KAPITOLA 1. NEURČITÝ INTEGRÁL Příklad 1.16 Vypočtěte integrál J j^r^^^dx. Zde jsou kořeny již naznačeny, příslušný rozklad bude x2 + 27x-32 A B C (x — ľ)2(x — 7 (x — l)2 x — 1 x — 7 Volba x = 1 dává —6 = — 6A =>- A = 1, volba x = 7 dává 136 = 36C =>- C = 3. Jelikož nám kořeny došly, zvolíme ještě x = 0 a s využitím již spočítaných koeficientu máme —32 = — 7 + 7B + 3; tedy B = —4. Je tedy -x2 + 27x - 32 , ľ f 1 4 3 ' dx = /--—---1--- dx (x — l)2(x — 7) J \{x — l)2 x — 1 x —7 -1 , (x-7)3 In i-b + c x — 1 (x S případem, kdy polynom má pouze reálne kořeny jsme se již popasovali, hodíme nyní čučku na případ, kdy ve jmenovateli je nerozložitelný (ireducibilní) polynom stupně dva. Budeme se pro začátek zabývat jen patřičným parciálním zlomkem. Příklad 1.17 Vypočtěte integrál J x2+QX3+10dx. Diskriminant je roven číslu-4, trojčlen nelze dál rozložit. Máme však štěstí, v čitateli je skoro derivace jmenovatele, jen je ten výraz poněkud malý. Tak si ho zvětšíme známým trikem, když ho vynásobíme jedničkou tvaru 2 • |. Máme tedy ľ ^ + 3 , 1 ľ 2x + 6 , 1, , 2 , , - / —--dx = - / —-dx =—mix + 6x + 10)+c = ln v x2 + 6x + 10+c J x2 + 6x + 10 2 J x2 + 6x + 10 2 K J Jistě víte, proč místo absolutní hodnoty tam vyskočily jen závorky a kde se nakonec vzala ta odmocnina. Příklad 1.18 Vypočtěte integrál j x2]^1Qdx. Tady nám zase vadí ta šestnáctka, kdyby tam byla jednička, tak je to jasný arkus tangens. Tož šiji tam podle Palackého věty vytvořme. 1 1 ľ 1 , 1 f Atdt 1 1 x —-dx = — / —~-dx = — / —- = - arctgŕ + c = - arctg —h c x2 + 16 16 J g + 1 16 J ť + 1 4 6 4 6 4 Příklad 1.19 Vypočtěte integrál j x2+^+5dx. Spojíme-li poznatky z předchozích dvou integrálů, můžeme psát x + 6_l[2x+12_l[2x + 4: f 4 x2 + Ax + 5 2 J x2 + 4x + 5 2 J x2 + 4x + 5 J x2 + 4x + 5 Jednotlivé sčítance budeme integrovat zvlášť. 1 f 2x + A , „ 4 ľ A :dx = /--—-dx = 4 arctg [x + 2) + c x2 + 4x + 5 7 + 2)2 + 1 Doufám, že jste si vzpomněli na úpravu zvanou doplnění na čtverec. Zadaný integrál je tedy :dx = ln (x2 + 4x + 5) + 4 arctg (rr + 2) + c x2 + Ax + 5 1.2. INTEGRACE RACIONÁLNÍ LOMENÉ FUNKCE 11 Na závěr ještě jeden příklad. Příklad 1.20 Vypočtěte integrál J (^i^^ä+io) ^x• Nejprve rozložíme na parciální zlomky 2x2 + 13a: - 25 A 5a: + C (x + l)(x2 + 2x+ 10) ~~ a:+ 1 + a:2 + 2a: + 10 Známou úpravou obdržíme 2x2 + 13a: - 25 = A(a:2 + 2a: + 10) + (Bx + C)(a; + 1) a po roznásobení 2x2 + 13a; - 25 = Ax2 + 2Ax + 1(L4 + 5a:2 + Bx + Cx + C Koeficienty spočítáme porovnáním součinitelů u odpovídajících si mocnin, čímž obdržíme soustavu tří rovnic o tří neznámých. 2 = 2A +C 13 = 2A +B +C -25 = 10A +C Řešením jsou čísla A = —A, B = 6, C = 15. Využitím výše uvedených poznatků máme 2x2 + 13a: - 25 , ľ -Adx 1 f 2x + 2 , f dx -dx = /--1— / —--dx (x + l)(x2 + 2x + 10) j i + l 3 J x2 + 2a: + 10 J (x + l)2 + 9 = -4 ln \x + 11 + i ln (x2 + 2a: + 10) + ^ arctg —y- + c Takto jsouce vybaveni, můžeme konečně přistoupit k výpočtu integrálu, bude hodně záležet na tom, kolik úloh spočítáte. Mám pro vás ještě dvě zprávy, jedna z nich je dobrá a druhá špatná, já je uvedu v opačném pořadí než jak je obvyklé u vtipů. Špatnou zprávou je to, že přes všechnu uvedenou teorii se nám nepodaří vždy najít analytické vyjádření neurčitého integrálu. Například j ^^dx =?. To je rozdíl oproti derivování, kdy umíme zderivovat každou elementární funkci. Dobrou zprávou je skutečnost, že existuje mnoho knížek, kde jsou různé integrály tabe-lovány, ty pořádné správočníky obsahují stovky vzorců a tím se vyhneme často mnoha složitým výpočtům. Dřinu strojům, pardon tabulkám. Závěrem této kapitolky vám ještě řeknu, že zatímco derivace představuje jednosměrnou ulici, k integrálu můžeme dojít mnoha cestami, jde o to vybrat tu nejpohodlnější. KAPITOLA 1. NEURČITÝ INTEGRÁL Kapitola 2 Určitý integrál a jeho užití 2.1 Definice a metody výpočtu Nyní se pustíme do integrálu určitého, nemohu jinak, než začít problémem z fyziky. Naším úkolem je vypočítat práci, kterou vykoná plyn při ději izotermickém, změní-li se jeho objem z hodnoty V\ = lm3 na hodnotu v2 = lm3. My sice máme k dispozici poměrně jednoduchý vzorec W = pAV, leč tento platí jen pro děj izobarický, kdy je tlak konstantní, kdežto při ději izotermickém se tlak mění v závislosti na objemu podle Boyle-Mariotova zákona pV = konst.. My pro jednoduchost položíme tuto konstantu rovnu jedné, prostě si vybereme tu teplotu, při níž to tak je. Je tedy p = y. Tato funkce je na intervalu [1; 2] spojitá, dokonce je zde klesající. Jelikož je to problém technický, nepotřebujeme výsledek přesný na miliony desetinných míst. Proto budeme uvažovat takto: Jelikož tlak během celého děje klesá, nemáme k dispozici žádný vzorec. Rozdělme tedy děj na několik fází, kdy se tlak sice změní, ale ne zase tak moc, abychom ho s přimhouřením obou očí nemohli považovat za konstantní. Pak můžeme použít vzorec pro děj izobarický a celkovou práci určíme tak, že sečteme jednotlivé dílčí výsledky. Řekněme, že budeme objem brát po 0,2, dostaneme následující hodnoty. V 1,0 1,2 1,4 1,6 1,8 2,0 p 1,00 0,83 0,71 0,63 0,56 0,50 Pesimista si řekne, že vezme nejmenší hodnotu z každého intervalu a vyjde mu W = (0, 83 + 0, 71 + 0, 63 + 0, 56 + 0, 5).0, 2=0,65. To optimista ví, že parní stroj je úžasné síly zdroj a vezme si naopak ty hodnoty největší, čímž získá W = (1 + 0, 83+0, 71+0, 63+0, 56).0, 2=0, 75. Je jasné, že pesimista to podcenil a že skutečně vykonaná práce bude větší, zatímco u optimisty je tomu naopak. Řekněme, že půjdeme zlatou střední cestou a prohlásíme, že plyn vykonal práci W = 0, 70 J. Také víme, že kdybych nebyl líný a rozdělil interval na více dílků, byl by výsledek přesnější. Také jsem mohl zvolit jiný způsob dělení a to takto. Zpočátku klesá tlak poměrně rychle, intervaly budou kratší. Ke konci pak mohu volit úseky delší, neboť tlak již tak prudce neklesá. K tomuto příkladu se ještě vrátíme. Zkusme ještě jeden příklad, a to výpočet obsahu obrazce omezeného křivkami y = 0, x = 0, x = 1, y = x2. Uděláme-li si obrázek, jedná se o lichoběžník, jehož jedna strana se nám poněkud pokřivila. Zopakujeme-li postup z předchozího příkladu, tak zjistíme, že dělíme-li po 0,2, obdrží pesimista hodnotu s = 0, 24 a 13 14 KAPITOLA 2. URČITÝ INTEGRÁL A JEHO UŽITÍ optimista S = 0,44. Tentokrát nebudeme líní a dělení zjemníme na polovinu, tedy délka úsečky na ose x bude 0,1. Pesimista získá hodnotu s = 0,285, optimista pak S = 0, 385. Již od dob Archimedových víme, že správná hodnota je |. K této hodnotě se oba přibližují, jeden odspodu a druhý od vrchu. Z hlediska konstrukce je zřejmé, že ji ani jeden z nich nemůže překročit. Nyní již můžeme definovat určitý integrál. Def. 2.1 Nechť je v intervalu [a,b] dána spojitá funkce f{x). Zvolme v tomto intervalu n — 1 bodů Xi,X2, ■ ■ ■, xn-i, které vyhovují nerovnostem a = XQ 0 existuje takové ô > 0, že nerovnost \Sn - I\ < e je splnéna pro každé dělení d daného intervalu, pro které platí v(d) < ô, a to nezávisle na volbě bodů Píšeme n k=l Zdůrazňuji, že tato limita je reálné číslo závislé na funkci a intervalu, nikoliv však na dělení intervalu či volbě bodů Konečně následuje pointa celého procesu. Věta 2.1 Nechť f(x) je spojitá v intervalu [a,b]. Pak limita integrálních součtů existuje a platí n „5 1 = iimny2f(Šk)Axk = / f(x)dx. k=l Ja Této větě se říká základní věta integrálního počtu. V dnešní době výkonných počítačů lze určitý integrál snadno spočítat například tak, že budeme zjemňovat dělení intervalu tak dlouho až se horní a dolní součet nebudou v rámci požadované přesnosti lišit. Integruje se vesele více než tři sta let a naši pradědové neměli co se týče výpočetní techniky takové možnosti jako my. Jak si tedy počínali? To naznačí následující definice. Def. 2.5 Nechť funkce f{x) je spojitá v intervalu I a nechť F{x) je libovolná funkce k ní primitivní. Pak je b f(x)dx = F(b) - F(á) Takto definovaný určitý integrál se nazývá Newtonův. Nyní máme integrály dva, matematici v rozmachu tvůrčí činnosti definovali integrály další (Lebesque, Stielts,...). Jestli si myslíte, že tím vnesli do matematiky pěkný nepořádek, tak jste na velkém omylu, neboť platí věta: Věta 2.2 Nechť k funkci f{x) existuje více určitých integrálů, pak jsou si rovny. 16 KAPITOLA 2. URČITÝ INTEGRÁL A JEHO UŽITÍ Tedy žádný zmatek, naopak výhoda, že výpočet nějakého integrálu můžeme nahradit výpočtem jiného, kdy je postup snazší. Vraťme se k úvodním příkladům. Práce plynu je počítána podle vzorce ľ2 dV r l9 / — = [\nV]l = ln2 - lni = ln 2=0, 69 Ji v V tomto případě jsme se skoro trefili, a to je dělení po 0,2 hodně hrubé. Podobně je \x ii _ }_ [3Jo" 3' Archimédes měl pravdu a my jsme se tím to postupem taky moc nezmýlili, vždyť kdybychom vzali průměrnou hodnot, tak bychom měli 0,335. Nyní si ukážeme výpočet určitého integrálu, jsme-li nuceni při hledání primitivní funkce použít metodu per partes. Postup je jednoduchý: pb u(x)v'(x)dx = [u(x)v(x)]ba — / u'{x)dx X (Í/X - i _ |q No a jeden příklad. Příklad 2.1 xsinxdx = [—xcosx]^ + / cosxdx = n o Jo S metodou substituční je to poněkud složitější, ale zvládneme to. Věta 2.3 Nechť je v uzavřeném intervalu I s krajními body a, b integrand tvaru f[g(x)]g'(x), kde funkce t = g{x) a její derivace g'{x) jsou spojité funkce via zároveň f(t) je spojitá funkce pro všechna t = g (x), kde x E I. Pak platí b rg(b) f[g(x)]g\x)dx = / f(t)dt. Jg(a) - ľ— fi Příklad 2.2 Vypočtěte integrál JQ2 sin xcosxdx. Samozřejmě že bychom se mohli po substituci navrátit k původním proměnným, ale byla by to zbytečná oklika. Nej-vodnější substitucí se jeví t = sin x, dt = cosxdx. Pretransformujeme meze x = 0=^ŕ = 0;:r=f=^ŕ=l. Potom je sin6 x cos xdx = / ŕdt o Jo n1 1 Věta 2.4 Nechť funkce f(x) je spojitá v uzavřeném intervalu I\ s krajními body a, b, funkce 0, přičemž rovnost platí pouze v případě, je-li funkce na tomto intervalu rovna nule. Věta 2.9 Nechť f{x) a g(x) jsou spojité na intervalu [a;b] a nechť pro všechny body tohoto intervalu platí f (x) < g{x). Pak je b pb f{x)dx < / g{x)dx. J a Věta 2.10 Nechť funkce f{x) je spojitá na intervalu [a; b] a nechť pro všechny body tohoto intervalu platí m <3 f{x) < M. Pak platí m(b - a) < í f(x)dx < M(b - a) J a Tato věta se nazývá první věta o střední hodnotě a používí se k odhadu hodnoty určitého integrálu v případě, že se nám nedaří najít primitivní funkci. Ukážeme si to na příkladě. 2.2. VLASTNOSTI URČITÉHO INTEGRÁLU 19 Příklad 2.5 Odhadněte hodnotu integrálu - ďr' o 10 + x3 - O, 5 cos10 x + V7x4 + 9 Funkci ve jmenovateli si označíme 0 v intervalu [a; b] je ja f(x)dx roven obshu křivočarého rovnoběžníka ohraničeného osou x, funkcí f (x) a přímkami o rovnicích x = a a x = b. Výpočet si ukážeme na příkladech. Je-li obrazec ohraničen křivkami f {x) a g (x), přičemž pro všechna x G [a; b] je f (x) > g (x), spočítáme jeho obsah podle vzorce P [f (x) - g (x)] dx. V tomto případě není nutné, aby funkční hodnoty obou funkcí byly nezáporné, vzpomeňte si na Cavalieriho princip. l+x2 Příklad 2.7 Vypočtěte obsah obrazce ohraničeného grafy funkcí f (x) g{x) = x2. Grafy obou funkcí mají dva společné body A[—l; 1] a B[l; 1]. Obsah tohoto obrazce je P x' 30 J doo x* 30 J doo 2 arctgrr--x & 3 7t- Příklad 2.8 Vypočtěte obsah rovinného obrazce ohraničeného křivkami y = 2x—x2 a y = —x. Grafem první funkce je parabola, která má vrchol V[l; 1] a je čumákem nahoru. Grafem druhé funkce je osa druhého a čtvrtého kvadrantu. Musíme spočítat x—ové souřadnice průsečíků těchto dvou křivek, budeme tedy řešit rovnici 2x — x2 = —x. Tato rovnice má dvě řešení, totiž x\ = 0 a x2 = 3, což jsou současně meze integrálu. P [2x — x2 — (—x)]dx -x x y 9 2' Objem tělesa, které vznikne rotací nezáporné funkce f{x) v intervalu [a; b] se spočítá podle vzorce V = ir f2(x)dx. Těleso jednoduše rozřežeme kráječem na salám na tenké plátky, které můžeme považovat za válec, jehož objem je -Kf2(xo)dx. Potom jednotlivé plátky zase složíme dohromady, to je ten integrál. Jako příklad si odvodíme vzorec pro objem rotačního kužele. Příklad 2.9 Odvoďte vzorec pro objem rotačního kužele. Toto těleso vznikne rotací přímky y = Lx kolem osy x. Je tedy V 7t v JI ~ 30 d/30 7t- 2 s r x -7ir2v A ještě jeden. 2.3. UŽITÍ URČITÉHO INTEGRÁLU 21 Příklad 2.10 Vypočtěte objem tělesa ohraničeného plochami, které vzniknou rotací parabol y = 1 — x2, y = x2 + 2 a přímek x = — 1 a x = 1 kolem osy x. Obrázky neumím, ale je zřejmé, že v daném intervalu je parabola číslo dva vždy nad parabolou číslo jedna, můžeme použít stejný trik jako při výpočtu obsahu křivočarého lichoběžníka. V = 7t j [(x2 + 2)2 - (1 - x2)2]dx = 2tt J (6x2 + 3)dx = 2tt [2x3 + 3x] J = IOtt. Pro výpočet délky křivky je důležitá následující věta. Věta 2.15 Nechí křivka je dána parametricky rovnicemi x = tp(t), y = ip(t), t G [«;/?]• Mají-li tyto funkce v tomto intervalu spojité derivace (v krajních bodech zprava či zleva), pak je tato křivka schopna rektifikace. Slova rektifikace se nelekejte, matematici tím jen vyjadřují skutečnost, že má konečnou délku. Tuto délku spočítáme na základě následujících vět. Věta 2.16 Jsou-li funkce tp'(t) a ip'(t) spojité v intervalu [a; (3], pak délka l křivky dané parametricky x = tp(t), y = ip(t), t G [a; (3] je dána vzorcem i= / ^cpf2(t) + 4>f2(t)dt Jestliže si uvědomíme, že i funkce y = f (x) představuje vlastně parametrické rovnice, kde x = t a y = f (t), lze snadno formulovat následující větu. Věta 2.17 Délka oblouku grafu funkce y = f (x), x G [a; b] je dána vzorcem ŕ l = / y/l + fv(x)dx. Příklad 2.11 Zkusme si nejprve odvodit vzorec pro délku kružnice. Zde nám přijdou vhod její parametrické rovnice, které jsou x = rcost, y = sinŕ, t G [0;27r]. Derivace jsou x' = — r siní a y' = cos ŕ. Obvod (délka) kružnice je tedy /»2-7t _ /»2-7t / = / y7t2 sin21 + r2 cos2 tdt = r dt = 2nr. Jo Jo A příklad na druhý vzorec. Příklad 2.12 Vypočtěte délku grafu funkce y = lncosrr, x G [0; |]. 7T 7T ľ3 /- ľ 3 1 ľ (x 7t\ l = y 1 + tg2 xdx = / -dx = ln tg ( —|— ) J o J o cos x L V 2 4/ =1,317. o Jednu fyzikální aplikaci jsme si uvedli v motivačním příkladu, podobně řešíme i další problémy, jako příklad mohu uvést výpočet polohy těžiště a statických momentů. Toto je však poměrně složité a jelikož jste většinou nefyzici, tak to ke kolokviu konkrétně vyžadovat nebudu. Závěrem vám odvodím vzorec pro potenciální energii pružiny. 22 KAPITOLA 2. URČITÝ INTEGRÁL A JEHO UŽITÍ Příklad 2.13 Je-li deformace pružná, je působící síla přímo úměrná výchylce, platí F = Kx, K je tuhost pružiny (Hookeův zákon). Koncový bod pružiny umístíme do počátku a pružinu potáhneme o d vpravo. Vypočítáme-li vynaloženou práci, máme i potenciální energii pružiny. d pd W= I Fdx= / Kxdx 'o Jo Kx2 d -. = -Kx2. o 2 Jak již bylo řečeno, v mnoha případech se nám nepodaří nalézt primitivní funkci, ačkoliv určitý integrál zcela jistě existuje. Jak bylo ukázáno, umíme aspoň odhadnout jeho hodnotu, ne vždy se s tím však můžeme spokojit. Z definice můžeme integrál spočítat alespoň přibližně, my si ukážeme tři metody. Metoda obdélníková spočívá v tom, že interval [a; b] rozdělíme na n stejných dílků, jejichž délka je Ax = ^-,to bude jedna strana obdélníka. Druhou bude funkční hodnota v bodě x,i. Vypočítáme obsah každého z obdélníků a výsledky sečteme. Obdržíme b b - a f(x)dx «-\f(xQ) + f(Xl) + ■■■ + /Orn_i)] n Za předpokladu, že funkce má v tomto intervalu spojitou derivaci, umíme odhadnout i chybu Rn. DAb — a)2 Rn < -, -Di = maxai- 1 jeho výpočtu limitu, je vše v pořádku, integrál konverguje k |. Výpočet limity je často velmi složitý, někdy nám však stačí znát, zda nevlastní integrál konverguje a pak ho se pokusíme vypočítat alespoň přibližně. Vět o konvergenci je více, my si uvedeme jen některé. Věta 2.19 Nechť funkce f{x) a g{x) jsou spojité v intervalu I a nechť zde platí 0 < f(x) < g{x). Konverguje-li integrál funkce jjg(x)dx, konverguje i jIf{x)dx. Diverguje-li jIf{x)dx, pak diverguje i jIg{x)dx. Interval I zahrnuje všechny výše uvedené případy. Příklad 2.16 Integrál JQ4 ^Jídx je konvergentní. V intervalu [0;4) zřejmě platí cos2 x 1 0 < , < V4 - x V4 - x Obě funkce jsou v tomto intervalu spojité. Dále je í4 1 / dx = [-2V4- x] = lim V4 - 4- 2 = 0 + 2 = 2. Jo ^/A-x 1 j0 *->4 Integrál větší funkce konverguje, konverguje tedy i zadaný integrál. Příklad 2.17 Integrál ^dx diverguje. V intervalu (0; 1] zřejmě platí 1 3" 0 < - < — a obě funkce jsou zde spojité. ľ 1 / — dx = ílnrrL = 0 — lim In x = +oo J0 x 10 x^O 'o Zadaný integrál musí tedy divergovat také. Věta 2.20 Nechť funkce f{x), g{x) a h{x) jsou spojité v intervalu I a nechť pro všechna x E I platí f (x) < g (x) < h(x). Konvergují-li integrály fľf(x)dx a jjh{x)dx, konverguje i jjg(x)dx. Interval I zahrnuje všechny možné případy. Závěrem několik slov o integrálu jako funkci horní meze. Zafixujeme-li v integrálu j f{x)dx dolní mez a horní mez budeme měnit, dostaneme pro každou hodnotu b jiné číslo. Každé hodnotě bi je přiřazeno číslo Jj. je tedy definována funkce U (x) = / f{t)dt = F (x) - F (a). J a Základní vlastnost tohot integrálu udává následující věta. 26 KAPITOLA 2. URČITÝ INTEGRÁL A JEHO UŽITÍ Věta 2.21 Je-li funkce f(t)dt spojitá v intervalu I, pak derivace určitého integrálu U (x) = f f(t)dt podle proměnné horní meze se v každém bodě x E I rovná hodnotě integrované funkce v tomto bodě, tedy je U'(x) = f (x). Je-li interval I někde uzavřen, jedná se o derivaci zprava či zleva. Kapitola 3 Funkce více proměnných V této kapitole se budeme věnovat funkcím více proměnných. Kdo se dobře učil v minulém semestru, bude mít úlohu značně usnadněnou, neboť řada věcí je stejných či alespoň hodně podobných. Některé věci jsou však znčně odlišné, tak si na to dávejte pozor čili bacha. Budu se snažit na to upozorňovat. Na druhé straně vám to zjednoduším tím, že se budeme skoro výhradně bavit o funkci dvou proměnných. 3.1 Limita a spojitost Def. 3.1 Reálná funkce dvou reálných proměnných je zobrazeni M —> R, kde M c IRxIR. Jinými slovy každé uspořádané dvojici [x; y] g M je přiřazeno právě jedno 2éR. Množina M se nazývá definiční obor, množina všech z, které jsou přiřazeny k nějaké uspořádané dvojici [x; y] se nazývá obor hodnot funkce. Píšeme z = f(x, y). Nebude-li hrozit nedorozuměni, budeme mluvit stručně o funkci dvou proměnných. Příklad 3.1 Určete definiční obor funkce z = arcsiny + ln (4 — x2 — y2). Budeme vycházet ze znalosti funkcí jedné proměnné a vzpomeneme si na definiční obory funkcí arkussinus a přirozený logaritmus. Obě mají jistá omezení, která musí platit současně, je tedy -1 < y < 1 n 4 - x2 - y2 > 0. Zatímco první podmínce vyhoví pás mezi přímkami y = — 1 a y = 1, druhé podmínce vyhoví všechny body uvnitř kruhu se středem v počátku a poloměrem r = 2, v prvním případě včetně hranice. Definiční obor je samozřejmě průnik obou oblastí, leč obrázek zatím neumím. Def. 3.2 Grafem funkce dvou proměnných nazýváme množinu uspořádaných trojic [x,y,z], kde body [x,y] patří do definičního oboru funkce. Jinými slovy je to plocha o rovnici z = f(x,y). Vrstevnice je křivka o rovnici f(x,y) = c. Takovým nejběžnějším příkladem grafu funkce dvou proměnných je obyčejná plastická mapa. Proměnné představují zeměpisná šířka a délka, funkční hodnotu pak nadmořská výška. Patří sem i ona původně normální mapa, která se nacházela v kanceláři 91. pěšího pluku a kterou učinil plastickou až kocour chovaný písaři. Jen připomínám, že prvním, kdo se o této změně dotykem přesvědčil byl oberst Schrôder a že to mělo pro písaře nepříjemné následky. Pojem vrstevnice je převzat z 27 28 KAPITOLA 3. FUNKCE VÍCE PROMĚNNÝCH geografie a má stejný význam—zlepšit představu o grafu funkce v dvourozměrném modelu. Uvedeme několik příkladů. 3.1. Z analytické geomtrie víte, že grafem funkce z = ax + by + c je rovina v IR3. 3.2. Grafem funkce z = a/9 — x2 — y2 je horní polovina kulové plochy se středem v počátku a poloměrem r = 3 nad podstavnou rovinou os x a y. 3.3. Grafem funkce z = x2 + y2 je rotační paraboloid s vrcholem v počátku, jehož osou je osa z. Vrstevnice tvoří soustředné kružnice o rovnicích x2 + y2 = c. Nyní přistoupíme k definici pojmů limita a spojitost. Začneme definicí okolí. Def. 3.3 Vzdálenost dvou bodů A[xi;yi] a B[x2;y2] rozumíme číslo d = a/(x2 - xi)2 + (y2 - yi)2 Vzdálenost bodů ve vícerozměrném prostoru si jistě odvodíte sami. Def. 3.4 ô-okolím bodu P nazýváme množinu všech bodů, jejichž vzdálenost od bodu P je menší než ô. Vyjmeme-li z této množiny samotný bod P, mluvíme o ryzím okolí. Def. 3.5 Řekneme, že funkce z = f(x,y) má v bodě M[xo;yo] limitu rovnou číslu L, jestliže ke každému e > 0 existuje ô > 0 takové, že pro pro všechny body z ryzího ô okolí bodu M platí \f(x,y) -L\[0;0] x y pokusíme se tedy spočítat limity pro různé cesty. Začněme přímkami y = kx. x + y x + kx 1 + k hm - = hm--— = hm-- [:r;j/H[0;0] x — y x^o x — k x x^o 1 — k Limity pro různé směry jsou různé, limita neexistuje. Tato situace paradoxně nenastane, pokud bychom se přibližovali po parabolách y = kx2. x + y x + kx2 1 + kx lim -= lim--—- = lim--— = 1 [x;j/]^[0;0] x — y x^o x — kx x^o 1 — kx Příklad 3.3 Zjistěte, zda existuje lim —+9^- • Zkusme se nejdřív přibližovat po přímkách y = kx. 2xy 2kx2 lim - = lim n [x;j/]^[0;0] xy + 2x — y x^o kx2 + 2x — kx ovšem s výjimkou k = 2, to je Ax lim x^o 2x + 2 - 2 Tato limita tedy neexistuje. K bodu Pq[xq; yo] se z bodu P[x; y] můžeme rovněž přibližovat po dvou kolmých přímkách x = p a y = q, kde p a g jsou konstanty, a to dvojím způsobem. Pak lze limitu funkce vypočítat postupným limitním přechodem funkce jedné proměnné, jak uvádí následující věta. 30 KAPITOLA 3. FUNKCE VÍCE PROMĚNNÝCH Věta 3.2 Označme lim y^yo lim f(x,y) Li lim X—>XQ lim f(x,y) y^yo Lo. Existuje-li limita lim f(x, y) = L, [x;y]->[x0;yo] pak platí L = Li = L2 Uvědomte si, že tato věta je implikací a představuje pouze podmínku nutnou, což značí, že bude sloužit k důkazu neexistence limity. Příklad 3.4 Rozhodněte, zda existuje limita ,. x-2y lim [*;*/]-> [0;0] 3x + y Určime postupné limity. lim lim x — 2y -2y y^O \_x^o 3x + y i- ľi- x~2y hm hm- x^O [y^o 3x + y lim y^O y lim — x^o 3y 1 3' Tato limita neexistuje. Podobně jako pro funkci jedné proměnné platí následující věta. Věta 3.3 Nechť pro všechny body x G O s výjimkou bodu M[xo;yo] platí f(x,y) g(x,y) a nechť je lim f(x,y) = L. Pak je i lim g(x,y) = L. [x;y]->[x0;yo] [x;y]->[x0;yo] Příklad 3.5 lim x — y lim (x — y){x2 + xy + y2) [x;y]"[í;-l] x4 - í/4 [x;y]"[í;-l] (x - y) (x + y) (x2 + y2) Tato funkce je v okolí bodu [1; — 1] shodná s funkcí x2 + xy + y2 z = - (x + y) (x2 + y2) Její limita je v tomto bodě rovna funkční hodnotě, tedy je x3 — y3 x2 + xy + y2 3 lim -= lim -= — [x;y]^[l;-l] xa - í/4 (x + y){x2 + y2 8 Při výpočtu limity můžeme použít i některé triky známé z funkce jedné proměnné, jeden příklad následuje. 3.2. PARCIÁLNÍ DERIVACE 31 Příklad 3.6 3(x2 + y2) lim [x;y]-+[0;0] ^x2 + + 4 - 2 Vynásobíme-li funkci jedničkou ve tvaru a/x2 + y2 + 4 + 2 v/rr2 + í/2 + 4 + 2 a upravíme-li, počítáme limitu lim 3(v/x2 + w2 + 4 + 2) = 12 [*»]->[0;0] Závěrem této části vám uvedu tři limity, které vám určitě něco připomenou. sin fix, y) lim M '.^ = 1 foj/Hbosž/o] /(z,?/) lim —-— = 1 [x;j,]^[x0;í/o] /(z,?/) ^ \ f(x,y) lim 1 + —-- I = e [x;y\-+[x0;y0\ \ f{X,y) 3.2 Parciální derivace Už problémy s limitou nám naznačují, že to s derivacemi vůbec nebude snadné. Pokud bychom chtěli udělat nějakou analogii s funkcí jedné proměnné, bylo by to značně obtížné. Proto půjdeme jinou cestou. Grafem funkce dvou proměnných je plocha. Pokud však plochu řízneme nějakou rovinou, tak řezem je křivka, křivku umíme popsat funkcí jedné proměnné—čajník je v kredenci. My budeme řezat rovinami kolmými k osám x a y. Def. 3.7 Nechi existuje limita f(x,yo) - f(x0,y0) lim >x0 X — Xq Pak tuto limitu nazveme parciální derivací podle x v bodě [x0; y0], značíme 9^qx^ či fx(xo,yo). Analogicky definujeme parciální derivaci podle y. Nechť existuje limita Hm f(x0,y) - f(x0,y0) y^y0 y-yo Pak tuto limitu nazveme parciální derivací podle y v bodě [xo;yo], značíme 9^v^ či fý(xo, yo). Má-li funkce parciální derivaci podle nějaké proměnné v každém bodě nějaké oblasti M, potom říkáme, že zde má parciální derivaci. Jinými slovy vznikne na této oblasti nová funkce, to je stejné jako u funkce jedné proměnné. 32 KAPITOLA 3. FUNKCE VÍCE PROMĚNNÝCH Protože jsme parciální derivace definovali jako derivace funkce jedné proměnné, platí pro ně všechna pravidla tak jak je znáte z kurzu MA1, nebudu je tedy uvádět. Stejně tak nebudu řešit derivace vyšších řádů, tam ovšem jeden problém přece jen vyvstane. Záleží na pořadí proměnných podle kterých derivujeme nebo nezáleží, to je oč tu běží. Odpověďnám dává Schwarzova věta. Věta 3.4 Nechť jsou derivace fxy(x,y) a fyx(x,y) jsou v bodě M[rr0;ž/o] spojité. Pak jsou si rovny. Obdobnou větu bychom mohli zformulovat i pro smíšené parciální derivace vyšších řádů. Jak vidíme, u spojitých funkcí je to s parciálními derivacemi jako s mušketýry, je jich o jednu víc než je jejich řád. Stejně jako tři mušketýři byli čtyři, tak i třetí parciální derivace jsou čtyři: fxxx, fxxy, fxyy a fý'ýy. Tak je tomu u parciálních derivací jakéhokoliv řádu. Příklad 3.7 Je-li z = u(x) + v(y), je z'x = u'(x), z' = v'(y), zxx = u"{x), z'ý = v"(y) a zxy = zyx = 0. Tak je-li z = 3x2 — 5y3 + 1, máme zx = 6x a z' = — 15y2. Pro druhé derivace vychází zxx = 6, z" = —30y a z" = z" = 0. Příklad 3.8 Je-li z = u{x)v{y), je z'x = u'(x)v(y), z' = u(x)v'(y), zxx = u"(x)v(x), z'ý = u(x)v"(y) a zxy = zyx = u'(x)v'(y). Tak je-li z = 5x2y4, máme z'x = 10xy4 a z'y = 20x2y3. Pro druhé derivace vychází zxx = 10y4, zyy = 60x2y2 a zxy = zyx = AOxy3. Příklad *í Q lp U 7 — 1P 7' - VĽäí 7> — A*W(y) 11 _ W'(x) „ _ „ _ r 1 iridii o.w je u z — je zx — v^ , zy — , zxx — ^ , zyx — zyx — —u ^2(y^ a z'ýy = —u{x)v ^V^l^xjv ^ ■ Konkrétní příklad dáme tento: z = Potom ie z' — — z' — — z" — D z" — — a z" — — ruuuiu je yy ' xxx x3 ' xxy x2 ' xyy yyy -luulíiiío diferenciál druhého řádu je d2z = z"xx(dx)2 + 2z"xydxdy + z'ýy(dy)2 = —^(dx)2 H—dxdy(dy)2 Totální diferenciál třetího řádu je pak 2y 3 d3z = 2///xxx(dx)3+34/^(^)2^+34/TO^(^)2+4ra(^)3 = ~^(dxf~^dx)2dy 3.3. TOTÁLNÍ DIFERENCIÁL 35 3.3.1 Extrémy funkce dvou proměnných Def. 3.11 Řekneme, že funkce f(x,y) má v bodě M[x0; y0] lokálni maximum (minimum), existuje-li okolí O bodu M takové, že pro všechna x E O platí f(xQ,y0) > f(x,y) (f(xQ,y0) < f(x,y). V případě ostrých nerovností mluvíme o ostrém lokálním maximu (minimu). Postup při stanovení extrémů je obdobný jako u funkce jedné proměnné. Věta 3.10 Nechť funkce f(x, y) má v bodě M[x0; y0] lokální extrém a nechť existují v bodě M parciální derivace prvního řádu. Pak je fx(xQ,y0) = fý(xQ,y0) = 0. Poznámka 3.4 Tak jako u funkce jedné proměnné budeme bod M nazývat bodem stacionárním. Stacionární (podezřelé z extrému) body budeme vyšetřovat pomocí následující věty. Věta 3.11 Nechť M je stacionární bod a nechť v jeho okolí existují spojité parciální derivace prvního a druhého řádu. Vypočtěme výraz D = fxx(X0,y0)fyy(X0,y0) ~ (fxy(X0,y0))2 Je-li D > 0, pak pro fxx > 0 je v bodě M lokální minimum a pro fxx(xo, y$) < 0 je v bodě M lokální minimum. Je-li D < 0, pak v bodě M extrém není, je-li D = 0, nemůžeme o extrému rozhodnout (extrém tady být může, ale nemusí). Poznámka 3.5 Označení D jsme nezvolili náhodou, D je de facto determinant druhého řádu, přičemž v hlavní diagonále jsou derivace podle xx ayy a ve vedlejší diagonále jsou derivace smíšené (ty jsou si samozřejmě rovny). Nyní ukážeme několik příkladů. Příklad 3.14 Nalezněte lokální extrémy funkce z = x3 + xy2 + 6x2 + y2. Nejprve spočítáme parciální derivace prvního řádu. f'x = 3x2 + y2 + 12x, f'y = 2xy + 2y . Položíme-li tyto derivace rovny nule, získáme čtyři stacionární body S\[—1;3], S2[—1; —3], S3[0;0] a S^[—4;0]. Spočteme tedy parciální derivace druhého řádu f xx = ^x 12, fxy = 2y, fyy = 2x + 2. Budeme postupně dosazovat jednotlivé stacionární body a počítat číslo D. V prvních dvou případech je D(Si) = —36, -D^) = —36, extrém nenastává. D(Ss) = 24 a protože je fxx(0, 0) = 12, je v počátku minimum. Naproti tomu je D(S4) = 24, ale tentokrát je f'x'x(—4,0) = —12, v bodě S4 je tedy maximum. 36 KAPITOLA 3. FUNKCE VÍCE PROMĚNNÝCH Příklad 3.15 Určete lokálni extrémy funkce z = x2 + Axy + 6y2 — 2x + 8y — 5. Začneme prvními parciálními derivacemi. z'x = 2x + Ay-2 z'y=Ax+12y + 8 Opět položíme obě derivace rovny nule, po vyřešení soustavy dvou lineárních rovnic získáme jediný stacionární bod S[7; —3]. Druhé parciální derivace jsou zxx = 2, zxy = 4 a z'ý = 12. Všechny jsou konstantní, není kam dosazovat a determinant má univerzální hodnotu D = 8 > 0. Jelikož je zxx = 2 > 0, je v bodě S lokální minimum. Jen pro zajímavost, jeho hodnota je z(7, —3) = —24. Příklad 3.16 Určete lokální extrémy funkce z = —3x4 — 5y4. Spočteme první derivace z'x = —12x3, z'y = —20y3. Jediným stacionárním bodem je počátek. Jdeme na druhé derivace. z"xx = — 36x2, z'ý = —60y2, zxy = 0. Zřejmě je D = 0, o extrému nemůžeme tímto způsobem rozhodnout. My si ale všimneme, že funkční hodnoty jsou mimo počátek záporné, je zřejmé, že v počátku bude maximum. Tak jako u funkce jedné proměnné můžeme určovat i extrémy absolutní, a to v případě, že je funkce definovaná na uzavřené oblasti. Ty pak mohou nastat buď v bodech lokálních extrémů nebo na hranici oblasti. Ukážeme si to na příkladu, nejdříve trochu teorie. Def. 3.12 Řekneme, že funkce f(x,y) má v bodě [xo;yo] absolutní maximum (minimum), jestliže pro všechny body [x;y] G M platí f (x, y) < f(xQ,yo) (f (x,y) > f(x0,yo)- Def. 3.13 Bod [xo;yo] nazveme vnitřním bodem množiny M, existuje-li okolí O tohoto bodu takové, že O G M. Množina, která obsahuje pouze vnitřní body, se nazývá otevřená. Bod [xo;yo] nazveme vnějším bodem množiny M, jestliže každé jeho obsahuje jak body množiny M, tak i body, které do ní nepatří. Množinu všech hraničních bodů nazýváme hranicí. Množina, která obsahuje všechny své hraniční body, se nazývá uzavřená. Def. 3.14 Množina M se nazývá omezená, existuje-li kruh K se středem v počátku tak, že M C K. Věta 3.12 Nechť f (x, y) je spojitá funkce definovaná na omezené uzavřené množině. Pak zde nabývá své nejmenší a největší hodnoty. Příklad 3.17 Stanovte absolutní extrémy funkce z = \/2x — x2 — Ay2. Po doplnění na čtverec zjistíme, že definičním oborem jsou vnitřní a hraniční body elipsy {x — l)2 + Ay2 < 1. Stacionární body určíme řešením soustavy rovnic ' - 2~2X _ g / _ -% _ q x ~ 2^2x - x2 - Ay2 ~ y ~ 2^2x - x2 - Ay2 ~ Existuje jediný stacionární bod S'fljO]. Je /(1,0) = 1. Stanovení extrémů na hranici je obecné velmi obtížné, vezmeme-li rozum do hrsti, tak vidím,e že funkce je na hranici rovna nule a jinak je kladná. Absolutní maximum je tedy ve středu elipsy a minimun na její hranici. 3.4. PŘEHLED LÁTKY 37 3.4 Přehled látky 3.1. Primitivní funkce a neurčitý integrál 3.2. Základní integrační metody 3.3. Integrace per partes 3.4. Substituční metoda 3.5. Integrace racionální lomené funkce 3.6. Určitý integrál (newtonův, Riemannův) 3.7. Vlastnosti určitého integrálu 3.8. Numerické metody výpočtu určitého integrálu 3.9. Užití určitého integrálu 3.10. Nevlastní integrál 3.11. Funkce dvou proměnných, definice, graf 3.12. Limita a spojitost funkce dvou proměnných 3.13. Parciální derivace 3.14. Totální diferenciál 3.15. Extrémy funkce dvou proměnných