(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 3.0, MathReader 3.0, or any compatible application. The data for the notebook starts with the line of stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. ***********************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 93621, 2486]*) (*NotebookOutlinePosition[ 94331, 2511]*) (* CellTagsIndexPosition[ 94287, 2507]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["Numerical Methods", "Title", TextAlignment->Center], Cell["\<\ Dennis Silverman Mathematical Physics 212 A Department of Physics and Astronomy U. C. Irvine \ \>", "Subtitle", TextAlignment->Center], Cell[BoxData[ \(Some\ of\ the\ examples\ are\ taken\ from\ "\", \nVersion\ 2.2, \ Wolfram\ Research\)], "Input"], Cell[CellGroupData[{ Cell["\<\ Interpolation \ \>", "Section"], Cell[BoxData[{ \(There\ are\ two\ ways\ to\ input\ \(data : \n\t\t \((1)\)\ \ If\ equally\ spaced\ xi\ integers\ starting\ at\ 1\), \ only\ need\ to\ give\ data\ points\ fi.\n\t\t\((2)\)\ \ If\ at\ various\ xi, \ have\ to\ give\ set\ of\ pairs\ - \n \t{{x1, f1}, \ {x2, f2}, \ ...}\), \(Ex.\ Least\ squares\ fit\ in\ x\ using\ given\ functions\ 1, x, \ x^2\n \)}], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(data\ = \ Table[n^2, \ {n, \ 1, \ 7}]\)], "Input"], Cell[BoxData[ \({1, 4, 9, 16, 25, 36, 49}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Fit[data, \ {1, x, x^2}, x]\)], "Input"], Cell[BoxData[ \(3.55271367880050092`*^-15 - 7.10542735760100185`*^-15\ x + 1.00000000000000133`\ x\^2\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(?? Fit\n\)\)], "Input"], Cell[BoxData[ \("Fit[data, funs, vars] finds a least-squares fit to a list of data as a \ linear combination of the functions funs of variables vars. The data can have \ the form {{x1, y1, ... , f1}, {x2, y2, ... , f2}, ... }, where the number of \ coordinates x, y, ... is equal to the number of variables in the list vars. \ The data can also be of the form {f1, f2, ... }, with a single coordinate \ assumed to take values 1, 2, ... . The argument funs can be any list of \ functions that depend only on the objects vars."\)], "Print"], Cell[BoxData[ InterpretationBox[ StyleBox[\(Attributes[Fit]\ = \ {Protected}\), ShowStringCharacters->True, NumberMarks->True], InputForm[ Definition[ Fit]], Editable->True, AutoDelete->True]], "Print"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(InterpolatingPolynomial[{{1, 4}, {2, 5}, {4, 6}}, x]\)], "Input"], Cell[BoxData[ \(4 + \((1 + \(2 - x\)\/6)\)\ \((\(-1\) + x)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(?? InterpolatingPolynomial\)], "Input"], Cell[BoxData[ \("InterpolatingPolynomial[data, var] gives a polynomial in the variable \ var which provides an exact fit to a list of data. The data can have the \ forms {{x1, f1}, {x2, f2}, ... } or {f1, f2, ... }, where in the second case, \ the xi are taken to have values 1, 2, ... . The fi can be replaced by {fi, \ dfi, ddfi, ... }, specifying derivatives at the points xi."\)], "Print"], Cell[BoxData[ InterpretationBox[ StyleBox[\(Attributes[InterpolatingPolynomial]\ = \ {Protected}\), ShowStringCharacters->True, NumberMarks->True], InputForm[ Definition[ InterpolatingPolynomial]], Editable->True, AutoDelete->True]], "Print"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(?? Interpolation\n\)\)], "Input"], Cell[BoxData[ \("Interpolation[data] constructs an InterpolatingFunction object which \ represents an approximate function that interpolates the data. The data can \ have the forms {{x1, f1}, {x2, f2}, ... } or {f1, f2, ... }, where in the \ second case, the xi are taken to have values 1, 2, ... ."\)], "Print"], Cell[BoxData[ InterpretationBox[ StyleBox[\(Attributes[Interpolation]\ = \ {Protected}\n\ \n Options[Interpolation]\ = \ {InterpolationOrder\ -> \ 3}\), ShowStringCharacters->True, NumberMarks->True], InputForm[ Definition[ Interpolation]], Editable->True, AutoDelete->True]], "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Rational Interpolation or Pade Approximant with numerator of degree m and \ denominator of degree k.\ \>", "Section"], Cell[BoxData[ \(<< NumericalMath`Approximations`\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?? RationalInterpolation\n\)\)], "Input"], Cell[BoxData[ \("RationalInterpolation[func, {x, m, k}, {x1, x2, ..., xn}, (opts)], (n \ = m+k+1), gives the rational interpolant to func (a function of the variable \ x), where m and k are the degrees of the numerator and denominator, \ respectively, and {x1, x2, ..., xn} is a list of m+k+1 abscissas of the \ interpolation points. An alternative form is \n RationalInterpolation[func, \ {x, m, k}, {x, x0, x1}, (opts)], which specifies the list of abscissas \ implicitly: the abscissas come from the interval (x0,x1). The function func \ must be Listable."\)], "Print"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ri1 = RationalInterpolation[Exp[x], {x, 2, 4}, {0, 1/3, 2/3, 1, 4/3, 5/3, 2}] \)], "Input"], Cell[BoxData[ \(\((1.`15.9546 + 0.379961505998214261578663758`15.9546\ x + 0.0469527572648758848274127968`15.9546\ x\^2)\)/ \((1 - 0.620028516690565809804294869`15.9546\ x + 0.1669139144430911226940924326`15.9546\ x\^2 - 0.02340576618306169162112492488`15.9546\ x\^3 + 0.001452790199322339652547522172`15.9546\ x\^4)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Plot[ri1 - Exp[x], \ {x, 0, 2}]\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.47619 0.222446 272629 [ [.2619 .20995 -9 -9 ] [.2619 .20995 9 0 ] [.5 .20995 -3 -9 ] [.5 .20995 3 0 ] [.7381 .20995 -9 -9 ] [.7381 .20995 9 0 ] [.97619 .20995 -3 -9 ] [.97619 .20995 3 0 ] [.01131 .08613 -39 -6.25 ] [.01131 .08613 0 6.25 ] [.01131 .35876 -33 -6.25 ] [.01131 .35876 0 6.25 ] [.01131 .49508 -33 -6.25 ] [.01131 .49508 0 6.25 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid .2619 .22245 m .2619 .2287 L s [(0.5)] .2619 .20995 0 1 Mshowa .5 .22245 m .5 .2287 L s [(1)] .5 .20995 0 1 Mshowa .7381 .22245 m .7381 .2287 L s [(1.5)] .7381 .20995 0 1 Mshowa .97619 .22245 m .97619 .2287 L s [(2)] .97619 .20995 0 1 Mshowa .125 Mabswid .07143 .22245 m .07143 .2262 L s .11905 .22245 m .11905 .2262 L s .16667 .22245 m .16667 .2262 L s .21429 .22245 m .21429 .2262 L s .30952 .22245 m .30952 .2262 L s .35714 .22245 m .35714 .2262 L s .40476 .22245 m .40476 .2262 L s .45238 .22245 m .45238 .2262 L s .54762 .22245 m .54762 .2262 L s .59524 .22245 m .59524 .2262 L s .64286 .22245 m .64286 .2262 L s .69048 .22245 m .69048 .2262 L s .78571 .22245 m .78571 .2262 L s .83333 .22245 m .83333 .2262 L s .88095 .22245 m .88095 .2262 L s .92857 .22245 m .92857 .2262 L s .25 Mabswid 0 .22245 m 1 .22245 L s .02381 .08613 m .03006 .08613 L s gsave .01131 .08613 -39 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (-5) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-7) show grestore .02381 .35876 m .03006 .35876 L s gsave .01131 .35876 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (5) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-7) show grestore .02381 .49508 m .03006 .49508 L s gsave .01131 .49508 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (1) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-6) show grestore .125 Mabswid .02381 .11339 m .02756 .11339 L s .02381 .14066 m .02756 .14066 L s .02381 .16792 m .02756 .16792 L s .02381 .19518 m .02756 .19518 L s .02381 .24971 m .02756 .24971 L s .02381 .27697 m .02756 .27697 L s .02381 .30424 m .02756 .30424 L s .02381 .3315 m .02756 .3315 L s .02381 .38602 m .02756 .38602 L s .02381 .41329 m .02756 .41329 L s .02381 .44055 m .02756 .44055 L s .02381 .46781 m .02756 .46781 L s .02381 .05887 m .02756 .05887 L s .02381 .03161 m .02756 .03161 L s .02381 .00434 m .02756 .00434 L s .02381 .52234 m .02756 .52234 L s .02381 .5496 m .02756 .5496 L s .02381 .57686 m .02756 .57686 L s .02381 .60413 m .02756 .60413 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02381 .22245 m .03279 .17688 L .04262 .13905 L .04758 .12427 L .05288 .11135 L .05785 .10176 L .06244 .0949 L .06503 .09182 L .06788 .08905 L .07075 .0869 L .07233 .08597 L .07379 .08527 L .07506 .08477 L .07645 .08436 L .07722 .08418 L .07792 .08406 L .0786 .08396 L .07932 .0839 L .07999 .08386 L .08072 .08386 L .08204 .08393 L .08333 .08409 L .08451 .08432 L .08567 .08463 L .08691 .08504 L .08916 .08598 L .09183 .08742 L .09427 .08903 L .10458 .09854 L .11449 .11103 L .12366 .12469 L .14443 .15972 L .16468 .19461 L .18369 .22409 L .20186 .24694 L .21207 .25698 L .22144 .26425 L .22618 .26721 L .23125 .26983 L .23562 .27166 L .24037 .27317 L .24327 .27386 L .24454 .27411 L .24594 .27435 L .24718 .27452 L .2485 .27467 L .24923 .27474 L .2499 .27479 L .25121 .27487 L .25247 .27491 L Mistroke .25363 .27492 L .25491 .2749 L .25562 .27488 L .25627 .27485 L .25746 .27477 L .25874 .27466 L .26107 .27437 L .26345 .27397 L .266 .27344 L .2706 .27218 L .27606 .27023 L .28104 .26804 L .30212 .25517 L .34202 .22179 L .36028 .20648 L .37066 .19867 L .3804 .19224 L .39074 .18657 L .40021 .18259 L .40515 .18102 L .41034 .17976 L .41304 .17928 L .41555 .17892 L .41699 .17877 L .4183 .17865 L .41902 .1786 L .41981 .17855 L .42055 .17852 L .42124 .17849 L .42249 .17847 L .42369 .17847 L .42476 .17849 L .42592 .17854 L .42719 .17861 L .42854 .17871 L .43095 .17898 L .43318 .17931 L .43527 .17969 L .44004 .18082 L .44522 .18247 L .44997 .18436 L .46069 .18989 L .47959 .2035 L .50028 .22273 L .52151 .24476 L .54081 .26417 L .55043 .27276 L .56087 .28072 L .56566 .28378 L .57073 .28656 L Mistroke .57539 .28865 L .57982 .29019 L .58218 .29082 L .58342 .2911 L .58474 .29136 L .58597 .29155 L .58708 .2917 L .58819 .29181 L .58924 .29189 L .59048 .29194 L .59178 .29195 L .59289 .29192 L .59411 .29185 L .59543 .29173 L .59683 .29155 L .59816 .29132 L .59937 .29107 L .60388 .28976 L .60642 .28875 L .60878 .28764 L .61768 .28187 L .62232 .27785 L .62732 .27274 L .63636 .26143 L .65697 .22605 L .69813 .12595 L .71867 .07339 L .72868 .05107 L .73777 .03426 L .74232 .02746 L .74726 .02154 L .74952 .01939 L .75192 .01753 L .75396 .0163 L .75504 .01579 L .75621 .01534 L .75752 .01498 L .75877 .01478 L .75999 .01472 L .76112 .01478 L .76181 .01488 L .76245 .01502 L .76367 .01538 L .76497 .01594 L .76638 .01673 L .7688 .01858 L .77134 .02121 L .7759 .02776 L .78088 .03781 L .78541 .04975 L Mistroke .79559 .08707 L .80502 .13574 L .81394 .19522 L .83391 .37862 L .85234 .60892 L Mfstroke .85234 .60892 m .85294 .61803 L s .97022 .61803 m .97619 .22245 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg?ooo`030000oooooooo02?oool00`000?ooooooo`04oooo000Zoooo00<0 003oooooool0a_ooo`030000oooooooo02?oool00`000?ooooooo`04oooo000Zoooo00<0003ooooo ool0a_ooo`030000oooooooo02?oool00`000?ooooooo`04oooo000Zoooo0P000"], ImageRangeCache->{{{82.875, 369.875}, {433.75, 256.813}} -> {-1.0466, 2.77925, 0.00837206, 1.46231}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[BoxData[ \(Instead, \ we\ can\ let\ Mathematica\ choose\ the\ interpolation\ \(points.\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(ri2 = RationalInterpolation[Exp[x], {x, 2, 4}, {x, 0, 2}]\)], "Input"], Cell[BoxData[ \(\((1.000000157557966967521263947`15.9546 + 0.37982664361058964954449162`15.9546\ x + 0.0468693215807399416119038629`15.9546\ x\^2)\)/ \((1 - 0.62016573039132743382724984`15.9546\ x + 0.1669778816332340770087054875`15.9546\ x\^2 - 0.0234118949793066443664426667`15.9546\ x\^3 + 0.001451916095874725742905961275`15.9546\ x\^4)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Plot[ri2 - Exp[x], {x, 0, 2}]\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.47619 0.328263 150181 [ [.2619 .31576 -9 -9 ] [.2619 .31576 9 0 ] [.5 .31576 -3 -9 ] [.5 .31576 3 0 ] [.7381 .31576 -9 -9 ] [.7381 .31576 9 0 ] [.97619 .31576 -3 -9 ] [.97619 .31576 3 0 ] [.01131 .0279 -39 -6.25 ] [.01131 .0279 0 6.25 ] [.01131 .17808 -39 -6.25 ] [.01131 .17808 0 6.25 ] [.01131 .47844 -33 -6.25 ] [.01131 .47844 0 6.25 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid .2619 .32826 m .2619 .33451 L s [(0.5)] .2619 .31576 0 1 Mshowa .5 .32826 m .5 .33451 L s [(1)] .5 .31576 0 1 Mshowa .7381 .32826 m .7381 .33451 L s [(1.5)] .7381 .31576 0 1 Mshowa .97619 .32826 m .97619 .33451 L s [(2)] .97619 .31576 0 1 Mshowa .125 Mabswid .07143 .32826 m .07143 .33201 L s .11905 .32826 m .11905 .33201 L s .16667 .32826 m .16667 .33201 L s .21429 .32826 m .21429 .33201 L s .30952 .32826 m .30952 .33201 L s .35714 .32826 m .35714 .33201 L s .40476 .32826 m .40476 .33201 L s .45238 .32826 m .45238 .33201 L s .54762 .32826 m .54762 .33201 L s .59524 .32826 m .59524 .33201 L s .64286 .32826 m .64286 .33201 L s .69048 .32826 m .69048 .33201 L s .78571 .32826 m .78571 .33201 L s .83333 .32826 m .83333 .33201 L s .88095 .32826 m .88095 .33201 L s .92857 .32826 m .92857 .33201 L s .25 Mabswid 0 .32826 m 1 .32826 L s .02381 .0279 m .03006 .0279 L s gsave .01131 .0279 -39 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (-2) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-6) show grestore .02381 .17808 m .03006 .17808 L s gsave .01131 .17808 -39 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (-1) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-6) show grestore .02381 .47844 m .03006 .47844 L s gsave .01131 .47844 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (1) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-6) show grestore .125 Mabswid .02381 .05794 m .02756 .05794 L s .02381 .08797 m .02756 .08797 L s .02381 .11801 m .02756 .11801 L s .02381 .14805 m .02756 .14805 L s .02381 .20812 m .02756 .20812 L s .02381 .23815 m .02756 .23815 L s .02381 .26819 m .02756 .26819 L s .02381 .29823 m .02756 .29823 L s .02381 .3583 m .02756 .3583 L s .02381 .38834 m .02756 .38834 L s .02381 .41837 m .02756 .41837 L s .02381 .44841 m .02756 .44841 L s .02381 .50848 m .02756 .50848 L s .02381 .53852 m .02756 .53852 L s .02381 .56855 m .02756 .56855 L s .02381 .59859 m .02756 .59859 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02381 .35193 m .03279 .3333 L .04262 .31852 L .04758 .31304 L .05288 .30849 L .05546 .30673 L .05789 .30532 L .06244 .30332 L .06487 .30257 L .06627 .30223 L .06754 .30197 L .06889 .30176 L .06958 .30168 L .07034 .3016 L .07103 .30155 L .07168 .30151 L .07295 .30148 L .07422 .30149 L .07541 .30154 L .07648 .30162 L .07764 .30174 L .07891 .30191 L .08025 .30214 L .08265 .30265 L .08806 .30428 L .0931 .30631 L .10458 .31241 L .14349 .33964 L .1633 .35226 L .17375 .3576 L .17947 .36003 L .18485 .36198 L .18998 .36351 L .19469 .36462 L .19991 .36553 L .20266 .36587 L .20404 .366 L .20555 .36611 L .20684 .36619 L .20806 .36624 L .20915 .36627 L .21034 .36628 L .21164 .36628 L .21229 .36627 L .21301 .36625 L .21428 .3662 L .21547 .36613 L .21779 .36595 L .21991 .36573 L .22469 .365 L Mistroke .22948 .36399 L .23457 .36261 L .24372 .35936 L .2546 .35429 L .26446 .34871 L .30457 .31939 L .34318 .28944 L .36315 .27698 L .37321 .27211 L .37887 .26986 L .38423 .26808 L .38914 .26678 L .39172 .26622 L .39446 .26573 L .39677 .26539 L .39799 .26524 L .39931 .26511 L .40051 .265 L .40162 .26493 L .40267 .26487 L .40378 .26483 L .4051 .2648 L .40629 .2648 L .40695 .26481 L .40767 .26482 L .40897 .26487 L .41018 .26494 L .41147 .26503 L .41378 .26526 L .41647 .26564 L .41899 .26608 L .42471 .26746 L .4297 .26907 L .43514 .27125 L .44504 .27637 L .46359 .2898 L .48334 .30898 L .50434 .33362 L .54205 .38237 L .563 .40735 L .57218 .41684 L .58222 .42581 L .59197 .4328 L .59664 .43547 L .60105 .43756 L .60345 .4385 L .60607 .43938 L .60824 .43998 L .61063 .4405 L .61188 .44072 L Mistroke .61322 .44091 L .61396 .44099 L .61464 .44105 L .61532 .4411 L .61596 .44114 L .61717 .44118 L .61848 .44118 L .61964 .44114 L .62088 .44106 L .62212 .44094 L .62349 .44075 L .62491 .4405 L .62624 .44022 L .63118 .43872 L .63399 .43754 L .63655 .43627 L .64231 .43271 L .65267 .42376 L .66248 .41231 L .68241 .3804 L .70083 .34182 L .74012 .24299 L .75983 .19514 L .76845 .17698 L .77789 .16017 L .78732 .1474 L .78981 .14482 L .79247 .14246 L .79496 .14064 L .79725 .13931 L .79854 .13871 L .79995 .13818 L .80137 .13779 L .80215 .13763 L .80287 .13752 L .80356 .13746 L .80421 .13742 L .80543 .13745 L .80675 .13759 L .80749 .13773 L .80818 .13789 L .80948 .13831 L .81071 .13881 L .81299 .14007 L .81545 .14188 L .81811 .1444 L .82305 .15063 L .82765 .15833 L .8381 .18278 L .84829 .21602 L Mistroke .85919 .26154 L .87905 .36652 L .90062 .494 L .91069 .54688 L .91621 .57069 L .9189 .58048 L .92137 .58814 L .92347 .59358 L .92576 .59826 L .92702 .60022 L .92819 .60163 L .9293 .6026 L .93049 .6032 L .93179 .60332 L .93253 .60313 L .9332 .60278 L .93445 .60171 L .93576 .59994 L .93694 .59776 L .93804 .5952 L .94054 .5874 L .94292 .57724 L .94512 .56527 L .95011 .52808 L .95451 .48223 L .95918 .41814 L .96727 .26351 L .97187 .14656 L .97619 .01472 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgcoool00`000?oo ooooo`04oooo0006oooo00<0003oooooool00oooo`030000oooo000000?oool01`000?oooooooooo o`000?ooo`00000Aoooo00<0003oooooool0k?ooo`030000oooooooo00Coool00004oooo00000000 00000oooo`030000oooooooo00?oool00`000?ooooooo`02oooo00L0003oooooooooool0003oool0 00004Oooo`030000oooooooo0>coool00`000?ooooooo`04oooo0007oooo00<0003oooooool00_oo o`030000oooo000000?oool01`000?ooooooooooo`000?ooo`00000Aoooo00<0003oooooool0k?oo o`030000oooooooo00Coool000Koool00`000?ooo`000002oooo00D0003oooooooooool000000_oo o`070000oooooooooooo0000oooo0000017oool20000kOooo`030000oooooooo00Coool000Ooool0 0`000?ooooooo`07oooo0P0000Coool00`000?ooooooo`0@oooo00<0003oooooool0k?ooo`030000 oooooooo00Coool0023oool300001oooo`030000oooooooo0>coool00`000?ooooooo`04oooo000P oooo00<0003oool000001oooo`030000oooooooo0>coool00`000?ooooooo`04oooo000Koooo0`00 00;oool300001oooo`030000oooooooo0>coool00`000?ooooooo`04oooo000Poooo00<0003ooooo ool01oooo`030000oooooooo0>_oool00`000?ooooooo`05oooo000Qoooo0P0000Ooool00`000?oo ooooo`3[oooo00<0003oooooool01Oooo`00:_ooo`80003/oooo00<0003oooooool01Oooo`00:_oo o`030000oooooooo0>_oool00`000?ooooooo`05oooo000Zoooo00<0003oooooool0joooo`030000 oooooooo00Goool002[oool00`000?ooooooo`3[oooo00<0003oooooool01Oooo`00:_ooo`030000 oooooooo0>_oool00`000?ooooooo`05oooo000Zoooo00<0003oooooool0joooo`030000oooooooo 00Goool002[oool00`000?ooooooo`3[oooo00<0003oooooool01Oooo`00:_ooo`030000oooooooo 0>_oool00`000?ooooooo`05oooo000Zoooo0P000>coool00`000?ooooooo`05oooo000Zoooo00<0 003oooooool0joooo`030000oooooooo00Goool002[oool00`000?ooooooo`3[oooo00<0003ooooo ool01Oooo`00:_ooo`030000oooooooo0>_oool00`000?ooooooo`05oooo000Zoooo00<0003ooooo ool0joooo`030000oooooooo00Goool002[oool00`000?ooooooo`3[oooo00<0003oooooool01Ooo o`00:_ooo`030000oooooooo0>[oool00`000?ooooooo`06oooo000Zoooo0P000>_oool00`000?oo ooooo`06oooo000Zoooo00<0003oooooool0j_ooo`030000oooooooo00Koool002[oool00`000?oo ooooo`3Zoooo00<0003oooooool01_ooo`00:_ooo`030000oooooooo0>[oool00`000?ooooooo`06 oooo000Zoooo00<0003oooooool0j_ooo`030000oooooooo00Koool002[oool00`000?ooooooo`2o oooo100002Ooool00`000?ooooooo`06oooo000Zoooo00<0003oooooool0_Oooo`800004oooo00<0 003oooooool09?ooo`030000oooooooo00Koool002[oool00`000?ooooooo`2loooo0P0000Koool0 0`000?ooooooo`0Soooo00<0003oooooool01_ooo`00:_ooo`80002moooo00<0003oooooool01_oo o`030000oooooooo02;oool00`000?ooooooo`06oooo000Zoooo00<0003oooooool0^oooo`030000 oooooooo00Soool00`000?ooooooo`0Qoooo00<0003oooooool01_ooo`00:_ooo`030000oooooooo 0;_oool00`000?ooooooo`08oooo00<0003oooooool08Oooo`030000oooooooo00Koool002[oool0 0`000?ooooooo`2joooo00<0003oooooool02_ooo`030000oooooooo023oool00`000?ooooooo`06 oooo0005oooo1@000003oooo0000oooo00;oool600000_ooo`030000oooooooo013oool00`000?oo ooooo`2ioooo00<0003oooooool02oooo`030000oooooooo023oool00`000?ooooooo`06oooo0007 oooo00<0003oooooool00_ooo`030000oooo000000?oool01`000?ooooooooooo`000?ooo`00000A oooo00<0003oooooool0^Oooo`030000oooooooo00coool00`000?ooooooo`0Ooooo00<0003ooooo ool01_ooo`0000Coool0000000000003oooo00<0003oooooool00oooo`030000oooooooo00;oool0 1`000?ooooooooooo`000?ooo`00000Aoooo00<0003oooooool0^?ooo`030000oooooooo00goool0 0`000?ooooooo`0Ooooo00<0003oooooool01_ooo`001oooo`030000oooooooo00;oool00`000?oo o`000003oooo00L0003oooooooooool0003oool000004Oooo`80002hoooo00<0003oooooool03ooo o`030000oooooooo01koool00`000?ooooooo`06oooo0007oooo00D0003oooooooooool000000ooo o`040000oooooooo00000oooo`030000oooo0000017oool00`000?ooooooo`2goooo00<0003ooooo ool03oooo`030000oooooooo01koool00`000?ooooooo`06oooo0006oooo0P0000Woool200001?oo o`030000oooooooo00Koool300001oooo`030000oooooooo0;Koool00`000?ooooooo`0@oooo00<0 003oooooool07_ooo`030000oooooooo00Koool0023oool00`000?ooo`000007oooo00<0003ooooo ool0]_ooo`030000oooooooo013oool00`000?ooooooo`0Noooo00<0003oooooool01_ooo`006ooo o`<00002oooo0`0000Ooool00`000?ooooooo`2eoooo00<0003oooooool04_ooo`030000oooooooo 01goool00`000?ooooooo`06oooo000Poooo00<0003oooooool01oooo`030000oooooooo0;Goool0 0`000?ooooooo`0Boooo00<0003oooooool07Oooo`030000oooooooo00Koool0027oool200001ooo o`030000oooooooo0;Coool00`000?ooooooo`0Coooo00<0003oooooool07Oooo`030000oooooooo 00Koool002[oool00`000?ooooooo`2doooo00<0003oooooool04oooo`030000oooooooo01goool0 0`000?ooooooo`06oooo000Zoooo0P000;Coool00`000?ooooooo`0Eoooo00<0003oooooool06ooo o`030000oooooooo00Ooool002[oool00`000?ooooooo`2coooo00<0003oooooool05Oooo`030000 oooooooo01_oool00`000?ooooooo`07oooo000Zoooo00<0003oooooool0/oooo`030000oooooooo 01Goool00`000?ooooooo`0Koooo00<0003oooooool01oooo`00:_ooo`030000oooooooo0;;oool0 0`000?ooooooo`0Foooo00<0003oooooool06oooo`030000oooooooo00Ooool002[oool00`000?oo ooooo`2boooo00<0003oooooool05oooo`030000oooooooo01[oool00`000?ooooooo`07oooo000Z oooo00<0003oooooool0/Oooo`030000oooooooo01Soool00`000?ooooooo`0Joooo00<0003ooooo ool01oooo`00:_ooo`030000oooooooo0;7oool00`000?ooooooo`0Hoooo00<0003oooooool06_oo o`030000oooooooo00Ooool002[oool20000/Oooo`030000oooooooo01Woool00`000?ooooooo`0J oooo00<0003oooooool01oooo`00:_ooo`030000oooooooo0;3oool00`000?ooooooo`0Joooo00<0 003oooooool06Oooo`030000oooooooo00Ooool002[oool00`000?ooooooo`2`oooo00<0003ooooo ool06_ooo`030000oooooooo01Woool00`000?ooooooo`07oooo000Zoooo00<0003oooooool0[ooo o`030000oooooooo01_oool00`000?ooooooo`0Ioooo00<0003oooooool01oooo`00:_ooo`030000 oooooooo0:ooool00`000?ooooooo`0Koooo00<0003oooooool06Oooo`030000oooooooo00Ooool0 02[oool00`000?ooooooo`2^oooo00<0003oooooool07Oooo`030000oooooooo01Soool00`000?oo ooooo`07oooo000Zoooo00<0003oooooool0[_ooo`030000oooooooo01goool00`000?ooooooo`0H oooo00<0003oooooool01oooo`00:_ooo`030000oooooooo05Soool90000COooo`030000oooooooo 01goool00`000?ooooooo`0Hoooo00<0003oooooool01oooo`00:_ooo`80001Foooo0`0000Woool3 0000BOooo`030000oooooooo01koool00`000?ooooooo`0Hoooo00<0003oooooool01oooo`00:_oo o`030000oooooooo05?oool200003oooo`800017oooo00<0003oooooool07oooo`030000oooooooo 01Ooool00`000?ooooooo`07oooo000Zoooo00<0003oooooool0D_ooo`030000oooooooo017oool0 0`000?ooooooo`13oooo00<0003oooooool08?ooo`030000oooooooo01Ooool00`000?ooooooo`07 oooo000Zoooo00<0003oooooool0D?ooo`80000Eoooo0P0004?oool00`000?ooooooo`0Poooo00<0 003oooooool05oooo`030000oooooooo00Ooool002[oool00`000?ooooooo`0coooo00<0003ooooo ool00_ooo`030000oooooooo00;oool00`000?ooooooo`0>oooo0P0001Woool200001oooo`D0000a oooo1@0000;oool00`000?ooooooo`02oooo00<0003oooooool06?ooo`030000oooooooo01Ooool4 00001_ooo`00:_ooo`030000oooooooo03;oool00`000?ooo`000007oooo00<0003oool000003_oo o`030000oooooooo01_oool00`000?ooooooo`06oooo00<0003oooooool0oooo`030000oooooooo02coool00`000?ooooooo`0Yoooo00<0003oooooool04oooo`030000 oooooooo00Soool002[oool00`000?ooooooo`0Noooo0P0001Soool20000?_ooo`030000oooooooo 02[oool00`000?ooooooo`0Zoooo00<0003oooooool04oooo`030000oooooooo00Soool002[oool0 0`000?ooooooo`0Poooo0P0001Goool00`000?ooooooo`0ooooo00<0003oooooool0:Oooo`030000 oooooooo02[oool00`000?ooooooo`0Coooo00<0003oooooool02?ooo`00:_ooo`80000Soooo0P00 017oool20000@oooo`030000oooooooo02Ooool00`000?ooooooo`0/oooo00<0003oooooool04_oo o`030000oooooooo00Soool002[oool00`000?ooooooo`0Toooo0P0000_oool40000A_ooo`030000 oooooooo02Koool00`000?ooooooo`0/oooo00<0003oooooool04_ooo`030000oooooooo00Soool0 02[oool00`000?ooooooo`0Voooo2`0004_oool00`000?ooooooo`0Uoooo00<0003oooooool0;?oo o`030000oooooooo01;oool00`000?ooooooo`08oooo000Zoooo00<0003oooooool0O?ooo`030000 oooooooo02Coool00`000?ooooooo`0]oooo00<0003oooooool04_ooo`030000oooooooo00Soool0 02[oool00`000?ooooooo`1moooo00<0003oooooool08oooo`030000oooooooo02goool00`000?oo ooooo`0Boooo00<0003oooooool02?ooo`00:_ooo`030000oooooooo07koool00`000?ooooooo`0Q oooo00<0003oooooool0;oooo`030000oooooooo017oool00`000?ooooooo`08oooo000Zoooo00<0 003oooooool0Ooooo`030000oooooooo023oool00`000?ooooooo`0_oooo00<0003oooooool04?oo o`030000oooooooo00Woool002[oool00`000?ooooooo`20oooo00<0003oooooool07_ooo`030000 oooooooo033oool00`000?ooooooo`0@oooo00<0003oooooool02Oooo`00:_ooo`800021oooo00<0 003oooooool07_ooo`030000oooooooo033oool00`000?ooooooo`0@oooo00<0003oooooool02Ooo o`00:_ooo`030000oooooooo087oool00`000?ooooooo`0Loooo00<0003oooooool0oooo00<0 003oooooool02Oooo`00:_ooo`030000oooooooo08Woool00`000?ooooooo`0>oooo00<0003ooooo ool0>Oooo`030000oooooooo00koool00`000?ooooooo`09oooo000Zoooo00<0003oooooool0R_oo o`030000oooooooo00coool00`000?ooooooo`0joooo00<0003oooooool03_ooo`030000oooooooo 00Woool002[oool00`000?ooooooo`2;oooo00<0003oooooool02Oooo`80000moooo00<0003ooooo ool03_ooo`030000oooooooo00Woool002[oool00`000?ooooooo`2 oooo2000047oool00`000?ooooooo`0=oooo00<0003oooooool02Oooo`00:_ooo`030000oooooooo 0=Ooool00`000?ooooooo`0=oooo00<0003oooooool02Oooo`00:_ooo`80003Hoooo00<0003ooooo ool03Oooo`030000oooooooo00Woool002[oool00`000?ooooooo`3Goooo00<0003oooooool03?oo o`030000oooooooo00[oool002[oool00`000?ooooooo`3Goooo00<0003oooooool03?ooo`030000 oooooooo00[oool002[oool00`000?ooooooo`3Hoooo00<0003oooooool02oooo`030000oooooooo 00[oool000Koool5000000?oool0003oool00_ooo`H00002oooo00<0003oooooool03oooo`030000 oooooooo0=Soool00`000?ooooooo`0;oooo00<0003oooooool02_ooo`002?ooo`030000oooooooo 00;oool00`000?ooo`000003oooo00L0003oooooooooool0003oool000004?ooo`030000oooooooo 0=Soool00`000?ooooooo`0;oooo00<0003oooooool02_ooo`002?ooo`030000oooooooo00?oool0 0`000?ooooooo`02oooo00L0003oooooooooool0003oool000004?ooo`030000oooooooo0=Soool0 0`000?ooooooo`0;oooo00<0003oooooool02_ooo`002?ooo`030000oooooooo00;oool00`000?oo o`000003oooo00L0003oooooooooool0003oool000004?ooo`030000oooooooo0=Soool00`000?oo ooooo`0;oooo00<0003oooooool02_ooo`002?ooo`050000oooooooooooo000000?oool010000?oo ooooo`0000?oool00`000?ooo`00000@oooo0P000=Woool00`000?ooooooo`0;oooo00<0003ooooo ool02_ooo`001oooo`800009oooo0P0000Coool00`000?ooooooo`0?oooo00<0003oooooool0fOoo o`030000oooooooo00[oool00`000?ooooooo`0:oooo000Qoooo0`0000Koool00`000?ooooooo`3I oooo00<0003oooooool02_ooo`030000oooooooo00[oool0027oool00`000?ooo`000006oooo00<0 003oooooool0fOooo`030000oooooooo00[oool00`000?ooooooo`0:oooo000Loooo0`0000;oool3 00001_ooo`030000oooooooo0=Woool00`000?ooooooo`0:oooo00<0003oooooool02_ooo`008Ooo o`030000oooooooo00Koool00`000?ooooooo`3Ioooo00<0003oooooool02_ooo`030000oooooooo 00[oool002;oool200001_ooo`030000oooooooo0=Woool00`000?ooooooo`0:oooo00<0003ooooo ool02_ooo`00:_ooo`80003Koooo00<0003oooooool02?ooo`030000oooooooo00_oool002[oool0 0`000?ooooooo`3Joooo00<0003oooooool02?ooo`030000oooooooo00_oool002[oool00`000?oo ooooo`3Joooo00<0003oooooool02?ooo`030000oooooooo00_oool002[oool00`000?ooooooo`3J oooo00<0003oooooool02?ooo`030000oooooooo00_oool002[oool00`000?ooooooo`3Joooo00<0 003oooooool02?ooo`030000oooooooo00_oool002[oool00`000?ooooooo`3Joooo00<0003ooooo ool02?ooo`030000oooooooo00_oool002[oool00`000?ooooooo`3Joooo00<0003oooooool02?oo o`030000oooooooo00_oool002[oool00`000?ooooooo`3Koooo00<0003oooooool01oooo`030000 oooooooo00_oool002[oool20000g?ooo`030000oooooooo00Ooool00`000?ooooooo`0;oooo000Z oooo00<0003oooooool0foooo`030000oooooooo00Ooool00`000?ooooooo`0;oooo000Zoooo00<0 003oooooool0foooo`030000oooooooo00Koool00`000?ooooooo`0"], ImageRangeCache->{{{82.875, 369.875}, {257.188, 80.25}} -> {-1.0466, -3.4634, 0.00837206, 2.65459}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[BoxData[ \(To\ minimize\ the\ maximum\ error\ over\ the\ entire\ range\ of\ interpolation, \ use\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?MiniMaxApproximation\)\)], "Input"], Cell[BoxData[ \("MiniMaxApproximation[func, {x, {x0, x1}, m, k}, (opts)] finds the \ mini-max approximation to func (a function of the variable x) on the interval \ (x0, x1), where m and k are the degrees of the numerator and denominator, \ respectively. The answer returned is {AbscissaList, {Approximation, \ MaxError}}, where AbscissaList is a list of the abscissas where the maximum \ error occurs, Approximation is the rational approximation desired, and \ MaxError is the value of the mini-max error. The function func must be \ Listable.\n MiniMaxApproximation[f, approx, {x, {x0, x1}, m, k}, (opts)] is a \ form that allows the user to start the iteration from a known approximation. \ Here approx must be in the form of an answer returned by \ MiniMaxApproximation."\)], "Print"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(mmlist\ = \ MiniMaxApproximation[Exp[x], {x, {0, 2}, 2, 4}]\)], "Input"], Cell[BoxData[ \({{0, 0.1063486487562546614033124614`15.9546, 0.400915308263480259506883759`15.9546, 0.816636147986608924532073332`15.9546, 1.262697793261275633137973799`15.9546, 1.649749456495247482297801461`15.9546, 1.909120496974659570292876186`15.9546, 2.`15.9546}, { \((1.000000206052106264524645031`15.9546 + 0.38088147329891583439907302`15.9546\ x + 0.0472394925260850347699151898`15.9546\ x\^2)\)/ \((1 - 0.619109229762559169793026569`15.9546\ x + 0.166282879427172251466502717`15.9546\ x\^2 - 0.02323044886220642823609590194`15.9546\ x\^3 + 0.00143324892348305835232108387`15.9546\ x\^4)\), \(-2.060521062645246450310702021`15.9546*^-7\)}}\)], "Output"] }, Open ]], Cell[BoxData[ \(To\ extract\ just\ the\ rational\ approximation\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(mmfunc\ = \ mmlist[\([2, 1]\)]\)], "Input"], Cell[BoxData[ \(\((1.000000206052106264524645031`15.9546 + 0.38088147329891583439907302`15.9546\ x + 0.0472394925260850347699151898`15.9546\ x\^2)\)/ \((1 - 0.619109229762559169793026569`15.9546\ x + 0.166282879427172251466502717`15.9546\ x\^2 - 0.02323044886220642823609590194`15.9546\ x\^3 + 0.00143324892348305835232108387`15.9546\ x\^4)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Plot[1 - mmfunc/Exp[x], {x, 0, 2}]\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.47619 0.309016 1.42829e+006 [ [.2619 .29652 -9 -9 ] [.2619 .29652 9 0 ] [.5 .29652 -3 -9 ] [.5 .29652 3 0 ] [.7381 .29652 -9 -9 ] [.7381 .29652 9 0 ] [.97619 .29652 -3 -9 ] [.97619 .29652 3 0 ] [.01131 .02336 -39 -6.25 ] [.01131 .02336 0 6.25 ] [.01131 .16619 -39 -6.25 ] [.01131 .16619 0 6.25 ] [.01131 .45185 -33 -6.25 ] [.01131 .45185 0 6.25 ] [.01131 .59467 -33 -6.25 ] [.01131 .59467 0 6.25 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid .2619 .30902 m .2619 .31527 L s [(0.5)] .2619 .29652 0 1 Mshowa .5 .30902 m .5 .31527 L s [(1)] .5 .29652 0 1 Mshowa .7381 .30902 m .7381 .31527 L s [(1.5)] .7381 .29652 0 1 Mshowa .97619 .30902 m .97619 .31527 L s [(2)] .97619 .29652 0 1 Mshowa .125 Mabswid .07143 .30902 m .07143 .31277 L s .11905 .30902 m .11905 .31277 L s .16667 .30902 m .16667 .31277 L s .21429 .30902 m .21429 .31277 L s .30952 .30902 m .30952 .31277 L s .35714 .30902 m .35714 .31277 L s .40476 .30902 m .40476 .31277 L s .45238 .30902 m .45238 .31277 L s .54762 .30902 m .54762 .31277 L s .59524 .30902 m .59524 .31277 L s .64286 .30902 m .64286 .31277 L s .69048 .30902 m .69048 .31277 L s .78571 .30902 m .78571 .31277 L s .83333 .30902 m .83333 .31277 L s .88095 .30902 m .88095 .31277 L s .92857 .30902 m .92857 .31277 L s .25 Mabswid 0 .30902 m 1 .30902 L s .02381 .02336 m .03006 .02336 L s gsave .01131 .02336 -39 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (-2) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-7) show grestore .02381 .16619 m .03006 .16619 L s gsave .01131 .16619 -39 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (-1) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-7) show grestore .02381 .45185 m .03006 .45185 L s gsave .01131 .45185 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (1) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-7) show grestore .02381 .59467 m .03006 .59467 L s gsave .01131 .59467 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (2) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-7) show grestore .125 Mabswid .02381 .05192 m .02756 .05192 L s .02381 .08049 m .02756 .08049 L s .02381 .10906 m .02756 .10906 L s .02381 .13762 m .02756 .13762 L s .02381 .19475 m .02756 .19475 L s .02381 .22332 m .02756 .22332 L s .02381 .25188 m .02756 .25188 L s .02381 .28045 m .02756 .28045 L s .02381 .33758 m .02756 .33758 L s .02381 .36615 m .02756 .36615 L s .02381 .39471 m .02756 .39471 L s .02381 .42328 m .02756 .42328 L s .02381 .48041 m .02756 .48041 L s .02381 .50898 m .02756 .50898 L s .02381 .53754 m .02756 .53754 L s .02381 .56611 m .02756 .56611 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02381 .01472 m .03279 .23369 L .04262 .4049 L .04758 .46786 L .05288 .52003 L .05789 .55634 L .06244 .57962 L .06507 .58918 L .0665 .5933 L .06785 .59649 L .06908 .59886 L .07022 .60057 L .07148 .60198 L .07217 .60254 L .07283 .60293 L .07354 .6032 L .07418 .60331 L .07493 .60329 L .07563 .60312 L .07687 .60247 L .07757 .60192 L .07822 .6013 L .07958 .59962 L .08108 .59723 L .0841 .59078 L .08922 .57539 L .09469 .55352 L .10458 .50325 L .14594 .23958 L .15598 .18139 L .16664 .12743 L .17651 .08605 L .18579 .05551 L .19455 .0344 L .19948 .02589 L .20401 .02021 L .20643 .01799 L .20868 .01645 L .20998 .01578 L .21116 .01531 L .21181 .01511 L .21252 .01494 L .2138 .01475 L .21504 .01472 L .21575 .01476 L .21641 .01485 L .2171 .01498 L .21774 .01514 L .21917 .01564 L .22038 .0162 L Mistroke .22171 .01697 L .22412 .01878 L .22668 .02126 L .22945 .02458 L .23446 .03221 L .24383 .05175 L .25476 .08234 L .26481 .1167 L .30435 .28717 L .34238 .4543 L .36193 .52288 L .3728 .55299 L .38285 .57492 L .38771 .58334 L .39305 .59095 L .39843 .59679 L .40137 .59921 L .40413 .60097 L .4066 .60213 L .40784 .60257 L .40855 .60277 L .40919 .60293 L .41036 .60315 L .4116 .60328 L .41275 .60332 L .41383 .60328 L .4151 .60313 L .41574 .60302 L .41644 .60287 L .41768 .60252 L .41884 .60211 L .42146 .60086 L .42427 .59904 L .42954 .59429 L .43442 .58837 L .44548 .56975 L .45523 .54768 L .4656 .51894 L .50331 .37978 L .54348 .20898 L .56194 .13808 L .5724 .10304 L .58214 .07476 L .59262 .04975 L .60216 .03247 L .6074 .02536 L .61235 .02027 L .61507 .01816 L .61758 .01666 L .61889 .01604 L Mistroke .62033 .0155 L .62106 .01528 L .62184 .01508 L .62257 .01493 L .62325 .01483 L .62454 .01472 L .62575 .01473 L .62684 .01482 L .62802 .01501 L .62931 .01533 L .63067 .01579 L .6331 .01694 L .63536 .01838 L .63746 .02004 L .64225 .02497 L .64748 .03218 L .65225 .04039 L .66298 .06438 L .68208 .12474 L .70284 .21134 L .74364 .41067 L .76391 .50129 L .77384 .53819 L .78293 .56591 L .79179 .58623 L .79624 .59364 L .79875 .59692 L .80111 .59942 L .80235 .60048 L .80368 .60145 L .8048 .6021 L .80602 .60266 L .8067 .6029 L .80735 .60307 L .80802 .60321 L .80875 .60329 L .81006 .60329 L .81126 .60312 L .81246 .60277 L .81358 .6023 L .8148 .60161 L .8161 .60068 L .81852 .59841 L .8207 .59574 L .8256 .58761 L .83076 .57579 L .84 .54632 L .84994 .50293 L .86084 .44257 L .88204 .29596 L Mistroke .90134 .15222 L .91094 .08933 L .916 .0616 L .92137 .03795 L .92363 .03014 L .92606 .02336 L .92745 .0203 L .92872 .01805 L .93 .01634 L .93122 .01527 L .93231 .01479 L .93347 .01478 L .93457 .01527 L .93557 .01615 L .93674 .01775 L .93801 .02018 L .93921 .02316 L .94032 .02652 L .94255 .03521 L .94458 .04545 L .94916 .07729 L .9534 .11882 L .95742 .17014 L .96209 .24607 L .96641 .33369 L .97619 .60332 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgoooo`030000oooooooo00Koool0 0`000?ooooooo`0;oooo000Zoooo0P0002Goool00`000?ooooooo`0@oooo00<0003oooooool0Cooo o`030000oooooooo01;oool00`000?ooooooo`0joooo00<0003oooooool01oooo`030000oooooooo 00_oool002[oool200009?ooo`030000oooooooo017oool00`000?ooooooo`1>oooo00<0003ooooo ool04oooo`030000oooooooo03[oool00`000?ooooooo`07oooo00<0003oooooool02oooo`00:_oo o`80000Toooo00<0003oooooool04Oooo`030000oooooooo04koool00`000?ooooooo`0Doooo00<0 003oooooool0>Oooo`030000oooooooo00Ooool00`000?ooooooo`0;oooo000Zoooo0P0002Coool0 0`000?ooooooo`0Aoooo00<0003oooooool0COooo`030000oooooooo01Goool00`000?ooooooo`0i oooo00<0003oooooool01oooo`030000oooooooo00_oool002[oool200008oooo`030000oooooooo 01?oool00`000?ooooooo`1?ooo`030000oooooooo00Ooool0 0`000?ooooooo`0;oooo000Zoooo0P0002?oool00`000?ooooooo`0Coooo00<0003oooooool0Booo o`030000oooooooo01Ooool00`000?ooooooo`0hoooo00<0003oooooool02?ooo`030000oooooooo 00[oool002[oool200008oooo`030000oooooooo01?oool00`000?ooooooo`1;oooo00<0003ooooo ool05oooo`030000oooooooo03Soool00`000?ooooooo`08oooo00<0003oooooool02_ooo`00:_oo o`80000Roooo00<0003oooooool05Oooo`030000oooooooo04[oool00`000?ooooooo`0Hoooo00<0 003oooooool0=_ooo`030000oooooooo00Woool00`000?ooooooo`0:oooo000Zoooo0P0002;oool0 0`000?ooooooo`0Eoooo00<0003oooooool0BOooo`030000oooooooo01Woool00`000?ooooooo`0f oooo00<0003oooooool02Oooo`030000oooooooo00[oool002[oool200008_ooo`030000oooooooo 01Goool00`000?ooooooo`19oooo00<0003oooooool06Oooo`030000oooooooo03Koool00`000?oo ooooo`09oooo00<0003oooooool02_ooo`00:_ooo`80000Qoooo00<0003oooooool05_ooo`030000 oooooooo04Woool00`000?ooooooo`0Joooo00<0003oooooool0=Oooo`030000oooooooo00Woool0 0`000?ooooooo`0:oooo000Zoooo0P00027oool00`000?ooooooo`0Foooo00<0003oooooool0BOoo o`030000oooooooo01[oool00`000?ooooooo`0eoooo00<0003oooooool02Oooo`030000oooooooo 00[oool002[oool200008Oooo`030000oooooooo01Ooool00`000?ooooooo`17oooo00<0003ooooo ool06oooo`030000oooooooo03Goool00`000?ooooooo`09oooo00<0003oooooool02_ooo`00:_oo o`80000Qoooo00<0003oooooool05oooo`030000oooooooo04Ooool00`000?ooooooo`0Koooo00<0 003oooooool0=Oooo`030000oooooooo00Woool00`000?ooooooo`0:oooo000Zoooo0P00027oool0 0`000?ooooooo`0Goooo00<0003oooooool0Aoooo`030000oooooooo01coool00`000?ooooooo`0d oooo00<0003oooooool02Oooo`030000oooooooo00[oool002[oool200008?ooo`030000oooooooo 01Soool00`000?ooooooo`17oooo00<0003oooooool07?ooo`030000oooooooo03?oool00`000?oo ooooo`0:oooo00<0003oooooool02_ooo`00:_ooo`80000Poooo00<0003oooooool06Oooo`030000 oooooooo04Goool00`000?ooooooo`0Moooo00<0003oooooool0oooo00<0 003oooooool02?ooo`00:_ooo`030000oooo000001coool00`000?ooooooo`0Poooo00<0003ooooo ool0??ooo`030000oooooooo02Koool00`000?ooooooo`0]oooo00<0003oooooool03_ooo`030000 oooooooo00Soool002[oool300007?ooo`030000oooooooo023oool00`000?ooooooo`0loooo00<0 003oooooool09oooo`030000oooooooo02coool00`000?ooooooo`0>oooo00<0003oooooool02?oo o`00:_ooo`030000oooo000001coool00`000?ooooooo`0Qoooo00<0003oooooool0>oooo`030000 oooooooo02Ooool00`000?ooooooo`0[oooo00<0003oooooool03oooo`030000oooooooo00Soool0 02[oool00`000?ooo`00000Koooo00<0003oooooool08_ooo`030000oooooooo03_oool00`000?oo ooooo`0Woooo00<0003oooooool0:oooo`030000oooooooo00ooool00`000?ooooooo`08oooo000Z oooo00<0003oool000006oooo`030000oooooooo02;oool00`000?ooooooo`0joooo00<0003ooooo ool0:?ooo`030000oooooooo02_oool00`000?ooooooo`0?oooo00<0003oooooool02?ooo`00:_oo o`030000oooo000001_oool00`000?ooooooo`0Roooo00<0003oooooool0>_ooo`030000oooooooo 02Soool00`000?ooooooo`0[oooo00<0003oooooool03oooo`030000oooooooo00Soool002[oool0 0`000?ooo`00000Koooo00<0003oooooool08_ooo`030000oooooooo03[oool00`000?ooooooo`0Y oooo00<0003oooooool0:_ooo`030000oooooooo00ooool00`000?ooooooo`08oooo000Zoooo00<0 003oool000006oooo`030000oooooooo02?oool00`000?ooooooo`0ioooo00<0003oooooool0:Ooo o`030000oooooooo02[oool00`000?ooooooo`0?oooo00<0003oooooool02?ooo`00:_ooo`<0000K oooo00<0003oooooool08oooo`030000oooooooo03Soool00`000?ooooooo`0Zoooo00<0003ooooo ool0:_ooo`030000oooooooo00ooool00`000?ooooooo`08oooo000Zoooo00<0003oool000006_oo o`030000oooooooo02Coool00`000?ooooooo`0hoooo00<0003oooooool0:_ooo`030000oooooooo 02Woool00`000?ooooooo`0@oooo00<0003oooooool02?ooo`00:_ooo`030000oooo000001[oool0 0`000?ooooooo`0Toooo00<0003oooooool0>?ooo`030000oooooooo02[oool00`000?ooooooo`0Y oooo00<0003oooooool04?ooo`030000oooooooo00Soool002[oool00`000?ooo`00000Joooo00<0 003oooooool05_ooo`030000oooooooo00;oool00`000?ooooooo`02oooo00<0003oooooool00_oo o`030000oooooooo02goool500001Oooo`030000oooooooo02Woool500000_ooo`030000oooooooo 00;oool00`000?ooooooo`0Noooo00<0003oooooool04?ooo`030000oooo000000800006oooo000Z oooo00<0003oool000006_ooo`030000oooooooo01Goool00`000?ooo`000007oooo00<0003oool0 00000oooo`030000oooooooo02ooool00`000?ooooooo`05oooo00<0003oooooool0:oooo`030000 oooooooo00Koool00`000?ooo`00000Ooooo00<0003oooooool04?ooo`030000oooo000000Soool0 02[oool00`000?ooo`00000Joooo00<0003oooooool05Oooo`030000oooo000000Woool01@000?oo ooooooooo`00000aoooo00<0003oooooool01?ooo`030000oooooooo02coool00`000?ooooooo`08 oooo00<0003oooooool07Oooo`030000oooooooo013oool010000?ooooooo`0000Ooool002[oool0 10000?ooooooo`0001Woool00`000?ooooooo`0Eoooo00<0003oool000001oooo`800004oooo00<0 003oooooool0;oooo`030000oooooooo00Coool00`000?ooooooo`0/oooo00<0003oooooool01_oo o`80000Poooo00<0003oooooool04?ooo`040000oooooooo00001oooo`00:_ooo`8000000oooo`00 0?ooo`0Goooo00<0003oooooool05_ooo`030000oooo000000Ooool00`000?ooooooo`04oooo00<0 003oooooool0;_ooo`030000oooooooo00Coool00`000?ooooooo`0/oooo00<0003oooooool01_oo o`030000oooooooo01ooool00`000?ooooooo`0@oooo00D0003oool0003oool000001_ooo`00:_oo o`040000oooooooo00006?ooo`030000oooooooo01Ooool00`000?ooooooo`06oooo0`0000Coool0 0`000?ooooooo`0]oooo0P0000Koool00`000?ooooooo`0[oooo0`0000Ooool300007_ooo`030000 oooooooo017oool010000?ooooooo`0000Ooool002[oool010000?ooooooo`0001Soool00`000?oo ooooo`0Woooo00<0003oooooool0=?ooo`030000oooooooo02koool00`000?ooooooo`0Voooo00<0 003oooooool04Oooo`030000oooooooo00Soool002[oool010000?ooooooo`0001Soool00`000?oo ooooo`0Woooo00<0003oooooool0=?ooo`030000oooooooo02koool00`000?ooooooo`0Voooo00<0 003oooooool04_ooo`030000oooooooo00Ooool002[oool010000?ooooooo`0001Soool00`000?oo ooooo`0Woooo00<0003oooooool0=?ooo`030000oooooooo02koool00`000?ooooooo`0Voooo00<0 003oooooool04_ooo`030000oooooooo00Ooool002[oool010000?ooooooo`0001Soool00`000?oo ooooo`0Xoooo00<0003oooooool0?oo o`030000oooooooo01koool00`000?ooooooo`0Foooo00<0003oooooool01_ooo`00:_ooo`050000 oooooooooooo000001?oool00`000?ooooooo`0aoooo00<0003oooooool0:?ooo`030000oooooooo 03Woool00`000?ooooooo`0Noooo00<0003oooooool05_ooo`030000oooooooo00Koool002[oool0 1@000?ooooooooooo`00000Coooo00<0003oooooool0Oooo`030000oooooooo01goool00`000?ooooooo`0Goooo00<0003ooooo ool01_ooo`00:_ooo`050000oooooooooooo000001?oool00`000?ooooooo`0boooo00<0003ooooo ool09_ooo`030000oooooooo03_oool00`000?ooooooo`0Loooo00<0003oooooool05oooo`030000 oooooooo00Koool002[oool01@000?ooooooooooo`00000Coooo00<0003oooooool0<_ooo`030000 oooooooo02Koool00`000?ooooooo`0koooo00<0003oooooool07?ooo`030000oooooooo01Ooool0 0`000?ooooooo`06oooo000Zoooo00D0003oooooooooool000004_ooo`030000oooooooo03Coool0 0`000?ooooooo`0Uoooo00<0003oooooool0>oooo`030000oooooooo01coool00`000?ooooooo`0G oooo00<0003oooooool01_ooo`00:_ooo`050000oooooooooooo000001;oool00`000?ooooooo`0d oooo00<0003oooooool09?ooo`030000oooooooo03coool00`000?ooooooo`0Loooo00<0003ooooo ool05oooo`030000oooooooo00Koool002[oool00`000?ooooooo`02oooo00<0003oooooool03ooo o`030000oooooooo03Coool00`000?ooooooo`0Toooo00<0003oooooool0??ooo`030000oooooooo 01coool00`000?ooooooo`0Goooo00<0003oooooool01_ooo`00:_ooo`800003oooo00<0003ooooo ool03oooo`030000oooooooo03Coool00`000?ooooooo`0Toooo00<0003oooooool0?Oooo`030000 oooooooo01_oool00`000?ooooooo`0Goooo00<0003oooooool01_ooo`00:_ooo`030000oooooooo 00;oool00`000?ooooooo`0?oooo00<0003oooooool0=Oooo`030000oooooooo02?oool00`000?oo ooooo`0moooo00<0003oooooool06oooo`030000oooooooo01Ooool00`000?ooooooo`06oooo000Z oooo00<0003oooooool00_ooo`030000oooooooo00ooool00`000?ooooooo`0eoooo00<0003ooooo ool08_ooo`030000oooooooo03koool00`000?ooooooo`0Joooo00<0003oooooool06?ooo`030000 oooooooo00Koool002[oool00`000?ooooooo`02oooo00<0003oooooool03oooo`030000oooooooo 03Goool00`000?ooooooo`0Roooo00<0003oooooool0?_ooo`030000oooooooo01[oool00`000?oo ooooo`0Hoooo00<0003oooooool01_ooo`001_ooo`D000000oooo`000?ooo`02oooo1P0000;oool0 0`000?ooooooo`0?oooo00<0003oooooool00_ooo`030000oooooooo00koool00`000?ooooooo`0f oooo00<0003oooooool08_ooo`030000oooooooo03ooool00`000?ooooooo`0Ioooo00<0003ooooo ool06?ooo`030000oooooooo00Koool000Soool00`000?ooooooo`02oooo00<0003oool000000ooo o`070000oooooooooooo0000oooo0000013oool00`000?ooooooo`02oooo00<0003oooooool03_oo o`030000oooooooo03Ooool00`000?ooooooo`0Poooo00<0003oooooool0@?ooo`030000oooooooo 01Woool00`000?ooooooo`0Hoooo00<0003oooooool01_ooo`002?ooo`030000oooooooo00?oool0 0`000?ooooooo`02oooo00L0003oooooooooool0003oool000004?ooo`030000oooooooo00;oool0 0`000?ooooooo`0>oooo00<0003oooooool0=oooo`030000oooooooo023oool00`000?ooooooo`10 oooo00<0003oooooool06Oooo`030000oooooooo01Soool00`000?ooooooo`06oooo0008oooo00<0 003oooooool00_ooo`030000oooo000000?oool01`000?ooooooooooo`000?ooo`00000@oooo0P00 00Coool00`000?ooooooo`0=oooo00<0003oooooool0=oooo`030000oooooooo023oool00`000?oo ooooo`10oooo00<0003oooooool06Oooo`030000oooooooo01Soool00`000?ooooooo`06oooo0008 oooo00D0003oooooooooool000000oooo`040000oooooooo00000oooo`030000oooo0000013oool0 0`000?ooooooo`03oooo00<0003oooooool03Oooo`030000oooooooo03Ooool00`000?ooooooo`0P oooo00<0003oooooool0@?ooo`030000oooooooo01Soool00`000?ooooooo`0Ioooo00<0003ooooo ool01_ooo`001oooo`800009oooo0P0000Coool00`000?ooooooo`07oooo00<0003oooooool01Ooo o`030000oooooooo00?oool00`000?ooooooo`0=oooo00<0003oooooool0=oooo`030000oooooooo 01ooool00`000?ooooooo`12oooo00<0003oooooool05oooo`030000oooooooo01Woool00`000?oo ooooo`06oooo000Roooo00<0003oooooool01Oooo`030000oooooooo00?oool00`000?ooooooo`0= oooo00<0003oooooool0>?ooo`030000oooooooo01koool00`000?ooooooo`12oooo00<0003ooooo ool05oooo`030000oooooooo01Woool00`000?ooooooo`06oooo000Loooo0`0000Coool00`000?oo ooooo`04oooo00<0003oooooool00oooo`030000oooooooo00coool00`000?ooooooo`0ioooo00<0 003oooooool07_ooo`030000oooooooo04;oool00`000?ooooooo`0Goooo00<0003oooooool06Ooo o`030000oooooooo00Koool002?oool00`000?ooooooo`04oooo00<0003oooooool00oooo`030000 oooooooo00coool00`000?ooooooo`0ioooo00<0003oooooool07Oooo`030000oooooooo04?oool0 0`000?ooooooo`0Goooo00<0003oooooool06_ooo`030000oooooooo00Goool0027oool300001_oo o`030000oooooooo00?oool00`000?ooooooo`0Oooo`030000oooooooo 01goool00`000?ooooooo`14oooo00<0003oooooool05Oooo`030000oooooooo01_oool00`000?oo ooooo`05oooo000Zoooo0P0000Coool00`000?ooooooo`0_ooo`030000 oooooooo01coool00`000?ooooooo`14oooo00<0003oooooool05Oooo`030000oooooooo01_oool0 0`000?ooooooo`05oooo000Zoooo00<0003oooooool00oooo`030000oooooooo00coool00`000?oo ooooo`0joooo00<0003oooooool07?ooo`030000oooooooo04Coool00`000?ooooooo`0Eoooo00<0 003oooooool06oooo`030000oooooooo00Goool002[oool00`000?ooooooo`03oooo00<0003ooooo ool03?ooo`030000oooooooo03[oool00`000?ooooooo`0Koooo00<0003oooooool0AOooo`030000 oooooooo01Goool00`000?ooooooo`0Koooo00<0003oooooool01Oooo`00:_ooo`030000oooooooo 00?oool00`000?ooooooo`0;oooo00<0003oooooool0>oooo`030000oooooooo01_oool00`000?oo ooooo`16oooo00<0003oooooool05?ooo`030000oooooooo01_oool00`000?ooooooo`05oooo000Z oooo00<0003oooooool01?ooo`030000oooooooo00[oool00`000?ooooooo`0koooo00<0003ooooo ool06oooo`030000oooooooo04Koool00`000?ooooooo`0Coooo00<0003oooooool07?ooo`030000 oooooooo00Goool002[oool00`000?ooooooo`04oooo00<0003oooooool02_ooo`030000oooooooo 03coool00`000?ooooooo`0Ioooo00<0003oooooool0Aoooo`030000oooooooo01?oool00`000?oo ooooo`0Loooo00<0003oooooool01Oooo`00:_ooo`030000oooooooo00Coool00`000?ooooooo`0: oooo00<0003oooooool0??ooo`030000oooooooo01Woool00`000?ooooooo`17oooo00<0003ooooo ool04oooo`030000oooooooo01coool00`000?ooooooo`05oooo000Zoooo0P0000Goool00`000?oo ooooo`0:oooo00<0003oooooool0??ooo`030000oooooooo01Woool00`000?ooooooo`18oooo00<0 003oooooool04_ooo`030000oooooooo01coool00`000?ooooooo`05oooo000Zoooo00<0003ooooo ool01?ooo`030000oooooooo00[oool00`000?ooooooo`0loooo00<0003oooooool06Oooo`030000 oooooooo04Soool00`000?ooooooo`0Aoooo00<0003oooooool07Oooo`030000oooooooo00Goool0 02[oool00`000?ooooooo`04oooo00<0003oooooool02Oooo`030000oooooooo03koool00`000?oo ooooo`0Goooo00<0003oooooool0BOooo`030000oooooooo017oool00`000?ooooooo`0Moooo00<0 003oooooool01Oooo`00:_ooo`030000oooooooo00Coool00`000?ooooooo`09oooo00<0003ooooo ool0?_ooo`030000oooooooo01Ooool00`000?ooooooo`19oooo00<0003oooooool04Oooo`030000 oooooooo01goool00`000?ooooooo`05oooo000Zoooo00<0003oooooool01?ooo`030000oooooooo 00Woool00`000?ooooooo`0noooo00<0003oooooool05oooo`030000oooooooo04[oool00`000?oo ooooo`0@oooo00<0003oooooool07Oooo`030000oooooooo00Goool002[oool00`000?ooooooo`04 oooo00<0003oooooool02Oooo`030000oooooooo03koool00`000?ooooooo`0Foooo00<0003ooooo ool0Boooo`030000oooooooo00ooool00`000?ooooooo`0Noooo00<0003oooooool01Oooo`00:_oo o`030000oooooooo00Coool00`000?ooooooo`08oooo00<0003oooooool0@?ooo`030000oooooooo 01Goool00`000?ooooooo`1;oooo00<0003oooooool03oooo`030000oooooooo01koool00`000?oo ooooo`05oooo000Zoooo0P0000Goool00`000?ooooooo`08oooo00<0003oooooool0@?ooo`030000 oooooooo01Goool00`000?ooooooo`1oooo00<0 003oooooool07_ooo`030000oooooooo00Goool002[oool00`000?ooooooo`05oooo00<0003ooooo ool01oooo`030000oooooooo047oool00`000?ooooooo`0Coooo00<0003oooooool0COooo`030000 oooooooo00goool00`000?ooooooo`0Ooooo00<0003oooooool01Oooo`00:_ooo`030000oooooooo 00Goool00`000?ooooooo`06oooo00<0003oooooool0@oooo`030000oooooooo01;oool00`000?oo ooooo`1>oooo00<0003oooooool03?ooo`030000oooooooo01ooool00`000?ooooooo`05oooo000Z oooo00<0003oooooool01Oooo`030000oooooooo00Koool00`000?ooooooo`13oooo00<0003ooooo ool04Oooo`030000oooooooo04ooool00`000?ooooooo`0"], ImageRangeCache->{{{82.875, 369.875}, {257.188, 80.25}} -> {-1.0466, -2.29412, 0.00837206, 2.79124}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "Roots", StyleBox["\n", "Input"] }], "Section"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?FindRoot\)\)], "Input"], Cell[BoxData[ \("FindRoot[lhs==rhs, {x, x0}] searches for a numerical solution to the \ equation lhs==rhs, starting with x=x0."\)], "Print"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(FindRoot[gam == Sinh[gam], {gam, 1}]\)], "Input"], Cell[BoxData[ \({gam \[Rule] 0.0178633039634514485`}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(FindRoot[gam == 0.5 Cosh[gam], {gam, 0}]\)], "Input"], Cell[BoxData[ \({gam \[Rule] 0.589387763462236335`}\)], "Output"] }, Open ]], Cell[BoxData[""], "Input"], Cell[BoxData[""], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["Integrating Numerical Data", "Section"], Cell[BoxData[ \(<< NumericalMath`ListIntegrate`\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?ListIntegrate\)\)], "Input"], Cell[BoxData[ \("ListIntegrate[{y0, y1, ..., yn}, h, k] uses an InterpolatingFunction \ object to give an approximation to the integral of a function with values \ equal to y0,...,yn at points equally spaced a distance h apart. If k is odd \ it is increased by 1. The default value of k is 4. ListIntegrate[{{x0,y0}, \ {x1,y1}, ..., {xn,yn}}, k] can be used for variable stepsize data. If the \ data are known to contain errors, you may be better off performing Integrate \ on the result of Fit applied to the data."\)], "Print"] }, Open ]], Cell[BoxData[ \(This\ integrates\ a\ collection\ of\ interpolating\ polynomials\ of\ degree\ \(k.\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(data\ = \ Table[n^2, \ {n, 0, 7}]\)], "Input"], Cell[BoxData[ \({0, 1, 4, 9, 16, 25, 36, 49}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ListIntegrate[data, 1]\)], "Input"], Cell[BoxData[ \(343\/3\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ListIntegrate[data, 1, 2]\)], "Input"], Cell[BoxData[ \(231\/2\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ListIntegrate[data, 1, 3]\)], "Input"], Cell[BoxData[ \(343\/3\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Integrate[x^2, \ {x, 0, \ 7}]\)], "Input"], Cell[BoxData[ \(343\/3\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Gaussian Quadrature", "Section"], Cell[BoxData[ \(<< NumericalMath`GaussianQuadrature`\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?GaussianQuadratureWeights\)\)], "Input"], Cell[BoxData[ \("GaussianQuadratureWeights[n, a, b, prec] gives a list of the pairs \ {abscissa, weight} to prec digits precision for the elementary n-point \ Gaussian quadrature formula for quadrature on the interval a to b. The \ argument prec is optional."\)], "Print"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(GaussianQuadratureWeights[2, \(-1\), 1]\)], "Input"], Cell[BoxData[ \({{\(-0.577350269189625731`\), 0.999999999999999822`}, { 0.577350269189625731`, 0.999999999999999822`}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(GaussianQuadratureError[2, f, \(-1\), 1]\)], "Input"], Cell[BoxData[ RowBox[{\(-0.00740740740740740921`\), " ", SuperscriptBox["f", TagBox[\((4)\), Derivative], MultilineFunction->None]}]], "Output"] }, Open ]] }, Open ]] }, Open ]] }, FrontEndVersion->"Microsoft Windows 3.0", ScreenRectangle->{{0, 640}, {0, 424}}, WindowSize->{461, 252}, WindowMargins->{{0, Automatic}, {Automatic, 5}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic} ] (*********************************************************************** Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. ***********************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1731, 51, 59, 1, 105, "Title"], Cell[1793, 54, 146, 7, 170, "Subtitle"], Cell[1942, 63, 169, 3, 70, "Input"], Cell[CellGroupData[{ Cell[2136, 70, 41, 3, 74, "Section"], Cell[2180, 75, 420, 8, 150, "Input"], Cell[CellGroupData[{ Cell[2625, 87, 71, 1, 30, "Input"], Cell[2699, 90, 59, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2795, 96, 60, 1, 30, "Input"], Cell[2858, 99, 126, 2, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3021, 106, 45, 1, 50, "Input"], Cell[3069, 109, 543, 7, 139, "Print"], Cell[3615, 118, 256, 8, 25, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[3908, 131, 85, 1, 30, "Input"], Cell[3996, 134, 78, 1, 42, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4111, 140, 59, 1, 30, "Input"], Cell[4173, 143, 397, 5, 101, "Print"], Cell[4573, 150, 296, 8, 25, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[4906, 163, 55, 1, 50, "Input"], Cell[4964, 166, 316, 4, 82, "Print"], Cell[5283, 172, 349, 9, 63, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[5681, 187, 128, 4, 95, "Section"], Cell[5812, 193, 65, 1, 30, "Input"], Cell[CellGroupData[{ Cell[5902, 198, 63, 1, 50, "Input"], Cell[5968, 201, 578, 8, 158, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[6583, 214, 126, 3, 30, "Input"], Cell[6712, 219, 386, 6, 68, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7135, 230, 64, 1, 30, "Input"], Cell[7202, 233, 22574, 605, 186, 5906, 395, "GraphicsData", "PostScript", "Graphics"], Cell[29779, 840, 130, 3, 29, "Output"] }, Open ]], Cell[29924, 846, 118, 3, 30, "Input"], Cell[CellGroupData[{ Cell[30067, 853, 90, 1, 30, "Input"], Cell[30160, 856, 421, 7, 68, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[30618, 868, 62, 1, 30, "Input"], Cell[30683, 871, 22195, 607, 186, 6119, 404, "GraphicsData", "PostScript", "Graphics"], Cell[52881, 1480, 130, 3, 29, "Output"] }, Open ]], Cell[53026, 1486, 125, 2, 30, "Input"], Cell[CellGroupData[{ Cell[53176, 1492, 58, 1, 30, "Input"], Cell[53237, 1495, 794, 11, 215, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[54068, 1511, 93, 1, 30, "Input"], Cell[54164, 1514, 812, 14, 144, "Output"] }, Open ]], Cell[54991, 1531, 80, 1, 30, "Input"], Cell[CellGroupData[{ Cell[55096, 1536, 64, 1, 30, "Input"], Cell[55163, 1539, 421, 7, 68, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[55621, 1551, 67, 1, 30, "Input"], Cell[55691, 1554, 34300, 767, 186, 6446, 419, "GraphicsData", "PostScript", "Graphics"], Cell[89994, 2323, 130, 3, 29, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[90173, 2332, 67, 3, 73, "Section"], Cell[CellGroupData[{ Cell[90265, 2339, 46, 1, 30, "Input"], Cell[90314, 2342, 144, 2, 44, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[90495, 2349, 69, 1, 30, "Input"], Cell[90567, 2352, 70, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[90674, 2358, 74, 1, 30, "Input"], Cell[90751, 2361, 69, 1, 29, "Output"] }, Open ]], Cell[90835, 2365, 26, 0, 30, "Input"], Cell[90864, 2367, 26, 0, 30, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[90927, 2372, 45, 0, 53, "Section"], Cell[90975, 2374, 64, 1, 30, "Input"], Cell[CellGroupData[{ Cell[91064, 2379, 51, 1, 30, "Input"], Cell[91118, 2382, 537, 7, 139, "Print"] }, Open ]], Cell[91670, 2392, 123, 2, 30, "Input"], Cell[CellGroupData[{ Cell[91818, 2398, 67, 1, 30, "Input"], Cell[91888, 2401, 62, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[91987, 2407, 55, 1, 30, "Input"], Cell[92045, 2410, 40, 1, 42, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[92122, 2416, 58, 1, 30, "Input"], Cell[92183, 2419, 40, 1, 42, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[92260, 2425, 58, 1, 30, "Input"], Cell[92321, 2428, 40, 1, 42, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[92398, 2434, 62, 1, 30, "Input"], Cell[92463, 2437, 40, 1, 42, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[92552, 2444, 38, 0, 53, "Section"], Cell[92593, 2446, 69, 1, 30, "Input"], Cell[CellGroupData[{ Cell[92687, 2451, 63, 1, 30, "Input"], Cell[92753, 2454, 277, 4, 82, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[93067, 2463, 72, 1, 30, "Input"], Cell[93142, 2466, 144, 2, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[93323, 2473, 73, 1, 30, "Input"], Cell[93399, 2476, 182, 5, 29, "Output"] }, Open ]] }, Open ]] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)