diff --git a/Lectures_my/NumMet/2016/Lecture7/Makefile b/Lectures_my/NumMet/2016/Lecture7/Makefile new file mode 100644 index 0000000..0bcbf45 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/Makefile @@ -0,0 +1,7 @@ +all: lecture7.pdf + +lecture7.pdf: lecture7.tex + xelatex lecture7.tex + xelatex lecture7.tex +clean: + rm *.aux *.bbl *.blg *.log *.out *.pdf *.sh *.xml *.snm *.toc *.nav diff --git a/Lectures_my/NumMet/2016/Lecture7/beamercolorthemeFlip.sty b/Lectures_my/NumMet/2016/Lecture7/beamercolorthemeFlip.sty new file mode 100755 index 0000000..d6b9c50 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/beamercolorthemeFlip.sty @@ -0,0 +1,179 @@ +% Color style file for the BEAMER FLIP THEME +% Copyright 2012 by Flip Tanedo +% This file may be distributed and/or modified +% 1. under the LaTeX Project Public License and/or +% 2. under the GNU Public License. + + + +\mode +\newif\ifbeamer@Flip@dark +\beamer@Flip@darkfalse +\DeclareOptionBeamer{dark}{\beamer@Flip@darktrue} +\ProcessOptionsBeamer + + +% These are a few colors that I like. + +\definecolor{crimsonred}{RGB}{153,0,0} % Neurtal red, good for dark or light bg +\definecolor{darkcharcoal}{RGB}{25,25,25} % Darker gray +\definecolor{charcoal}{RGB}{51,51,51} % Darker gray +\definecolor{ash}{RGB}{100,100,100} % medium gray +\definecolor{paleblue}{RGB}{0,102,102} % More of an `ocean' color +\definecolor{turtlegreen}{RGB}{51,153,0} % A more neutral green +\definecolor{paleale}{RGB}{204,204,102} % Only for dark BG +\definecolor{lager}{RGB}{140,110,10} % Use instead of pale ale for white BG +\definecolor{regal}{RGB}{90,0,120} % A more neutral purple +\definecolor{jeans}{RGB}{20,30,150} % A more neutral blue + + +% Define some "invariant colors" so that a color that doen't work with a +% light/dark background will automatically be substituted with one that does. +% These will be rewritten in beamercolorthemeFlipDark + +\ifbeamer@Flip@dark + \definecolor{FlipSand}{RGB}{204,204,102} % paleale, Only for dark BG + \definecolor{FlipGreen}{RGB}{0,255,0} % Neon-like green, only for dark BG +\else + \definecolor{FlipSand}{RGB}{140,110,10} % lager, Only for light BG + \definecolor{FlipGreen}{RGB}{51,153,0} % turtlegreen, works for both +\fi + +%% Future revision: It would be nice to define these in terms of the palette. + + +\ifbeamer@Flip@dark + \setbeamercolor{normal text}{fg=white, bg=black} + \setbeamercolor{structure}{fg=gray, bg=black} % fg controls bullets + \setbeamercolor*{footerlike}{fg=white,bg=black} % Controls the footerline +\else + % Default to black on white + \setbeamercolor*{footerlike}{fg=white,bg=charcoal} % for charcoal footer + \setbeamercolor{structure}{fg=charcoal, bg=white} % fg controls bullets + % \setbeamercolor{normal text}{fg=white, bg=black} +\fi + + + + + +% These defined colors can be used in \textcolor{colorname}{text in that color}. + +% Colors for Keynote background: DARK +\definecolor{keynotebottom}{rgb}{0.32,0.3,0.38} +\definecolor{keynotemiddle}{rgb}{0.08,0.08,0.16} +\definecolor{keynotetop}{rgb}{0,0,0} +%%% Example of usage: +%% \setbeamertemplate{background canvas}[vertical shading] +%% [bottom=keynotebottom, middle=keynotemiddle, top=keynotetop] + + + + +% ALERT AND COMMENT COLORS (just copies of the above colors) +% (I don't know how ot do this more elegantly) +\ifbeamer@Flip@dark + \definecolor{ALERT}{RGB}{153,0,0} % crimsonred + \definecolor{Alert}{RGB}{51,153,0} % turtlegreen + \definecolor{alert}{RGB}{204,204,102} % paleale + \definecolor{charcoal}{RGB}{51,51,51} % charcoal + \definecolor{comment}{RGB}{80,80,80} % light charcoal + \definecolor{Comment}{RGB}{100,100,100} % ash + \definecolor{COMMENT}{RGB}{80,20,120} % regal +\else +% \definecolor{ALERT}{RGB}{153,0,0} % crimsonred + \definecolor{ALERT}{RGB}{28,134,238} % DodgerBlue2 +% \definecolor{Alert}{RGB}{51,153,0} % turtlegreen + \definecolor{Alert}{RGB}{255,127,36} % chocolate1 +% \definecolor{alert}{RGB}{140,110,10} % lager + \definecolor{alert}{RGB}{0,205,102} % SpringGreen3 + \definecolor{charcoal}{RGB}{80,80,80} % light charcoal + \definecolor{comment}{RGB}{51,51,51} % charcoal + \definecolor{Comment}{RGB}{100,100,100} % ash + \definecolor{COMMENT}{RGB}{80,20,120} % regal +\fi + + +% \definecolor{crimsonred}{RGB}{153,0,0} % Neurtal red, good for dark or light bg +% \definecolor{charcoal}{RGB}{51,51,51} % Darker gray +% \definecolor{ash}{RGB}{100,100,100} % medium gray +% \definecolor{paleblue}{RGB}{0,102,102} % More of an `ocean' color +% \definecolor{turtlegreen}{RGB}{51,153,0} % A more neutral green +% \definecolor{paleale}{RGB}{204,204,102} % Only for dark BG +% \definecolor{lager}{RGB}{140,110,10} % Use instead of pale ale for white BG +% \definecolor{regal}{RGB}{90,0,120} % A more neutral purple +% \definecolor{jeans}{RGB}{20,30,150} % A more neutral blue + + +% A few useful colors +\setbeamercolor{alerted text}{fg=alert} + +%% I would like ot use these colors, but there seems to be a bug in Beamer +%% where the footnote color infects the normal text color +%% See: http://tex.stackexchange.com/questions/8264/beamer-changing-the-footnote-color-also-changes-the-normal-text-color/8268#8268 + +% \setbeamercolor{footnote}{parent=normal text, fg=FlipSand} +% \setbeamercolor{footnote mark}{parent=normal text, fg=FlipSand} + +% Blocks +\setbeamercolor{block title}{fg=gray!30!white,bg=charcoal} +\setbeamercolor{block body}{parent=normal text,bg=gray, fg=white} + +\setbeamercolor{block title example}{fg=turtlegreen!30!white,bg=turtlegreen!40!black} +\setbeamercolor{block body example}{parent=normal text,bg=turtlegreen, fg=white} +% +\setbeamercolor{block title alerted}{fg=crimsonred!30!white,bg=crimsonred!40!black} +\setbeamercolor{block body alerted}{parent=normal text,bg=crimsonred, fg=white} + + + +% These "setbeamercolor" commands are more standardized. + + +%% From the Beamer User Guide: +% Outer themes base the color of navigational elements and, possibly, also of other elements, on the +% four palette colors. The “primary” palette should be used for the most important navigational elements, +% which are usually the ones that change most often and hence require the most attention by the audience. +% The “secondary” and “tertiary” are less important, the “quaternary” one is least important. +% By default, the palette colors do not have a background and the foreground ranges from structure.fg to black. +% For the sidebar, there is an extra set of palette colors, see palette sidebar primary. + +\setbeamercolor*{palette primary}{fg=crimsonred,bg=charcoal} +\setbeamercolor*{palette secondary}{fg=paleale,bg=lager} +\setbeamercolor*{palette tertiary}{fg=turtlegreen,bg=jeans} +\setbeamercolor*{palette quaternary}{fg=paleblue,bg=regal} + +% \setbeamercolor*{palette primary}{fg=crimsonred,bg=white} +% \setbeamercolor*{palette secondary}{fg=charcoal,bg=white} +% \setbeamercolor*{palette tertiary}{fg=paleblue,bg=white} +% \setbeamercolor*{palette quaternary}{fg=turtlegreen,bg=white} + + + +%% These are some structures that I define for use in the outer theme. + +\setbeamercolor*{crimsontown}{fg=paleblue,bg=crimsonred} +\setbeamercolor*{bluetown}{fg=white,bg=DodgerBlue2} +%\setbeamercolor*{titlelike}{fg=crimsonred} +\setbeamercolor*{titlelike}{fg=DodgerBlue2} +\setbeamercolor*{topbarlike}{fg=white, bg=charcoal} + + + + + +%% Some old definitions that I might want to get back to +% +% \setbeamercolor*{frametitle}{bg=white,fg=charcoal} +% \setbeamercolor*{part title}{bg=white,fg=charcoal} +% \setbeamercolor*{item}{fg=charcoal, bg=white} +% +% \setbeamercolor*{separation line}{} +% \setbeamercolor*{fine separation line}{} + + + + + +\mode + diff --git a/Lectures_my/NumMet/2016/Lecture7/beamerinnerthemeFlip.sty b/Lectures_my/NumMet/2016/Lecture7/beamerinnerthemeFlip.sty new file mode 100755 index 0000000..8dc12a0 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/beamerinnerthemeFlip.sty @@ -0,0 +1,95 @@ +% Inner style file for the BEAMER FLIP THEME +% Copyright 2010 by Flip Tanedo +% This file may be distributed and/or modified +% 1. under the LaTeX Project Public License and/or +% 2. under the GNU Public License. +% +% Based on: beamerinnerthemeFlip by Marco Barisione + +\mode + +% % Use alternative title page style. +% \DeclareOptionBeamer{alternativetitlepage}[true]{\def\beamer@Flip@alternativetitlepage{#1}} + +% Logo to use in the alternative title page. +% \def\beamer@Flip@titlepagelogo{} +% \DeclareOptionBeamer{titlepagelogo}{\def\beamer@Flip@titlepagelogo{#1}} + +% Bullet shape. +\DeclareOptionBeamer{bullet}{\def\beamer@Flip@bullet{#1}} +\DeclareOptionBeamer{shadow}[true]{\def\beamer@themerounded@shadow{#1}} +\ExecuteOptionsBeamer{shadow=false} + + +% \ExecuteOptionsBeamer{alternativetitlepage=false,bullet=square} +\ProcessOptionsBeamer + + + +% Margins. +\newlength{\beamer@Flip@normalmargin} +\setlength{\beamer@Flip@normalmargin}{.06\paperwidth} +\setbeamersize{text margin left=\beamer@Flip@normalmargin} +\setbeamersize{text margin right=\beamer@Flip@normalmargin} +\setlength\leftmargini{.6\beamer@Flip@normalmargin} +\setlength\leftmarginii{.6\beamer@Flip@normalmargin} +\setlength\leftmarginiii{.6\beamer@Flip@normalmargin} + +% Normal title page. +\defbeamertemplate*{title page normal}{Flip theme}[1][] +{ + \vbox{} + \vfill + \begin{centering} + \begin{beamercolorbox}[wd=\paperwidth,sep=8pt,center,#1]{title page header} + \usebeamerfont{title}\inserttitle\par% + \ifx\insertsubtitle\@empty% + \else% + \vskip0.25em% + {\usebeamerfont{subtitle}\usebeamercolor[fg]{subtitle}\insertsubtitle\par}% + \fi% + \end{beamercolorbox}% + %\vskip0.5em\par + \begin{beamercolorbox}[sep=8pt,center,#1]{author} + \usebeamerfont{author}\insertauthor + \end{beamercolorbox} + \begin{beamercolorbox}[sep=8pt,center,#1]{institute} + \usebeamerfont{institute}\insertinstitute + \end{beamercolorbox} + \begin{beamercolorbox}[sep=8pt,center,#1]{date} + \usebeamerfont{date}\insertdate + \end{beamercolorbox}\vskip0.5em + {\usebeamercolor[fg]{titlegraphic}\inserttitlegraphic\par} + \end{centering} + \vfill +} + + +% Items. +\defbeamertemplate{itemize item}{squarealt}% +{\tiny\raise.5ex\hbox{\donotcoloroutermaths$\blacksquare$}} +\defbeamertemplate{itemize subitem}{squarealt}% +{\tiny\raise.4ex\hbox{\donotcoloroutermaths$\square$}} +\defbeamertemplate{itemize subsubitem}{squarealt}% +{\tiny\raise.3ex\hbox{\donotcoloroutermaths$\blacksquare$}} +% +\defbeamertemplate{itemize item}{circlealt}% +{\small\raise.2ex\hbox{\donotcoloroutermaths$\bullet$}} +\defbeamertemplate{itemize subitem}{circlealt}% +{\small\raise.1ex\hbox{\donotcoloroutermaths$\circ$}} +\defbeamertemplate{itemize subsubitem}{circlealt}% +{\scriptsize\raise.1ex\hbox{\donotcoloroutermaths$\bullet$}} +\setbeamertemplate{items}[circlealt] +% +\def\circletext{circle} +\ifx\beamer@Flip@bullet\circletext + \setbeamertemplate{items}[circlealt] +\else + \setbeamertemplate{items}[squarealt] +\fi + +\setbeamertemplate{blocks}[rounded][shadow=\beamer@themerounded@shadow] + + +\mode + diff --git a/Lectures_my/NumMet/2016/Lecture7/beamerouterthemeFlip.sty b/Lectures_my/NumMet/2016/Lecture7/beamerouterthemeFlip.sty new file mode 100755 index 0000000..d12f6d6 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/beamerouterthemeFlip.sty @@ -0,0 +1,189 @@ +% Outer style file for the BEAMER FLIP THEME 212 +% Copyright 2012 by Flip Tanedo +% This file may be distributed and/or modified +% 1. under the LaTeX Project Public License and/or +% 2. under the GNU Public License. +% +% Based on: beamerouterthemesplit by Till Tantau + + +\mode + +% Watermark +\DeclareOptionBeamer{watermark}{\def\beamer@Flip@watermark{#1}} +\ExecuteOptionsBeamer{watermark=} % Default value is empty + + +\newif\ifbeamer@Flip@bigpagenumber +\beamer@Flip@bigpagenumberfalse +\DeclareOptionBeamer{bigpagenumber}{\beamer@Flip@bigpagenumbertrue} + + +\def\beamer@Flip@truetext{true} % Kludge for true/false so I can use ifx + +\DeclareOptionBeamer{topline}[true]{\def\beamer@Flip@topline{#1}} +\ExecuteOptionsBeamer{topline=false} % Default + + + +% \ExecuteOptionsBeamer{topline=} % Default value is empty + + +\ProcessOptionsBeamer + + +\usepackage{tikz} % For fancy decorations + + +% I should fix these up to be more robust against color changes +\setbeamercolor{section in head/foot}{parent=footerlike} +\setbeamercolor{subsection in head/foot}{parent=footerlike} +\setbeamercolor{author in head/foot}{parent=footerlike} +\setbeamercolor{title in head/foot}{parent=footerlike} + + +% \setbeamercolor{author in head/foot}{parent=section in head/foot} +% \setbeamercolor{title in head/foot}{parent=subsection in head/foot} + +%\setbeamercolor{topbar}{parent=crimsontown} +\setbeamercolor{topbar}{parent=bluetown} + + + + +\usesectionheadtemplate + {\hfill\insertsectionhead} + {\hfill\color{fg!50!bg}\insertsectionhead} + + +%%%%%%%%%%%%%%%%%% +% The footer bar % +%%%%%%%%%%%%%%%%%% + +\defbeamertemplate*{footline}{Flip theme} +{% + \leavevmode% +% OBSERVATION: the "%" symbols inside hbox are all very important here. +% The \hbox enviroment will insert spaces whenever there's whitespace +% Adding a % at the END of each line ensures that any hard returns aren't +% interpreted as white space. This allows the color boxes to be flush against +% one another. If a faint white line appears one could probably do something like +% hspace{-1px} and add an additional pixel-or-two wide beamercolorbox. + +\hbox{% +\begin{beamercolorbox}[wd=.4\paperwidth,ht=2.5ex,dp=1.125ex,leftskip=.3cm, rightskip=.3cm plus1fil]{author in head/foot}% + \usebeamerfont{author in head/foot}\insertshortauthor \end{beamercolorbox}% +% +\begin{beamercolorbox}[wd=.6\paperwidth,ht=2.5ex,dp=1.125ex,leftskip=.3cm,rightskip=.3cm plus1fil]{title in head/foot}% + \usebeamerfont{title in head/foot}{\textit \insertshorttitle}% + \hskip2ex plus1fill% + % \insertpagenumber\,/\,\insertpresentationendpage %% Want FRAME, not SLIDE number + \insertframenumber/\inserttotalframenumber + \end{beamercolorbox}% +}% + +% This is a complete kludge for the faint white line which appears +% between the beamer color boxes above. (It doesnt appear in presentation mode +% but sometimes looks funny when the pdfs are viewed in a window) +\begin{tikzpicture}[overlay] + \usebeamercolor[bg]{author in head/foot} + \draw[thick] (0,3.625ex) -- (\paperwidth,3.625ex); + \draw (.4\paperwidth,0) -- (.4\paperwidth,3.625ex); + \draw[thick] (0,0) -- (\paperwidth,0); +\end{tikzpicture} + + + +\ifbeamer@Flip@bigpagenumber + \begin{tikzpicture}[overlay] + \usebeamercolor[bg]{author in head/foot} + \draw[fill] (\paperwidth-5ex,3.625ex) circle (6ex); + \usebeamercolor[fg]{author in head/foot} + \draw[fill] (\paperwidth-5ex,3.625ex) circle (5ex); + \usebeamercolor[bg]{author in head/foot} + \node at (\paperwidth-5ex,4ex) {\small$^\text{\insertframenumber}/_\text{\inserttotalframenumber}$}; + \end{tikzpicture} +\fi +} + + + +\defbeamertemplate*{frametitle}{Flip theme}[1][left] +{ +% \ifbeamercolorempty[bg]{frametitle}{}{\nointerlineskip}% +\@tempdima=\textwidth% +\advance\@tempdima by\beamer@leftmargin% +\advance\@tempdima by\beamer@rightmargin% + +\vbox{}\vskip-3.5ex% +% +\begin{beamercolorbox}[sep=0.3cm,#1,wd=\the\@tempdima]{titlelike} + \usebeamerfont{frametitle}% + \vbox{}\vskip-1ex% + \if@tempswa\else\csname beamer@fte#1\endcsname\fi% + % \strut{\bf\insertframetitle}\strut\par% %% Include the \bf + \strut{\insertframetitle}\strut\par% + {% + \ifx\insertframesubtitle\@empty% + \else + {\usebeamerfont{framesubtitle} + \usebeamercolor[fg]{framesubtitle} + \insertframesubtitle\strut\par}% + \fi + }% + % + \vskip-1ex% + \if@tempswa\else\vskip-.3cm\fi% set inside beamercolorbox... evil here... +\end{beamercolorbox}% +%% + +% \def\beamer@fteright{\vskip0.35cm\advance\leftskip by 1.7cm\advance\rightskip by1.7cm} +} + + + + +% Sidebar right is important: anything placed here will be +% visible *behind* the main text. This is where to put any +% interesting watermarks. +% +% +\defbeamertemplate*{sidebar right}{Flip theme} +{ +% \llap{\includegraphics[width=50px]{footdecoration}} +% \llap{\includegraphics[width=\paperwidth,height=\paperheight]{upperleft_watermark_transp}} +% +% THIS WOULD BE A NICE WAY TO PUT THE PARTICLE BUBBLE CHAMBER BG +% ... make a png of the background with TRANSPARENT BG +% ... then put the image here (use \hfill to get it on the left) +% Now one can FURTHER apply a slight gradient in the BG +% +% +% % Top Line +\ifx\beamer@Flip@topline\beamer@Flip@truetext + \begin{tikzpicture}[overlay] + \usebeamercolor[fg]{titlelike} % alternately, topbarlike would be a good beamercolor +% \usebeamercolor[fg]{topbarlike} % alternately, topbarlike would be a good beamercolor + \draw[fill] (-\paperwidth,0) rectangle (0,-3pt); + \end{tikzpicture} +\fi +% +% + \vskip.2\beamer@leftmargin% + %\llap{\insertlogo\hskip.5\beamer@leftmargin}% + \llap{\insertlogo\hskip.08\beamer@leftmargin}% + \vfill% + \if \beamer@Flip@watermark\@empty + % This is kind of a kludge since if I put anything in here, I automatically get errors + % Somehow LaTeX wants to read the \else part! + \else% + \llap{\includegraphics[width=\paperwidth,height=\paperheight]{\beamer@Flip@watermark}} + \fi +} + + + +\mode +\endinput + + diff --git a/Lectures_my/NumMet/2016/Lecture7/beamerouterthemedecolines.sty b/Lectures_my/NumMet/2016/Lecture7/beamerouterthemedecolines.sty new file mode 100755 index 0000000..669e493 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/beamerouterthemedecolines.sty @@ -0,0 +1,168 @@ +% Copyright 2007 by Marco Barisione +% +% This file may be distributed and/or modified +% +% 1. under the LaTeX Project Public License and/or +% 2. under the GNU Public License. + +\mode + +% String used between the current page and the total page count. +\def\beamer@decolines@pageofpages{/} +\DeclareOptionBeamer{pageofpages}{\def\beamer@decolines@pageofpages{#1}} + +% Show a line below the frame title. +\DeclareOptionBeamer{titleline}[true]{\def\beamer@decolines@titleline{#1}} + +% Image used for the watermark. +\def\beamer@decolines@watermarkorig{} +\DeclareOptionBeamer{watermark}{\def\beamer@decolines@watermarkorig{#1}} + +% Height of the watermark. +\def\beamer@decolines@watermarkheight{100px} +\DeclareOptionBeamer{watermarkheight}{\def\beamer@decolines@watermarkheight{#1}} + +% The original image height is watermarkheightmult * watermarkheight. +\def\beamer@decolines@watermarkheightmult{1} +\DeclareOptionBeamer{watermarkheightmult}{\def\beamer@decolines@watermarkheightmult{#1}} + +\ExecuteOptionsBeamer{titleline=false} +\ProcessOptionsBeamer + +% Enable/disable the watermark. +\def\watermarkon{% + \def\beamer@decolines@watermark{\beamer@decolines@watermarkorig}% +} +\def\watermarkoff{\def\beamer@decolines@watermark{}} + +% Initially enable the watermark. +\watermarkon + +% Colors. +\setbeamercolor*{lineup}{parent=palette primary} +\setbeamercolor*{linemid}{parent=palette secondary} +\setbeamercolor*{linebottom}{parent=palette tertiary} +\setbeamercolor*{page header}{parent=titlelike} + +% Lengths +\newlength{\headerheight} +\setlength{\headerheight}{.045\paperheight} +\newlength{\beamer@decolines@lineup} +\setlength{\beamer@decolines@lineup}{.025\paperheight} +\newlength{\beamer@decolines@linemid} +\setlength{\beamer@decolines@linemid}{.015\paperheight} +\newlength{\beamer@decolines@linebottom} +\setlength{\beamer@decolines@linebottom}{.01\paperheight} + +% The height of the watermark part below the 3 bottom lines. +\newlength{\beamer@decolines@watermarkheightbottom} +\addtolength{\beamer@decolines@watermarkheightbottom}{\beamer@decolines@lineup} +\addtolength{\beamer@decolines@watermarkheightbottom}{\beamer@decolines@linemid} +\addtolength{\beamer@decolines@watermarkheightbottom}{\beamer@decolines@linebottom} + +% The height of the watermark part over the 3 bottom lines before shrinking. +\newlength{\beamer@decolines@watermarkheightupperorig} +\setlength{\beamer@decolines@watermarkheightupperorig}{\beamer@decolines@watermarkheight} +\addtolength{\beamer@decolines@watermarkheightupperorig}{-\beamer@decolines@watermarkheightbottom} +\multiply\beamer@decolines@watermarkheightupperorig by \beamer@decolines@watermarkheightmult + +% % Footer. +% \defbeamertemplate*{footline}{decolines theme} +% { +% \leavevmode% +% % Page number. +% \hbox{% +% \begin{beamercolorbox}[wd=.2\paperwidth,ht=0ex,dp=0ex,center]{}% +% \usebeamerfont{palette primary}\insertframenumber{} \beamer@decolines@pageofpages{} \inserttotalframenumber% +% \end{beamercolorbox}% +% \begin{beamercolorbox}[wd=.8\paperwidth,ht=0ex,dp=0ex]{}% +% \end{beamercolorbox}% +% } % +% % First line. +% \hbox{% +% \begin{beamercolorbox}[wd=.2\paperwidth,ht=\beamer@decolines@lineup,dp=0pt]{}% +% \end{beamercolorbox}% +% \begin{beamercolorbox}[wd=.8\paperwidth,ht=\beamer@decolines@lineup,dp=0pt]{lineup}% +% \end{beamercolorbox}% +% } % +% % Second line. +% \hbox{% +% \begin{beamercolorbox}[wd=\paperwidth,ht=\beamer@decolines@linemid,dp=0pt]{linemid}% +% \end{beamercolorbox}% +% } % +% % Third line. +% \hbox{% +% \begin{beamercolorbox}[wd=.1\paperwidth,ht=\beamer@decolines@linebottom,dp=0pt]{}% +% \end{beamercolorbox}% +% \begin{beamercolorbox}[wd=.9\paperwidth,ht=\beamer@decolines@linebottom,dp=0pt]{linebottom}% +% \end{beamercolorbox}% +% }% +% % This seems to fix some alignment problems with the watermark. It has to be +% % always applied if you do not want to see the footer moving up and down when +% % moving from a page with watermark to a page without or vice versa. +% \vskip-.5px% +% % Watermark. +% \if\beamer@decolines@watermark\@empty\else% +% \vskip-\beamer@decolines@watermarkheightbottom% +% \llap{\includegraphics[height=\beamer@decolines@watermarkheightbottom,clip=true,% +% trim=0pt 0pt 0pt \beamer@decolines@watermarkheightupperorig]{\beamer@decolines@watermark}\hskip-\paperwidth}% +% \fi% +% } + +\defbeamertemplate*{headline}{decolines theme} +{ + \leavevmode% + \hbox{% + \begin{beamercolorbox}[wd=\paperwidth,ht=\headerheight,dp=0pt]{page header}% + \end{beamercolorbox}% + } % + \vskip0pt% +} + +\defbeamertemplate*{frametitle}{decolines theme}[1][left] +{ + \ifbeamercolorempty[bg]{frametitle}{}{\nointerlineskip}% + \@tempdima=\textwidth% + \advance\@tempdima by\beamer@leftmargin% + \advance\@tempdima by\beamer@rightmargin% + \vbox{}\vskip-.5\beamer@leftmargin% + \begin{beamercolorbox}[sep=\beamer@leftmargin,#1,wd=\the\@tempdima]{} + \usebeamerfont{frametitle}\usebeamercolor[bg]{framesubtitle}% + \vbox{}\vskip0ex% + \if@tempswa\else\csname beamer@fte#1\endcsname\fi% + \strut\insertframetitle\strut\par% + {% + \ifx\insertframesubtitle\@empty% + \else% + {\usebeamerfont{framesubtitle}\usebeamercolor[bg]{framesubtitle}\insertframesubtitle\strut\par}% + \fi + }% + \vskip-1ex% + \if@tempswa\else\vskip-\beamer@leftmargin\fi + \end{beamercolorbox}% + \def\beamer@decolines@truetext{true}% + \ifx\beamer@decolines@titleline\beamer@decolines@truetext% + \vskip-.5\beamer@leftmargin% + \begin{beamercolorbox}[wd=\textwidth,ht=.1ex,dp=0ex]{linemid}% + \end{beamercolorbox}% + \fi +} + +% Frame title continuations, default +\defbeamertemplate*{frametitle continuation}{decolines theme}{(\insertcontinuationcount)} + + + +\defbeamertemplate*{sidebar right}{decolines theme} +{ + \vskip.1\beamer@leftmargin% + \llap{\insertlogo\hskip.5\beamer@leftmargin}% + \vfill% + \if\beamer@decolines@watermark\@empty\else% + \llap{\includegraphics[height=\beamer@decolines@watermarkheight]{\beamer@decolines@watermark}}% + \vskip-\beamer@decolines@watermarkheightbottom% + \fi +} + +\mode + diff --git a/Lectures_my/NumMet/2016/Lecture7/beamerthemeFlip.sty b/Lectures_my/NumMet/2016/Lecture7/beamerthemeFlip.sty new file mode 100755 index 0000000..a615ff5 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/beamerthemeFlip.sty @@ -0,0 +1,59 @@ +% Main style file for the BEAMER FLIP THEME 2012 +% Copyright 2012 by Flip Tanedo +% This file may be distributed and/or modified +% 1. under the LaTeX Project Public License and/or +% 2. under the GNU Public License. +% +% Comments: This is still a work in progress. +% A good way to present: http://www.cs.hmc.edu/~oneill/freesoftware/pdftokeynote.html + +% \usepackage{beamerthemesplit} +\mode + + +%% These are the options that are fed in through the driver file +%% Some of them get passed on to the other theme files +\DeclareOptionBeamer{bullet}{\PassOptionsToPackage{bullet=#1}{beamerinnerthemeFlip}} +\DeclareOptionBeamer{bigpagenumber}{\PassOptionsToPackage{bigpagenumber}{beamerouterthemeFlip}} +\DeclareOptionBeamer{topline}[true]{\PassOptionsToPackage{topline=#1}{beamerouterthemeFlip}} +\DeclareOptionBeamer{shadow}[false]{\PassOptionsToPackage{shadow=#1}{beamerinnerthemeFlip}} + +\DeclareOptionBeamer{watermark}{\PassOptionsToPackage{watermark=#1}{beamerouterthemeFlip}} + +\ProcessOptionsBeamer + +\useinnertheme{Flip} % Calls beamerinnerthemeFlip.sty +\useoutertheme{Flip} % Calls beamerouterthemeFlip.sty +\usecolortheme{Flip} % Calls beamercolorthemeFlip.sty + + + +%%%%%%%%%%%%%%%%%%%%%%%% +% Background Gradient % +%%%%%%%%%%%%%%%%%%%%%%%% + + +%% Usually I would leave this on +%% However, there seems to be problems with XeLaTeX +%% http://tex.stackexchange.com/questions/29497/xelatex-preventing-beamer-from-using-different-backgrounds +\setbeamertemplate{background canvas}[vertical shading][bottom=blue!.1, middle=white, top=white] + +%% HOWEVER: it seems to make it impossible to change this later on in the document + + +%% For a dark background: +% \setbeamertemplate{background canvas}[vertical shading][bottom=keynotebottom, middle=keynotemiddle, top=keynotetop] + + +%% An alternate way: \beamertemplateshadingbackground{blue!.1}{red!2} +%% In general: you want the gradient to be subtle! + + +%%%%%%%%%%%%%%%%%%%%%%% +% Navigation symbols % +%%%%%%%%%%%%%%%%%%%%%%% + +\setbeamertemplate{navigation symbols}{} % Turns off PDF navigation symbols + +\mode +\endinput \ No newline at end of file diff --git a/Lectures_my/NumMet/2016/Lecture7/images/BG_lower.png b/Lectures_my/NumMet/2016/Lecture7/images/BG_lower.png new file mode 100755 index 0000000..073034a --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/images/BG_lower.png Binary files differ diff --git a/Lectures_my/NumMet/2016/Lecture7/images/bubble.gif b/Lectures_my/NumMet/2016/Lecture7/images/bubble.gif new file mode 100644 index 0000000..ecff282 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/images/bubble.gif Binary files differ diff --git a/Lectures_my/NumMet/2016/Lecture7/images/bubble.jpg b/Lectures_my/NumMet/2016/Lecture7/images/bubble.jpg new file mode 100644 index 0000000..ceb584d --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/images/bubble.jpg Binary files differ diff --git a/Lectures_my/NumMet/2016/Lecture7/images/bubble2.jpg b/Lectures_my/NumMet/2016/Lecture7/images/bubble2.jpg new file mode 100644 index 0000000..0fd03bb --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/images/bubble2.jpg Binary files differ diff --git a/Lectures_my/NumMet/2016/Lecture7/images/runge-int.pdf b/Lectures_my/NumMet/2016/Lecture7/images/runge-int.pdf new file mode 100644 index 0000000..1c3c0b4 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/images/runge-int.pdf Binary files differ diff --git a/Lectures_my/NumMet/2016/Lecture7/images/uzh-transp.pdf b/Lectures_my/NumMet/2016/Lecture7/images/uzh-transp.pdf new file mode 100644 index 0000000..977da67 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/images/uzh-transp.pdf Binary files differ diff --git a/Lectures_my/NumMet/2016/Lecture7/lecture7.tex b/Lectures_my/NumMet/2016/Lecture7/lecture7.tex new file mode 100644 index 0000000..612d9fd --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/lecture7.tex @@ -0,0 +1,315 @@ +% vim: set sts=4 et : +\documentclass[11 pt,xcolor={dvipsnames,svgnames,x11names,table}]{beamer} + +\usetheme[ + bullet=circle, % Other option: square + bigpagenumber, % circled page number on lower right + topline=true, % colored bar at the top of the frame + shadow=false, % Shading for beamer blocks + watermark=BG_lower, % png file for the watermark + ]{Flip} + +\usepackage{bm} +\usepackage{tikz,pgfplots} +\usepgfplotslibrary{fillbetween} +\usepackage{graphicx} + +% SOME COMMANDS THAT I FIND HANDY +% \renewcommand{\tilde}{\widetilde} % dinky tildes look silly, dosn't work with fontspec +\newcommand{\comment}[1]{\textcolor{comment}{\footnotesize{#1}\normalsize}} % comment mild +\newcommand{\Comment}[1]{\textcolor{Comment}{\footnotesize{#1}\normalsize}} % comment bold +\newcommand{\COMMENT}[1]{\textcolor{COMMENT}{\footnotesize{#1}\normalsize}} % comment crazy bold +\newcommand{\Alert}[1]{\textcolor{Alert}{#1}} % louder alert +\newcommand{\ALERT}[1]{\textcolor{ALERT}{#1}} % loudest alert +%% "\alert" is already a beamer pre-defined +\newcommand*{\Scale}[2][4]{\scalebox{#1}{$#2$}}% + +\graphicspath{{images/}} % Put all images in this directory. Avoids clutter. + +% suppress frame numbering for backup slides +% you always need the appendix for this! +\newcommand{\backupbegin}{ + \newcounter{framenumberappendix} + \setcounter{framenumberappendix}{\value{framenumber}} +} +\newcommand{\backupend}{ + \addtocounter{framenumberappendix}{-\value{framenumber}} + \addtocounter{framenumber}{\value{framenumberappendix}} +} + +\begin{document} +\tikzstyle{every picture}+=[remember picture] + +{ +\setbeamertemplate{sidebar right}{\llap{\includegraphics[width=\paperwidth,height=\paperheight]{bubble2}}} +\begin{frame}[c]%{\phantom{title page}} +\begin{center} +\begin{center} + \begin{columns} + \begin{column}{0.9\textwidth} + \flushright%\fontspec{Trebuchet MS} + \bfseries \Huge {Numerical Integration} + \end{column} + \begin{column}{0.2\textwidth} + %\includegraphics[width=\textwidth]{SHiP-2} + \end{column} + \end{columns} +\end{center} + \quad + \vspace{3em} +\begin{columns} +\begin{column}{0.6\textwidth} +\flushright \vspace{-1.8em} {%\fontspec{Trebuchet MS} + \large Marcin Chrzaszcz, Danny van Dyk\\\vspace{-0.1em}\small \href{mailto:mchrzasz@cern.ch}{mchrzasz@cern.ch}, \href{mailto:danny.van.dyk@gmail.com}{danny.van.dyk@gmail.com}} + +\end{column} +\begin{column}{0.4\textwidth} +\includegraphics[height=1.3cm]{uzh-transp} +\end{column} +\end{columns} + +\vspace{1em} +\vspace{0.5em} + \textcolor{normal text.fg!50!Comment}{Numerical Methods, \\ 26. September, 2016} +\end{center} +\end{frame} +} + +\begin{frame}{Plan for today} + \begin{itemize} + \item \alert{General problem of estimating the integral of a continuous function $f(x)$ on a finite support}\\ + \vfill + \item \alert{Specific problems in which properties of the integrand can be used to our advantage}\\ + What properties of $f(x)$ can we use make our live easier? + \end{itemize} +\end{frame} + +\begin{frame}{Setup} +\end{frame} + +\begin{frame}{Constant approximation} +\end{frame} + +\begin{frame}{Linear approximation} +\end{frame} + +\begin{frame}{Wait a minute, I know this!} + Integrate the interpolating polynomial! + + Newton-Coates! + + \begin{equation} + \sum_k \omega_k = 1\,. + \end{equation} +\end{frame} + +\begin{frame}{Pathological example 1} + \resizebox{.9\textwidth}{!}{ + \begin{tikzpicture} + \begin{axis}[% + axis x line=center, + axis y line=left, + ymin=0.0,ymax=1.05, + xmin=-1,xmax=+1, + samples=160 + ] + \addplot[thick,black,domain=-1:0] { + 1 / (1 + 25 * x * x) + }; + \addplot[thick,black,domain=0:+1] { + 1 / (1 + 25 * x * x) + }; + \addplot[thick,red, domain=-1:0] { + 1.0 - 16.8552 * x^2 + 123.36 * x^4 - 381.434 * x^6 + 494.91 * x^8 - 220.942 * x^(10) + }; + \addplot[thick,red, domain=0:+1] { + 1.0 - 16.8552 * x^2 + 123.36 * x^4 - 381.434 * x^6 + 494.91 * x^8 - 220.942 * x^(10) + }; + \addplot[thick,red,only marks,mark=*,mark size=3pt] coordinates { + (-1.0, 0.0384615) + (-0.8, 0.0588235) + (-0.6, 0.1) + (-0.4, 0.2) + (-0.2, 0.5) + ( 0.0, 1.0) + ( 0.2, 0.5) + ( 0.4, 0.2) + ( 0.6, 0.1) + ( 0.8, 0.0588235) + ( 1.0, 0.0384615) + }; + \end{axis} + \end{tikzpicture} + } +\end{frame} + +\begin{frame}{Conclusion 1: We know how to fix it!} + \only<1>{ + When interpolating a function, we saw that using splines fixes + the problem of oscillating interpolating polynomials.\\ + \vfill + + We can now use this to stabilise our interpolation formula as well.\\ + \vfill + + \begin{itemize} + \item Let $M$ be the number of intervals in which we want to integrate. + \item Let $N$ be the degree of the interpolating polynomial in each interval. + \item Enforce that the interpolating function is continous, but not differentiable + at the interval boundaries. + \end{itemize} + \vfill + The integral can then be approximated as + \begin{equation} + I \approx \sum_{k=0}^{N \cdot M + 1} \omega_k f(x_k) + \end{equation} + with $x_k = a + (b - a) \frac{k}{N\cdot M}$\,. + } + \only<2->{ + \resizebox{.9\textwidth}{!}{ + \begin{tikzpicture} + \begin{axis}[% + axis x line=center, + axis y line=left, + ymin=0.0,ymax=1.05, + xmin=-1,xmax=+1, + samples=160 + ] + \addplot[thick,black,domain=-1:0] { + 1 / (1 + 25 * x * x) + }; + \addplot[thick,black,domain=0:+1] { + 1 / (1 + 25 * x * x) + }; + \only<3>{ + \addplot+[smooth,red,name path=A1,domain=-1:-0.6,mark=none] { + 0.348416 + 0.570136 * x + 0.260181 * x^2 + }; + \addplot+[draw=none,mark=none,name path=B1] {0}; + \addplot+[gray] fill between[of=A1 and B1,soft clip={domain=-1:-0.6}]; + } + \only<4>{ + \addplot+[smooth,red,name path=A2,domain=-0.6:-0.2,mark=none] { + 1. + 3. * x + 2.5 * x^2 + }; + \addplot+[draw=none,mark=none,name path=B2] {0}; + \addplot+[gray] fill between[of=A2 and B2,soft clip={domain=-0.6:-0.2}]; + } + \only<5>{ + \addplot+[smooth,red,name path=A3,domain=-0.2:+0.2,mark=none] { + 1. - 12.5 * x^2 + }; + \addplot+[draw=none,mark=none,name path=B3] {0}; + \addplot+[gray] fill between[of=A3 and B3,soft clip={domain=-0.2:+0.2}]; + } + \only<6>{ + \addplot+[smooth,red,name path=A4,domain=+0.2:+0.6,mark=none] { + 1. - 3. * x + 2.5 * x^2 + }; + \addplot+[draw=none,mark=none,name path=B4] {0}; + \addplot+[gray] fill between[of=A4 and B4,soft clip={domain=+0.2:0.6}]; + } + \only<7>{ + \addplot+[smooth,red,name path=A5,domain=+0.6:+1,mark=none] { + 0.348416 - 0.570136 * x + 0.260181 * x^2 + }; + \addplot+[draw=none,mark=none,name path=B5,domain=+0.6:+1] {0}; + \addplot+[gray] fill between[of=A5 and B5,soft clip={domain=0.6:1}]; + } + \addplot[thick,red,only marks,mark=*,mark size=3pt] coordinates { + (-1.0, 0.0384615) + (-0.8, 0.0588235) + (-0.6, 0.1) + (-0.4, 0.2) + (-0.2, 0.5) + ( 0.0, 1.0) + ( 0.2, 0.5) + ( 0.4, 0.2) + ( 0.6, 0.1) + ( 0.8, 0.0588235) + ( 1.0, 0.0384615) + }; + \end{axis} + \end{tikzpicture} + } + } +\end{frame} +\begin{frame}{} + \begin{equation} + \omega_k = \begin{cases} + \frac{1}{6} & k = 0 \lor k = N \cdot M + 1\\ + \frac{4}{6} & k\,\text{even}\\ + \frac{2}{6} & k\,\text{odd} + \end{cases} + \end{equation} + +\end{frame} + +\newcommand{\eps}{\varepsilon} +\begin{frame}{Pathological example 2} + \begin{equation*} + I \equiv \int_{-3}^{+4} \mathrm{d}x\, \left[1 - \frac{\eps}{x^2 + \eps}\right] + \, \quad \text{with }\eps = \frac{1}{9} + \end{equation*} + \begin{columns} + \begin{column}{.5\textwidth} + \resizebox{\textwidth}{!}{ + \begin{tikzpicture} + \begin{axis}[% + axis x line=center, + axis y line=left, + ymin=0.0,ymax=1.05, + xmin=-3,xmax=+4, + samples=160 + ] + \addplot[thick,black,domain=-3:0] { + 1 - 1/9 / (x^2 + 1/9) + }; + \addplot[thick,black,domain=0:+4] { + 1 - 1/9 / (x^2 + 1/9) + }; + \addplot[thick,blue,only marks,mark=+,mark size=5pt] coordinates { + ( 1, 0.54030 ) + }; + \end{axis} + \end{tikzpicture} + } + \begin{equation*} + \end{equation*} + \end{column} + \begin{column}{.5\textwidth} + \begin{block}{Exact result} + ~\\ + \vspace{-2\bigskipamount} + \begin{align*} + I & = 5.57396 + \end{align*} + \end{block} + \begin{block}{Entire range, 7 data points} + ~\\ + \vspace{-2\bigskipamount} + \begin{align*} + I_7 & = 6.25746 & + \delta_7 & = 12.26\% + \end{align*} + \end{block} + \begin{block}{Two integration, 5 data points} + ~\\ + \vspace{-2\bigskipamount} + \begin{align*} + I_5 & = 5.65645 & + \delta_5 & = 1.47\% + \end{align*} + \end{block} + \end{column} + \end{columns} +\end{frame} + +\backupbegin + +\begin{frame}\frametitle{Backup} + +\end{frame} + +\backupend +\end{document} diff --git a/Lectures_my/NumMet/2016/Lecture7/runge-int.nb b/Lectures_my/NumMet/2016/Lecture7/runge-int.nb new file mode 100644 index 0000000..200574f --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/runge-int.nb @@ -0,0 +1,523 @@ +(* Content-type: application/vnd.wolfram.mathematica *) + +(*** Wolfram Notebook File ***) +(* http://www.wolfram.com/nb *) + +(* CreatedBy='Mathematica 10.3' *) + +(*CacheID: 234*) +(* Internal cache information: +NotebookFileLineBreakTest +NotebookFileLineBreakTest +NotebookDataPosition[ 158, 7] +NotebookDataLength[ 20735, 514] +NotebookOptionsPosition[ 19551, 469] +NotebookOutlinePosition[ 19888, 484] +CellTagsIndexPosition[ 19845, 481] +WindowFrame->Normal*) + +(* Beginning of Notebook Content *) +Notebook[{ + +Cell[CellGroupData[{ +Cell[BoxData[ + RowBox[{ + RowBox[{"runge", "[", "x_", "]"}], " ", "=", " ", + RowBox[{"1", " ", "/", " ", + RowBox[{"(", + RowBox[{"1", " ", "+", " ", + RowBox[{"25", " ", + RowBox[{"x", "^", "2"}]}]}], ")"}]}]}]], "Input", + CellChangeTimes->{{3.6855202639388638`*^9, 3.6855202996308117`*^9}}], + +Cell[BoxData[ + FractionBox["1", + RowBox[{"1", "+", + RowBox[{"25", " ", + SuperscriptBox["x", "2"]}]}]]], "Output", + CellChangeTimes->{{3.6855202971715717`*^9, 3.685520299998887*^9}, + 3.685529753291757*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"rungeData11", " ", "=", " ", + RowBox[{"Table", "[", "\[IndentingNewLine]", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"runge", "[", "x", "]"}]}], "}"}], ",", "\[IndentingNewLine]", + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"-", "1"}], ",", " ", + RowBox[{"+", "1"}], ",", " ", "0.2"}], "}"}]}], "\[IndentingNewLine]", + "]"}]}]], "Input", + CellChangeTimes->{{3.685520443605824*^9, 3.6855204506374273`*^9}}], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"-", "1.`"}], ",", "0.038461538461538464`"}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.8`"}], ",", "0.05882352941176469`"}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.6`"}], ",", "0.1`"}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.3999999999999999`"}], ",", "0.20000000000000007`"}], + "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.19999999999999996`"}], ",", "0.5000000000000001`"}], + "}"}], ",", + RowBox[{"{", + RowBox[{"0.`", ",", "1.`"}], "}"}], ",", + RowBox[{"{", + RowBox[{"0.20000000000000018`", ",", "0.49999999999999956`"}], "}"}], ",", + RowBox[{"{", + RowBox[{"0.40000000000000013`", ",", "0.1999999999999999`"}], "}"}], ",", + + RowBox[{"{", + RowBox[{"0.6000000000000001`", ",", "0.09999999999999998`"}], "}"}], ",", + + RowBox[{"{", + RowBox[{"0.8`", ",", "0.05882352941176469`"}], "}"}], ",", + RowBox[{"{", + RowBox[{"1.`", ",", "0.038461538461538464`"}], "}"}]}], "}"}]], "Output", + CellChangeTimes->{{3.685520443884881*^9, 3.6855204591362753`*^9}, + 3.685529753470796*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"rungeApprox11", " ", "=", " ", + RowBox[{ + RowBox[{"InterpolatingPolynomial", "[", + RowBox[{"rungeData11", ",", " ", "x"}], "]"}], " ", "//", " ", + "Simplify"}]}]], "Input", + CellChangeTimes->{{3.6855204739097023`*^9, 3.6855204741090593`*^9}, { + 3.685529748047958*^9, 3.685529757823937*^9}}], + +Cell[BoxData[ + RowBox[{"1.`", "\[VeryThinSpace]", "+", + RowBox[{"1.3322676295501878`*^-15", " ", "x"}], "-", + RowBox[{"16.855203619909503`", " ", + SuperscriptBox["x", "2"]}], "-", + RowBox[{"3.019806626980426`*^-14", " ", + SuperscriptBox["x", "3"]}], "+", + RowBox[{"123.35972850678728`", " ", + SuperscriptBox["x", "4"]}], "+", + RowBox[{"1.4210854715202004`*^-13", " ", + SuperscriptBox["x", "5"]}], "-", + RowBox[{"381.43382352941154`", " ", + SuperscriptBox["x", "6"]}], "-", + RowBox[{"2.2737367544323206`*^-13", " ", + SuperscriptBox["x", "7"]}], "+", + RowBox[{"494.9095022624431`", " ", + SuperscriptBox["x", "8"]}], "+", + RowBox[{"1.1368683772161603`*^-13", " ", + SuperscriptBox["x", "9"]}], "-", + RowBox[{"220.94174208144779`", " ", + SuperscriptBox["x", "10"]}]}]], "Output", + CellChangeTimes->{3.685530508856052*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"rungePlot", " ", "=", " ", + RowBox[{"Show", "[", "\[IndentingNewLine]", + RowBox[{ + RowBox[{"Plot", "[", + RowBox[{ + RowBox[{"runge", "[", "x", "]"}], ",", " ", + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"-", "1"}], ",", " ", + RowBox[{"+", "1"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", + RowBox[{"ListPlot", "[", + RowBox[{"rungeData11", ",", " ", + RowBox[{"PlotStyle", " ", "\[Rule]", " ", "Red"}]}], "]"}], ",", + "\[IndentingNewLine]", + RowBox[{"Plot", "[", + RowBox[{"rungeApprox11", ",", " ", + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"-", "1"}], ",", " ", + RowBox[{"+", "1"}]}], "}"}], ",", " ", + RowBox[{"PlotStyle", " ", "\[Rule]", " ", "Red"}]}], "]"}]}], + "\[IndentingNewLine]", "]"}]}]], "Input", + CellChangeTimes->{{3.6855203042714443`*^9, 3.685520320302579*^9}, { + 3.685520372318378*^9, 3.68552039602993*^9}, {3.685520460981647*^9, + 3.6855204931812057`*^9}}], + +Cell[BoxData[ + RowBox[{ + StyleBox[ + RowBox[{"InterpolatingPolynomial", "::", "list"}], "MessageName"], + RowBox[{ + ":", " "}], "\<\"List expected at position \[NoBreak]\\!\\(\\*RowBox[{\\\"1\ +\\\"}]\\)\[NoBreak] in \[NoBreak]\\!\\(\\*RowBox[{\\\"InterpolatingPolynomial\ +\\\", \\\"[\\\", RowBox[{\\\"rungeData\\\", \\\",\\\", RowBox[{\\\"-\\\", \ +\\\"0.9999591428571428`\\\"}]}], \\\"]\\\"}]\\)\[NoBreak]. \ +\\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ +ButtonFrame->None, ButtonData:>\\\"paclet:ref/message/General/list\\\", \ +ButtonNote -> \\\"InterpolatingPolynomial::list\\\"]\\)\"\>"}]], "Message", \ +"MSG", + CellChangeTimes->{3.6855297537585983`*^9}], + +Cell[BoxData[ + RowBox[{ + StyleBox[ + RowBox[{"InterpolatingPolynomial", "::", "list"}], "MessageName"], + RowBox[{ + ":", " "}], "\<\"List expected at position \[NoBreak]\\!\\(\\*RowBox[{\\\"1\ +\\\"}]\\)\[NoBreak] in \[NoBreak]\\!\\(\\*RowBox[{\\\"InterpolatingPolynomial\ +\\\", \\\"[\\\", RowBox[{\\\"rungeData\\\", \\\",\\\", RowBox[{\\\"-\\\", \ +\\\"0.9591428163265305`\\\"}]}], \\\"]\\\"}]\\)\[NoBreak]. \ +\\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ +ButtonFrame->None, ButtonData:>\\\"paclet:ref/message/General/list\\\", \ +ButtonNote -> \\\"InterpolatingPolynomial::list\\\"]\\)\"\>"}]], "Message", \ +"MSG", + CellChangeTimes->{3.685529753818915*^9}], + +Cell[BoxData[ + RowBox[{ + StyleBox[ + RowBox[{"InterpolatingPolynomial", "::", "list"}], "MessageName"], + RowBox[{ + ":", " "}], "\<\"List expected at position \[NoBreak]\\!\\(\\*RowBox[{\\\"1\ +\\\"}]\\)\[NoBreak] in \[NoBreak]\\!\\(\\*RowBox[{\\\"InterpolatingPolynomial\ +\\\", \\\"[\\\", RowBox[{\\\"rungeData\\\", \\\",\\\", RowBox[{\\\"-\\\", \ +\\\"0.9183264897959182`\\\"}]}], \\\"]\\\"}]\\)\[NoBreak]. \ +\\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ +ButtonFrame->None, ButtonData:>\\\"paclet:ref/message/General/list\\\", \ +ButtonNote -> \\\"InterpolatingPolynomial::list\\\"]\\)\"\>"}]], "Message", \ +"MSG", + CellChangeTimes->{3.685529753874078*^9}], + +Cell[BoxData[ + RowBox[{ + StyleBox[ + RowBox[{"General", "::", "stop"}], "MessageName"], + RowBox[{ + ":", " "}], "\<\"Further output of \ +\[NoBreak]\\!\\(\\*StyleBox[RowBox[{\\\"InterpolatingPolynomial\\\", \\\"::\\\ +\", \\\"list\\\"}], \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed \ +during this calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ +ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ +ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \ +\\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG", + CellChangeTimes->{3.685529753881908*^9}], + +Cell[BoxData[ + GraphicsBox[{{{}, {}, + {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[ + 1.], LineBox[CompressedData[" +1:eJw1mHk01F/YwJEkZUklW4WsKdrs8fySrCmJEipkrUQqImuWKFtIdmmGSEmy +h2snS4xhjG3M2LLO114k3nnPed+/7vmc+5xz733OXT7PFbVxMbZjYWJismRm +YvrfVg74pre2MMSxmX5c7SsOTO5Kha6tY2jRo7pfrgQHXm+UJZdWMOQ31XVS +ohIHTTPmNhMzGNqa2Dcm8AMHVgnp/W29GCrGQfARGg5eL0i3Jn7BUJNe+OGy +PXhYxql9On0bQ+5KNcJrj/FQzm79wKGKjl6b8YTKamWCXztppsZ6Domf3HNd +9uQHmMoUOKAkPIs+cqTfuqeWA0eHdXdyD02jd0yKV+l+uYDfi1vuez6F8NE3 +pyqM8qDOoiUCwSSSCNp1npibD8LpwiIJtAnEd/ar84nYAjA4b+FJDBtHrbm2 +XcwNhXD4r72/jN4YUnSXUXt3rARmexzw5efHkKtA8U26UgmU5js2G2iMob2S +tgFa50vAyP4uz4PTYyju5uUpdvMS8Cc4vys4OIa4OAx3Lr0ogZEPj6vVFkdR +5/LPvVfHSyDL5DnTpeRRJH7Dy+xpVikcz0v1d5sbQRf6SVPiWuWgW2nKbTU5 +glKO111qNSoHmzbONMPRERR45KOf181ySJj2+S7dN4L+YV9sl9zLgVXq1iql +fgR18NjwG30sh/70Q/cNUkZQ2X/i0hW83yHkdfp1cYMRJMRMqkiY/g7Uxxly +pGwaUjK6fceotBJ8FZ7oDuFoqPXSUrFIUyUIr+rajKXRkCtnJ+5vTyWYecy/ +WYyjoZ/Wwy6EpUogeMI/rgAaSsFWqP/kq6DWd6hV5wYN1ZqpN9p8rAJ8qIBT +2U4aOt8FLbgPCJxSXuNS7lIRl6bGtxOj1bDMrh+mZEdF1YJPT/bOVoP/ExaX +rttU1Az5r4JXqyH+0iNVdlMqGjyvu7y2swZqmUwJj4CKkjwcPDVP1oCgHf+W +/j4q+jdUGHEjoAbeWGgULT4bRvUPf1cOSNWCme+NI5/EKKitzUXgVHQd7KpW +rk/ipyD+0IJ7esl1UMXMbxfGRUHvI7ZSnLLqQCKYlGW/PoRa3EZQfUUdLL66 +elS0awgpegpsTUzVQXii4cn4gCGkZRCqPXqhHqoLz4E/dRCdZ8fJbbA2gPSM +jIVJ2gAqC+nEDWY0QlFQV+zd2AGUaO3+PCW/ETQPPWvzDx1AbHciGp1QI1he +aVP/7DaAptTPlR8baoTokgcibDoDyJzIvmgq0ARrgd/GS7B+FCKI5jXjmqBF +WN1V6Fw/Sv4e4LAS1wzXi8ezTyr2IyWb23/9Mpth7HIkTUe2Hyno4xcEipth +8/mw8eP9/SjgEJ++d28znJoKUGif6kOts7NRG4I/IKmocd03pg9914kZWsf/ +gHuXjIJHRsnozdZ5rQcNLbCpf6Jbto+Mmi8Klef3tkCMDs+RJz/JSFoiRY11 +ugXK/+us3lFGRtF9VZ8oXK3Aftpo41gkGVGiJczv3GiFrANGjzyUyWj/DdOK +ioVWGKFett4d0YucTX79Uj3RDuaPLqsrKpJQdpaZH/tUB7wip1orHyOh/B1v +3NI3O6BCfTZYVYyE7CgCi+IcnXCIPeynBhcJ4Tqx5tnDnUBLqb2tM9GDrraP +iy0ZdIJjo8Jzs/gedLHjrdVoZie4Cwg3ea12I+knsolJ1gT44Ht3xnu2G0Vn +jRG/OhOAPFrK7TfSjbbztAp1exJANc/ULPBnNxr8K7mlEkOAf5pRU6+yutFF +Xxn98DoCBN1n2Z16rRs1/4FKE+kuiEWTRqiYiBR3kkUb1rsgzZuJcyqXiKg8 +qTJeO4mQrcL/gzeDiC7Z6Y2q8xOh6pvOOftXRMRpF56zoECEyQ9ZJzmtiEiU +Rf3EbjciaETb8lrsJKI0n8OlgnNEmLKiElctutBqiPv8Lqwblg/+iRK50oXc +n6j6FG11w2Y/90V97S50rkimyIWnB/aZQn3qiS40Om7twHeqBzT00gq1tnch +LqWtZDn3Hog7afkmJo+AhBzfhyQwkwBYyNfkmAloZ7taUb1EL0xI3Zrk+d2J +vP7zJFgp9UKk4Zjn0mwnGj+eocyp1wuUhPmUUnInkr2hL/z2fi/4yu8c1czv +RNmJEtmuhb1QYanqcv12J9p1tkZJW4cMjx+ocDM/60BFj62quj374C651jVz +qQ3Jfi7KF9k/CM5E67qyoTZ0+MjZgJ3HB8H1JxPfz6Y2ZLnHdID9wiB41Gt8 +X01uQwdMmc0MngxCUH75dl2tNrTgcC6JShqEtNCCpOm4VnTAdq3OJG0IiMrv +6+UVW9DvsoDKJI1hOJsYwF/+tAkdvKwln2VDA0PzqLdn7jShwI2gs2oPaXBL +KJXvi2ETcupoMhzwo4Ffauk+vFgTkjUXv3M2lQa1GRhPZFsj6hMcSwgg00Dn +482dNiKNqPFm3ak9V0bA+LvK+s7meqSdLCu3oTcKjkMLAzcO1KI7xgIGPtfH +4TL12Odm8XJUktb+oNN6Ci7aSlnBWhm6dszaV+/RFOhPiu4tbi9DGm8eQGXQ +FGjN83ni3MtQJ3Iqif0wBcpMzBd8m0vRxnxbyPLsFIiJ9Ayevl+CPJjOYefd +p2Hlts/u9G+F6Ne2NDnMfwYWRz3Q/tBClK9xoFv69QzMO7i5hVsWomdC8TJX +M2Zg2sWB7Lm9EBWHvQx+WjMDFL8reJPr31BU/H7TDaZZaE6XOMvx9ytS4RJj +SvKeBUWU3DsdlIeidF+059+ZA8HdrFNpp/OQXf4//Xuuc/DP7P668chnJKss +HcPlMwf1i2cPfofPiDPyuSlX/BxclaTYvFrPRaWxOe+mG+fANUKELuuag+Qq +FjxJR+hg0h+6NXwoBzku7ND9K0cHJakFnrj2bAQ5ODNWVTpsVtec3jiajaYO +sh2tvUyH8GUbz7bxLKQgI8M07kmHHItMVmcLPLrvvvdleDOjP5uTT5QDj2zm +lHJ7CHRwXXki1VOKQzsqv01s9tNBOUpbX50Phzav4IvmZ+nQWPsripOQgVy4 +L/AZcGPwkdsoo8YvA/Fq/HCMPIBBhGVpwRO5DPTSuP3Pl8MYmKyG9gy9eod2 +qevMh8tj0HjVRNEyIg0tV9v8l23IiJ8Tv1E/nYo+Dqn49Jow4l+sPDumm4qi +HMr1Ri0wGC2Pr9lgSUFRdUYtUU4YMIn1XUx9mogIUb2SSoEYNH3PcdlOSkCi +BdIJF8IwiDT1inE+nYAo/Tai8lEYCIcJkdXp8Wg58xghJgkDlfmbdyh34lDG +RMzphC8Y3C/HT3nUxyLKv1tlNoUYpAXNuOyRiEUPB2QS2Msw2Cbg6a018Rod +cpWijNdgoDiKWCjar1HduifbwUYMHD+zhXl8iEZ3znFzyLdgkBgoLWtKiESS +q1wnRgkYbB4oVgqNiUB8b32FgnswsPl0Xuv71XBkLwmPmfswONZz66YoKQxh +R6gcYcMYvHaadTJ5G4o4fv0SiBvBYPWfp/sLsxcoMjhG2XMcg2rJN1Fz/UGo +9vorY8I0BhLfxVJEUgIRr9x1Ka05DF5ezs++evM52kauq4rGMMBG1YtCDgUg +C1VF9pIFDD7Mxwc+Qn6oJEN8uXgJA1eBgpS+cz4oqf/A4+gVDFQ124ug3gs9 +Ufzkdf43Bqz3Jn9maj9FIT4xfzv+YPAzdtvkrh9P0O326BmFdQwSKg4xuxk8 +QvttA/978pexvnEVQXK7K+q9fnMhYgOD68If7b3f3UdMajKrfv8wmHamJuhm +OiCLTzt1DTYx8EV8rfs+2qDvRBUaxmAxsaT1mpsWqLCFs9J5CwO9Brbb5g5X +kNXCYm89g6sz2jlb69SRpbyr/AqDdXFb1Vmy4pApIlm7yeCg9/a/Sn/rQrhc +RxSNwUJ52txabtfgSxwkpjLYa7tzj+kZKxi8Ykc9xeCJb4q4MHk72L1d1zGD +Mf5VG6aHlUfvwrVidoUpxnyreVo0FiRc4NcrmgEHg+OH9UJieNzg+D2Jb6yM +9bUVyzkfi3wM20LAjczIB0vkXpPGXR4g1V3wOngNA2W7P6pWYZ4w/VRAaDcj +fy5nh0TX2bxBc7+4wP1VDLL21rLHBfmCClvi28xlDB48If7mC/CHP1FXcaWL +GOAvllLyN59DrbdT8j06BlB6tMtHPQjeWg1k7ZrFoP9IaoO+dzC45lX6BE9h +wPM34NPY2gvgiCrsZRvDINd+Of2rchggCcvnXDQMtLvsY309XsIPx7c680MY ++OQYePGvhEMut6O3Zi8G/PurnMdPR0JiNlXoGxGDQv8T1gVuUbDEkq+w2YHB +Si6HuOSfaPh22OC8eDMGsoEYR/WTWJgRl54pKMbAUm0v+2exOFjrfHYu5isG +4UuK25M64mCnnWm49icM6Hd8tx7JxMODmaZ8uQwMCrQ4l6WGEoD0ZtyEhXF+ +RzdOLOx7mQh90ab2kc8x2FdkQmdWSoJ/OzOw+WcYPJFImRyITga5/CZN1QeM +/LLJDkVppYFTvLrQmBEGTlWX+r0X0iDn1W/eo3oYJLm79TqlpcOwh7yx7jnG +flsanBT1fgcyD/JxXCcxGL4jPvq3JwMcpXkLf3Bh8F7rW++XMDxkymwrVqij +Q7nIhWuSVDyEUuNqDpfRgbhB6k5RzITwiduPJ/LowFq0TggbzYTS7G7WXUl0 +sJPQbLuj/gHyuOLvlz2ggzQbofrAQg5Irb7VXNxDhy9N9Bw/sy+wOhsVI6E5 +B5RbnD7yOV8gJkt765TCHHCtyhoNr32B3jdyw0el50Ax5syzdMN8eO9zX3iA +cw7iqh+KbtHz4ZVNV9gF8ixcOTjrXHWyAMQz5AWeO81CC4nKpl5SCIr0Qsc3 +fjNQqdeirFpbBhV8vh3XwqZg5uN7p0tmtbDEMSL4WnwcVvsdWRqca4F0Scav +de84MO2ST1YLrIW+wUZWzm3jwHfve5t0Xi2QhTInM2ljoCnbLb+NtQ5+Pg2T +PJU+Bgm521dLvtSBlTjfo1ahMdD65PhcdEcDtJ9sFH0nNAqpn+WSl4uaYOCd +ymF7FRo82M0ftvajCS6+ZKduyNJA4x6zx+ZQE9QV34TUQzSgSHcbs7M1w6Ew +z8n1bTQQwXtyCF9rhrdaM+eOdlABn1TveX61Gfw6I20tHKjw6YW5WYxiC6Qb +7vEQY/hLhXXI/hMlbVBpt3dR9/og9PEMEcfb2sDnY0ncIxiEVXQ6JnmkDZq/ +95sXSg3CiUM0rh1c7RCiFCH59M8A4PrV2Cm27ZCjwBG0I3kAXl5dWH/F+xN4 +Q9ButbF+MLtgSf3l3AF5JKJxckgfeO1L7LnMS4Crlubf+dh74d62a8MPhQkw +/HbYOG2VBJaLvFOxkgTQK1bfqzjO8MfO8I1eVQLEupv0vaklAWu4/xGrOwTg +f8bKtt+XBJHbHB+6FhLgJoidYfvTA7hFRc7X17rgxhVO+2u/u6Gts1urK4nh +yeuex3+IEKHY0ly5Bk+E2e1J6sp7iZAxOSybn8dgm/m8b9uJ4M40wxtZSwR/ +zTMhHdNdIHKCmaY3TYSGET2UWNQFT6KO+9SodENRW9qZB4ZdcOhSSFE+uRsc +0mS8m4IJ8LBVSTKKjwQ6bWtnyXyd8HxrKeCaKAkODVx4L7yrE2JO5w8dPEaC +GZH+XjKj/ihMln7z6RwJtBTOrlh0dsCfewKsLfdJkHyh3czlcQf47v5LY60j +wXYKi39R9U94aViV6uXSC5tpZDFBp3bAd2jtt28mAze/9orjcgvwJ8ckP+4i +g43xSpjsaAu8sqeKBg6SYb556fg/Qgu4/Xsm926eDPwuC3btX1rgnOw3nX7+ +PrjSwjx5+H4LUEJEvS459UGcUuua368fwK+xSVHk6Ic0oZnyhslmCP9Ukr3j +4gB4CDO/t+VtAkmRs2/p1wbAmWOexWhbE1THVgf3WA/AgO3qPr3lRljy+mGD +8xiAUx7iPLakRjDXGxAG3ACsctBr1JIbQWZiM9p9fQC+7rDVspNqhMbD2p7j +2YOghnvcra7XAMyx3Xp12ylgrGy5KZZfB5up0zYL3BRo7W3WHMDXwXo2k/dh +QQp0p1LfpSTWwXLVsbxnchTQ8SMEqwbWwcR0EO+Z6xTgLD/xm3atDto0FQbw +2RSoORZY0btVCwmLb+6HGAxDu9PPFTHGOT+Qk9uE16XCs2ovfw+RGuD9J6j3 +3pAKZ0xvKdw5UAO7r7xsSTemQsPwtlpTrhpgWXNsT7KkgvLjBt3LG9VA15Uk +RrtSIejollksuRoaJ99RvBOo4MlNZjaPqQZ36fgVk0kqDPGQa4o5qoH0wf/I +9lAabAzm2fcLVjHq6119DyJoEL9+wdKTowoWf7+JJMfQ4JWUxH6J9Urg0s1d +y2XUDRvPPdPe91WC9mRPh/E3Glxz6i44k1gJJdKy3u8oNJBn6ZL9KFgJCTk9 +3WoKI5C4pOH3QrICLHKPhriNjEByYp9u+u1yEDy82en3awSGYp95/jUuh/4Y +glDE7AhYLyfaWWmXg7nX068fVkegrS7iitbxcjDTaxwc5BiFphjK+0d/y8Dk +l/VpndOj8CPwbqJTYhlcPJJIFQoahfThb/1+A6VwNnmHWoPEGKjzkhN73Etg +NKpsVTxoHIL1pdf4Ogqh8IjfjOjgBLh6RPvdTy4AKs7jFOvBSVB0OjXeSc6H +LukL/ktPp+AF1/D7Gfs8OJx0QaS5chpCyqtLpypyIf7jjTWVjRlIeK9TYeGW +A0Sk2VdqOgcugpjPjfMf4F7NrV0HE+gQW1AYYeCRCTbpffEVqXTQtI06UuWS +CeY+JmKWODqQAzWyzjhmgq6Kvkoy452qOubgcPxGJoh/VXAQaKDDNvlnkYZq +mTD4bnfd/kU6VBQUHtFmyoSL/uVeXBcxGLLqesYWgYfztxTY8q5gUML6I2sh +GA+qZ/NfG17H4IujoDPVFw/SfzKzw20wKOWSqW12xQOr62sShycGCkPj+vUm +eKi87XhqRxbDQ8sa620P4qFQfbQyKxeDENUzoaF8ePgkdEtPm+EJO9gMa75w +4yGp96pVcAUGG57XX7Gz4MH9MuNKZHiI+Qo/jE7gwPl4uQCOzHj3q4Vdxag4 +sN2lkKlJwUDuXw6rXR8OjJuOVgQwPMipwu3ycisO9DIzdUQZHlxUocZ8rgEH +/wWKEKsZnlUolcP5ugoHStbJt6wY3iatFuI2VoKD//uvhf//r/0fJxxYHw== + + "]]}}, {{}, {{}, + {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[ + 1.6], + PointBox[{{-1., 0.038461538461538464`}, {-0.8, + 0.05882352941176469}, {-0.6, 0.1}, {-0.3999999999999999, + 0.20000000000000007`}, {-0.19999999999999996`, 0.5000000000000001}, { + 0., 1.}, {0.20000000000000018`, 0.49999999999999956`}, { + 0.40000000000000013`, 0.1999999999999999}, {0.6000000000000001, + 0.09999999999999998}, {0.8, 0.05882352941176469}, {1., + 0.038461538461538464`}}]}, {}}, {}}, {}}, + AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], + Axes->{True, True}, + AxesLabel->{None, None}, + AxesOrigin->{0, 0}, + DisplayFunction->Identity, + Frame->{{False, False}, {False, False}}, + FrameLabel->{{None, None}, {None, None}}, + FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, + GridLines->{None, None}, + GridLinesStyle->Directive[ + GrayLevel[0.5, 0.4]], + Method->{ + "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> + AbsolutePointSize[6], "ScalingFunctions" -> None}, + PlotRange->{{-1, 1}, {0., 0.9999994079562567}}, + PlotRangeClipping->True, + PlotRangePadding->{{ + Scaled[0.02], + Scaled[0.02]}, { + Scaled[0.05], + Scaled[0.05]}}, + Ticks->{Automatic, Automatic}]], "Output", + CellChangeTimes->{ + 3.6855203207892113`*^9, {3.685520384995545*^9, 3.685520396339651*^9}, + 3.685520461663477*^9, 3.6855205001903887`*^9, 3.6855297538885403`*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"rungeInt", " ", "=", " ", + RowBox[{"Integrate", "[", + RowBox[{ + RowBox[{"runge", "[", "x", "]"}], ",", " ", + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"-", "1"}], ",", " ", + RowBox[{"+", "1"}]}], "}"}]}], "]"}]}]], "Input", + CellChangeTimes->{{3.6855204254572487`*^9, 3.685520430517417*^9}, { + 3.685520516108758*^9, 3.685520560228278*^9}}], + +Cell[BoxData[ + FractionBox[ + RowBox[{"2", " ", + RowBox[{"ArcTan", "[", "5", "]"}]}], "5"]], "Output", + CellChangeTimes->{{3.685520518937872*^9, 3.6855205284217854`*^9}, + 3.685520568516642*^9, 3.685529754132867*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"rungeIntApprox11", " ", "=", " ", + RowBox[{"Integrate", "[", + RowBox[{"rungeApprox11", ",", " ", + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"-", "1"}], ",", " ", + RowBox[{"+", "1"}]}], "}"}]}], "]"}]}]], "Input", + CellChangeTimes->{{3.685520529848077*^9, 3.685520566340179*^9}}], + +Cell[BoxData[ + RowBox[{ + StyleBox[ + RowBox[{"InterpolatingPolynomial", "::", "list"}], "MessageName"], + RowBox[{ + ":", " "}], "\<\"List expected at position \[NoBreak]\\!\\(\\*RowBox[{\\\"1\ +\\\"}]\\)\[NoBreak] in \[NoBreak]\\!\\(\\*RowBox[{\\\"InterpolatingPolynomial\ +\\\", \\\"[\\\", RowBox[{\\\"rungeData\\\", \\\",\\\", \\\"x\\\"}], \ +\\\"]\\\"}]\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ +ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ +ButtonData:>\\\"paclet:ref/message/General/list\\\", ButtonNote -> \ +\\\"InterpolatingPolynomial::list\\\"]\\)\"\>"}]], "Message", "MSG", + CellChangeTimes->{3.68552975425229*^9}], + +Cell[BoxData[ + RowBox[{ + SubsuperscriptBox["\[Integral]", + RowBox[{"-", "1"}], "1"], + RowBox[{ + RowBox[{"InterpolatingPolynomial", "[", + RowBox[{"rungeData", ",", "x"}], "]"}], + RowBox[{"\[DifferentialD]", "x"}]}]}]], "Output", + CellChangeTimes->{3.685520537852358*^9, 3.685520569087985*^9, + 3.685529754312867*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{ + RowBox[{"(", + RowBox[{"rungeIntApprox11", " ", "-", " ", "rungeInt"}], ")"}], " ", "/", + " ", "rungeInt"}]], "Input", + CellChangeTimes->{{3.6855205445248613`*^9, 3.68552058114006*^9}}], + +Cell[BoxData[ + FractionBox[ + RowBox[{"5", " ", + RowBox[{"(", + RowBox[{ + RowBox[{"-", + FractionBox[ + RowBox[{"2", " ", + RowBox[{"ArcTan", "[", "5", "]"}]}], "5"]}], "+", + RowBox[{ + SubsuperscriptBox["\[Integral]", + RowBox[{"-", "1"}], "1"], + RowBox[{ + RowBox[{"InterpolatingPolynomial", "[", + RowBox[{"rungeData", ",", "x"}], "]"}], + RowBox[{"\[DifferentialD]", "x"}]}]}]}], ")"}]}], + RowBox[{"2", " ", + RowBox[{"ArcTan", "[", "5", "]"}]}]]], "Output", + CellChangeTimes->{3.6855205817013617`*^9, 3.685529754579514*^9}] +}, Open ]] +}, +WindowSize->{1054, 1179}, +WindowMargins->{{0, Automatic}, {Automatic, 19}}, +FrontEndVersion->"10.3 for Linux x86 (64-bit) (October 9, 2015)", +StyleDefinitions->"Default.nb" +] +(* End of Notebook Content *) + +(* Internal cache information *) +(*CellTagsOutline +CellTagsIndex->{} +*) +(*CellTagsIndex +CellTagsIndex->{} +*) +(*NotebookFileOutline +Notebook[{ +Cell[CellGroupData[{ +Cell[580, 22, 308, 8, 32, "Input"], +Cell[891, 32, 216, 6, 51, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[1144, 43, 490, 12, 99, "Input"], +Cell[1637, 57, 1203, 35, 55, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[2877, 97, 326, 7, 32, "Input"], +Cell[3206, 106, 863, 21, 58, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[4106, 132, 1021, 25, 121, "Input"], +Cell[5130, 159, 689, 13, 23, "Message"], +Cell[5822, 174, 687, 13, 23, "Message"], +Cell[6512, 189, 687, 13, 23, "Message"], +Cell[7202, 204, 578, 12, 23, "Message"], +Cell[7783, 218, 8890, 158, 238, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[16710, 381, 397, 10, 32, "Input"], +Cell[17110, 393, 221, 5, 49, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[17368, 403, 332, 8, 32, "Input"], +Cell[17703, 413, 645, 12, 23, "Message"], +Cell[18351, 427, 331, 9, 49, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[18719, 441, 215, 5, 32, "Input"], +Cell[18937, 448, 598, 18, 63, "Output"] +}, Open ]] +} +] +*) + +(* End of internal cache information *) diff --git a/Lectures_my/NumMet/2016/Lecture7/subdivision.nb b/Lectures_my/NumMet/2016/Lecture7/subdivision.nb new file mode 100644 index 0000000..1388964 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture7/subdivision.nb @@ -0,0 +1,537 @@ +(* Content-type: application/vnd.wolfram.mathematica *) + +(*** Wolfram Notebook File ***) +(* http://www.wolfram.com/nb *) + +(* CreatedBy='Mathematica 10.3' *) + +(*CacheID: 234*) +(* Internal cache information: +NotebookFileLineBreakTest +NotebookFileLineBreakTest +NotebookDataPosition[ 158, 7] +NotebookDataLength[ 22366, 528] +NotebookOptionsPosition[ 21255, 486] +NotebookOutlinePosition[ 21592, 501] +CellTagsIndexPosition[ 21549, 498] +WindowFrame->Normal*) + +(* Beginning of Notebook Content *) +Notebook[{ +Cell[BoxData[ + RowBox[{ + RowBox[{"f", "[", + RowBox[{"x_", ",", " ", "\[Epsilon]_"}], "]"}], " ", ":=", " ", + RowBox[{"1", " ", "-", " ", + RowBox[{"\[Epsilon]", "/", + RowBox[{"(", + RowBox[{ + RowBox[{"x", "^", "2"}], " ", "+", " ", "\[Epsilon]"}], + ")"}]}]}]}]], "Input", + CellChangeTimes->{{3.685522664678437*^9, 3.68552268302462*^9}, { + 3.685522742927903*^9, 3.6855227663596888`*^9}, {3.685522810599539*^9, + 3.6855228111512327`*^9}, {3.6855228661831017`*^9, 3.6855228728946533`*^9}, { + 3.685523930821341*^9, 3.6855239310050707`*^9}}], + +Cell[BoxData[ + RowBox[{ + RowBox[{"fDataGeneric7", " ", "=", " ", + RowBox[{"Table", "[", "\[IndentingNewLine]", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"f", "[", + RowBox[{"x", ",", " ", + RowBox[{"1", "/", "9"}]}], "]"}]}], "}"}], ",", + "\[IndentingNewLine]", + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"-", "3"}], ",", " ", + RowBox[{"+", "4"}], ",", " ", + RowBox[{"7", "/", "6"}]}], "}"}]}], "\[IndentingNewLine]", "]"}]}], + ";"}]], "Input", + CellChangeTimes->{{3.685522896958269*^9, 3.68552304398955*^9}, { + 3.685523098901132*^9, 3.685523192412067*^9}, {3.685523871750024*^9, + 3.68552387185387*^9}, {3.685523911277623*^9, 3.685523911461322*^9}}], + +Cell[BoxData[{ + RowBox[{ + RowBox[{"fDataSpecificA", " ", "=", " ", + RowBox[{"Table", "[", "\[IndentingNewLine]", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"f", "[", + RowBox[{"x", ",", " ", + RowBox[{"1", "/", "9"}]}], "]"}]}], "}"}], ",", + "\[IndentingNewLine]", + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"-", "3"}], ",", " ", "0", ",", " ", + RowBox[{"3", "/", "2"}]}], "}"}]}], "\[IndentingNewLine]", "]"}]}], + ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"fDataSpecificB", " ", "=", " ", + RowBox[{"Table", "[", "\[IndentingNewLine]", + RowBox[{ + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"f", "[", + RowBox[{"x", ",", " ", + RowBox[{"1", "/", "9"}]}], "]"}]}], "}"}], ",", + "\[IndentingNewLine]", + RowBox[{"{", + RowBox[{"x", ",", " ", "0", ",", " ", + RowBox[{"+", "4"}], ",", " ", "2"}], "}"}]}], "\[IndentingNewLine]", + "]"}]}], ";"}]}], "Input", + CellChangeTimes->{{3.685523648623685*^9, 3.685523662839806*^9}, { + 3.685523875534165*^9, 3.685523877101618*^9}, {3.685523914581822*^9, + 3.685523916509296*^9}}], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"Show", "[", "\[IndentingNewLine]", + RowBox[{ + RowBox[{"Plot", "[", "\[IndentingNewLine]", + RowBox[{ + RowBox[{"f", "[", + RowBox[{"x", ",", " ", + RowBox[{"1", "/", "9"}]}], "]"}], ",", "\[IndentingNewLine]", + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"-", "5"}], ",", " ", + RowBox[{"+", "4"}]}], "}"}], ",", "\[IndentingNewLine]", + RowBox[{"PlotRange", " ", "\[Rule]", " ", "All"}]}], + "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", + RowBox[{"ListPlot", "[", + RowBox[{"fDataGeneric7", ",", " ", + RowBox[{"PlotStyle", " ", "\[Rule]", " ", "Black"}]}], "]"}], ",", + "\[IndentingNewLine]", + RowBox[{"ListPlot", "[", + RowBox[{ + RowBox[{"Join", "[", + RowBox[{"fDataSpecificA", ",", " ", "fDataSpecificB"}], "]"}], ",", " ", + RowBox[{"PlotStyle", " ", "\[Rule]", " ", "Red"}]}], "]"}]}], + "\[IndentingNewLine]", "]"}]], "Input", + CellChangeTimes->{{3.6855226872609577`*^9, 3.685522730751862*^9}, { + 3.685522771991954*^9, 3.685522799287492*^9}, {3.685522861703998*^9, + 3.685522862270746*^9}, {3.685523049709237*^9, 3.685523089260763*^9}, + 3.685523196308086*^9, {3.685523674303648*^9, 3.6855236969673*^9}, { + 3.685523880558132*^9, 3.685523880637704*^9}, {3.685523919717774*^9, + 3.6855239198052998`*^9}}], + +Cell[BoxData[ + GraphicsBox[{{{}, {}, + {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[ + 1.], LineBox[CompressedData[" +1:eJwV13k4Fd8bAPBrp3ATCamskYQoCV/vJKVSKm2ShIRkJxRKkYpIpSgUEZKK +rFGdsVfWLFH2/eLeGfu+/Ob31zyfZ57zzrzvOe+ZM9LWLiaX2Wk02gAbjfb/ +a0hg+sTKijCeg/2OwFtJqNv23Fx5SRjXqHI4n0t5XeOtMtN5YVzMNEc6nXK8 +3MmozElhnPn9ZEIU5dyyOe1LQ8K4aIibozPlHh7DO+X1wniUD5/aBsp7Hnav +CXsrjOddFDrv9o8ERqSoioSRMJ4YkfmB/peEqIZUkaXgtfj+ax3xv5pIMHFa +XuyoFcLvizcZmtWTELykdc1UUgg/YmyWEV5DwqrCmBcvz63B9VPlPbp/krA7 +JfBWniEdbzxw+XJfKQkL6lJyld8E8Bsm30bTvpOQ8T6YIaDDj5/1PWwUlkfC +rX8eGcb1q/C9u2ea6j6RQF7KrJMy4sP5rMRrx96S4JN4bXK0nQfn9TaLXH5J +glJb4JhpCw9uZbWM11JuE3lC4PU8+J/s7dXxlLHgj4yICh68qzs3UJ8y75XB +NrVMHvxfo7xt0AsqP5VzZa5BPPgcZ1Q/WzRVry+6z0eVeHBvwQujPZEkTNRy +aI16c+Mu0XnFu8JJWHPY40aHGzfOG/Fs9VwYCcqlPV+rrnLj33fzyX2lbJNf +tPfdRW68XL2veC/lP/G3jKwNufE59a+fDz8kocBtwaJRlBv39n6vfDKEhADh +ibv5OVx4hpIn73/BJAie7WkImODEFeqmztj5k8DBvCjYzOLEy8cii6Qpzwa0 +H9zO4MTDfjVNtvqR0J3W8u1vGyf+wsEk+xjl7OWaFPVyTlx+/nvETl8SzicX ++vZFc+KSBq9DJnxISJ2KlD2ox4m/SvCvNPAkYV+koQf9AQdeNTaULO9AwhuH +LRMjdzjwYpdenp9XSGDfy+Xxw48DV2M4CjtSxlnF7rddOXBlzjDFTHsS9Az1 +3CfOcuBid67CbjsStOZ3ujVv4cC1KhpxdRsSVCxlXF6XsuNtW+tmey+QEK7J +Rvp+Z8cPDD5q9KXM4u9yNs1nx+tkAo8IU07/EucslM6Oi8rw0veak7BVWNw5 +8Ck7vjvYry7KjAS5crqTnTU7zh898nvHWRLElZccVGnsuN2SC5/IcWo+lM6n +KM6z4fu1csLij1Hzr/ilV3qSDa+Tn4pXpmwrd81chMGGj0+Ku+8zJuGzJOvo +XC0bHmqT8/HqERIOC7Srlbxmw92sOPmTD5LgR3ydPgNseGev8tw3jIRHTAmN +43vY8NXtE9t1KL8Z9nE5pMGGey5fa88DEn4M7GToKLDhXqat7z7pUf3Zmf5X +is6Gj1nu6ovWJeFjbezX4Q4arlVzW+eEFgldGX53bt6k4af9TomcUiWh6icj +t9OHhj/SHNvzU4WEvJ5TI5gHDe/OWWn6j3K4yPZTHPY0vP/La2f57STo+rTL +3z9Bw//CiAChREI06P18LEfD/ZUrLW22UPtB7Mrmw89W0H/+vjOdkiR41zTf +5w9ZQWbbIrz0KRvTMsZqbq6g7ZsUTyZtIGHxkmXJSfsVtH7+irStBAmm24ps +L+isoEiJ2Lnu9dT6L7jzwbV7GeW0ztZ8WUut32Yunajty0hdYrtjLw9V/4fx +KrbSy2hO+uQzKcotmI7MrnXLyO8sd6s5Nwlqaa58DYtLyKbleW49J5W/X2uz +YNUS2nnG4l4WGzVfMpmewQ5LKIoZYXxwgQCak/nHaymLqNhVuSmLRUC69EyC +Qcwi4rS++Pk3kwDTP4+fCT9aRA8Ofb1BjBCQARV+mV6LaI5zskxumADLtepG +rP2LSNHyuUjQAAF4Hg/Dpm8BvZWOmV3dScBt9izpU9ILCH26t1ejloAO5hWu +iXULKOVl2G3dGgK0m6WGHq9aQB8dLrwwqCZg8n14Ru3EPEpz4zp0opIA+9MO +YFQ+jy6Ym3qaVxBgnCp9Qf/qPIpTZxzVQgRsMI6IVsuZQ+mtwU0zHwjw0TL0 +q303h4rktfe0pBPQKLN80fnVHLol8kQ97z0B4TNXFdLvzSEB3Z9tru8IYE8w +zFUwm0MzT7yr6pMIYEwsN2xamUUnaU/HjGMIyHnhRBc4NIsaN6pziN4jYIB7 +0D9fdxZpztu0pdwlQMzTknlJbRb1TuSk7A4iwM/45M+C9bOo444S14nbBBhw +7Am0H5hBO969d3TwJaDekXOmJHAGmai+O7fThYBRiOm48X0avaocnOU6TYDM +B5GjWz5PI+vra1qOniTglER44e+306hOYEL86QkC8icDohXDppHoda7MdcYE +3Hpna/Ln/DT65usYNX+AAEFh9fIdc1Mo02LSTU+LgG0DPz4wNKaQf0ByY5AY +AaVvsGfVClMoLgR98BYlwOJivt/nDVPo9oVnB+xECIhoSTHy45hCJ7bcjNVd +Q8D0r+AhesMkyubIFEM8BBR/MpDf7T6JvqV5C2+cZoHpdTz2bsYEShBgbuWp +Y8H4Lq0gh6QJVKKuGRpYzYLQ8U9Xj0VPIA/hzk8Lv1jw/eprbfGACSTBom/q +LGOBrMXNlvRjE8g0+0egWyELCH1dkUZiHAW/5a8pfMuCIP78UJnt4yhktZZd +jBcLapaNLjVLj6PdH8xmMj1YIDbWqf1QdBwNaVxNLnFlwYcmnuGplTE0pz6i +2eLAgj+vzh78WT+G5O+eHy+1YIHCjhlOV58x9KVg58Kb/Sz4dVrz5reSURRy +RE5Jks4CkYO/Trvnj6ImL6zu1moWWGhbbFf4MIr4pzDudh4WTGwKbot4Poo6 +/30fCKSxQJLRpG1zZRRdj7Gtvz3GBKcb12ZW0UdRDe1EuFo9EzbrXenOyCbQ +tLXGbHcYEywURZM1YglkdNjk0tR9JsSuLXHIDSTQppfHUpYDmSA2KDlZaEIg +RdXfxuPXmbA2oo67YpSFZMIXN5y/zATunt3K7cosNL1ryyNxbSYYVPWNXhRh +IY3GEkP/nUwIzH2c07PARP1KpZxNKkyghY7oMX4x0fDxCj5LWSbMabw+MWHP +RBMiSy+6VjOBeZfbm+/tCPIsLNFz+jcCSq5ZOg8fjiATiYqNbI0jcMXMkkb3 +HEE+v/Dy4OoRGFApeCCybwQdO1/fbIuPQFezU+zm7mEUIPztzrW3I9Cwtalo +18Zh5H664E2K4wgUVCfyWz9joH39nZtFiWEwk32QGXSTgVY5BT3cODAM8z7O +Z1LsGCjSsfmEcMcwaMvtiWfuYSDda+vlf1YPQ8GNGg3vjkH0TXNfoGk6ZYU5 +s4cKg0g/Vr6q5TIVz79j5eOaQXQ6O+jfB3MqXn1J0u+5AXSkUnjK6SQV72Y4 +IVo1gLYWTth8xKjxjbJ3ElwH0GrbucStEpRvG7/L/dKPHGoY3fU/hsCsWcP4 +75t+lMou/lvn+xDMK4tPLIT2o+HJ5rrQrCHQbunV3WvRjzRJDUZn3BAUqFyv +q2TvR1ND1U+Pu1FuTZzpMupDyXfeca4Vprxrbv/qjh4Uck6e0QwM2KrU2hZf +0YO2Kmocrd/BgOhN3zw0M3uQsF5DSposA7x4bidYBfUg2Zb5sgluBmj85V3K +U+pB+Y7rlNsrByHdTyzbxrsbBTW/rc8xHoQNbvOH5y270Q1tg0QrvUEIudzW +/ehwNxJ781Wqf/sg2BvH0ws3dqOq3gpXN/5BkJNSvCpU2oWUSbfzTysGIK5k +t8z3NV3I59iA5YtdAyDmquGeMtmJStQMA0JlBmDadOZ1fUUnWn5WvdWYPgBZ +SrcWFJ07kX1S6Z59jH5Qrgn/3FTQgaJmWjk6o/ph47qPUqqn2lHXXvY7bqw+ +oCUyF7uD/yG/E2x9Mc69oB5VcBrO/ENq11qGb5j0gk3o/Y+x8v/Q86InAms0 +e+Gnh5zl2dK/qNF87RWzxR54YmBeUsX2F5XKFptY3u0BuYHKB3l+zchVc2xD +8f1uKtSLHpEjzWi1UmdbmV033K+x03Hf0IyY9tH7fQ50AyuXg6Vc8Ac58in+ +GWTvhtx7OsffzDShP4YKA8PeXXBo63vRMPdGtHbCW6ksqQP8Nl53Gd7biLRn +0919TDvgo5DhD0OhRhQ9w1cdvLoD1s73XGfPaEDK6/bT81zaoe2XRLs3qx4l +68cb+Km1gbNjSKK1/W/kVxP5OjvyL8hdbLB0/VqFJss6QnJnG4A+dfxpmlUV +ChqjW/5IaID5kJqyPu4qFKHLxxVzuAHqcn8qnTtRifBVxRuiYurBTxCf3Dvw +E5k2Ha0P1PkNDd/S7wuvqUAMYfV72FgVfD+pVHg0pxxFVtA2iOytgtShFNY9 +s3IkFXFURPtxJdxcl2iymFSGIgvlryD1X6Dk9EKyf08pKnVevBfpVwF3JIMz +ci8VoR8zB3SQTjE4feboHeUtQl9Prfqm+6IITA/eXrftI46yBc6r683iYM6Y +ndQbQ8i++f6fA14IsqP676/Z8A3Fp3jWFH4ogJCSLlbQ5a+IMfXhOH73C1iS +rSZznwrRrLXLZhXrfOA3rJfsMShAUjf32c0p5sKl6e8ZWc55KGB4X+SuiUzQ +kilYp/glF2HEYJzDRAYIGufciOXIRaNnHTT/jX+CgrfvDe5GZyOakXWtgGo6 +RPxOfjfXk4Wk5u4YDlxMA9ulBEHn7VkIe/pN9ENcKqw9Hd1ypjgTjQoeXGO0 +9y0wAp7+V8mfiTLWV07gKYnwPT38DZzNQBGfks2tJd9AI++qv8dbPiBsuWU0 +wfoVBFjSEvc4paP4qF+rzifFgHL+tKMM+3tE8/8uksr5AoLsemmT21IRXp8c +wyMbCXXiSrLOb5IRbX34H5Wsx7CxynU/Q+wtCjjhHSPbEA4ON/PsrB8lIqz2 +aotXWAjkqi0/aON6gwKGjjwyfx4M7L0G6Wf84hHtVWfKD+IOvFJaH5fQH4cC +5BIdtVVvgdajpO5n9S9RwG7VoZOy1+H3xI4tISgKBTwPDz1b6AEOpsjhZnok +omU00XeMOQHHtyOf3F88RrQzl2zu5l8GrX3Zg9X8YYgWsH2DXdwFuMpf65o7 +EIxoQQeVex1MIK5paO4VHoBoi8s7Hubvg9f/8XHTd11DNOy70504VeCJ91fD +f1tR48OeHDjJDTjTPN1LwghotCM3jRr+6tmJpgcvWLgADZMUEzgiDbU+Upnp +C75Ak3oscS9IB3Z5XKtlyQYBjXlso+CtI8CRZnaA63AI0M5Inm/TNwX70dcW +gloRQEt9EmPiaQ21oWslJIyeAk0mObNT1AE0Fe42yVs8B9ri3h0vAt0grngm +YofbCwgoPe4cmuAFnBYOR/4LigVaUM0Yp40fOMy18RyKeg20TOHnNTa3wf5X +ipXS6QTAQl0LrRlBIClmBRkdbyBglxF7f+F9qLsssVHTPglo97Pcr3Q8hMCs +hvnCsbeAn/pkeD0lAoaND+RWcL4DS3qXVV/IM4iLXXl6NDwNLHdsPdO3MxqO +D+e7NaxPB9pZdYmXHDGQd3fb9k6lT4BxhdgYBMdDzhX3u+pmGRChZjyXIfMG +Hq055jm1KROkxF1ExOsSQd+C18Q39TNESDQlnTdPAUmufkzPKQssP2x96aL9 +DqbfF6myqWcDrdvdd1zhPbybuyFwrzAHLAuc1eJVP8GaZ6wfT2ryQS2RNtAx +ngXDOr/yTj/9Aox48sk38Rwo7UlOFjMtgIziiIDYA7ngrWYZ9Kq7EBjKngte +mfnQUVkPaZPfIaPcK/pC2Tc4w2beahlWDKP5+g/2mZVC8N/mkGlWMXDWFl2X +EiiDvMyT2g+NS+B+5fpKzaIyELM2epFHL4XZ10trMpQr4F+x9lnBJ2UQ7aFw +4rjgL7AIEm8ofP4D+jj2Lr8dqgU7nuZfovE1UCUoODVf3ATRnSY30mm1UFPr +8/TCxj/wI69mq75VLURcSG455f0H1u28VKopVAfuvaaeo1ubwavPTNa6vg6+ +Pvyz1e1BC2jtP9ydf6oeGDwuE10qrfCVe+sF23NN4PnU05efowuEJL083YKb +QM36x+c0uS6w3VES6pfVBEMqUy/ZDnTBGvMLBU8E/kB2lt3nv/e64HLm4/Wo +6A8IDvC37eLtBkHz+d/rlVrgCPNubsVCN1xyMxySOdsCAZmt5znEeyA/OJKm +EtQCm4TTeRd29YB1poqqQUcLSD5v79vs0gO53Jceuj75C+nRb+FwRw9czKw8 +8HP+HziyhDbdzOiF7HIxi0aFVmjZzC2b9rMXVrVdvtZ5qhUyy/1Z0T3UfW5a +0tTHVjjBLVaAC/cBr/lONplLbSAmuSrzgEcfZHDHFd6obAcH0fu/Xij0A5vK +bjFR4y4wGBG7Xms/AOxZpuF3bLvg06ilWYMfZa0bnOTNLsheKnyQEEFZ/9to +xacuiJXx2/Y6bwA4z+z7eV2oGw50kumqHIPAe/PYjfbGbqCXdHjLPqbM4TZ6 +iNUNnuOXjq95Mwh895/Y5nL1gN+JI6frPg/CqidNJo80e+Clx/xoZcMg8Cef +37Y3ugfmpkYCkDAD1tTYtSWd7wWfB6EnOR9SPvnARMizF3RoseUWLxkg1JL2 +w/9hL/jWL3BEpjJgbQ8r+8y3XjASeBodUcoAkWmPMN7NfZBZvLDSPs8AsU0B +eo49fTBk23FLzmoIpJyj4zWuDEAaf+zcrOgwoHu3NWkBA3BOuSmhUGoYLBIc +qqqeU+ea5tvzdkrDENuoO3u5dAAWLz76fPO/YViv030sevMg4OsOVqZZD4Mg +r9LKYtMg5NS116i+G4YP0sLPfjIHQZ1D7ndq5jAc0VlUes7BgHJ+x1X0gmEI +ca45o0qd88Rnd7MSfw0Dd5P7Rysqz9xRwe7J4WFYSii8UK4/BP+UyvgvKI5A +bGHSxJNzQ5Deo7sJVx0B7aawBxddh2C0o/aE8O4R8OG1zJmlzp1hL5wrQvaP +wKQzp8C2uSHQX73LodhqBJg6R78++kTFDfAnvaJG4F9Th8Q5SSb81ZVYPzY7 +AmnHK79/VmdCe/PBr0HLI3CjMs+a/xATMsQMLAU4mSBeFJGGrjGpPgwXWhZk +gmm6vvaWWiYMpo2zW8oxoflO8rmJ2yx4qqnnrXuUCY2qztFhgwQcPyuyYTGK +CUlpZrqDS9R/Y/Xp1RdjmeApb9iFiZDAloG3fYlngvAGKcVJjIQMr9MuZ98x +wYS7Ps/sJQk4XnthZwET7ImXI0+3jEIWRy3/3X/Uc0uXV1n9HYUiAcnrU+tY +cOOWusnt/lEoUe0dLBJnAb+27YuEsVFIVytQvreRBeqfqhV6+MZg6Pt7/0U5 +Ftx8EWtwSXsMXnGcXxepwQJRZ+2bl2PG4KbF4QN3j7HAYL3X2JUL4+AjbcLR +cocFzb/faYVcGQdbe70j7XdZ4PCw/VbatXEI3FwM7fdZEMG2X3AkbBzmg/jO +VoezoHVEeJvj93GwiMp94vKSBe54po3z5gmY7XuwuzyDBfFXWM3u3RNQUsU1 +y/WPBXcITtUS5gQ0KcjKH2ljgY2HZLDw7AT0JaVVh3ewQNHfaGe24CQw30h/ +4OtlQWbEu4gpnUkgNCTt/oywoDjv8iGf55Ow8Ko99vciC3q5Ogr8jaaA73p4 +1A1JAspDpoRqzkyBk1JimvkmAlLpAlc2WU+Bt4fQko4UAY7iuuuRzxTYb+72 +HJclYFL5pQcteQpGf88CbCOA89QZ5UDaNHj+OBFjvoeALW9q4u7nTMNmcYv+ +8pMEaCTJOlQVTYOt/MMF39MEYMk+mmtqpsHjHccdlbMEmKXJ1EYNTIP1aP1E +qBkB4Z+9aCnrZ0An23NS0YqAmeJNl8quz0DsUe0I0omAH73Oiux6syBZba9e +c5eApv6Syf2HZ+E2UlXXvkdAz6BY0YMz1LvVNGxPuk/A4kjxOSGXWdDLSVBw +DyVAbVI0VCp+Fhwt2R9ORxDwghOx9NjnQIHzoXpSDAEO8vQs3/I5eFIjr5P2 +iarPmOEj8vcc/Gl4LfsjgwCZ7wFXL7XPwf6tRm/7MgloPTMuazQ5B189r/qL +ZxNw5H5TpITMPMBDQ9w9nwCVkRivL/7zYF9h7t+PEzCeoag9s3MB7Dl7Lvyr +I8DY30r0KrYAYRyKEbW/CUg79HK8w2gBwnnV4krqCbjUs/p9ufUC7HmWY5za +SECjMLnhecQCfEyW3mvTQkCuV87STuYCxFU19qV0EuCrq1/knrgI23zTQ6+w +CHD94if34eMiJN6VUN1DEGCjmXdv8MsiMALpwzwk9Xw1ZWPzukW4l9j//s0o +AVJyoq37lxYhNmhTUPUElf/q4Umxs0vgmTP2dm6OgLWtjxXRqmUouBunkc1F +Ao9pVejcumV4JsitcYWbhMVGblJDehm6WI6uG3lI6K/xzUndvQxZnr7vAnlJ +yC++rP/EZhlm3CTVDVeTYJG2x/zy92VwtAoJfL+GhLTr3Y9Wu6+AfnTnEilB +wi+Pc/Z7bq7AyPsx0dgNJAw5/cbsQlZA5dwOV0NJEhSti8aK36xAIs/rwZiN +JKQcTjh5o2EFnK1CeXWlSKgPywj1lKBhE7vWrJyTJ6HOwvVYry8NS7X4WSSs +SkLX+XSd+ds07O7uXR3vKI+aMhSE7tMwXPa6GaiR1PfFkgaRNGw0TcDLfgcJ +xwyPZ75Mp2GvnLNOZGmQUKOqJmzSTsPYRy+ZKmuRUEkb/YPrsWHTOdUGmdS+ +xd67sdLFgA0TEx9UU95LgnaZEdp0mA3D9selJFNOfZCS4neaDeO+qLr7pT4J +wWsv+mg5smFTyn2ErwEJ+nI14hkv2DCRb37/KR0kocAw3fz1JBumsfRrreYx +any4fY9/GjvmtN5kQt6chGuqux4nZ7BjVUGQ94iyfh0bVpvLjqWTz9fNUW4T +iomTKmHHXkncS/x1gcrvWbVpSSs7FqliuvbKRRK8Y9RreQU4MB72J6lPrUk4 +kLpY+NSFAzu4+fLMG3sS+ooinr3byYm9TtvGKPIgYUm7416ONifmhU/tE/Ek +QTR7240ijBPbWpZ11pbyweRyi79HODEuwjSV7xoJH0IWFfguc2IcY4dnjniR +4HXSvsD+GSemGWv/pNyHBL5+6FSY5sSeul9tvedPggovqZicx4W9t448ZRJM +goutiuzhb1yYdKjS7QjKGaVOG4liLszZ2fpmDeUdt5lCmjVc2Cdvr5FD90jY +OT80X9bHhTF4zJ/p3afqz+yr6hfixiQS8eObQ6j86v65yjtyY3rZXq9RGAkX +o8vzk6R4sMEF0a9BkSRo5iQ5FW/hwYbU1uJfKfPX35HpUubBLo7eip6knL8a +Hm7Yw4Oxu9tWXHpGgtCtvItPTvBgnNnaMthzEkoup3IH3OHBYlKkx0ejSFDQ +CDl1vo8Ha9vZJrUzhgSy+uioUCovZtpgY3AogQRm7cYPFR94sfJdZY89KQ// +Zl3xz+LFvNOiil5THmgK62V858UKvi79mKLc0V79BzXxYlFPs6/HvyGhinX0 +mxMHH+aQH3t9NJHqJwHjkJ8WfJhAx2yibzIJtkeN5W6vW4UNfFyu3feehJHO +IK0i39XYGqn0R9FZVH+svzrpQhPAvKwDtpsUkhDz4le61FVBrOagbszdIhKO +N2Rv9XQTxKok4cJrytwCr5MrvAUxvMWK/oWye4BHvHOQIHZxSduASfmwnWTk +1zhBbL4iaexEMQnzO519TesEsYRt7rYSJSSY1gkdjthFxwa8dQTjS0kQXLVY +3qtDx271ezDzKJfuG9i3W5+OuZSd/FpLWS2v4L8OYzr2M2X3jhXKPK8u7VC2 +p2Ptq2pcL5SRkHs1R+zHSzqGabOdWF9OguPb1882JNCxw3YzxiqUZTofrHVJ +oWNxWwxhP+VwEwt+0Ww6ts9u04w75ct7eGg21XRMaPRPaTVlCY8xv7wGOlY0 +V3qmj3Jdeuvcqn907EyMVes8ZV2pjInMATqmYlP3VaGChPFzL125WHRs+nie +uB7l1KdBLNMJOvbD0uvKKcoW1c4O6XN0TLp9fboDZRGec4MrK3TsBldydwDl +/wF7ZZ/N + "]]}}, {{}, {{}, + {GrayLevel[0], PointSize[0.012833333333333334`], AbsoluteThickness[1.6], + PointBox[{{-3., 0.9878048780487805}, {-1.8333333333333333`, + 0.968}, {-0.6666666666666666, 0.8}, {0.5, 0.6923076923076923}, { + 1.6666666666666667`, 0.9615384615384616}, {2.8333333333333335`, + 0.9863481228668942}, {4., 0.993103448275862}}]}, {}}, {}}, {{}, {{}, + {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[ + 1.6], PointBox[{{-3., 0.9878048780487805}, {-1.5, 0.9529411764705882}, { + 0., 0.}, {0., 0.}, {2., 0.972972972972973}, {4., + 0.993103448275862}}]}, {}}, {}}}, + AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], + Axes->{True, True}, + AxesLabel->{None, None}, + AxesOrigin->{0, 0}, + DisplayFunction->Identity, + Frame->{{False, False}, {False, False}}, + FrameLabel->{{None, None}, {None, None}}, + FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, + GridLines->{None, None}, + GridLinesStyle->Directive[ + GrayLevel[0.5, 0.4]], + Method->{ + "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> + AbsolutePointSize[6], "ScalingFunctions" -> None}, + PlotRange->{All, All}, + PlotRangeClipping->True, + PlotRangePadding->{{ + Scaled[0.02], + Scaled[0.02]}, { + Scaled[0.05], + Scaled[0.05]}}, + Ticks->{Automatic, Automatic}]], "Output", + CellChangeTimes->{{3.685522731198388*^9, 3.6855228125078278`*^9}, { + 3.685522862545142*^9, 3.685522874429077*^9}, {3.685523073648672*^9, + 3.685523196670637*^9}, 3.685523697644681*^9, 3.685523881702211*^9, { + 3.685523920312171*^9, 3.6855239436288424`*^9}}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[{ + RowBox[{"fInt", " ", "=", " ", + RowBox[{"Integrate", "[", + RowBox[{ + RowBox[{"f", "[", + RowBox[{"x", ",", " ", + RowBox[{"1", "/", "4"}]}], "]"}], ",", " ", + RowBox[{"{", + RowBox[{"x", ",", " ", + RowBox[{"-", "3"}], ",", " ", + RowBox[{"+", "4"}]}], "}"}]}], "]"}]}], "\[IndentingNewLine]", + RowBox[{"N", "[", "%", "]"}]}], "Input", + CellChangeTimes->{{3.685523491048286*^9, 3.6855235220325212`*^9}, { + 3.685523889437829*^9, 3.68552388953345*^9}}], + +Cell[BoxData[ + RowBox[{"7", "-", + FractionBox[ + RowBox[{"ArcTan", "[", "6", "]"}], "2"], "-", + FractionBox[ + RowBox[{"ArcTan", "[", "8", "]"}], "2"]}]], "Output", + CellChangeTimes->{{3.685523511725808*^9, 3.685523522484685*^9}, + 3.68552389016968*^9, 3.685523946035872*^9}], + +Cell[BoxData["5.573955509185797`"], "Output", + CellChangeTimes->{{3.685523511725808*^9, 3.685523522484685*^9}, + 3.68552389016968*^9, 3.68552394603734*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[{ + RowBox[{ + RowBox[{"\[Omega]3", " ", "=", " ", + RowBox[{ + RowBox[{ + RowBox[{"{", + RowBox[{ + "1", ",", " ", "4", ",", " ", "2", ",", " ", "4", ",", " ", "2", ",", + " ", "4", ",", " ", "1"}], "}"}], " ", "/", " ", "6"}], " ", "/", " ", + "3"}]}], ";"}], "\n", + RowBox[{"fIntNC3", " ", "=", " ", + RowBox[{ + RowBox[{ + RowBox[{"Sum", "[", + RowBox[{ + RowBox[{ + RowBox[{"\[Omega]3", "[", + RowBox[{"[", "k", "]"}], "]"}], " ", "*", " ", + RowBox[{"fDataGeneric7", "[", + RowBox[{"[", + RowBox[{"k", ",", " ", "2"}], "]"}], "]"}]}], ",", " ", + RowBox[{"{", + RowBox[{"k", ",", " ", "1", ",", " ", "7"}], "}"}]}], "]"}], " ", + "*", " ", + RowBox[{"(", + RowBox[{"4", " ", "-", " ", + RowBox[{"(", + RowBox[{"-", "3"}], ")"}]}], ")"}]}], " ", "//", " ", + "N"}]}]}], "Input", + CellChangeTimes->{{3.685523370174515*^9, 3.685523405937571*^9}, { + 3.685523787278356*^9, 3.68552381490996*^9}, 3.6855238975432453`*^9}], + +Cell[BoxData["6.2574588650388545`"], "Output", + CellChangeTimes->{{3.685523481304573*^9, 3.685523485045783*^9}, + 3.685523531318762*^9, {3.685523804103651*^9, 3.685523819289815*^9}, + 3.6855238991612787`*^9, 3.685523947746043*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{ + RowBox[{"(", + RowBox[{"fIntNC3", " ", "-", " ", "fInt"}], ")"}], " ", "/", " ", + "fInt"}]], "Input", + CellChangeTimes->{{3.685523996844348*^9, 3.685523998276311*^9}}], + +Cell[BoxData["0.12262447282305239`"], "Output", + CellChangeTimes->{3.685523998952448*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[{ + RowBox[{ + RowBox[{"\[Omega]AB", " ", "=", " ", + RowBox[{ + RowBox[{"{", + RowBox[{"1", ",", " ", "4", ",", " ", "1"}], "}"}], " ", "/", " ", + "6"}]}], ";"}], "\[IndentingNewLine]", + RowBox[{"fIntA", " ", "=", " ", + RowBox[{ + RowBox[{ + RowBox[{"Sum", "[", + RowBox[{ + RowBox[{ + RowBox[{"\[Omega]AB", "[", + RowBox[{"[", "k", "]"}], "]"}], " ", + RowBox[{"fDataSpecificA", "[", + RowBox[{"[", + RowBox[{"k", ",", " ", "2"}], "]"}], "]"}]}], ",", " ", + RowBox[{"{", + RowBox[{"k", ",", " ", "1", ",", " ", "3"}], "}"}]}], "]"}], " ", "*", + " ", "3"}], " ", "//", " ", "N"}]}], "\[IndentingNewLine]", + RowBox[{"fIntB", " ", "=", " ", + RowBox[{ + RowBox[{ + RowBox[{"Sum", "[", + RowBox[{ + RowBox[{ + RowBox[{"\[Omega]AB", "[", + RowBox[{"[", "k", "]"}], "]"}], " ", + RowBox[{"fDataSpecificB", "[", + RowBox[{"[", + RowBox[{"k", ",", " ", "2"}], "]"}], "]"}]}], ",", " ", + RowBox[{"{", + RowBox[{"k", ",", " ", "1", ",", " ", "3"}], "}"}]}], "]"}], " ", "*", + " ", "4"}], " ", "//", " ", "N"}]}], "\[IndentingNewLine]", + RowBox[{"fIntAB", " ", "=", " ", + RowBox[{"fIntA", " ", "+", " ", "fIntB"}]}]}], "Input", + CellChangeTimes->{{3.685523706930175*^9, 3.685523783583617*^9}, { + 3.685523822798047*^9, 3.685523856469739*^9}}], + +Cell[BoxData["2.399784791965567`"], "Output", + CellChangeTimes->{ + 3.685523772217701*^9, {3.685523833015215*^9, 3.6855238572254143`*^9}, + 3.685523902925892*^9, 3.6855239501498938`*^9}], + +Cell[BoxData["3.256663560111836`"], "Output", + CellChangeTimes->{ + 3.685523772217701*^9, {3.685523833015215*^9, 3.6855238572254143`*^9}, + 3.685523902925892*^9, 3.685523950151155*^9}], + +Cell[BoxData["5.656448352077403`"], "Output", + CellChangeTimes->{ + 3.685523772217701*^9, {3.685523833015215*^9, 3.6855238572254143`*^9}, + 3.685523902925892*^9, 3.685523950152231*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{ + RowBox[{"(", + RowBox[{"fIntAB", " ", "-", " ", "fInt"}], ")"}], " ", "/", " ", + "fInt"}]], "Input", + CellChangeTimes->{{3.68552397854193*^9, 3.685523989548272*^9}}], + +Cell[BoxData["0.014799695253336416`"], "Output", + CellChangeTimes->{3.685523989813211*^9}] +}, Open ]] +}, +WindowSize->{1054, 1179}, +WindowMargins->{{0, Automatic}, {Automatic, 19}}, +FrontEndVersion->"10.3 for Linux x86 (64-bit) (October 9, 2015)", +StyleDefinitions->"Default.nb" +] +(* End of Notebook Content *) + +(* Internal cache information *) +(*CellTagsOutline +CellTagsIndex->{} +*) +(*CellTagsIndex +CellTagsIndex->{} +*) +(*NotebookFileOutline +Notebook[{ +Cell[558, 20, 565, 13, 32, "Input"], +Cell[1126, 35, 750, 19, 99, "Input"], +Cell[1879, 56, 1181, 32, 187, "Input"], +Cell[CellGroupData[{ +Cell[3085, 92, 1343, 29, 209, "Input"], +Cell[4431, 123, 11863, 205, 241, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[16331, 333, 503, 13, 55, "Input"], +Cell[16837, 348, 283, 7, 49, "Output"], +Cell[17123, 357, 157, 2, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[17317, 364, 1041, 30, 55, "Input"], +Cell[18361, 396, 234, 3, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[18632, 404, 197, 5, 32, "Input"], +Cell[18832, 411, 89, 1, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[18958, 417, 1387, 38, 99, "Input"], +Cell[20348, 457, 188, 3, 32, "Output"], +Cell[20539, 462, 186, 3, 32, "Output"], +Cell[20728, 467, 186, 3, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[20951, 475, 195, 5, 32, "Input"], +Cell[21149, 482, 90, 1, 32, "Output"] +}, Open ]] +} +] +*) + +(* End of internal cache information *)