Analýza: derivácie, integrály a limity

Príklad 1. Definujeme polynomickú funkciu premennej x a nájdeme jej prvú a tretiu deriváciu.

f=x^5-2x^3-1
D[f,x]
D[f,{x,3}]

-1 - 2 x^3 + x^5

-6 x^2 + 5 x^4

-12 + 60 x^2

Tie isté príkazy sa dajú zapísať pomocou palety matematických symbolov. Zápis f' vykoná substitúciu. Výpočet derivácie môžeme zabezpečiť symbolmi  ∂_     и    ∂_ (, ) . V symbole pre tretiu deriváciu pridáme príslušný počet x :

f = x^5 - 2x^3 - 1 f^′ f1 = ∂_x f f3 = ∂_ (x, x, x) f

-1 - 2 x^3 + x^5

(-1 - 2 x^3 + x^5)^′

-6 x^2 + 5 x^4

-12 + 60 x^2

Týmto spôsobom môžeme deklarovať ľubovoľnú funkciu:

h[x_] := x^5 - 2x^3 - 1 ∂_x h[x] h[1]

-6 x^2 + 5 x^4

-2

Príklad 2. Zobrazíme grafy funkcií z príkladu 1. Priebeh funkcie je zaujímavý na intervale [-2, 2], preto grafy derivácií zobrazíme iba na tomto intervale.

g1 = Plot[f, {x, -5, 5}] g2 = Plot[f1, {x, -2, 2}] g3 = Plot[f3, {x, -2, 2}]

[Graphics:HTMLFiles/index_14.gif]

⁃Graphics⁃

[Graphics:HTMLFiles/index_16.gif]

⁃Graphics⁃

[Graphics:HTMLFiles/index_18.gif]

⁃Graphics⁃

Príklad 3. Definujeme funkciu dvoch premenných a nájdeme jej parciálne derivácie podľa premennej x, potom podľa premennej y a napokon druhú zmiešanú parciálnu deriváciu ∂_ (x, y).

g[x_, y_] := 7x^2 - y^3 - Log[x] ∂_x g[x, y] ∂_y g[x, y] ∂_ (x, y) g[x, y]

-1/x + 14 x

-3 y^2

0

Príklad 4. Vynulujeme premennú f. Definujeme symbolickú deriváciu.

f = .
D[f[x],x]
D[5* x *f[x], x]

f^′[x]

5 f[x] + 5 x f^′[x]

Príklad 5. Premennej f priradíme novú hodnotu, pričom stará sa vymaže automaticky. Vypočítame neurčitý integrál z f - pomocou symbolov z palety, aj zápisom príkazovej funkcie, a napokon vypočítame určitý integrál funkcie f na danom intervale a vyjadríme jeho číselnú hodnotu.

f = (1 - 2x^3)/(2 + 3x) ∫fxIntegrate[f, x]    &nb ...         (* Here the symbol % takes the preceding result *)

(1 - 2 x^3)/(2 + 3 x)

1/81 (-6 x (4 - 3 x + 3 x^2) + 43 Log[2 + 3 x])

1/81 (-6 x (4 - 3 x + 3 x^2) + 43 Log[2 + 3 x])

1/81 (-24 + 43 Log[5/2])

0.19013

Príklad 6. Definujeme funkciu dvoch premenných s názvom newf. Vypočítame jej deriváciu podľa premennej x a dvojný integrál na danej oblasti.

newf[x_, y_] := x^4 + 3y d = ∂_x newf[x, y] ∫_0^1∫_1^2newf[x, y] yx %//N

4 x^3

47/10

4.7

Príklad 7. Definujeme trigonometrickú funkciu. Chceme vypočítať určitý integrál na danej oblasti, ale systém má problémy a výpočet je veľmi pomalý, funkcia cos[ x/y] nie je definovaná v bode y = 0. Zdĺhavé výpočty môžeme prerušiť súčasným stlačením klávesníc ALT + , alebo pomocou Kernel/Abort Evaluation menu.

f=. f[x_, y_] := Cos[( x)/y]  ∫_0^π∫_0^πf[x, y] xy  ... , π}]     (* Approximate integration by numerical method *)   

-1/4 π^2 (π - 2 (Cos[1] + Sin[1] + SinIntegral[1]))

3.73594

3.73594

NIntegrate :: slwcon : Numerical integration converging too slowly; suspect one of the followi ... r integrand is oscillatory try using the option Method->Oscillatory in NIntegrate.  More…

NIntegrate :: ncvb : NIntegrate failed to converge to prescribed accuracy after 13 recursive bisections in y near {x, y} = {3.000902059754438`, 0.0019899063331735444`} .  More…

3.73532

Príklad 8. Zobrazíme graf funkcie z predchádzajúceho príkladu, aby sme zistili, kde nastávajú singularity.

Plot3D[f[x,y],{x,0,1},{y,0,1}]

                                      1 Power :: infy : Infinite expression  --- encountered. More…                                      0.`

∞ :: indet : Indeterminate expression 0.` ComplexInfinity encountered. More…

                                      1 Power :: infy : Infinite expression  --- encountered. More…                                      0.`

∞ :: indet : Indeterminate expression 0.` ComplexInfinity encountered. More…

                                      1 Power :: infy : Infinite expression  --- encountered. More…                                      0.`

General :: stop : Further output of Power :: infy will be suppressed during this calculation. More…

∞ :: indet : Indeterminate expression 0.` ComplexInfinity encountered. More…

General :: stop : Further output of ∞ :: indet will be suppressed during this calculation. More…

Plot3D :: plnc : f[x, y] is neither a machine-size real number at {x, y} = {0.`, 0.`} nor a list of a real number and a valid color directive. More…

Plot3D :: plnc : f[x, y] is neither a machine-size real number at {x, y} = {0.041666666666666664`, 0.`} nor a list of a real number and a valid color directive. More…

Plot3D :: plnc : f[x, y] is neither a machine-size real number at {x, y} = {0.08333333333333333`, 0.`} nor a list of a real number and a valid color directive. More…

General :: stop : Further output of Plot3D :: plnc will be suppressed during this calculation. More…

Plot3D :: gval : Function value Indeterminate at grid point xi = 1, yi = 1 is not a real number. More…

Plot3D :: gval : Function value Indeterminate at grid point xi = 2, yi = 1 is not a real number. More…

Plot3D :: gval : Function value Indeterminate at grid point xi = 3, yi = 1 is not a real number. More…

General :: stop : Further output of Plot3D :: gval will be suppressed during this calculation. More…

[Graphics:HTMLFiles/index_60.gif]

⁃SurfaceGraphics⁃

Príklad 9. Definujeme novú funkciu a vypočítame jej dvojný integrál na oblasti, ktorej hranice vzhľadom na y sú určené ako funkcie premennej x.

f=. ;     f = x^2 + 1/(1 + y^2) ∫_0^1∫_0^x^(1/2) fyx N[%]

x^2 + 1/(1 + y^2)

-5/7 + π/2

0.856511

Príklad 10. Príklad trojného integrálu - najprv určíme jeho exaktnú hodnotu, potom vypočítame približnú numerickú hodnotu výsledku s dvojnásobnou presnosťou.

∫_0^1∫_x/2^x∫_ (x + y)/2^(x + y) 1/(1 + x + y + z) zyx %//N

1/432 (-1408 Log[4] + 675 Log[5] - 216 Log[8] + 676 Log[13] - 72 Log[256])

0.0462077

RowBox[{0.0462077, }]

Predchádzajúci výpočet sa dá zapísať aj v štandarnej forme, ako:

f= 1/(1+x+y+z)
Integrate[f,{x,0,1},{y,x/2,x},{z,(x+y)/2,x+y}]
N[%]

1/(1 + x + y + z)

1/432 (-1408 Log[4] + 675 Log[5] - 216 Log[8] + 676 Log[13] - 72 Log[256])

0.0462077

Príklad 11. Vyšetrujme teraz inú funkciu. Systém zobrazí jej graf na intervale [0, 10] napriek tomu, že funkcia nie je v bode x = 0 definovaná. Pri zobrazení na väčšom intervale však už graf nie je nakreslený v okolí bodu 0.

f = Sin[x]/x Plot[f, {x, 0, 10}] Plot[f, {x, 0, 100}]

Sin[x]/x

[Graphics:HTMLFiles/index_75.gif]

⁃Graphics⁃

[Graphics:HTMLFiles/index_77.gif]

⁃Graphics⁃

Príklad 12. Vypočítajme limity a hodnotu nevlastného integrálu funkcie f.

Limit[f, x0] Limit[f, xInfinity] ∫_0^∞fx

1

0

π/2

Príklad 13. Príklad ďalšej funkcie s bodmi nespojitosti.

h = (x + 1)/(1 + x^2 - x^4) Plot[h, {x, -10, 10}]

(1 + x)/(1 + x^2 - x^4)

[Graphics:HTMLFiles/index_85.gif]

⁃Graphics⁃

Príklad 14. Ukážky grafu funkcie z predchádzajúceho príkladu na vhodných intervaloch.

Plot[h,{x,-1, 1}]
Plot[h,{x,-1.5, -1.2}]  

[Graphics:HTMLFiles/index_87.gif]

-Graphics-

[Graphics:HTMLFiles/index_88.gif]

-Graphics-

Created by Mathematica  (October 10, 2007)