(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     23742,        520]
NotebookOptionsPosition[     22316,        469]
NotebookOutlinePosition[     22911,        491]
CellTagsIndexPosition[     22868,        488]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.498398015903*^9, 3.498398020979*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"$Assumptions", "=", 
  RowBox[{"{", 
   RowBox[{"n", "\[Element]", " ", "Integers"}], "}"}]}]], "Input",
 CellChangeTimes->{{3.498398027125*^9, 3.4983980388269997`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"n", "\[Element]", "Integers"}], "}"}]], "Output",
 CellChangeTimes->{3.4983980392980003`*^9, 3.5300137449131403`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"Clear", "[", "\[Lambda]", "]"}]], "Input",
 CellChangeTimes->{{3.498398087329*^9, 3.498398091085*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"k", "[", "n_", "]"}], ":=", 
  RowBox[{"n", " ", 
   FractionBox[
    RowBox[{"2", "\[Pi]"}], "\[Lambda]"]}]}]], "Input",
 CellChangeTimes->{{3.4983980410629997`*^9, 3.4983980510620003`*^9}, 
   3.530013748119461*^9, 3.530013780574706*^9, 3.530013815735222*^9}],

Cell[BoxData[
 RowBox[{"Clear", "[", "B", "]"}]], "Input",
 CellChangeTimes->{{3.498398128776*^9, 3.49839813073*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"y", "[", "x_", "]"}], ":=", 
  RowBox[{"x", " ", 
   FractionBox["B", "\[Lambda]"]}]}]], "Input",
 CellChangeTimes->{{3.498398105564*^9, 3.498398120698*^9}, 
   3.5300137513797865`*^9, 3.5300137832399726`*^9, 3.530013818367485*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"a0", "=", 
  RowBox[{
   FractionBox["2", "\[Lambda]"], 
   RowBox[{
    SubsuperscriptBox["\[Integral]", "0", "\[Lambda]"], 
    RowBox[{"1", 
     RowBox[{"y", "[", "x", "]"}], 
     RowBox[{"\[DifferentialD]", "x"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.4983984414709997`*^9, 3.49839851867*^9}}],

Cell[BoxData["B"], "Output",
 CellChangeTimes->{{3.498398504078*^9, 3.498398519307*^9}, 
   3.530013836541302*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Simplify", "[", 
  RowBox[{
   RowBox[{"a", "[", "n_", "]"}], "=", 
   RowBox[{
    FractionBox["2", "\[Lambda]"], 
    RowBox[{
     SubsuperscriptBox["\[Integral]", "0", "\[Lambda]"], 
     RowBox[{
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"k", "[", "n", "]"}], " ", "x"}], "]"}], 
      RowBox[{"y", "[", "x", "]"}], 
      RowBox[{"\[DifferentialD]", "x"}]}]}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.498398522092*^9, 3.4983985359*^9}, {
  3.5300138054391923`*^9, 3.5300138096066093`*^9}}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{
  3.498398539738*^9, {3.5300137690825567`*^9, 3.5300138281014585`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Simplify", "[", 
  RowBox[{
   RowBox[{"b", "[", "n_", "]"}], "=", 
   RowBox[{
    FractionBox["2", "\[Lambda]"], 
    RowBox[{
     SubsuperscriptBox["\[Integral]", "0", "\[Lambda]"], 
     RowBox[{
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"k", "[", "n", "]"}], " ", "x"}], "]"}], 
      RowBox[{"y", "[", "x", "]"}], 
      RowBox[{"\[DifferentialD]", "x"}]}]}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.498398543981*^9, 3.498398557657*^9}, {
  3.530013849210569*^9, 3.53001385362201*^9}}],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox["B", 
   RowBox[{"n", " ", "\[Pi]"}]]}]], "Output",
 CellChangeTimes->{
  3.4983985594890003`*^9, {3.530013846153263*^9, 3.530013854678116*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"y", "[", "x_", "]"}], ":=", 
   RowBox[{
    FractionBox["a0", "2"], "+", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"n", "=", "1"}], "300"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"b", "[", "n", "]"}], " ", 
       RowBox[{"Sin", "[", 
        RowBox[{
         RowBox[{"k", "[", "n", "]"}], " ", "x"}], "]"}]}], ")"}]}]}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.498398569955*^9, 3.498398625111*^9}, {
   3.4983987457980003`*^9, 3.498398749387*^9}, 3.5300138626939173`*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"B", "=", "1"}], "\[IndentingNewLine]", 
 RowBox[{"\[Lambda]", "=", "1"}]}], "Input",
 CellChangeTimes->{{3.498398652013*^9, 3.4983986572200003`*^9}}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{3.498398657758*^9, 3.5300138676964173`*^9}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{3.498398657758*^9, 3.5300138677074184`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"y", "[", "x", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "0", ",", 
     RowBox[{"2", "\[Lambda]"}]}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4983986275810003`*^9, 3.4983986499110003`*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwVl3k4V10Xhk2VBimR9CpkCEmSkqhHiSQVQoZeKUKmUMpQhqJMUUhE5nlO
ZEgoXpWUQimZkmTIOT8hleF8+/vrXOusvdd61jpr7+vcYqfPGZzhYGNjK2Rn
Y/v/U+fM0Lu64TN7kqW1teOYXnBprHXpFdoPpWUKWtr7G7AtKTbhoZAhLL8V
l5z/nl07ICzhf0nKAkOOIi9c/imrPfS6o+S+kDX8Rp5bFUSH10qErctgRM6C
zct3m6liFmZ1rO6eknKG9MiF8DV3MjAS0N5wKsYVbEHeLc73jtQet55+Eil0
AX2SPgWXHRJqGzTWltcnXEQd505Oj1chtQriu4snRTyR7JgkMF4UhQQOyxzJ
NG+o77cIE+JPBnf/1VRjKR/UcdWsU9K8iaq+1tf1hn7w41Z+ai3kV+uiZ61z
IsYfbA5LXQT8rtZK1U01TnZchWiTmV6KcgC6t9zQuCkUAPWxg/bCFvcQnbSm
TtI8EHWy/ZvtD8dAhzdXrSbhOtR3WBgFylwEm69qpXHPDRKPyeFPs6l9RDVv
p0WCoR6Z6OAOLzhaWJTcOBWCt5rqAvzzd7DhDS0vmhYKhVv7LKJm4/Bxt39e
xUAY2NiKZB6fjUB4AZ+0vlQ42Ow7WL+un4fmuvT0YdsIqJv7iBV88YJP470V
Tw1vQT3bunfmfSRMfb4euTl1C7fCyoZ3bY+H0g65m6Yxt6F3/ZHva9m74KUu
vJJUjoR6XbTUh1/BGMl4svhnB7H9fsbPbLqG//5dqF3jEQXLnjs5EYq3kSxw
9HqIUDT6hCnaiCse3q/vNhhXRUNPIXSvkFoCjAP7OMTN74DVZV39oCAaW3fL
7KVn7qDP5cABn6dhWDbl6vs4IQbJX0/f9DC8je/5VU9u7L4L9d2jxnP69/DM
mnP2WM9dsNxSM25qJeG+sO4uUd9Y1M1lVc6ZxcOjPdrjh0gc1HVOJnbpR+NY
WPejiro46LkHVC44EAX5/VJTAafuwa9bnT1/IA6LZ5236XPEg9WTX1Npk4yB
h+Wu69KILeoe8Uk/GbUObMXDGgnw+yN+vj3jHu6JH6TKBhKgVyW2h70jBhc/
35a7GngfoqWmm0IN7kE/qtP+iFQiknva1vhXJEPukHjO2ueJeHvM7dtm5VRw
czp+H7RNgsIH7pfPjifD7d/OSFOuZLCtFPvhGB+Pk5XzD2oNk1Gn4sJ9cyIe
ugLi7yQzkiFq/3XPoR/J2OV6gBU6lQw7XbXd6m5pkH7twPtTMwXFlzm5PH+l
QkDmlrxJTArqHk2qtAQkgyOw9HDNYApYU4aqzZ8Twer76ChBdKwJGN1brpuC
HrW50JAbqRBdOGZ2UDEdr2LF8lgdqRiSXXOWbWkGKiY1m4yl02AZs/Vm4IE0
ZOrZD1d7pGHFIPteE5cUROeHc4u/TEOFZ7Oi6YpU+HM/3BgslA7W1n1+fWXp
kD1wO6jEJR0v+hbNyMdnojXw3PDn5+mwzHo2Fnc1A94Nh3UWiGSg2LOnfBOJ
K8Eplyd/MQOW5wSLnGvT0Lx3yTKT1xnI1u3YZWOeCXe/IUd/iUxI1/8sCJLI
xvraxte53plYsa2b21E2C41z6fLtrZlYw5Z1P083A+fUrkXMyWRB+nkTxReZ
gTXep1hS/lkQdR8/N8qRjbpK6Ot9zILCaqnmnw9zcPb3uhLPLeR9slPhfFs2
+JRn+dKuZ0NdX7IrVzgLVe6d55u7s7Ezbg11wTILVqUV7VNKOfCLqP28ND8H
yyZitouE5eBt8byMsnYeyra6x2h/Jf40y6hfPrmwcDk27borF9pFN8S5h7Ox
sGirSfztXHz8ReXZbspB0RhvZcMQsc2GX3ob58FEjhKikAeP9+ornPryyT3R
7CV4Nw/FolcSeDflIycn97M6lQfR9g37HUtyYTAUpGavmQ/p/l2coRO5mJGy
vR+VkI+gylPs7BwFSD+jOV89kY9YJRfvfu9C6KaLnxzUKYBe+KpJgY8FmOpn
r+NNLYAoV1SIqWc+EsX6RFX+EPshu5ZZZT60LGv8T+sVwkR/N5tzYSHoxIT+
0KxCRDetvLSQpxix3V4aZfOF6PMTdv8SWAQT05Qi1ngpTI58bDx9uRRLGf8T
B9bXQu3jxrsXn9VAUp8aqg9tAN9eh/SDEQ249EYFsZMvodWR63uEtwmCpSt7
2NlbUBXN87BkYQse9Ws0t6i1YselmtW7t7fimf39f0x/tGN6YP/pH+Q5JbWb
t8CpA00hm0XlL3WA40/ein8kO+FZya0eS547GMtHxwq7oH5JMWgwpwvNAUpv
GmV7gf56R2mpXojcGv2TQ/WhIFG/l2uiD6PTvT3vFfpR1MOj+0qxH+kjj5uW
h39FR+SP5qmQr1i+LoseGBiAwCTnIYnBAYSdLbNJlhhEldwnHR+ZQfxSzxuw
8/2OI1ST7M7L33FBsaF0df0Qlmhv1r5dPYTewobQd2IjkLAMPen+zwhm7Nvj
ZbxGsaB3xEcoYBSb9i6at3n8A5bLHgkV1P3AgVQx2weCFAy+WldNcVJQil0u
1ShMQXjE6dpNUQqi4TMDnWIUbs7G6v3ZSYEnYCiNlqJQrH93cu92Cn8935/m
kqPQaZthkytC4fu5Z2JCWyky75KCdWsotJ8p6tu8g4LCjrbfGVvIe/OEpH2q
FN7H7i/hOEqhQD/Y4rg6hTvh6ioRxhTuHbi4zlGTgrJpv2LZPgrXd1t1+elQ
qPbv1kwjcc9v04u/Q/YtkjrPvNxDwVJmt1muIYUSv7jpAUsKh0VkhWpNKYSt
vb4x8gKFXQKCH9ssKEjH5uk62VKQWsp1d8iK1MXXFbKNxFnFPm40Z0fsa0vr
uonde7Ju/8EqEq/4VKVHHIVmowJOi1pS31yTtc4TChWH7j1za6CwrNRKMn+A
QsbeG/43XlIoxfh6qW8UIpUvqCe8ocAVmc0pU0/Bd/MppriNwmXViPcChRQc
xY/U/PeR+B8KXe9qpGAipHqls5vCq6UyzrG/KGjySqvR/eR7vHPZt3IxDcUF
AjOcQxRMlZ/aRv2gsH6GvWrNGAVDZpeiHMm7dJzy2PyTQm73qOHMWwq/Bz8r
75umMCpQd2EPB41vXS9+Gc9S0PF4sTpkI43W1rIyB3YaL0pqXTTW0Kh5kXrB
byGNv6ufp+aPkpuhJmLbnaU05KZsP/n2UIgtvfwzZwUNL9cO3zZuGgG5Zx/U
CNCIiZoeCdtDwzXZ2KVtLY1J+UfFLho0LGI0tgyJ0OAvXtwXRtYdClOgZiVo
VGy3/ryX6Nx5dV3BSlkaPSOW7iHLaEh4LHGU2kKj8eU1nnZdGiudp2VVlWgo
GXEp9drTmLcaGD6qQsO972p4gSaNEdN32dYk7yq6pe7pchodR2tsPUnebLmN
lY5EZ4NmnlS4No1ilxwrHgMaD1Rjv6UepmHMtLXJhdJI3BqYXk7ey1oY3fni
SSNko5tV83FSx7E7ywa307i07uSGLydo/F7suE6M9Oe0noJx1ykaGpIqLiBx
jlxjD+mwIfVvuXYqKZHGrketT1odaPTNe3sveEBDajht/LULjdmO1qT/ztHg
E3aXfOlOQ2VOvOAJqXvuiJZpgxcNmXVlxpvUaAz5C96s9aUxWjdfdCaeRnvp
UF1VANFX/qvlWy+Nuu+Vk2XBNO5uej59rpJG3tpQ6QfhNARS9Z21DImOwydO
5EfRSBv2UEnbTOOq3+ZbWbE0HFvi+g/503B+OF+fep/oPp5XPvyHhulgy/T9
VBqnAv0vf1vAwn6hlE1xWWS/ncyVRJJfQdftZHQ+DepYvUKUHA1hX42oCFKX
+pJPbzbq0VhUwv885BGpy1ZpYv0EjZ8D3/4GPiZ9SgpzGNzCQo9gubx/HZkn
rWuULInfpBN0+vJ/NDwK/GuuaNEou2Iac6mJxqMuevFHMp/JxbJNbi00uISM
BhNI3LCvM3NO7TQuTPbE6Ziy4LH69dazn2iwHvIWXrZkwepg4hnrHhpvVNJb
DEjdRy6fizv5lUaThff65WfJeS9Sf202RCO2c7jZiybnvX8lu/EYjeWdizdL
x7OwUuCrkv5P4r/x2LzQdhyzB0rtdKdpCPp+Z1sux8KQV2DCgVkyL0plAm0W
Qxgx3mSXwsbCtOaTnO/rOmHJWsUhwslCdIjS8lsp5fgQPBufQOrUO7pB1U+i
pPaQ+Lfta7lZ+PiOo5TDvaS2rvp1y90lLLBVRPEb/LLADuNHZwV4WPh97KRK
al4O8ulEziheknf73e9vR1KwIfjG/RV8ZL2agUrltH9t7AYX5XB+FuoaTM1q
ShNql1ebvFsqSPwqq9NOiPnWBhjtdQgWIv6EigTX73fxh5JZsEiYhRVGK7ns
VibBOYgvKWA9C37ulnw+6YEYEJvZySHGgnpOh0texrVas8dfW33FiX/Hja11
//jUvjVsdpyTZEFUZUWF3LVgaFGlC72lWej75WSkSf5rq2/cT/4tywIrGl/U
2qKhKHZ918XNZD/U0szNXJBV5dw+QeaAbZ1+89mLZ2qFDY87uyoS/yv55IcS
vogcAzetxEJx5w5L9cYYcN+QTnVUZkFhqsScvyQOPqIr1UZUSH779V4vucMx
Wfnnva0a0Tu75Yh2shvsj/Wf+7aH1KtgTL8VvIy+H02LrfYSPQJK1S+komB0
/WFanwaJtyq6M2tPPJpEEnZbaLFwy1XKnFpD/vsrAzo+a7NgeaDwGDMTjEcG
Tq6mh1hI3qdvO7r+GuR+GC3tOEz6p3rjGRtzCxWGurUhekT/8/i6D6fiofFk
3/k9x0g/3rk/HrgTjzeSKht/GpH11StkG/siYRa+5XOGCdEnvzvE/mQYBn9J
Rpias+BSzivImEXD7aSwBo8F8f+8xTj53Mf8c77pOsv/z8+VE827ExGssDjv
ghXRH9Ew3/YrBvxxjIW0DdHXnLhRJDkSyey/+LrsSLw6uwi6/C7k7H80RjiQ
fr4wNUv/noSK1n4vDWcyX++fu7MFJWO/6if5aRfSz921fVwD99CS1tKfe57E
cxa8HPQ+BmbLGmMsLpL+xWoYaK2Ox+CFah0+T9IfSma5wL4UuHWXzP/nzcKL
9QaNa5kUzGnmlHj6kPWZNkqDpYkIKkyy2exP+qteqR/Qfw/8gjFrv1wj8f/l
bOO/SLjFN+xN9HViF/JHWqemQm7o6lXtYDIfnCpSlC/5z9fz3DEbSuKt1duX
vjMFGpXnRorCSX2vgk5HtybijZhNotVtsn52gWrmaDLMQk4YCEYTPVoZr4d4
0zH402DhqxhSr/acWsOiDLiZH6zyiWNhp+EWWoI3DfP1cFZMYMFOYGFVgnEK
guV2bBhMZCHI9k39SEIq+O/IfYhLYeHtVLqo7Y0MJM9tCDmcTvTadPK+jc6E
nI3QHvYs8r0/9w8/7E1H+Rven6U5JN91VflG1TTsU16YaZdP6ldzzX11JB1v
kmZNhYvIfteMO2kLs2DGPcHz9gE53+XOPO2S2Rh0GX56rZTca87HT4TlZcLt
U6+7cjmZrxcSnbfWZmBu7weZ0UoyHyfLz+1mz0RQbnN3YjULsVsmVi4Mywb/
qvrbBrXE/keigqrMQbJ3pebCZ0Tfq8hfkeezITdQ9KeygeSLEnx2YDITFbqZ
BU7PyflxPNy4vTQLGmUJp8SayHx9rWiwWp2LN+uiBN43s2ByIeyrt14ezK4H
vwxqIf3j0S933ZmLQcr3ilorqdfELLP3ZTZcj1/cymon+bw6QvNsczBb6/gt
rYMF6Sx/ju70PARJW8Ud7yT+lXV+7MP54L9tenhpN/l+3TY5OX/zkPTnKHtt
L5n36RPiKxNyIXtaq8ytn6zPXZIW808eHjWpnZX6Rvy7jKSf7yjAvm3b1nV+
J9/HSXK7+LVCvI6XeXdzhOjptIhOrCqACZdo4N4xsn/cq++sQz4GHFerTNHE
No96Trfkw+X9srHsnyyssWjtvP+qELO7OVNOTJH+GosvPcxbjKDMP4YrfpPz
s0fxe6JHEVbxsrgb/pJ82QbyJjsLwVUbXySeMI6AK+a0dXMlHJQef1p86Cda
m8bfPLj8DHkHTyfX/PoJnxkt6/8ynmNqdO2KwZgJnB6LrFRObkbZsQOc0fsn
EV0/KvJuzTv4agcmlQ5OIn+plU1vehvEunU+pgVOwYO9V1mK/wN2KQdECyn+
gpNAyLKTzz7ixPqWP+xtv6CWb8DprPcZU003+138p+G8sGpMc3kPYrksLhhI
/IaKsMk5S90+rOopS+Vo+A0+nj1qCwK/YMEVf0MV1z/waxfOXNXXD7Yci6J+
vr/waqlK0t44ALd7RTf5nv3FAzzcNRn0DU82BXGvdJ6BzeJhdqdPg1ip4yZ3
dfksTCrfBqkLDqFr89s2z+pZZDSUS6zyHkZxIiQGLOaw+Runj3rDCP4t2rFy
3fwcrK1MTj5kRuEXmOLflTuPI/aiOvKaY+C5tEjnVME8BJc+vnxLfwz37Jz4
Bovm0b/+U9S2s2OQMmvrtC+Zh3C5deMj9zE8PKSSRpfOw+NM/vE+pzGo705y
uFA+D4NjvU1zlmNoll+g9KdyHhKXeBLLrMdgKuowe6V6Hr/nPF598RrD4Mp3
DRy189CzrTvMFTEGN07lmzeezmPehTXeGD4GZjLBaFnDPDKti1P3+I0hbJBj
/e3GeYgyomvHyP41H+0GBV7OQyksJv966BgyXr4pvPdqHqkOrYxF+hgUHytd
EnlD9MZs6LxdMIaa/HtIfzsPP8phdWvKGA4lsnHLtM1j39HQBok7Y1hdfjLc
fWYetEPud9PBMaRl/2f8d47oc/i3XXBoDAr3Non4sjGIsTmZ9J2dcE/o7e+c
nAxO2Ece4/qHwsEr00VBCxgsED8xuoDYH5z/9eDhZuAZfVMzeBGF05b16pFL
GDQn83ZNz42B1pdZLMjD4N2lvgVaKwiHaES8i+dloGej1/NNnsLi7VNxonwM
uC7kl87tohAjZX46g5+Bx6n9fnc2Udiw5qmsrCAD+VN5VROE54oWb5woFGIg
oXE/r4XwntpM2ONtwgzUPX/w+hEee/Hj57WK9QyuFMzPGBOeM+ox0d0txsAn
4qf/YT0K/S01/M/EyX6eB1MKKhTOPZXo1pJisFbl+mSPIoXZkpCMV9IM7NyS
dW8eohCUznLS28TgVlf4wFJHCgIxxjveb2ag6ypQJ3yeQuqN6nlTBQaGej0c
+f9SkPfc8LxHkcGv1TVeB8j+x/ZBEVbbGSy3nU0LNaGgfYI6PqTMoGHC55Lq
ZcKnhw1FnXYxqOVw7ou4TeEUqobG1Yieheac+gEUKAXRBxfB4Hrdnep+8n/p
teG658xeBnySw6e1CTcu4v+x128/gxUxq1r2BVGIXmCwZMEBBo/8dHh7cwhX
Tpe3Bh9kYNPKfvNoLuHdoXXxy3UZfOi8lqgZTuq7H2jZc4RB6VJtgydXKFgb
jEkW6jNwqRrN/xpPQX2R0egVQwbCoctEhv4jXF5dXXz4OIP0oWf39xN+/O0i
cXGdGYOPS2KbLdMotEmGqY6dYGBaZ2S/OpZCYecE25OTDFYnyUGnhkJIhHlj
2GmiZ8nzabvvFGz214eeOEP0iv7t+kp4b98fWX05Owbv7flWD5UT3iyMXD1r
z+D+yzuyYQ8p/Dn99/MrJwaXe0fNJ76Q/gmeTol3YXCRvSrfh/BhcfNLG4fz
DGbk3pwdnyIc7r9VTvUiqbfq6XPWKwp2O+LGl3iSfvHovjlBeHX/KFt5pzcD
53aD/HHCoSLJdpdzfRi0/7NRf0KacITh271e/gyGnq3pUCJc+mHxzkU6AQxk
E5SfhPcSTq9Jaha6waBlU/XpNopC+PlFkcPBpN6trpZWhIfspc8drwxjkG+5
czqIcKFm9wfh4AhSb/bRMcv1NEQj9/SbRDJIcHrM8Y7w6axWZpb0HWIbhJkZ
8NL4OMPj9PsuOX8TnA+uGNEoLXZXfHGPnJeMid4JwoG3znRP373PIFuxvp6N
cKzDWs0ntskk/sGj9tULaGi15F9VTmNgsSk40k6RxoYAfu1FmQz2xRrcSSXc
N7/zMk9HNgMBe/va19dofBr72pqZx6BvRfBXfcKHpamHYi8WMqgz+rR4H+HR
W8cf/qv1gNT37G77nAnhvGX/iK8uZeC60Uz3JuE27adXh749Iue74vYBD8Kn
EhdHCsoqGYQcGbUe30+DkTU4H1jNoGAurkWO8Ghnb+VOo1oGFbJeL//PnY+i
xeYlnjFY//OVRUU3jciDwfWTDQwUCxeceVdGeHKeFdTwnMGc56Lk5YTndB6a
HIluYqDB19b3ZxcNSbu6VdavyfzlJwyLEl5lWyf9adtbBjqj/BY7lrLQ9S4i
kbONwWSIoLnucxoV16et2t4zWPJtl/EywvnRqidl0j4yqF+Zaul9kIYLq5Fy
+8yghim9osDQOJQhX7qvh8E52kQjiXCGlFmMJ98Xcp7P634xfEWDg3d+T/9X
Br/HFQV1ltDoqT/DVTJI+s2p/kbwGI0qj9cv/YcZmPlEetuqshCzeXuE/g8G
B/+oBm8k/9Gu/QmGYjS5//JVtpsRbta9y7V2fJyBW5eYzCNyv0jrOvbWTZL7
sL9r/y1zwq3s7em3phmUCUS+uLhwHH1lqvaWf8l9EeulWSQ+jmr7tC0Kcww4
csovS3H8wF2RpVMMw0B57LP5NrY+/A+6YJn8
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  PlotRange->{{0, 2}, {-0.08244995961961232, 1.0878127654363876`}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{3.498398661237*^9, 3.498398804244*^9, 
  3.5300138822958775`*^9},
 ImageCache->GraphicsData["CompressedBitmap", "\<\
eJzFXAtcncWV/xIukGhaY0CJ5iaCFSMqGhJQSMAHjYmiRoOIgIBGVDSYkIQg
kphMaq1ufXWtJRoVamvV+n6sYouKr2qrRo3bFpQoxPej9iVuUpPdpXPOnLnf
/8797tef++vu9td4hzNzHv8z55xvZr7HyUtWX3DesiWrL2xakr1w1ZIVF1zY
1Ja9oHWVJqWM87zxUc8bt3Cyp9uXeinbvP3G1Jj+/ww15nl76/+keJ7Sf3uG
nCrkiCeESPi4r09uXRdIzh5t6wkg52S0dQeTe7xd1K7us3rSXLad1G7riY6a
/nRXm+nvzslM7Af+7uwAftbfptuT1M752zzP+2uLHsFtai2OeuM1MWX2CBu3
XP+bTET6H9No4JbpGdJi9uUszR0UIG0ZmVC7gccRfZI36Ux/MEiMuIDI4NTa
DTPmqrGlM8dt9CaypCsXnqV/W6mtts9KX3iWR4wTjOP80ThgN2Aku9OWzjTG
rd88eyTV0WLN0CG1UjfTNR4yeExkTTByCzYGGDGRZyd+/HLXDuQVm9IC9KS4
3lhFg9Zv9r02e8QaM0vTgozJHnXGBxgDvDvnf36mGOPwOcb0GGOWzvQHzZhr
jCEXJvFMhjM+wRjhXcH0xdOpnRqgJ9gzzqB0iToSviLBGB0rIcZMBN5Vpk3h
wrxhxnCWBk1ThJGYQNdzzHGVEDdr7ks6VWkOPyVGyvrN7K3IP54u8ZARFguw
VCOEka0Kip22rrjxaFA68K5grwyw+jHjHWOUoy/eqIrahPROg4AkD0U0Ko8m
weS3DqCKGn/4chyRDqwrOKE4fDTFG7cxLUmCG0Pyqowhfh1idKkaiSd/m2Be
HI2rNSUDPgcOmOAwk52pbKf5O81RhtZkj+a3+2W+oCuxjOdkmP7WdTw2aX9z
aU5GflVcf4T6o6PVfV/RH82l0dHcjOpNeCGKkAU5mSW9MiInMyejbAeOSCUZ
ekzGok7vb/R3cyn9FR2tG/LsJW+Svvrpq7HnTSHIxJVmNBNnczRT+VfbvZO0
s/5JY0S5UZ2T0Vz6z1L+dcdPjJvknExvvJ6olPyqth5Nm/J/7JF0PR2ZFCrZ
owVdHf1exv+3U9JtVNFEVRR7l8Cgtabf9maPVtQk9CcITDWp4Hu6pZ4HtDuM
FOM8igfp6CzVxIth0JVW4kQ/hOk3N6N1nfcl0SmPfOPzCs2qcFzQfKdSKyej
oH3tsEdrtDSTYWBBUe9aqqj/bbl3N46pqGkuXdSZV0g6ynZo7i90324m9RYU
t9RX1ORVkXbOVUj2vMKW+vEmS/PbSZUuul9Y2KSLO9u6C7rIispNvmKuJTrH
+cK1dqSoi1Q37vL+qri2UBUxfR39JQOkWvf9VxxmLY2qRqQ8ixS3NPzFcOZk
6krBHWU7SGnruv+M0ylsWqdmzMido69dfzac2aN1W6VvuDwrOppXpft2AXMu
+Z6qb2RBcU5mfvvakT+JzgytM91w6i49ATrod6rYKs90jVTU5OqSq737R8vW
Um+qObFRwdd9X4HCfAonXrhU1FDVtIwlg5ZO2sqz4pj0ZUPPKzMVM9zPDVNR
r6WTrgXFXF5jTNQ5yWgq2qit/4NossIiZGJuYUXtDlfVN6jduKuoV3vrM1A1
yXB19FfUVG7aDlzc+03F15iSwdYG71NAtQe1y7Maye9fuqq4s627ZKC51PsE
VDF97UjZQFM0nok6J1NbR9Fg0zTvY0A12TCV9OpEG3U17cmdw2U7NNNHoGlP
I6ygS4P9wsU0BZg+BExMb+uhfOD4jtOUYZgY0wegiekt9cz0FxdTpsWkvfc+
tZumle2wdMH3Z1fTXmD7e9Su21qeZel1QwuK9e+fhClFWLWzB/Xwd2EI0ys3
LerUv390dWSZaS3q1SG0zTCxjiwTd3VD+vdzF81UhtpQMqiZRkAT08unciz8
wdW0j4FaNqhzZhjQTIVJQCZ26r7WV9qp74CmfQxTUa/252eupmkmwFn624Bp
X5mEAZ3Fn7qYoiB9K7SngVc/cTVNp/aizkpy7RBgihomTrGPXUwzwLVvgabp
4NWPXE37UbtsB7v2TcA0A7z6oYsp2zDxNWwQNO0HXv3A1ZRt0qyoS6fZIGDK
BgHvu5hyhKlXF9wBGMj06r6KGv37nqtp/5iXerzfA6b9wavvqlhR3t+Gqrb6
d6DjWzJ8k/7d5uo4QPlZ+VtAcwDoGHHR5IpTKUP/HTTlgqZhV9OBgOANaB8I
mt5xZ2gmSN8C7Zmg6W1X00FqbPyizi2A5iDQsdVFk0ftilo6ZvJeAx15oGPI
1XEwB+lUDtJXAc3BoOktF80hJt6aKGc2g6ZDQNObribulPzfDJgOZbOL2WyK
Q95GBJ1V1W6I0NkSnyGN8Z7SImcRnPAN3itgT77yr1kDLDoNz4loSRF3HhQB
c/ONzCK6Cr0MjjnMyKwj9/+eZbrHN7yZheMd67XDTBqUDOhi8RJYebjys/d3
RqKzqTfb49gZjV7j5htZvIj/DTiTZekL7vbmqJUVdDoCRyzWh7M46xpKBrTE
X4N1BdQu6m1pkBxJlMiHP3BOYj1YYFKSJb4IHixQtqavsxKDzrbgkMN6cLbN
V10EXwAbZyupSlriGzIn8ScQ5vwm8aTC2jpH2Uo57P0K/DlHoqBLl64tiRGU
qpwDB/RoIch8HqxlOq06tMzXrUw/4HmhC6cGqWBmkRHJi5fnwKVFJhyKKBxe
c7P1CND9LLSZLnX2VTdbj5SBdNV+BhxypJkCvSIe4SyOq0DFIP1paBfbedO/
r7iaSqgty5Z+wFQC7ZddTHOVX2meAk1Ml+vQS66meVxPnwI0hvKSi4MyUa+3
uSQ9AdKZLlf0X6vYgqxMmZJA5D6wuswkNpefF10ERzHT1EaqI78EHWWSwVQm
X3QRHA0SfwE4jlJjsRXPCy6aY5RfXx4HTUcrv8L9ytV0jJhBK/fHAdMxyi9i
z7tM5aDpMdB0rGEqoP3Lc8DEncwk5eAxwMR0KR/PAhN3ftuYx3uER0ET06WK
PYOaKI7ng6Z/A0zzlZ+oT6MmYjqObe8poWx9BDQdB0z9LqYFyqZ4t/cwYGJ6
o9H6lItpIdu+roTy8SHQtBDaT7qYjlc2iYe9BwHT8SYRmOkJFVtbMVky8QGQ
ewLkU5+LpoKZojzF9wOaCs6ePhfHiQxyF3vvPtBxIkT8L1wcJ8HAewHHSZBb
j7tzc7LyS8c9IOBk0NTrolmkOLkZ6t2Ahuk63qdJ8MZhOkX5mf9z0HQKRPyj
LqZTlZ+tdwGmU5SfW8hUHmOSZdJdoGkxtWkzqn8fcTEttpGt8/9OwBSj64h/
xMVUqcZie9U7QFOl8nPrYRfTaWostuP5GWA6DTQ95GKqMky8t7odNFUpqUA6
tx70g7TKxrRO4p8CmtMhqx5w0VRTW3ZVPwEd1ZZJ0+930Zwh6UC7qtsAzRmQ
Vfe5aOjaYndVPwZNNcrP33vdGaqFgT2AqdZYwNfTe1xMdcDUDe065Wfx3S6m
M5Wf+bcCJqZLFv/cxVSv/Py/BTQxXXL5LhdTgzBR+t0MmBqgfaeLqVFx+lUT
wybQxPTqmxeRQ+5wMdFS2i66bwJMTK/s5AJA0ejRkVjQ3sFsF2R5mQrAzzaS
OUdvBHPOFjOput2uEhZ+ExLFMdsSMHQjuIHpUj9+qoJ2DekBBhLnOcqvMF1g
4BLl1wmRGHDr0d59T6F/ZLymjecdBAOUnUIXuPQckPsTFbTeJ7lp5g6g0YT+
PBfMvQHMbVJ+sbktUKy5GU/bpkS3nmuM5bOKG2DOz1VjsSXMjwOdwFLVzuOu
b0l07Xkilbh/CAnAdClXPSpoO8HH2Grx7CMAPgfq+cZ7RXSUdz1fHpuNMN6p
dScJJLpn6MSlYTObhH8F25guZe5WK86P+GRxeYEaix3B/AAcyHQpgLe4qXqh
saGgS9em68AGpksBvNlN1aXKL53XsgdapIZkScbHFZ6W2PAR7xrQ0QKxc5Ny
Cs9FbEC0bLtmuhrQMF3WAze6aJYpv2heBZqWQXuji2a50cRF8/uMphUyuUvF
rlk8sNEk+7+AxFbl15cfuThWQN5dCTiYXradlyI3uDhWAsgrQNNKyLwfujhW
se1X8H9X2bDXiwUcyPPBzxfIaczlID3GpHPlehfHavaJOUT6LuBoM9lQ0CVB
HIejHTxzGWhaDZH/AxdHu/Kr02WM5mKI7OtcNBcbq3n4d0BHjEkH0LUumg41
FttMbQA0HRDZ17hoLokxNXgKNF1imfTv1S6aTmMGb6bWM5o1ENlXuWjWwPB1
oGMNRPb3/ZBksmygLgUcfN9Up+92/cv06j5+0g3vVPOFgjtlM7UW9F0Ksb4W
0iFBwDo1FtsjrWF86yHi1yj/WpnAut6yWjeJ7vUQ/Z1hupXRxH6+BNArtqNT
dAeybjBOqYgFg+jeALo7wnR/h3V08H85tPVCIhZ2yRDbgX02zEXrZaCpPUzr
dwHTakDMdMnO1WG4L1d+RraBBZeD1rYwC3igbGzaGP33lJ+vbWHouZJJ9VoJ
ur+n/AxeFab7CuVXtZWA/gqwYGUY+iuNAN72rAALrowJbjB0xwKdaWaI2fy0
Mm5KQ1tDWsNw80A5XlgOWkUAi1wehvsqNRbbFC0D3IZuisWyMNxXWwE6yy8C
C64GAReFWcCXbamLLYyeSqqtPC1h6K81rLxlWgq6Dd1UnqVhuq8zmriQXwjo
md48r4TWNReGoefFkNTOC8ACpkvpuiDMAl6cNUfLaAPXzOjp8qjXOyyyOQy9
DORieD7oZrqsfc53ddOd7sAnefkEWp4STAUtvLKVbch54CKmyz7xPNdFicfc
6aSITuTNAtuenrOKGwRIbAMgQG6A9rkJQNwdkH2c0dBgDyNafgTGN7Gned8i
RzlNrqcD90MBexm9FzJbClPYzwGbWb6cI52TaH/Avsin6clw7N8IkpcAlo3g
/CXKnYhALeSl2QN8J8N6zmq5Eew8GzQyXQrs2YlYgvaO1C64Ju6J7hTQdBNI
OYvngylS/89OnI/gp1A5YNNB7iZjKa9+GwGBSOcCfVYCgriAZXHinDRXtFwM
GmEKhM7VtjEoF+J3WBF++F0eT7Wyb47JaIidjGyydDnkagirJLco/5pRz/7k
4xup7fVhleRW0bFRDzwTdMcE6OJ0ZphuPl2SYl8HjmG61PY61zEooMcI4EeU
asECpreYB19qwyzgkzS5ZtQwej6Qk8paE4b+NmEd1APPgJRlutTXMwJ06ys3
D5EyXw1DzAniVv9kMRluPp88xxT7081AXj8xXUrK6WG4E85DZaV2O8xCVRh6
PoXVq1Xy7mmA/meg9bQwC+6Q6SJNlTDwDmhXhvmAT5xl5boYfHAHZFtlmAV3
Ab5TwQd3ioDB5nlWcLAP+EheDsdPAR/cJSlFp06nhlnAAuQ0LOGsX3LylDAf
8O0Hqh/K3n4wPrgbLFgUZsE9yi9BJ4MP7gELmC43FBIE3Kv8OnMS+MCndxu6
+ClBAN+rkZ30iWCq0Lk2nGhgBgvge0Oyc68AH9wvAug8jeklg3ECdB7yEF1n
6CHMEwD9A6D7hDDdD6qx2DnS8YD+QeXXpePD0D8UEzBs7rgJeqZLXVoYZsHD
YMECQM90ObFbEIA+JuCRmA9GzN1F8cEjYMFxYRbwPU1Zy84HHzBdzuvmh/mA
76TK6cO3wQePxgQMG3oyCx4DC8rBB49ZAZaezAe9YMGx4AOm493kZBbwvWo5
mzgGfIA3vo8J84G5rW5OKI4GHzBdCvrRYRb8Eiw4CnzAdLkqHhXmgz6++h0F
6PuUX6HLAnTrBOoD80oB9xNgQWkY7idZKz7h8KTyK+68MMRPAdNc0Md0WXLO
DUPcD1hLoM0C9IqPHsGdG2bB02BBMaDvNwK4qJaEoWcBcqkqBmFMl+d/isMs
eMZoYgFHgg+YLjXxyDAfPGss0HW+1DyZIz5gutSlI8IseE6NxU4gi8AHz1kB
Sp4TSuaD5wUrOasQfMB0yclC14KQ7ah5QTS20hct/HSVVPQ54CimN+46lkrE
HNdRQdtRvcEigmwgcCX+ghqL3aOeDa6UR8Y4M2cnAAncjpo3u+JoVgs/0iZn
FAXg7xchEwtcfwdvSnX7Y9JSu8Fi0ZtStlbK/iyYDn44T+rIrEQUQXfV1M7j
rpm1LWEq+MFBKeuHw1T8BrQdnjAVQRoiJgRmj+Bm6yUQfxhMA9NlFXNYIoCg
G3jaMfwAqH9XizXIc5nwrKeZgpcgGQ9LnIKgZwMtLSCc7BOldLpxKLjmZZOx
fOKYn4Aj+ADFEBIOUDbDnB4CM/EK5PShQUnhPDaYJOdeFavpIoEP/G6GXDwk
rLa8BjblgaNftZmmY/TgsNrymswInQzkgQ9fgyDPC7PgdSOAN/8HgTWvmyzk
8niQ6yEUsMW4kl+ymAk+2AJROjPMAn7WW+6GHgg+eAMEHBjmA37CXCp5LviA
6fJQWG6YBfxcu6yuDgAf/BbS4IAAH+jlAQ+RM8lvAXr7fP1US0+mm5/blwcr
9wf0ll5s6cnQD0CI5wD6AV537B+mexCGZwPuQRCZHTb3b4KA/QD9m6w7O0z3
W8A6A3C/Be0ZYbiJMaVskE+zpoOwIQOE29PDLOCBsjSZDui3ggXRMPRbQUAU
0L8NFkwLs+BtK0Avw6aBBe+ABfuG+YAHyh7WDqQV5TBYsE+YBcMigHbRMrCF
Bo6ABVPDfDACFkwFU7eBBVkBFujUYVaqcVo3DykxQ94F3XuHoX8XdNuBhP49
0L1XGPr3AP1egP59sCAzDP37YEEmmPoBWJARZsEHYEEG+OBDEDYlzAcfggV2
4DxLr9uak6GvC1PCLPjIWqAF7Ak+IAERik16x3vPMCd8bE3QaTAZbCXJkY7+
sh15hU3RyWE2fAKJtAd44WMjwewz9ghzw6cg4ZvgBqGTYkNPZsJnIOAb4Aam
N+6it2oNPZkX+PXR6k382u4kMJUERNq6cwtz57TUc0d+e1t3MgmRtp78qpzM
5tLdwQ3c0VJPLyu3NOweJoJeNIw0l+YV5mbWDU00+ZXf3jyPX1uMNO7KLcwr
bC7ljqZp0VF9xd5NRGgJMqa5lAopP+/UuDV7VO+a6BXJSN2QhqCZuWNRZ04m
MMf08xvJZTv4wJJXrHVDFD78ZmZEM2WQE9KNBP6yR7CE6r6czOo+I4F8ryUQ
PULvBNNfaWJDRjIbIk1R80UFXqxV6lnRPqCOVFq80kvSdUNmIVfZGeeGOCn8
xjS9SF4+1Swqy3awP+iN1jTaPZM1C4o7+vnmDb0hHSSJRu9Oz6mRRQVddL1q
XRflTyfwC7W7t64r4HfB89sp4hqH6PX0INfIG+r57ead97zTNQ52B73Lq5cv
eVX2pfjyrJJBAhlkDg2e2Dwvt9B//Z+9w2/e1/UhNdgzNJBeB8eX/DlVvrTz
ar8QQK0gCf+h/6UTZqtMe4JoEyj/fAPaepiTbo3x8ps/dBT3VSRYpae68tt6
7GcAWP52mxkxqymDSCnpoJf28Tss+GkTlExS0uq2+m6jZV8aRaL1126uuKDv
EEVckfZ7AeRJFkmrWWN9MpGJ2129tSXkqZRz4tW/0d/08jvGVIKogBuh1joS
QJ8VYCeyNDryxdj62tIWTDWu4w+iEMrkbgu6d4qiqvsojFtKvzJ/hLkrYGtq
RRF3ik4ecthOaqOzaPrtF3FCQ+MrU5d9KSQx3knmYzYsBD5WkyCEL7nzjJD8
9uAk5IFtPVyOeCCVxuQDu32z6KsPQdVFJPoDW+qDCwgP7OiP/7xa0EDuJGVg
RqCN4nS+pMnAOBv1WvErGYKyAq1LlPW/bV3CwK9lI7fJQhVwOUcbuV25iT9m
EjxQJI7Ff78owbSw/lD+vKqO/rB+Mi1Jv/439g+/DJji/c8/LGia/EFDb9zf
AaM+QCM=\
\>"]]
}, Open  ]]
},
WindowSize->{1904, 1079},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{32000, 32000},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{600, 780}},
Magnification:>FEPrivate`If[
  FEPrivate`Equal[FEPrivate`$VersionNumber, 6.], 1.5, 1.5 Inherited],
FrontEndVersion->"8.0 for Microsoft Windows (32-bit) (February 23, 2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 86, 1, 43, "Input"],
Cell[CellGroupData[{
Cell[656, 25, 193, 4, 43, "Input"],
Cell[852, 31, 157, 3, 42, "Output"]
}, Open  ]],
Cell[1024, 37, 126, 2, 43, "Input"],
Cell[1153, 41, 296, 7, 70, "Input"],
Cell[1452, 50, 117, 2, 43, "Input"],
Cell[1572, 54, 266, 6, 69, "Input"],
Cell[CellGroupData[{
Cell[1863, 64, 325, 9, 70, "Input"],
Cell[2191, 75, 114, 2, 42, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2342, 82, 539, 15, 70, "Input"],
Cell[2884, 99, 120, 2, 42, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3041, 106, 536, 15, 70, "Input"],
Cell[3580, 123, 187, 5, 62, "Output"]
}, Open  ]],
Cell[3782, 131, 557, 17, 81, "Input"],
Cell[CellGroupData[{
Cell[4364, 152, 174, 3, 71, "Input"],
Cell[4541, 157, 91, 1, 42, "Output"],
Cell[4635, 160, 91, 1, 42, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4763, 166, 267, 7, 43, "Input"],
Cell[5033, 175, 17267, 291, 349, 9124, 155, "CachedBoxData", "BoxData", \
"Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)