\input /home/pasky/school/fastex/lib.tex

\chapter{Taylorův polynom}{

\scpart{Opakování}{
	{\bf Polynom} je každá funkce tvaru
	$$a_nx^n + a_{n-1}x^{n-1} + \cdots + a_1x + a_0$$

	$a_0,\ldots,a_n \in \real$ pevné --- koeficienty. $x \in \real$ je proměnná.

	\notes{
	\list{
	\listItem{
		Všechny polynomy jsou definovány a spojité na $\real$.
	}
	\listItem{
		Jestliže $a_n \ne 0$, pak $n = $ stupeň polynomu.
	}
	\listItem{
		Je-li $\cP$ polynom, stupeň $\cP = n$, pokud $\cP'$ je polynom, stupeň $\cP' = n-1$.
	}
	}

	\examples{
		$$\eqalign{ \deg \cP = 0 & \Rightarrow \cP \textbox{konstantní},\ \cP(x) \equiv a_0 \cr
		            \deg \cP = 1 & \Rightarrow \cP(x) = a_0x + a_0 \textbox{lineární} \cr
		            \deg \cP = 2 & \Rightarrow \textbox{kvadratická funkce, atd.} }$$
	}
	}

}

\scpart{Motivace}{
	Aproximace funkcí pomocí polynomů.

	Mějme $f$, $a \in \real$, $\existss f'(a) \in \real$. Tečna $t(x) = f(a) + f'(a)(x-a)$.
	Potom
	$$\eqalign{f(a) &= t(a)\cr
		   f'(a)&= t'(a)}$$

	Navíc
	$$\eqalign{ \lim_{x\to a} (f(x)-t(x)) &= 0 \textbox{(triviální), ale i}\cr
		   \lim_{x\to a} {f(x)-t(x)\over x-a} \becauseof{L'H}{=} \lim_{x\to a} (f'(x)-t'(x)) &= 0 \textbox{(aproximace je lepší než lineární)}}$$

	Chceme aproximace vyšších řádů.
}

\scpart{Idea}{
	Mějme $f$, $a \in \real$, $\existss f'(a) \in \real$. Chceme funkci $f$
	aproximovat polynomem v bodě $a$.

	\scpart{Aproximace polynomem stupně 1}{
		$$P_1(x) = \a x + \beta$$

		Polynom by měl procházet bodem $a$:
		$$P_1(a) = f(a) \Rightarrow P_1(x) = c(x-a) + f(a)$$
		(tedy hodnota funkce a aproximujícího polynomu splývají).

		Zároveň chci:
		$$ P_1'(a) = f'(a) \Rightarrow P_1(x) = f(a) + f'(a)(x-a)$$
		Nejlepší aproximace polynomem stupně $\le 1$ je totiž
		{\bf tečna} $t(x) = f(a) + f'(a)(x-a)$.
		Potom $$\lim_{x\to a} {f(x)-t(x) \over x-a} = 0$$

		Ukazuje se, že tyto dva požadavky jsou ekvivalentní.
	}

	\scpart{Aproximace polynomem stupně 2}{
		Příklad: $\cos x$ --- lineární aproximace nic moc, kvadratická aproximace
		je už lepší. $P_2(x) = f(a) + f'(a)(x-a) + c(x-a)^2$, $c = ?$.

		Chci: $$\eqalign{P_2(a) &= f(a)\cr
		                 P_2'(a) &= f'(a)} \et P_2''(a) = f''(a)$$

		$$\Rightarrow P_2'(x) = f'(a) + 2c(x-a) \Rightarrow P_2'(a) = f'(a) $$
		(první derivace splývá)
		$$\Rightarrow P_2''(x) = 2c \Rightarrow c = {f''(a) \over 2}$$
		$$\Rightarrow P_2(x) = f(a) + f'(a)(x-a) + {f''(a) \over 2}(x-a)^2$$
		To je nejlepší kvadratická aproximace $f$ v bodě $a$.

		Navíc platí: $$\lim_{x\to a} {f(x)-P_2(x) \over (x-a)^2} = 0$$
		(Důkaz: \DCV)
	}
}


\definition{
	Nechť $f$ je reálná funkce, $a \in \real$, $n \in \nat$.
	Pak funkci
	$$ T_n^{f,a}(x) = f(a) + f'(a)(x-a) + \cdots + {f^{(n)}(a)\over n!}(x-a)^n $$
	nazýváme {\bf Taylorovým polynomem funkce $f$ řádu $n$ v bodě $a$}.
}

\notes{
\list{
	\listItem{$T_n^{f,a}$ je polynom, $\deg T_n^{f,a} \le n$.}
	\listItem{Příklad: $f(x) = \sin x$, $a = 0$, potom:
		$$ \eqalign{T_1^{\kern0.8pt \sin,0}(x) &= x \cr
		            T_2^{\kern0.8pt \sin,0}(x) &= x\qquad \textnote{neboť $\sin''(0)=0$} \cr
		            T_3^{\kern0.8pt \sin,0}(x) &= x - {x^3\over6}\qquad \textnote{cvičení} \cr
		            T_4^{\kern0.8pt \sin,0}(x) &= x - {x^3\over6}} $$

		Poučení: $\deg T_n^{f,a}$ nemusí být roven $n$.
	}
	\listItem{Derivace Taylorova polynomu:
		$$ \left(T^{f,a}_n\right)'(x) = 0 + f'(a) + f''(a)(x-a) + {f'''(a)\over2}(x-a)^2 + $$
			$$ + \cdots + {f^{(n)}(a)\over(n-1)!}(x-a)^{n-1} = T^{f',a}_{n-1} $$

		Tedy $\left(T^{f,a}_n\right)'(x) = T^{f',a}_{n-1}$. (!)
	}
	\listItem{Splývání polynomu a funkce:
		$$ \sesac{\eqalign{T^{f,a}_n(a) &= f(a)\cr
		            \left(T^{f,a}_n\right)'(a) &= f'(a)\cr
			    &\vdots\cr
		            \left(T^{f,a}_n\right)^{(n)}(a) &= f^{(n)}(a)}}{(n+1) \textbox{požadavků.}} $$
	}
}
}

\scpart{Otázka}{
	$$ \lim_{x\to a}{f(x) - T^{f,a}_n(x) \over (x-a)^n} \hypothetically{=} 0 $$
}


\theorem{33}{o aproximaci funkce Taylorovým polynomem}{
	Nechť $f$ je reálná funkce, $a \in \real$ a nechť existuje vlastní $f^{(n)}(a)$.
	Nechť $\cP$ je polynom, $\deg \cP \le n$. Pak platí:
	$$ \lim_{x\to a} {f(x)-\cP(x) \over (x-a)^n} = 0 \Longleftrightarrow \cP = T^{f,a}_n $$

	\note{
		\proofleftimpl{
			$$ \lim_{x\to a}{f(x) - T^{f,a}_n(x) \over (x-a)^n} = 0 $$
			(odpověď na otázku je kladná)
		}
		\proofrightimpl{
			Tuto vlastnost má mezi polynomy stupně $\le n$ jen Taylorův.
		}
	}

	\proof{
		\proofleftimpl{
			Indukcí:
			\list{
			\namedListItem{$n=1$}{
				$$ \lim_{x\to a}{f(x) - T^{f,a}_1 \over x-a} = \lim_{x\to a}{f(x)-f(a)-f'(a)(x-a)\over x-a} $$
				$$ \lim_{x\to a}\left({f(x)-f(a)\over x-a} - f'(a)\right) = 0 $$
			}
			\namedListItem{$n-1 \derived n$}{
				$$ \lim_{x\to a}{f(x) - T^{f,a}_n(x) \over (x-a)^n} \becauseof{L'H ``${0\over0}$''}{=}
					\lim_{x\to a}{f'(x) - \left(T^{f,a}_n\right)'(x) \over n(x-a)^{n-1}} = $$
				$$ = {1\over n}\lim_{x\to a}{f'(x)-T^{f',a}_{n-1}(x) \over (x-a)^{n-1}} = 0
					\textnote{dle indukčního předp.} $$
			}
			}
		}

		\proofrightimpl{
			\list{
			\namedListItem{1. krok}{
				$$ \lim_{x\to a}{\cP(x) - T^{f,a}_n(x) \over (x-a)^n}
					= \undernote{\lim_{x\to a}{\cP(x) - f(x) \over (x-a)^n}}{0 \textnote{předp.}}
					+ \undernote{\lim_{x\to a}{f(x) - T^{f,a}_n(x) \over (x-a)^n}}{0 \textnote{právě dok.}}
				$$
			}
			\lemma{}{
				$Q$ polynom, $\deg Q \le n$, $a \in \real$ a platí
				$$\lim_{x\to a}{Q(x) \over (x-a)^n} = 0 $$
				pak $Q \equiv 0$.

				\proof{
					Indukcí:
					\list{
					\nlistItem{$n=1$:

						$Q$ je lineární, $Q(a) = 0$ (viz níže)
						$\Rightarrow$ $Q(x) = c(x-a)$, $c \in \real$.

						Víme:
						$$ \eqalign{\lim_{x\to a}{Q(x)\over x-a} &= 0 \cr
						 \Rightarrow 0 = \lim_{x\to a}{c(x-a)\over x-a} &= c
						   \Rightarrow c = 0 \Rightarrow Q(x) \equiv 0 } $$
					}
					\nlistItem{$n-1 \derived n$:
						$$ \sesac{\eqalign{\deg Q &\le n\cr
						                   \lim_{x\to a}{Q(x) \over (x-a)^n} &= 0}
							}{\Rightarrow Q(a) = 0} $$
						$$\Rightarrow \textbox{$a$ je kořen polynomu $Q$}$$
						$$\Rightarrow Q(x) = (x-a)R(x), \textbox{$R$ je polynom, $\deg R \le n-1$}$$

						Potom
						$$ 0 = \lim_{x\to a}{Q(x)\over(x-a)^n} = \lim_{x\to a}{(x-a)R(x)\over(x-a)^n}
							= \lim_{x\to a}{R(x)\over (x-a)^{n-1}} $$
						$$ \becauseof{ind. předp.}{\Longrightarrow} R \equiv 0
							\Rightarrow Q(x) = (x-a)\cdot0 \Rightarrow Q \equiv 0 $$
					}
					}

					\sqed\  (důkaz lemmatu)
				}
			}
			\penalty-1000
			\namedListItem{2. krok}{

				Víme:
				$$ \lim_{x\to a}{\cP(x) - T^{f,a}_n(x) \over (x-a)^n} = 0 $$
				($\cP(x) - T^{f,a}_n(x)$ je polynom, $\deg \le n$)
				$$ \becauseof{lemma}{\Longrightarrow} \cP(x) - T^{f,a}_n(x) \equiv 0 \Rightarrow \cP = T^{f,a}_n $$
			}
			}
			\qed
		}
	}

	\example{
		$$ \sesac{\eqalign{f &= \sin\cr x &= 0}}{\Rightarrow \lim_{x\to 0}{x-\sin x\over x^3} = {1\over6}} $$
		(plyne z implikace ``$\Leftarrow$'')
	}
}


\theorem{34}{obecný tvar zbytku}{

	Nechť $f$ má vlastní $(n+1)$-ní derivaci v intervalu $[a,x]$, $a,x\in\real$, $x>a$.

	Nechť $\varphi$ je spojitá funkce na $[a,x]$, která má na $(a,x)$ vlastní a nenulovou derivaci.

	Pak existuje $\xi \in (a,x)$ takové, že
		$$ R^{f,a}_n(x) \defined f(x) - T^{f,a}_n(x) = {1\over n!}{\varphi(x)-\varphi(a)\over \varphi'(\xi)}
			f^{(n+1)}(\xi)(x-\xi)^n $$
	(obrázek)

	\scpart{Důsledek 1 (Lagrangeův tvar zbytku):}{
		$$R^{f,a}_n(x) = {1\over (n+1)!}f^{(n+1)}(\xi)(x-a)^{n+1}\qquad (\existss\xi \in (a,x))$$
		\proof{
			Volím $\varphi(t) = (x-t)^{n+1}$. \DCV
		}
	}

	\scpart{Důsledek 2 (Cauchyův tvar zbytku):}{
		$$R^{f,a}_n(x) = {1\over n!}f^{(n+1)}(\xi)(x-\xi)^n(x-a)\qquad (\existss\xi \in (a,x))$$
		\proof{
			Volím $\varphi(t) = t$. \DCV
		}
	}

	\proof{
		Pro $t \in [a,x]$ definujeme
		$$ F(t) := f(x) - T^{f,t}_n(x) = f(x) - \left(f(t) + {f'(t)\over1!}(x-t) + \cdots + {f^{(n)}(t)\over n!}(x-t)^n\right) $$

		\penalty-100
		Potom:
		\list{
			\listItem{$F$ je spojitá na $[a,x]$ a má vlastní derivaci na $(a,x)$.}
		        \listItem{$F(a) = f(x) - T^{f,a}_n(x) = R^{f,a}_n(x)$}
			\listItem{$F(x) = 0$}
			\listItem{$\varphi$ je spojitá, $\existss \varphi' \ne 0$ vlastní na $(a,x)$}
		}

		Tedy podle \ias{}{23}{} (Cauchy o střední hodnotě):
		$$ \existss \xi \in (a,x):\qquad {F(x)-F(a) \over \varphi(x) - \varphi(a)} = {F'(\xi) \over \varphi'(\xi)} $$

		Přitom:
		$$ F'(t) = -\left(f'(t) + f''(t)(x-t) + f'(t)(-1) + {f'''(t)\over2}(x-t)^2 + {f''(t)\over2}2(x-t)(-1) +\right.$$
			$$ \left. + \cdots + {f^{(n)}\over n!}n(x-t)^{n-1}(-1) + {f^{(n+1)}(t)\over n!}(x-t)^n\right) $$
		$$ \becauseof{požírání}{\Longrightarrow} F'(t) = -{f^{(n+1)}(t)\over n!}(x-t)^n\qquad \forall t\in(a,x) $$
		$$ \Rightarrow F'(t) = -{f^{(n+1)}(\xi)\over n!}(x-\xi)^n $$

		Platí tedy:
		% Customized to make the middle equal signs centered
		$$ \vcenter{\ialign{\strut\hfil #&# \hfil\cr
		             $\displaystyle{{F(x) - F(a)\over \varphi(x) - \varphi(a)}}$
				     &$\displaystyle{= {0-R^{f,a}_n(x)\over \varphi(x)-\varphi(a)}}$\cr
			     \hfill $=$ \hfill & \hfill $=$ \hfill\cr
			     $\displaystyle{{F'(\xi)\over \varphi'(\xi)}}$
				     &$\displaystyle{= {-{f^{(n+1)}(\xi)\over n!}(x-\xi)^n\over \varphi'(\xi)}}$\cr
			   }} $$
		
		To znamená:
		$$ R^{f,a}_n(n) = {1\over n!}{\varphi(x)-\varphi(a)\over \varphi'(\xi)}f^{(n+1)}(\xi)(x-\xi)^n $$

		\qed
	}

	\shortnote{\ias{}{32}{} platí i pro $x < a$.}
}

\example{
	Rozvoj funkce $e^x$:
	$$ f'(x) = e^x,\ f(x) = e^x, \ldots ,\ f^{(n)}(x) = e^x $$
	$$ a = 0 \Rightarrow f(a) = f'(a) = \cdots = f^{(n)}(a) = 1 $$

	Zaveďme si:
	$$ T^{\hskip1pt\exp,0}_{\hskip1pt n} = \sum_{j=0}^n {x^j \over j!} $$

	Pak:
	$$ e^x = \sum_{j=0}^n {x^j\over j!} + R^{\exp,0}_n(x) = \sum_{j=0}^n {x^j\over j!} + e{\xi n\over(n+1)!}x^{n+1} $$

	Pro každé pevné $x$ je $\xi n \in (a,x)\quad \forall n \in N$, tedy $e^{\xi n} \le e^x$.  Tudíž
	$$ |R^{\exp,0}_n(x)| \le {e^x|x|^{n+1}\over(n+1)!} $$
	a platí
	$$ \lim_{n\to \infty} |R^{\exp,0}_n(x)| = 0 $$
	$$ \Rightarrow \lim_{n\to\infty}\sum_{j=0}^n {x^j\over j!} = e^x
		\Rightarrow e^x = \sum_{j=0}^\infty {x^j\over j!}\qquad \forall x \in \real$$
}

\impl{
	$$ \sum_{j=0}^\infty {x^j\over j!} = \lim_{n\to\infty}\left(1+{x\over n}\right)^n $$
	a speciálně pro $x = 1$:
	$$ e = \sum_{j=0}^\infty {1\over j!} = \lim_{n\to\infty}\left(1 + {1\over n}\right)^n $$
	$$ e = 1 + 1 + {1\over2} + {1\over6} + {1\over24} + {1\over120} + \cdots $$

	Chceme-li spočítat $e$ s přesností $0,001$:
	$$ e = \sum_{j=0}^n {1\over j!} + \sum_{j=n+1}^\infty {1\over j!}
		= \sum_{j=0}^n {1\over j!} + \sum_{j=0}^\infty {1\over (n+1+j)!} \le $$
	$$ \le \sum_{j=0}^n {1\over j!} + {1 \over (n+1)!}\sum_{j=0}^\infty {1\over(n+1)^j}
		= \sum_{j=0}^n {1\over j!} + {1 \over (n+1)!}{1 \over 1-{1\over n+1}}
		= \sum_{j=0}^n {1\over j!} + {1 \over n\cdot n!}$$

	Tedy:
	$$ \sum_{j=0}^n {1\over j!} < e < \sum_{j=0}^n {1\over j!} + {1 \over n\cdot n!}\qquad \forall n \in \nat $$
	(na přenost $0,001$ stačí $n = 6$)
}

\scpart{Důkaz iracionality e}{
	Nechť $e = {p / q}$, kde $p,q\in\nat$. Potom pro
	$$m = \sum_{j=0}^n {1\over j!}$$
	máme
	$$m < {p \over q} < m + {1\over n\cdot n!}$$
	$$qmn! < pn! < qmn! + {q\over n}$$

	Nyní stačí zvolit $n > q$.

	\XXX
}

\scpart{Rozvoj funkce sin}{
	Mějme $\sin$, $a = 0$:

	$$ \eqalign{\eqalign{\sin'(x) &= \cos(x)\cr
	                     \sin''(x) &= -\sin(x)\cr
			     \sin'''(x) &= -\cos(x)\cr
			     \sin''''(x) &= \sin(x)}
		  & \eqalign{\sin'(0) &= 1\cr
		             \sin''(0) &= 0\cr
			     \sin'''(0) &= -1\cr
			     \sin''''(0) &= 0}} $$

	$$ \Rightarrow T^{\hskip1pt\sin,0}_{\hskip1pt 2n+1}(x) = \sum_{j=1}^n (-1)^{j+1} {x^{2j+1}\over(2j+1)!} $$
	$$ R\hskip1pt^{\sin,0}_{2n+1}(x) = {1\over(2n+2)!}\sin^{(2n+2)}(\xi n)\cdot x^{2n+2} $$
	$$ |R\hskip1pt^{\sin,0}_{2n+1}(x)| \le {|x|^{2n+2}\over(2n+2)!} \to 0\qquad \forall x,\ n\to\infty $$
	$$ \Rightarrow \sin x = \sum_{j=1}^\infty (-1)^{j+1} {x^{2j+1}\over(2j+1)!},\qquad x \in \real $$

	Obdobně
	$$ \cos x = \sum_{j=0}^\infty (-1)^j {x^{2j}\over(2j)!},\qquad x \in \real $$
}

\scpart{Rozvoj funkce log(x+1)}{
	Mějme $f(x)\colon y=\log(x+1)$. Pro $x=0$ pak $f(0)=0$.

	$$ \eqalign{f'(x) &= {1\over1+x}\cr
	            f''(x) &= -{1\over(1+x)^2}\cr
		    f'''(x) &= {2\over(1+x)^3}\cr
		    &\ \vdots\cr
		    f^{(n)}(x) &= {(n-1)!\over(1+x)^n}(-1)^{n+1}} $$
	$$ f^{(n)}(0) = (-1)^{n-1}(n-1)! $$
	$$ \Rightarrow T^{\log(x+1),0}_n(x) = \sum_{j=1}^n {1\over j!}(-1)^{j-1}(j-1)!\cdot x^j
		= \sum_{j=0}^n(-1)^{j-1}{x^j\over j} $$
	$$ R^{\log(x+1),0}_n(x) = {1\over(n+1)!} {(-1)^n\cdot n!\over(1+\xi)^{n+1}}x^{n+1}
		= {(-1)^n\cdot x^{n+1}\over n(1+\xi)^{n+1}} $$
	To nelze odhadnout pro $|x|\ge1$, pro $|x|<1$ ano:
	$$ \Rightarrow \log(x+1) = \sum_{j=1}^\infty (-1)^{j-1} {x^j\over j},\qquad x \in (-1,1] $$
}

\definition{
	Řekneme, že funkce $f$ je v bodě $a$ zprava {\bf malá řádu $n$}, jestliže
	$$ \lim_{x\to a_+}{f(x)\over(x-a)^n} = 0. $$
	Píšeme
	$$ f(x) = \sigma\left((x-a)^n\right),\qquad x \to a_+ $$
	(Laudaův symbol).

	\shortdef{Pozor}{Jde pouze o symboliku, {\it neplatí} žadná rovnost!}
}

\examples{
\list{
	\listItem{$e^x = \sum_{j=0}^n {x^j\over j!} + \sigma(x^n),\quad x\to0_+$}
	\listItem{$x-\sin x = \sigma(x^2),\quad x\to0$}
	\listItem{$\cos x = 1 - {x^2\over2} + \sigma(x^3),\quad x\to0$}
}
}

\example{
	Mějme sumu
	$$ \sum_{n=1}^\infty {n! e^n\over n^{n+p}} $$
	Pro která $p$ konverguje?

	Raabe:
	$$ \lim_{n\to\infty} n\left({a_n\over a_{n+1}}-1\right)
		= \lim_{n\to\infty} n\left({1\over e}\left(1+{1\over n}\right)^{n+p} - 1\right) = $$
	$$ = \lim_{n\to\infty} n\left({1\over e}\cdot e^{(n+p)\log\left(1+{1\over n}\right)}-1\right)
		= \lim_{n\to\infty} n\left({1\over e}\cdot e^{(n+p)\left({1\over n}-{1\over2n^2}+\sigma\left(1\over n^2\right)\right)}-1\right) = $$
	$$ = p - {1\over2} - 1 + \sigma\undernote{\left(1\over n\right)}{\to0},$$
	tedy konverguje pro $p > {3\over2}$.
}



}


\bend