(* 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[     81495,       1535]
NotebookOptionsPosition[     80956,       1514]
NotebookOutlinePosition[     81458,       1534]
CellTagsIndexPosition[     81415,       1531]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"clamped", "=", 
   RowBox[{"Function", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"a0", ",", "Lr", ",", "\[Alpha]"}], "}"}], ",", 
     RowBox[{"{", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"r0", "=", "1.0"}], ";", "\[IndentingNewLine]", 
       RowBox[{"l", "=", 
        RowBox[{"Lr", "*", "r0"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"rm", "=", 
        RowBox[{"3", "*", "l"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"lm", "=", 
        RowBox[{"2", "*", "l"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"nd", "=", "2"}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{"\[Mu]0", "=", "10000"}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{"m", "=", "0"}], ";", "\[IndentingNewLine]", 
       RowBox[{"\[Nu]", "=", "0.3"}], ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"\[Mu]", "[", "z_", "]"}], "=", 
        RowBox[{"\[Mu]0", "+", 
         RowBox[{"m", "*", "z"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"\[Lambda]", "[", "z_", "]"}], "=", 
        RowBox[{"2", "*", "\[Nu]", "*", 
         RowBox[{
          RowBox[{"\[Mu]", "[", "z", "]"}], "/", 
          RowBox[{"(", 
           RowBox[{"1", "-", 
            RowBox[{"2", "*", "\[Nu]"}]}], ")"}]}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"\[Lambda]\[Mu]", "=", 
        RowBox[{"2", "*", 
         RowBox[{"\[Nu]", "/", 
          RowBox[{"(", 
           RowBox[{"1", "-", 
            RowBox[{"2", "*", "\[Nu]"}]}], ")"}]}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"\[Rho]s", "=", 
        RowBox[{"0.7", "*", "20000"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"\[Rho]p", "=", "20000"}], ";", "\[IndentingNewLine]", 
       RowBox[{"\[Omega]", "=", 
        RowBox[{"a0", "*", 
         RowBox[{
          SqrtBox[
           RowBox[{
            RowBox[{"\[Mu]", "[", "l", "]"}], "/", "\[Rho]s"}]], "/", 
          "r0"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"P0", "=", "10000"}], ";", "\[IndentingNewLine]", 
       RowBox[{"M0", "=", "0"}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{"Ip", "=", 
        RowBox[{
         RowBox[{"\[Pi]", "/", "4"}], "*", 
         RowBox[{"r0", "^", "4"}]}]}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{"\[Gamma]1new", "=", "0.3"}], ";", "\[IndentingNewLine]", 
       RowBox[{"\[Gamma]1old", "=", "0.001"}], ";", "\[IndentingNewLine]", 
       RowBox[{"\[Gamma]4new", "=", "0.3"}], ";", "\[IndentingNewLine]", 
       RowBox[{"\[Gamma]4old", "=", "0.001"}], ";", "\[IndentingNewLine]", 
       RowBox[{"Ep", "=", 
        RowBox[{"1000", "*", 
         RowBox[{"\[Mu]", "[", "l", "]"}], "*", "2", "*", 
         RowBox[{"(", 
          RowBox[{"1", "+", "\[Nu]"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"\[Mu]0", "=", 
        RowBox[{"10000", "*", 
         RowBox[{"(", 
          RowBox[{"1", "+", 
           RowBox[{"0.0", "\[ImaginaryI]"}]}], ")"}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"m", "=", "0"}], ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"\[Mu]", "[", "z_", "]"}], "=", 
        RowBox[{"\[Mu]0", "+", 
         RowBox[{"m", "*", "z"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"\[Lambda]", "[", "z_", "]"}], "=", 
        RowBox[{"2", "*", "\[Nu]", "*", 
         RowBox[{
          RowBox[{"\[Mu]", "[", "z", "]"}], "/", 
          RowBox[{"(", 
           RowBox[{"1", "-", 
            RowBox[{"2", "*", "\[Nu]"}]}], ")"}]}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"While", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"Abs", "[", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"\[Gamma]1old", "-", "\[Gamma]1new"}], ")"}], "/", 
             "\[Gamma]1old"}], "]"}], ">", "0.001"}], "&&", 
          RowBox[{
           RowBox[{"Abs", "[", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"\[Gamma]4old", "-", "\[Gamma]4new"}], ")"}], "/", 
             "\[Gamma]4old"}], "]"}], ">", "0.001"}]}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"\[Gamma]1", "=", "\[Gamma]1new"}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"\[Gamma]2", "=", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"\[Lambda]\[Mu]", "+", "1"}], ")"}], "/", 
            RowBox[{"(", 
             RowBox[{"\[Lambda]\[Mu]", "+", "2"}], ")"}]}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"\[Gamma]3", "=", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"\[Lambda]\[Mu]", "+", "3"}], ")"}], "/", 
            RowBox[{"(", 
             RowBox[{"\[Lambda]\[Mu]", "+", "2"}], ")"}]}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"\[Gamma]4", "=", "\[Gamma]4new"}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"\[Gamma]5", "=", 
           RowBox[{"(", 
            RowBox[{"\[Lambda]\[Mu]", "+", "1"}], ")"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"\[Gamma]6", "=", 
           RowBox[{"(", 
            RowBox[{"\[Lambda]\[Mu]", "+", "3"}], ")"}]}], ";", 
          "\[IndentingNewLine]", "\[IndentingNewLine]", 
          RowBox[{"sq", "=", "10"}], ";", "\[IndentingNewLine]", 
          RowBox[{"While", "[", 
           RowBox[{
            RowBox[{"sq", ">", "0"}], ",", "\[IndentingNewLine]", 
            RowBox[{
             RowBox[{"s", "=", 
              RowBox[{"First", "[", 
               RowBox[{"NDSolve", "[", 
                RowBox[{
                 RowBox[{"{", 
                  RowBox[{
                   RowBox[{
                    RowBox[{
                    RowBox[{"r", "*", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "''"}], "[", "r", "]"}]}], "+", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "'"}], "[", "r", "]"}], "-", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{"\[Gamma]1", "/", "r0"}], ")"}], "^", "2"}], "*", 
                    "r", "*", 
                    RowBox[{"\[Phi]1", "[", "r", "]"}]}], "-", 
                    RowBox[{"\[Gamma]2", "*", 
                    RowBox[{
                    RowBox[{"\[Phi]2", "'"}], "[", "r", "]"}]}], "-", 
                    RowBox[{"\[Gamma]3", "*", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "[", "r", "]"}], "-", 
                    RowBox[{"\[Phi]2", "[", "r", "]"}]}], ")"}], "/", 
                    "r"}]}]}], "\[Equal]", "0"}], ",", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"r", "*", 
                    RowBox[{
                    RowBox[{"\[Phi]2", "''"}], "[", "r", "]"}]}], "+", 
                    RowBox[{
                    RowBox[{"\[Phi]2", "'"}], "[", "r", "]"}], "-", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{"\[Gamma]4", "/", "r0"}], ")"}], "^", "2"}], "*", 
                    "r", "*", 
                    RowBox[{"\[Phi]2", "[", "r", "]"}]}], "+", 
                    RowBox[{"\[Gamma]5", "*", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "'"}], "[", "r", "]"}]}], "+", 
                    RowBox[{"\[Gamma]6", "*", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "[", "r", "]"}], "-", 
                    RowBox[{"\[Phi]2", "[", "r", "]"}]}], ")"}], "/", 
                    "r"}]}]}], "\[Equal]", "0"}], ",", 
                   RowBox[{
                    RowBox[{"\[Phi]1", "[", "r0", "]"}], "\[Equal]", "1"}], 
                   ",", 
                   RowBox[{
                    RowBox[{"\[Phi]2", "[", "r0", "]"}], "\[Equal]", "1"}], 
                   ",", 
                   RowBox[{
                    RowBox[{"\[Phi]1", "[", "rm", "]"}], "\[Equal]", "0"}], 
                   ",", 
                   RowBox[{
                    RowBox[{"\[Phi]2", "[", "rm", "]"}], "\[Equal]", "0"}]}], 
                  "}"}], ",", 
                 RowBox[{"{", 
                  RowBox[{"\[Phi]1", ",", "\[Phi]2"}], "}"}], ",", 
                 RowBox[{"{", 
                  RowBox[{"r", ",", "r0", ",", "rm"}], "}"}], ",", 
                 RowBox[{"MaxSteps", "\[Rule]", "\[Infinity]"}], ",", 
                 RowBox[{"MaxStepFraction", "\[Rule]", 
                  RowBox[{"1", "/", "100"}]}], ",", 
                 RowBox[{"Method", "\[Rule]", 
                  RowBox[{"{", 
                   RowBox[{"\"\<ExplicitRungeKutta\>\"", ",", " ", 
                    RowBox[{"\"\<DifferenceOrder\>\"", "\[Rule]", "9"}]}], 
                   "}"}]}], ",", " ", 
                 RowBox[{"InterpolationOrder", "\[Rule]", "All"}]}], "]"}], 
               "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
             RowBox[{"If", "[", 
              RowBox[{
               RowBox[{
                RowBox[{
                 RowBox[{"Abs", "[", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"\[Phi]1", "[", "rm", "]"}], "/.", "s"}], ")"}], 
                  " ", "]"}], ">", "0.0001"}], "||", "  ", 
                RowBox[{
                 RowBox[{"Abs", "[", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"\[Phi]2", "[", "rm", "]"}], "/.", "s"}], ")"}], 
                  " ", "]"}], ">", "0.0001"}]}], " ", ",", 
               RowBox[{
                RowBox[{"rm", "=", 
                 RowBox[{"rm", "*", "0.8"}]}], ";", 
                RowBox[{"sq", "=", "1"}]}], ",", 
               RowBox[{"sq", "=", "0"}]}], "]"}]}]}], "]"}], ";", 
          "\[IndentingNewLine]", "\[IndentingNewLine]", 
          RowBox[{"m1", "=", 
           RowBox[{
            RowBox[{
             RowBox[{"\[Pi]", "*", "\[Rho]s", "*", 
              RowBox[{"Integrate", "[", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "[", "r", "]"}], "/.", "s"}], ")"}], 
                    "^", "2"}], "+", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]2", "[", "r", "]"}], "/.", "s"}], ")"}], 
                    "^", "2"}]}], ")"}], "*", "r"}], ")"}], ",", 
                RowBox[{"{", 
                 RowBox[{"r", ",", "r0", ",", "rm"}], "}"}], ",", 
                RowBox[{"Assumptions", "\[Rule]", 
                 RowBox[{"r0", ">", "0"}]}]}], "]"}]}], "+", 
             RowBox[{"\[Pi]", "*", 
              RowBox[{"r0", "^", "2"}], "*", "\[Rho]p"}]}], "//", "N"}]}], 
          ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
          RowBox[{"c", "=", 
           RowBox[{
            RowBox[{"Pi", "*", 
             RowBox[{"Integrate", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "[", "r", "]"}], "/.", "s"}], ")"}], 
                   "^", "2"}], "+", 
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]2", "[", "r", "]"}], "/.", "s"}], ")"}], 
                   "^", "2"}]}], ")"}], "*", "r"}], ",", 
               RowBox[{"{", 
                RowBox[{"r", ",", "r0", ",", "rm"}], "}"}], ",", 
               RowBox[{"Assumptions", "\[Rule]", 
                RowBox[{"r0", ">", "0"}]}]}], "]"}]}], "//", "N"}]}], ";", 
          "\[IndentingNewLine]", "\[IndentingNewLine]", 
          RowBox[{"k", "=", 
           RowBox[{
            RowBox[{
             RowBox[{"\[Pi]", "*", 
              RowBox[{"(", 
               RowBox[{"\[Lambda]\[Mu]", "+", "2"}], ")"}], "*", 
              RowBox[{"Integrate", "[", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    SubscriptBox["\[PartialD]", "r"], " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "[", "r", "]"}], "/.", "s"}], ")"}]}], 
                    ")"}], "^", "2"}], "*", "r"}], ")"}], ",", 
                RowBox[{"{", 
                 RowBox[{"r", ",", "r0", ",", "rm"}], "}"}], ",", 
                RowBox[{"Assumptions", "\[Rule]", 
                 RowBox[{"r0", ">", "0"}]}]}], "]"}]}], "+", 
             RowBox[{"\[Pi]", "*", 
              RowBox[{"(", 
               RowBox[{"\[Lambda]\[Mu]", "+", "3"}], ")"}], "*", 
              RowBox[{"Integrate", "[", 
               RowBox[{
                RowBox[{
                 RowBox[{"r", "*", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"(", 
                    FractionBox[
                    RowBox[{" ", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "[", "r", "]"}], "-", " ", 
                    RowBox[{"\[Phi]2", "[", "r", "]"}]}]}], "r"], ")"}], "^", 
                    "2"}], ")"}]}], "/.", "s"}], ",", 
                RowBox[{"{", 
                 RowBox[{"r", ",", "r0", ",", "rm"}], "}"}], ",", 
                RowBox[{"Assumptions", "\[Rule]", 
                 RowBox[{"r0", ">", "0"}]}]}], "]"}]}], "+", 
             "\[IndentingNewLine]", 
             RowBox[{"\[Pi]", "*", 
              RowBox[{"Integrate", "[", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    SubscriptBox["\[PartialD]", "r"], 
                    RowBox[{"(", " ", 
                    RowBox[{
                    RowBox[{"\[Phi]2", "[", "r", "]"}], "/.", "s"}], ")"}]}], 
                    ")"}], "^", "2"}], "*", "r"}], ")"}], ",", 
                RowBox[{"{", 
                 RowBox[{"r", ",", "r0", ",", "rm"}], "}"}], ",", 
                RowBox[{"Assumptions", "\[Rule]", 
                 RowBox[{"r0", ">", "0"}]}]}], "]"}]}], "+", 
             RowBox[{"2", "*", "\[Pi]", "*", 
              RowBox[{"(", "\[Lambda]\[Mu]", ")"}], "*", 
              RowBox[{"Integrate", "[", 
               RowBox[{
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "[", "r", "]"}], "-", " ", 
                    RowBox[{"\[Phi]2", "[", "r", "]"}]}], ")"}], "*", 
                   RowBox[{
                    SubscriptBox["\[PartialD]", "r"], " ", 
                    RowBox[{"\[Phi]1", "[", "r", "]"}]}]}], ")"}], "/.", 
                 "s"}], ",", 
                RowBox[{"{", 
                 RowBox[{"r", ",", "r0", ",", "rm"}], "}"}], ",", 
                RowBox[{"Assumptions", "\[Rule]", 
                 RowBox[{"r0", ">", "0"}]}]}], "]"}]}], "+", 
             RowBox[{"2", "*", "\[Pi]", "*", 
              RowBox[{"Integrate", "[", 
               RowBox[{
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "[", "r", "]"}], "-", " ", 
                    RowBox[{"\[Phi]2", "[", "r", "]"}]}], ")"}], "*", 
                   RowBox[{
                    SubscriptBox["\[PartialD]", "r"], " ", 
                    RowBox[{"\[Phi]2", "[", "r", "]"}]}]}], ")"}], "/.", 
                 "s"}], ",", 
                RowBox[{"{", 
                 RowBox[{"r", ",", "r0", ",", "rm"}], "}"}], ",", 
                RowBox[{"Assumptions", "\[Rule]", 
                 RowBox[{"r0", ">", "0"}]}]}], "]"}]}]}], "//", "N"}]}], ";", 
          "\[IndentingNewLine]", "\[IndentingNewLine]", 
          RowBox[{"m2", "=", 
           RowBox[{
            RowBox[{
             RowBox[{"\[Pi]", "*", 
              RowBox[{"r0", "^", "2"}], "*", "\[Rho]s"}], "+", 
             RowBox[{"\[Pi]", "*", "\[Rho]s", "*", 
              RowBox[{"Integrate", "[", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "[", "r", "]"}], "/.", "s"}], ")"}], 
                    "^", "2"}], "+", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]2", "[", "r", "]"}], "/.", "s"}], ")"}], 
                    "^", "2"}]}], ")"}], "*", "r"}], ")"}], ",", 
                RowBox[{"{", 
                 RowBox[{"r", ",", "r0", ",", "rm"}], "}"}], ",", 
                RowBox[{"Assumptions", "\[Rule]", 
                 RowBox[{"r0", ">", "0"}]}]}], "]"}]}]}], "//", "N"}]}], ";", 
          "\[IndentingNewLine]", "\[IndentingNewLine]", 
          RowBox[{"cs", "=", 
           RowBox[{
            RowBox[{
             RowBox[{"Pi", "*", 
              RowBox[{"Integrate", "[", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]1", "[", "r", "]"}], "/.", "s"}], ")"}], 
                    "^", "2"}], "+", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Phi]2", "[", "r", "]"}], "/.", "s"}], ")"}], 
                    "^", "2"}]}], ")"}], "*", "r"}], ")"}], ",", 
                RowBox[{"{", 
                 RowBox[{"r", ",", "r0", ",", "rm"}], "}"}], ",", 
                RowBox[{"Assumptions", "\[Rule]", 
                 RowBox[{"r0", ">", "0"}]}]}], "]"}]}], "+", 
             RowBox[{"\[Pi]", "*", 
              RowBox[{"r0", "^", "2"}]}]}], "//", "N"}]}], ";", 
          "\[IndentingNewLine]", "\[IndentingNewLine]", 
          RowBox[{"\[Chi]1", "=", 
           SqrtBox[
            RowBox[{
             FractionBox[
              RowBox[{"c", " ", "\[Mu]0"}], 
              RowBox[{"2", " ", "Ep", " ", "Ip"}]], "+", 
             FractionBox[
              SqrtBox[
               RowBox[{
                RowBox[{
                 RowBox[{"-", "4"}], " ", "Ep", " ", "Ip", " ", "k", " ", 
                 "\[Mu]0"}], "+", 
                RowBox[{
                 SuperscriptBox["c", "2"], " ", 
                 SuperscriptBox["\[Mu]0", "2"]}], "+", 
                RowBox[{"4", " ", "Ep", " ", "Ip", " ", "m1", " ", 
                 SuperscriptBox["\[Omega]", "2"]}]}]], 
              RowBox[{"2", " ", "Ep", " ", "Ip"}]]}]]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"\[Chi]2", "=", 
           SqrtBox[
            RowBox[{
             FractionBox[
              RowBox[{"c", " ", "\[Mu]0"}], 
              RowBox[{"2", " ", "Ep", " ", "Ip"}]], "-", 
             FractionBox[
              SqrtBox[
               RowBox[{
                RowBox[{
                 RowBox[{"-", "4"}], " ", "Ep", " ", "Ip", " ", "k", " ", 
                 "\[Mu]0"}], "+", 
                RowBox[{
                 SuperscriptBox["c", "2"], " ", 
                 SuperscriptBox["\[Mu]0", "2"]}], "+", 
                RowBox[{"4", " ", "Ep", " ", "Ip", " ", "m1", " ", 
                 SuperscriptBox["\[Omega]", "2"]}]}]], 
              RowBox[{"2", " ", "Ep", " ", "Ip"}]]}]]}], ";", 
          "\[IndentingNewLine]", "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"u", "[", "z_", "]"}], "=", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              SuperscriptBox["\[ExponentialE]", 
               RowBox[{
                RowBox[{"-", "z"}], " ", 
                RowBox[{"(", 
                 RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]], " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"c", " ", "M0", " ", "\[Mu]0", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"z", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", 
                    RowBox[{"2", " ", "\[Chi]2"}]}], ")"}]}]], " ", "\[Chi]1",
                     " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "-", "\[Chi]2"}], ")"}]}], "+", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"z", " ", "\[Chi]1"}], "+", 
                    RowBox[{"2", " ", "l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]}]], " ", 
                    "\[Chi]1", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "-", "\[Chi]2"}], ")"}]}], "+", 
                   RowBox[{"2", " ", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"z", " ", "\[Chi]1"}], "+", 
                    RowBox[{"l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]}]], " ", 
                    "\[Chi]1", " ", "\[Chi]2"}], "+", 
                   RowBox[{"2", " ", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"z", " ", "\[Chi]2"}], "+", 
                    RowBox[{"l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]}]], " ", 
                    "\[Chi]1", " ", "\[Chi]2"}], "+", 
                   RowBox[{"2", " ", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}], "+", 
                    RowBox[{"z", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"2", " ", "\[Chi]1"}], "+", "\[Chi]2"}], 
                    ")"}]}]}]], " ", "\[Chi]1", " ", "\[Chi]2"}], "+", 
                   RowBox[{"2", " ", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}], "+", 
                    RowBox[{"z", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", 
                    RowBox[{"2", " ", "\[Chi]2"}]}], ")"}]}]}]], " ", 
                    "\[Chi]1", " ", "\[Chi]2"}], "-", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"z", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"2", " ", "\[Chi]1"}], "+", "\[Chi]2"}], ")"}]}]],
                     " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "-", "\[Chi]2"}], ")"}], " ", 
                    "\[Chi]2"}], "+", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"z", " ", "\[Chi]2"}], "+", 
                    RowBox[{"2", " ", "l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]}]], " ", 
                    "\[Chi]2", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "\[Chi]1"}], "+", "\[Chi]2"}], ")"}]}], "-", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"z", " ", "\[Chi]1"}], "+", 
                    RowBox[{"2", " ", "l", " ", "\[Chi]2"}]}]], " ", 
                    "\[Chi]1", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}], "-", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"2", " ", "l", " ", "\[Chi]1"}], "+", 
                    RowBox[{"z", " ", "\[Chi]1"}], "+", 
                    RowBox[{"2", " ", "z", " ", "\[Chi]2"}]}]], " ", 
                    "\[Chi]1", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}], "-", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"2", " ", "l", " ", "\[Chi]1"}], "+", 
                    RowBox[{"z", " ", "\[Chi]2"}]}]], " ", "\[Chi]2", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}], "-", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"2", " ", "z", " ", "\[Chi]1"}], "+", 
                    RowBox[{"2", " ", "l", " ", "\[Chi]2"}], "+", 
                    RowBox[{"z", " ", "\[Chi]2"}]}]], " ", "\[Chi]2", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]}], ")"}]}], 
                "+", 
                RowBox[{"Ep", " ", "Ip", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"-", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"z", " ", "\[Chi]1"}], "+", 
                    RowBox[{"2", " ", "l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]}]]}], " ", 
                    SuperscriptBox["\[Chi]1", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "P0"}], "+", 
                    RowBox[{"M0", " ", "\[Chi]1"}]}], ")"}], " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "-", "\[Chi]2"}], ")"}]}], "-", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"z", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", 
                    RowBox[{"2", " ", "\[Chi]2"}]}], ")"}]}]], " ", 
                    SuperscriptBox["\[Chi]1", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{"P0", "+", 
                    RowBox[{"M0", " ", "\[Chi]1"}]}], ")"}], " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "-", "\[Chi]2"}], ")"}]}], "-", 
                   RowBox[{"2", " ", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}], "+", 
                    RowBox[{"z", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"2", " ", "\[Chi]1"}], "+", "\[Chi]2"}], 
                    ")"}]}]}]], " ", 
                    SuperscriptBox["\[Chi]1", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "P0"}], "+", 
                    RowBox[{"M0", " ", "\[Chi]1"}]}], ")"}], " ", "\[Chi]2"}],
                    "-", 
                   RowBox[{"2", " ", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"z", " ", "\[Chi]2"}], "+", 
                    RowBox[{"l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]}]], " ", 
                    SuperscriptBox["\[Chi]1", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{"P0", "+", 
                    RowBox[{"M0", " ", "\[Chi]1"}]}], ")"}], " ", "\[Chi]2"}],
                    "+", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"2", " ", "l", " ", "\[Chi]1"}], "+", 
                    RowBox[{"z", " ", "\[Chi]1"}], "+", 
                    RowBox[{"2", " ", "z", " ", "\[Chi]2"}]}]], " ", 
                    SuperscriptBox["\[Chi]1", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "P0"}], "+", 
                    RowBox[{"M0", " ", "\[Chi]1"}]}], ")"}], " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}], "+", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"z", " ", "\[Chi]1"}], "+", 
                    RowBox[{"2", " ", "l", " ", "\[Chi]2"}]}]], " ", 
                    SuperscriptBox["\[Chi]1", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{"P0", "+", 
                    RowBox[{"M0", " ", "\[Chi]1"}]}], ")"}], " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}], "-", 
                   RowBox[{"2", " ", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}], "+", 
                    RowBox[{"z", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", 
                    RowBox[{"2", " ", "\[Chi]2"}]}], ")"}]}]}]], " ", 
                    "\[Chi]1", " ", 
                    SuperscriptBox["\[Chi]2", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "P0"}], "+", 
                    RowBox[{"M0", " ", "\[Chi]2"}]}], ")"}]}], "-", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"z", " ", "\[Chi]2"}], "+", 
                    RowBox[{"2", " ", "l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]}]], " ", 
                    SuperscriptBox["\[Chi]2", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "\[Chi]1"}], "+", "\[Chi]2"}], ")"}], " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "P0"}], "+", 
                    RowBox[{"M0", " ", "\[Chi]2"}]}], ")"}]}], "+", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"2", " ", "z", " ", "\[Chi]1"}], "+", 
                    RowBox[{"2", " ", "l", " ", "\[Chi]2"}], "+", 
                    RowBox[{"z", " ", "\[Chi]2"}]}]], " ", 
                    SuperscriptBox["\[Chi]2", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}], " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "P0"}], "+", 
                    RowBox[{"M0", " ", "\[Chi]2"}]}], ")"}]}], "-", 
                   RowBox[{"2", " ", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"z", " ", "\[Chi]1"}], "+", 
                    RowBox[{"l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]}]], " ", 
                    "\[Chi]1", " ", 
                    SuperscriptBox["\[Chi]2", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{"P0", "+", 
                    RowBox[{"M0", " ", "\[Chi]2"}]}], ")"}]}], "+", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"z", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"2", " ", "\[Chi]1"}], "+", "\[Chi]2"}], ")"}]}]],
                     " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "-", "\[Chi]2"}], ")"}], " ", 
                    SuperscriptBox["\[Chi]2", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{"P0", "+", 
                    RowBox[{"M0", " ", "\[Chi]2"}]}], ")"}]}], "+", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"2", " ", "l", " ", "\[Chi]1"}], "+", 
                    RowBox[{"z", " ", "\[Chi]2"}]}]], " ", 
                    SuperscriptBox["\[Chi]2", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}], " ", 
                    RowBox[{"(", 
                    RowBox[{"P0", "+", 
                    RowBox[{"M0", " ", "\[Chi]2"}]}], ")"}]}]}], ")"}]}]}], 
               ")"}]}], ")"}], "/", 
            RowBox[{"(", 
             RowBox[{"Ep", " ", "Ip", " ", "\[Chi]1", " ", "\[Chi]2", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"c", " ", "\[Mu]0", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"2", " ", "l", " ", "\[Chi]1"}]], "+", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"2", " ", "l", " ", "\[Chi]2"}]], "-", 
                    RowBox[{"4", " ", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]]}], "+", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"2", " ", "l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]]}], ")"}], 
                    " ", 
                    SuperscriptBox["\[Chi]1", "2"]}], "-", 
                   RowBox[{"2", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "1"}], "+", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"2", " ", "l", " ", "\[Chi]1"}]]}], ")"}], " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "1"}], "+", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"2", " ", "l", " ", "\[Chi]2"}]]}], ")"}], " ", 
                    "\[Chi]1", " ", "\[Chi]2"}], "+", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"2", " ", "l", " ", "\[Chi]1"}]], "+", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"2", " ", "l", " ", "\[Chi]2"}]], "-", 
                    RowBox[{"4", " ", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]]}], "+", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"2", " ", "l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]]}], ")"}], 
                    " ", 
                    SuperscriptBox["\[Chi]2", "2"]}]}], ")"}]}], "+", 
                RowBox[{"Ep", " ", "Ip", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"\[Chi]1", " ", 
                    SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "-", "\[Chi]2"}], ")"}], "2"], " ", 
                    "\[Chi]2"}], "+", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"2", " ", "l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]], " ", 
                    "\[Chi]1", " ", 
                    SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "-", "\[Chi]2"}], ")"}], "2"], " ", 
                    "\[Chi]2"}], "-", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"2", " ", "l", " ", "\[Chi]1"}]], " ", "\[Chi]1", 
                    " ", "\[Chi]2", " ", 
                    SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}], "2"]}], "-", 
                   RowBox[{
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"2", " ", "l", " ", "\[Chi]2"}]], " ", "\[Chi]1", 
                    " ", "\[Chi]2", " ", 
                    SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}], "2"]}], "+", 
                   RowBox[{"4", " ", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{"l", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Chi]1", "+", "\[Chi]2"}], ")"}]}]], " ", 
                    RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["\[Chi]1", "4"], "+", 
                    SuperscriptBox["\[Chi]2", "4"]}], ")"}]}]}], ")"}]}]}], 
               ")"}]}], ")"}]}]}], ";", "\[IndentingNewLine]", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"u2", "[", "z_", "]"}], "=", "0"}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"du2", "[", "z_", "]"}], "=", "0"}], ";", 
          "\[IndentingNewLine]", "\[IndentingNewLine]", 
          RowBox[{"n1", "=", 
           RowBox[{
            RowBox[{
             RowBox[{"Integrate", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"\[Mu]", "[", "z", "]"}], "*", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   SubscriptBox["\[PartialD]", "z"], " ", 
                   RowBox[{"(", 
                    RowBox[{"u", "[", "z", "]"}], ")"}]}], ")"}], "^", 
                 "2"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"z", ",", "0", ",", "l"}], "}"}]}], "]"}], "+", 
             "\[IndentingNewLine]", 
             RowBox[{"Sum", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"\[Mu]", "[", 
                 RowBox[{"l", "+", 
                  RowBox[{
                   RowBox[{"lm", "/", "nd"}], "*", "j"}], "-", 
                  RowBox[{
                   RowBox[{"lm", "/", "nd"}], "/", "2"}]}], "]"}], "*", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"du2", "[", 
                   RowBox[{"l", "+", 
                    RowBox[{
                    RowBox[{"lm", "/", "nd"}], "*", "j"}], "-", 
                    RowBox[{
                    RowBox[{"lm", "/", "nd"}], "/", "2"}]}], "]"}], ")"}], 
                 "^", "2"}], "*", 
                RowBox[{"lm", "/", "nd"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"j", ",", "1", ",", "nd"}], "}"}]}], "]"}]}], "//", 
            "N"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
          RowBox[{"n2", "=", 
           RowBox[{
            RowBox[{
             RowBox[{"Integrate", "[", 
              RowBox[{
               RowBox[{"\[Rho]s", "*", 
                RowBox[{"\[Omega]", "^", "2"}], "*", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"u", "[", "z", "]"}], ")"}], "^", "2"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"z", ",", "0", ",", "l"}], "}"}]}], "]"}], "+", 
             RowBox[{"Sum", "[", 
              RowBox[{
               RowBox[{"\[Rho]s", "*", 
                RowBox[{"\[Omega]", "^", "2"}], "*", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"u2", "[", 
                   RowBox[{"l", "+", 
                    RowBox[{
                    RowBox[{"lm", "/", "nd"}], "*", "j"}], "-", 
                    RowBox[{
                    RowBox[{"lm", "/", "nd"}], "/", "2"}]}], "]"}], ")"}], 
                 "^", "2"}], "*", 
                RowBox[{"lm", "/", "nd"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"j", ",", "1", ",", "nd"}], "}"}]}], "]"}]}], "//", 
            "N"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
          RowBox[{"n3", "=", 
           RowBox[{
            RowBox[{
             RowBox[{"Integrate", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"\[Lambda]", "[", "z", "]"}], "+", 
                  RowBox[{"2", "*", 
                   RowBox[{"\[Mu]", "[", "z", "]"}]}]}], ")"}], "*", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"u", "[", "z", "]"}], ")"}], "^", "2"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"z", ",", "0", ",", "l"}], "}"}]}], "]"}], "+", 
             "\[IndentingNewLine]", 
             RowBox[{"Sum", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"\[Lambda]", "[", 
                   RowBox[{"l", "+", 
                    RowBox[{
                    RowBox[{"lm", "/", "nd"}], "*", "j"}], "-", 
                    RowBox[{
                    RowBox[{"lm", "/", "nd"}], "/", "2"}]}], "]"}], "+", 
                  RowBox[{"2", "*", 
                   RowBox[{"\[Mu]", "[", 
                    RowBox[{"l", "+", 
                    RowBox[{
                    RowBox[{"lm", "/", "nd"}], "*", "j"}], "-", 
                    RowBox[{
                    RowBox[{"lm", "/", "nd"}], "/", "2"}]}], "]"}]}]}], ")"}],
                 "*", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"u2", "[", 
                   RowBox[{"l", "+", 
                    RowBox[{
                    RowBox[{"lm", "/", "nd"}], "*", "j"}], "-", 
                    RowBox[{
                    RowBox[{"lm", "/", "nd"}], "/", "2"}]}], "]"}], ")"}], 
                 "^", "2"}], "*", 
                RowBox[{"lm", "/", "nd"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"j", ",", "1", ",", "nd"}], "}"}]}], "]"}]}], "//", 
            "N"}]}], ";", "\[IndentingNewLine]", 
          RowBox[{"n4", "=", 
           RowBox[{"n3", "/", 
            RowBox[{"(", 
             RowBox[{"(", 
              RowBox[{"\[Lambda]\[Mu]", "+", "2"}], ")"}], ")"}]}]}], ";", 
          "\[IndentingNewLine]", Cell[GraphicsData["Metafile", "\<\
CF5dJ6E]HGAYHf4PEfU^I6mgLb15CDHPAVmbKF5d0@00032@0@0006`00000000000000385000C0P00
0000000000055P00b`P00215CDH00040T30009d0000700000000000000000000Wa80090J003:0000
8040000000000000000003dD0`1kI0@0AP0002`0000P0000ADe6:`500@0L00004000008@`=/00000
F08005P200160000G00005000015CDH[8T0400`0000000007T0900`000000000940100`000000000
<4020100000400000020?b501`0<0000000000A0000<0000000001H0000<0000600000X0000@0000
0000000000090000400003<5000C0P00DP00070100010000Y?ooo`00000000000000090100000000
14008T<0H@1/06T0HP1b06T000000000000000000000000000000000000000000000000000000000
00000000000001/0Z=HK0100000<fQ/0S=LK029@>Eh<fQ/01=LK0100001df1/0l=TK0?A?>Eh<fQ/0
1=LK0200001iME5P1=LK00cJ6`0P0000oooooc`7ZP7dME5Poooooooo0H3oo`60;@81P?ooool00000
00P00008001300000@000000001H0P009@0003L^T040000000000000003_0P2PNb00@00000000000
W`00000000130640K01Y0680LP00000000000<cG6`2o8E5P:7hbHBcK6`0hea/0Y3I:H0H000010000
M=LK07CG6`2`NDQP1P0009cG6`0l1jX1I7H02000000U000030000040000U000030000040000U0000
30000040000B000030000040000H000030000000009D0000E000000000000000=@0006l000010000
DYV7@8R7Qd000000E`000040001<00001000000000000000=0D001@2001@0000800003H000160000
:00001`00017A4U30P000?oooooooooo=0D001@200000000AP000>07003D1`00ADe6:bY0000T0000
60000000P3l000200000P000P3l000200000P2501@0<0000000000Q000Dl1`00<0L0008@`=/20000
0P0000H7003GcLJJ0000000001@0200900000152000100T000>30`00100L0000000500002@800000
1@0000810@0000D000010_ooo`050000;P4H00001@0000/2000000D0000<0P0801@B00009PH?01X0
ooooo`0040000<3ooon`oooo`1<00;07000;00009PH?00`0CF5dJ5AiL6D00=0120000?X2000@0000
000000@0000]0@001@0001@240CF00D0000C0Q0410880000nP8000P000000000100002d10@050000
508T0Z471@0001<290:E2@@0000]0@001@0001@240B@10D0000C0Q04/a<L0000n`;0o@0000000901
00000P@2011CNFeRKf`0Mj8a2W;PoQ/0T?0B00cH17O0H0Mge11VKP@0000]0@802000038:i`KG1`40
003bN@P0000b2Ph3=@d10000lWT80000<PX>0k040@000?9i70000?/2H?l00000002@0@0000840P0@
DgU]HVm/07L/5`[[0?lK093`4P0<f0Ag`607Mm@@IVh40000;@4300@0003`0@802000038:]@CZ1`40
002UN@P0000b2]`0B0d10000YGT80000<P[L0<<40@000:Ei70000?/2`?h00000002@0@0000840P0@
DgU]HVm/07NR<@Ych?hK093`4P0<f0Ag`607Mm@@IVh40000;@4200@0003`0@<02000038:C@KQ2P40
000[N@P0000b2W@2=P`10000:gT80000<P[P0[d90@000?Mi2000038:MP:m2@40003gN@P0000b2Zd3
_@T10000n7T80000<PZY0Kd90@000?Ii2000038:h0;l1P40003WN@P0000b2WH2o0H10000igT80000
<PZ]0o`60@000>Qi2000038:Z@7l1P40003VN@P0000b2V04Q@<10000?GT80000<PYe12`20@000?Mi
2000038:QPD/0P40003hN@P0000b2ZP3;0810000mWT80000<PYe13400@000>Mi2000038:QPDa0040
003XN@P0000b2ZP3<@010000iWTL0000n`9Po`0000000901000000@2011DJFeULb1>IGLPDVm]HFh0
T?0B00cH17O0H0Mge11VKP@0000]0@<010000?010P080000<PZm1od70@00031i2000038:R`GZ3P40
000bN@P0000b2[d5R@h10000:WT80000<P[T0e/=0@00031i2000038:/P4A4P40000bN@P0000b2^@1
/1410000:WT80000<P[T0Il@0@00039i2000038:i0?F1040000`N@P0000b2/X0B`X10000<WT80000
<PXL0A890@0002Yi2000038:b@:j0P40000bN@P0000b2^H3QP410000<GTL0000n`;0oP0000000901
000000@2011DJFeULb1>IGLPDVm]HFh0T?0B00cH17O0H0Mge11VKP@0000]0@8010000?010`080000
<PY=1VX=0@0002Ui2000038:C@KK2`40000bN@P0000b2Td6J`T10000:7TL0000n`;0oP0000000901
0@0000@2011DJFeULb1>IGLPDVm]HFh0T?0B00cH17O0H0Mge11VKP@0000]0@<010000?010P080000
<PY=1Tl?0P0006Aj2000038:C@KC3@40001eNPP0000b2Td6MP`10000AgX80000<PYd0WHB0P0006Aj
2000038:M0;j4040001eNPP0000b2W@2b@X20000I7X80000<PY]0`X80P0006Aj2000038:[06l1`80
001TM@P0000b2W@23@H10000AgD80000<PYI1A410@00059e70000?/2H?l00000002@0@4000040P0@
E6U]IG<PCVEg859_KF5^093`4P0<f0Ag`607Mm@@IVh40000;@4200@0003`0@<02000038:a09I3`40
001cMA`0003k0/3n00000000T0410002108045=iKF9_K01g;1L:kP3o6`2@l1803=P4Ml1P1gOD46I^
100002d10`040000l04200P0000b2Td6f@T10000K7D80000<PY=1XX80@00071e2000038:M0:N3`40
001gM@P0000b2W@2j0d20000L7880000<PYd0V<50@00071b2000038:UP?R0040001WLPX0000V1Pl0
2P3ooooo0@00000070000?/240070000002l0P0000010P8RDgUcM6E]06kD46I^000:04D0RP400000
0P000:cb4P040000;@4200@0003`0@<00`0000000008@048900001P000024<3K0@0000<000000000
00000000000K@000@00003@0000100000P00000000000000002PA@0004D30000oinZ?OooocjYJZI4
oooo?_nOZSgnk`A48@0000P0001R000030000040000Q000020000680000<00000@000140000<0000
200000X0000@000000000000000900004000000D000020003000010000000000000000/0000@0000
<`D001@2000Q0000200002D0000<00001`00P2D0000<00000000P1H0000<0000000001P0000<0000
000001T0000<0000oooo01@0000<00003@000580001`0@000P000:Coool0000000000000002@0@00
000000M002930640K01Y0680LP1Y0000000000000000000000000000000000000000000000000000
0000000000000000002`0Me2NgKH9@T0000901000000000071B1008000000000=QB10040001e2000
0`0001Rf6`0SfTUPm_l0032f6`080000000003Bf6`1YID]P:VI;H4bf6`1oI4]PO2@90>6X]1Ah`a/0
SLi:H22kG@L50000[`XQhA00T04004YP0>R/0LQW=64kWZYfm5M;H03X[040000000000?ooool00000
000000000000000000000000000100000000000000000000000008l20000000000000000003L]Q/0
3CakM]bf6`14?7]f16I22hSDG0OmJ29gWo=>MFAf00P000009@0000`000020000600000`000000000
4P0000`0000100006@0000`0003oool05P0000`0000H00002P00010000000000000000T0000@0000
01@00008000V0000700000<00000000040000000000000009@0000`0000300006`000100003F0000
40@003H0000@000010800104000V0000700000@00000000020000000000000009@0000`000040000
6`000100002Q1`00908003H0000@0000U@T002@2000U0000300000<0000K000040000904000@1000
=P000100002c4`0040@00580001`0@001@000<3mool0000000000000002@0@0000000PL2011C07T0
K@1R06l0K00000000000000000000000000000000000000000000000000000000000000000000000
001kMTRg6`3D@g]f0P000:l:8N4I0;01gD9kMP0000000000000009F[]1A4]A/0@li:H0002@1l]A/0
`61nM[6FNWJ<UWYfO;DK06@1000000000000030OO`D40000000001@M2Z<T]a/00YMjMPbGNWH00000
00000000000000008;]M1`00002_2R7Q51d:X`80002Toooo0000000000000000T04000000007@00R
@`1Q06`0J@1R07800000000000000000000000000000000000000000000000000000000000000000
00000000/07M@W]ff2D906Af00P000009@0000`000050000E00005@0000:0P00=04003<2003[0@00
0@0006h0SCoBdX`oe`L00>L600010000C00000000000000000000?ooooooooooD0000?;`002N0000
E00005@0001_0`00=00009P3003[00000@0006h0SCoBdX`o=@d000h300010000C000000000000000
00000?ooooooooooD0000?;`002N0000E00005@0000h0@00=0000641003[00000@0006h0SCoBdX`o
/0@000h300010000C00000000000000000000?ooooooooooD0000?;`002N0000DP00070100060000
H?ooo`000000000000000901000000021`8045<0N@1]0680K`1/0000000000000000000000000000
00000000000000000000000000000000000000000000000000000000000007]fB;LK0=A3NgH20000
[`XQhAT0/07M@W]f0000000000000000UJ^d54Be6`000000:P000000001Z0K`4[`XQh@0000000000
l_00000000000000<1mo1@@00000000051d:XbBg6`02UgYf39MjMP000000000000000000000P^ed7
00000:l:8N4D7@ZS0P000:Coool0000000000000002@0@00000000M002800640K01Y0680LP000000
0000000000000000000000000000000000000000000000000000000000000000I7H02000000U0000
300000H0000X0000300000D0001D0000E00000l2000?0@00;080044100010000KP2=?m;BS3oZ1`00
]@@00040001<00000000000000000000oooooooooom@0000YO0007<0001D0000E00007@3000?0000
T@<0044000010000KP2=?m;BS3m83@00g0000040001<00000000000000000000oooooooooom@0000
YO0007<0001D0000E00003d1000?0000FP40044000010000KP2=?m;BS3o31000g0000040001<0000
0000000000000000oooooooooom@0000YO0007<0001B0000L04000D00030o_oo0000000000000000
T040000000870P0@D`1i06d0HP1_06`0000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000NgI8]a/0e4=kMP80002_2R7Q
6@2`0Me2NgH0000000000000000Z000000000>81_0B_2R7Q00000000002Ul27Q00000000003bl000
00000000000`7gl510000000000D7@ZS000K00:GNWH<UgYf00000000000000000000022kG@L00000
[`XQh@00000002/0^6iI1`0000000000100^0000I@0000001d008P00H@20KET70@000;1^F@Oh]A/0
jVHRMkQ^F@L:8PYc0000042f6`2Q[JYf000[0000001TMP08000002D0000<00001@0002P0000<0000
1P0005@0001D0000e080050100010`00]0400040001^08dod]:<?n4:001=1P000@0004`000000000
00000000003ooooooooooe00000[l000/@0005@0001D0000;@<00500001J0`00]0000040001^08do
d]:<?cH<001d0P000@0004`00000000000000000003ooooooooooe00000[l000/@0005@0001D0000
R08006`0002W0P00d0000040001^08dod]:<?kd9003P0P000@0004`00000000000000000003ooooo
oooooe00003gl000N`0005@0001D0000R0800540002W0P00]@000040001^08dod]:<?kd9001f0P00
0@0004`00000000000000000003ooooooooooe00003gl000N`0005@0001D0000Q@800:40002W0P00
1@400040001^08dod]:<?kd9002]0`000@0004`00000000000000000003ooooooooooe00003hl000
N`0005@0001D0000Q@8001/0002W0P00O`000040001^08dod]:<?kd9002Y0@000@0004`000000000
00000000003ooooooooooe00003fl000N`0005@0001D0000d@4006`0003`0@00d0000040001^08do
d]:<?o`6003P0P000@0004`00000000000000000003ooooooooooe00003Wl000N`0005@0001D0000
d@400540003`0@00]@000040001^08dod]:<?o`6001f0P000@0004`00000000000000000003ooooo
oooooe00003Wl000N`0005@0001D0000d@400:40003d0@001@400040001^08dod]:<?o`6002]0`00
0@0004`00000000000000000003ooooooooooe00003Xl000N`0005@0001D0000d@4001/0003d0@00
O`000040001^08dod]:<?o`6002Y0@000@0004`00000000000000000003ooooooooooe00003Vl000
N`0005@0001D0000jP000=00000G0@00=0400040001^08dod]:<?hD3001P10000@0004`000000000
00000000003ooooooooooe00000ml000/@0005@0001D0000T@000=D0002`0000>@400040001^08do
d]:<?b`2001e10000@0004`00000000000000000003ooooooooooe00003gl000N`0005@0001D0000
SP0001`1002`0000P0400040001^08dod]:<?b`200261@000@0004`00000000000000000003ooooo
oooooe00003hl000N`0005@0001D0000SP000:00002`000010400040001^08dod]:<?b`2002X0`00
0@0004`00000000000000000003ooooooooooe00003fl000N`0005@0001D00003@000=D0000/0000
>@400040001^08dod]:<?c40001e10000@0004`00000000000000000003ooooooooooe00003Wl000
N`0005@0001D00003@0001`1000`0000P0400040001^08dod]:<?c4000261@000@0004`000000000
00000000003ooooooooooe00003Xl000N`0005@0001D00003@000:00000`000010400040001^08do
d]:<?c40002X0`000@0004`00000000000000000003ooooooooooe00003Vl000N`000580001`0@00
1P00063oool0000000000000002@0@00000000L2011D06T0K@1U07<0801>06D0M`0P0580K`1]0640
KP000000000000000000000000000000000000000000000000000000000000000000000000000000
0000000000000000001kMTRg6`3D@g]f0P000:l:8N4I0;01000002X000000000>P>l1:l:8N400000
00000>K`8N40000000000:G`8N40000000000?;`000000000000030OO`D00000000001@M2Z<001/0
0YMjMPbGNWH00000000000000000000000000000:`3HKUT700000000:`0402h0001U0000000402h0
001U081]F@L10000d6iI1oRe6`3ZIR9gf6iI1f<C2TT00000@;HK0:6]ZWH002/0000006Af00P00000
9@0000`000060000:00000`000050000E00005@0000D0P00g@4002P2000<0P000@0006h0SCoBdX`o
o@L00;d700010000C00000000000000000000?ooooooooooD0000300001A0000E00005@0003Q0`00
B`400?D3001j0@000@0006h0SCoBdX`ojPh008/500010000C00000000000000000000?oooooooooo
D0000380001A0000E00005@000370`00F0400=/300270@000@0006h0SCoBdX`oR@h00;d500010000
C00000000000000000000?ooooooooooD00002X0001A0000E00005@0001i0`00g@0008d3000<0@00
0@0006h0SCoBdX`oF`d00>@300010000C00000000000000000000?ooooooooooD0000300001A0000
E00005@0002b1000B`000<H4001j00000@0006h0SCoBdX`o4A800;8100010000C000000000000000
00000?ooooooooooD0000380001A0000E00005@0002I1000F0000:d4002700000@0006h0SCoBdX`o
/1400>@100010000C00000000000000000000?ooooooooooD00002X0001A0000E00005@0001B1000
F00006H4002700000@0006h0SCoBdX`oWa000>@100010000C00000000000000000000?oooooooooo
D0000380001A0000E00005@000120@00g@0005H1000<0@000@0006h0SCoBdX`oeP@00>@300010000
C00000000000000000000?ooooooooooD0000300001A0000E00005@0002]0P003`000<42000n0000
0@0006h0SCoBdX`oB`X00<X000010000C00000000000000000000?ooooooooooD0000380001A0000
E00005@0001L0P0090000702001C00000@0006h0SCoBdX`o4PT001`100010000C000000000000000
00000?ooooooooooD00002X0001A0000E00005@0002e0000T`000<T0003200000@0006h0SCoBdX`o
^P800<T200010000C00000000000000000000?ooooooooooD0000380001A0000E00005@0001U0000
g@0007T0000<0@000@0006h0SCoBdX`oQP400>H300010000C00000000000000000000?oooooooooo
D0000340001A0000DP00070100050000`?koo`000000000000000901000000001`8045@0J@1]06D0
L`0P04h0I@1g0200DP1_06d0H@1^0000000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000000007]fB;LK0=A3NgH00000
i0@0000000121;`4[`XQh@0000000000<@0Qh@0000000000i_0Qh@0000000000YO0Qh@0000000000
l_00000000000000<1mo1@000000000051d:X`006`02UgYf39MjMP0000000000000[0;Q^F@L00000
000[00@0;P0006D0000[00@0;P0006D0P6eI1`40002`KUT7n;DK0>YV8WNhKUT72R8:M@000010]Q/0
XJfZMP00:`000000I7H02000000U0000300000D0000X0000300000H0001D0000E00007d3001I0@00
V0<00;H100010000KP2=?m;BS3mZ3@00C@H00040001<00000000000000000000oooooooooom@0000
:@0006`0001D0000E00001D3001I0@00?P<00;H100010000KP2=?m;BS3oK2`00C@H00040001<0000
0000000000000000oooooooooom@0000<P000:80001D0000E00007<2001I0@00SP800;H100010000
KP2=?m;BS3m[2@00C@H00040001<00000000000000000000oooooooooom@0000:00006`0001B0000
L04000H00030o_oo0000000000000000T040004000070P0@E01Y06d0I@1c0200CP1U07L0801B06l0
K@1Q06h0000000000000000000000000000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000003T100000000;X4_0B_2R7Q
00000000000X027Q00000000000a027Q00000000003Vl27Q00000000002Ul27Q00000000003bl000
00000000000`7gl500000000000D7@ZS000K0000000002/0f6iI1`00000002/0100^0000I@0002/0
100^0000I@20KET70@000=1^F@Oh]A/0jVHRMmQ^F@MS4`Y;0000042f6`2Q[JYf000[0000001TMP08
000002D0000<00001P0002P0000<00001@0005@0001H0000n`<005T100141000]P400040001^08do
d]:<?dl?001=1P000P0004`00000000000000000003ooooooooooe00001T07X0XP0007/0001D0000
E00009P3001I0@00`@<00;H100010000KP2=?m;BS3oC3@00C@H00040001<00000000000000000000
oooooooooom@0000M@000:80001D0000E00003d3001I0@00N0<00;H100010000KP2=?m;BS3mf3000
C@H00040001<00000000000000000000oooooooooom@0000A`000>L0001D0000F0000<d4001I0000
5PD00;H000010000KP2=?m;BS3mf4P00M0800080001<00000000000000000000oooooooooom@0000
I01j0:80001k0000E00005@0001Z1000F@0009<4002f00000@0006h0SCoBdX`onQ0007@200010000
C00000000000000000000?ooooooooooD00007D0002R0000E00005P0003>0P00F@0001L3002f0000
0@0006h0SCoBdX`ob@X007@200020000C00000000000000000000?ooooooooooD00006@0NP2R0000
N`0005@0001H00005`8009X0001P0P00m`000040001^08dod]:<?`X8001]0`000P0004`000000000
00000000003ooooooooooe00001T07X0XP0007/0001D0000F00000<2000U0000EP80088000010000
KP2=?m;BS3nl1`00[0400080001<00000000000000000000oooooooooom@0000I01e0:80002Q0000
E00005@0002C0@00F@000<h1002f00000@0006h0SCoBdX`o3@H007@200010000C000000000000000
00000?ooooooooooD00004L0003W0000E00005@0001600006P4007T0001g0@000@0006h0SCoBdX`o
4@4005T500010000C00000000000000000000?ooooooooooD000058000340000DP00070100050000
H?ooo`000000000000000901000100001`8045@0J@1]06D0L`0P04h0I@1g0200DP1_06d0H@1^0000
00000000000000000000000000000000000000000000000000000000000000000000000000000000
000000000000000000000000000000000000000000000000i0@00000002R1K`4[`XQh@0000000000
DP1e000000000000:00Qh@0000000000<@0Qh@0000000000i_0Qh@0000000000YO0Qh@0000000000
l_00000000000000<1mo1@0000000000000[0;Q^F@L00000000[00@0;P0006D0000[00@0;P0006D0
P6eI1`40002`KUT7n;DK0>YV8WNhKUT72R8:M`000010]Q/0XJfZMP00:`000000I7H02000000U0000
300000D0000X0000300000H0001D0000E0000?d3002B000030@00<4000010000KP2=?m;BS3mI3`00
a0800040001<00000000000000000000oooooooooom@0000L`0003h0001B0000L04000H00030o_oo
0000000000000000T040004000870P0@D`1i06d0HP1_06`000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000000000000000000000000000
0000000000000000000000000000000000000000003T100000000?X5_0B_2R7Q00000000001c027Q
00000000001B07D000000000000X027Q00000000000a027Q00000000003Vl27Q00000000002Ul27Q
00000000003bl00000000000000002/0f6iI1`00000002/0100^0000I@0002/0100^0000I@20KET7
0@000=1^F@Oh]A/0jVHRMmQ^F@MS4`Y=0000042f6`2Q[JYf000[0000001TMP08000002D0000<0000
1P0002P0000<00001@0005@0001D0000S`800501002n0P00]0400040001^08dod]:<?mT9001=1P00
0@0004`00000000000000000003ooooooooooe00001/l000/@0005@0001D0000>0800501001]0P00
]0400040001^08dod]:<?hX8001=1P000@0004`00000000000000000003ooooooooooe00001`l000
Z@0005@0001D00003`@0050000171000]0000040001^08dod]:<?ih?001d0P000@0004`000000000
00000000003ooooooooooe00001gl000b00005@0001H0000W@<00500003k0`00]0000040001^08do
d]:<?nP=001d0P000P0004`00000000000000000003ooooooooooe00001`l7;`Z@000:d0001D0000
E00006L1001@0000W0400;@000010000KP2=?m;BS3mS1@00M0800040001<00000000000000000000
oooooooooom@0000L?2ZMZT0001D0000E00003/0002K0000J@000?l000010000KP2=?m;BS3oR0000
UP<00040001<00000000000000000000oooooooooom@0000Io1P1h<0001B0000L04000D0000@0000
1`00000000000000_080000000070P8RD`1i07<0M01U06d000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000000000000000000000000000
0000000000000000000000000000000000000000000Z000000000:86_0B_2R7Q00000000001Wl7;`
00000000001c027Q00000000001B07D000000000000X027Q00000000000a027Q00000000003Vl27Q
00000000002Ul27Q00000000000002/0^6iI1`00000002/0100^0000I@0002/0100^0000I@20KET7
0@000;1^F@Oh]A/0jVHRMkQ^F@L:8PYi0000042f6`2Q[JYf000[0000001TMP08000002D0000<0000
1@0002P0000<00001P0002D0000<00001`00P2D0000<00000000P300000<00003`00P2D0000<0000
0@0004/0000@0000000000D0000X0000300000<0000X0000300000@0000X0000300000D0000U0000
30000000080U0000300000L0080R000030000?oooolU0000300000d0080X000030000080000R0000
30000?oooolR000030000?oooom60000L00006@00015CDH[;40002@0000H0000oan5?P0008000020
oln4?_nOZSgoooln:T0002@0000H00000020?`00080000200020?`00080000208D0700`000000000
140000`000000000AP0001@000080000AdA9@`<0000U0000300000h0080U0000300000h008160000
=00002P00015CDH[:T0002@0000H00000020?`00080000200020?`00080000208@0000P0001R0000
30000040001<0000I000000000000000<PD001<200000000000003<5000D0P00:@2Z000000000000
0020?`00000000000020?`0000000000000000000000000000000000000000008P0000`0003ooooo
AP0001`0000@0000ADe6:`90000<0000000000h0000D000000000100000D0000
\>"], "Graphics",
           GeneratedCell->False,
           CellAutoOverwrite->False,
           ImageSize->{213, 85},
           ImageMargins->0,
           ImageRegion->{{0., 1.}, {0., 1.}}], ";", "\[IndentingNewLine]", 
          Cell[GraphicsData["Metafile", "\<\
CF5dJ6E]HGAYHf4PEfU^I6mgLb15CDHPAVmbKF5d0@0002b40@0006`000000000000003l5000C0P00
00000000000g5P00b`P00215CDH00040Q2`009H0000700000000000000000000Wa80090J003:0000
8040000000000000000003dD0`1kI0@0AP0002`0000P0000ADe6:`500@0L00004000008@`=/00000
F08005P200160000G00005000015CDH[8T0400`0000000007T0900`000000000940100`000000000
<4020100000400000020?b501`0<0000000000A0000<0000000001H0000<0000600000X0000@0000
000000000009000040000405000C0P009@0000`0000>00209@0000`0000>00204P0000`000010000
DP00070100010000Y?ooo`0000000000000009010000000014008T<0H@1/06T0HP1b06T000000000
00000000000000000000000000000000000000000000000000000000000001/0Z=HK0100000<fQ/0
S=LK029@>Eh<fQ/01=LK0100001df1/0l=TK0?A?>Eh<fQ/01=LK0200001iME5P1=LK00cJ6`0P0000
oooooc`7ZP7dME5Poooooooo0H3oo`60;@81P?ooool0000000P00008001300000@000000001H0P00
9@0003L^T040000000000000003_0P2PNb00@00000000000W`00000000130640K01Y0680LP000000
00000<cG6`2o8E5P:7hbHBcK6`0hea/0Y3I:H0H000010000M=LK07CG6`2`NDQP1P0009cG6`0l1jX1
I7H02000000U000030000040000U000030000040000U000030000040000H000030000000009D0000
E000000000000000=@0006l000010000DYV7@8R7Qd000000E`000040001<00001000000000000000
?`D001@2001@0000800003H000160000:00001`00017A4U30P000?oooooooooo@0D001@200000000
AP0003P7000/1`00ADe6:bY0000T000060000000P3l000200000P000P3l000200000P2501@0<0000
000000Q000FD1P00R0H0008@`=/200000P0005h6003GcLJJ0000000081@0200900000352000100T0
00<_0`00100L0000000500002@8000001@0000810@0000D000010_ooo`050000;P4H00001@0000/2
000000D0000<0P0881@B00009PH?01X0ooooo`0040000<3ooon`ooooh1<00;07000;00009PH?00`0
CF5dJ5AiL6D00=0120000?X2000@0000000000@0000]0@001@0001@240CF00D0000C0Q048P880000
nP8000P000000000100002d10@050000508T0[l71@0001<290:c2@@0000]0@001@0001@240B^10D0
000C0Q04dA<L0000n`;0o@000000090100000P@2011CNFeRKf`0MaTF2S98C2H0j?0B00cH17O0H0Mg
B@mVCP@0000]0@802000038:i`KB2@40003bN@P0000b2Ph3D`d10000lWT80000<PX>0lh40@000?9i
70000?/2H?l00000002@0@0000840P0@DgU]HVm/07O><PX724`V0>S`4P0<f0Ag`607MdT?ITh40000
;@4300@0003`0@802000038:]@CU2@40002UN@P0000b2]`0IPd10000YGT80000<P[L0>440@000:Ei
70000?/2`?h00000002@0@0000840P0@DgU]HVm/07LI5PXcB4`V0>S`4P0<f0Ag`607MdT?ITh40000
;@4200@0003`0@<02000038:M09D3040000[N@P0000b2^02f`T10000mgT80000<PYf0]/90@000?Mi
2000038:[@?K2@40003hN@P0000b2ZT1f`T10000mWT80000<P[P0QX70@000>Mi2000038:MP8J1`40
003WN@P0000b2Zd36PL10000j7T80000<PZY0AX70@000>Ii2000038:H0BS0`40000mN@P0000b2WD4
BP810000mgT80000<PZ61DX20@000?Qi2000038:Z0=:0P40003fN@P0000b2WD4<@010000igT80000
<PZ61C400@000>Qi2000038:Z0<a0040003VNA`0003k0V3o00000000T0400000108045AYKFEc84iU
Mb1BKfeQKP3Xl1803=P4Ml1P1gM93fI>100002d10`040000l04200P0000b2[d7n0T10000<7T80000
<PZ;1B/=0@00039i2000038:_@G:3040000ZN@P0000b2^@3N@d10000<7T80000<PZb0BlB0@00039i
2000038:i07>4@40000ZN@P0000b2^@1_A010000<WT80000<P[T0o@40@00031i2000038:bP1Y2P40
000bN@P0000b2Q`1<0T10000:WT80000<P[90]P20@00039i2000038:iP>H0@40000dNA`0003k0/3n
00000000T0410000108045AYKFEc84iUMb1BKfeQKP3Xl1803=P4Ml1P1gM93fI>100002d10P040000
l04300P0000b2Td6T0d20000I7X80000<PY=1Rl;0P0004Me2000038:M0:D4P80001TNPP0000b2W@2
61410000MGX80000<PYd0^L:0P0006Aj2000038:K@<X2080001TNPP0000b2Z`1fPL20000I7D80000
<PYd0R/60@0004Me2000038:F@DP0@40001BMA`0003k0V3o00000000T0410000108045AYKFEc84iU
Mb1BKfeQKP3Xl1803=P4Ml1P1gM93fI>100002d10`040000l04200P0000b2/@2M`l10000LgDL0000
n`;0oP00000009010@000P@2011CNFeRKf`0MaTF2SE8C2H0j?0B00cH17O0H0MgB@mVCP@0000]0@80
10000?010`080000<PY=1XD:0@00071e2000038:M0:l3`40001gM@P0000b2W@21Ph20000L7880000
<PYd0X450@00071b2000038:UP?R0040001WLPX0000V1Pl02P3ooooo0@00000070000?/240070000
002l0P0000010P8RDgUcM6E]04i93fI>000:04D0RP4000000`0000Cc4P040000;@4300@0003`0@80
0`0000000008@048900001P000024<3K0@0000<00000000000000000000K@000@00003@000010000
0P00000000000000002QA@0004D30000oinZ?OooocjXnZM4oooo?_nOZSgnk`A48@0000P0001R0000
30000040000Q000020000680000<00000@000140000<0000200000X0000@00000000000000090000
4000020D000020003000010000000000000000/0000@0000@0D001@2000Q0000200002D0000<0000
1`00P2D0000<00000000P1H0000<0000000001P0000<0000000001T0000<0000oooo01@0000<0000
3@000580001`0@000P000:Coool0000000000000002@0@00000000M002930640K01Y0680LP1Y0000
0000000000000000000000000000000000000000000000000000000000000000002`0Me2NgKH9@T0
000901000000000071B1008000000000=QB10040001e20000`0001Rf6`0SfTUPm_l0032f6`080000
000003Bf6`1YID]P:VI;H4bf6`1oI4]PO2@90>6X]1Ah`a/0SLi:H12lG@L50000AQHQQa00T04004YP
0>R/0LQW=64kWZYfm5M;H03X[040000000000?ooool0000000000000000000000000000000010000
0000000000000000000008l20000000000000000003L]Q/03CakM]bf6`14?7]fM6E22hSDG0OmJ29g
Wo=>MFAf00P000009@0000`000020000600000`0000000004P0000`0000100006@0000`0003oool0
5P0000`0000H00002P00010000000000000000T0000@000081@00008000V0000700000<000000000
40000000000000009@0000`0000300006`000100003F000040@003H0000@00008P800104000V0000
700000@00000000020000000000000009@0000`0000400006`000100002o1`00908003H0000@0000
/`T002@2000U0000300000<0000K000040000:h4000@1000=P000100003A4`0040@00580001`0@00
1@000<3mool0000000000000002@0@0000000PL2011C07T0K@1R06l0K00000000000000000000000
000000000000000000000000000000000000000000000000001kMTRg6`3D@g]f0P0004HF8HLI0;01
gD9kMP0000000000000009F[]1A4]A/0@li:H0002@1l]A/0`61nM[6FNWJ<UWYfO;DK06@100000000
0000060j>P040000000002dM2TTT]a/00YMjMPbGNWH0000000000000000000004;aM1`0000165R67
;Ad:B@80002Toooo0000000000000000T04000000007@00R@`1Q06`0J@1R07800000000000000000
00000000000000000000000000000000000000000000000000000000/07M@W]ff2D906Af00P00000
9@0000`000050000E00005@0002@0P00=0400;T2003[0@000@0009QnSCoBdX`odPT00>L600010000
C00000000000000000000?ooooooooooD0000?;`002M0000E00005@0001j0`00=0000:<3003[0000
0@0009QnSCoBdX`oD`d000h300010000C00000000000000000000?ooooooooooD0000?;`002M0000
E00005@000110@00=00006X1003[00000@0009QnSCoBdX`ocP@000h300010000C000000000000000
00000?ooooooooooD0000?;`002M0000DP00070100060000H?ooo`00000000000000090100000002
1`8045<0N@1]0680K`1/000000000000000000000000000000000000000000000000000000000000
000000000000000000000000000007]fB;LK0=A3NgH20000AQHQQaT0/07M@W]f0000000000000000
UJ^d54Be6`000000:P000000001Z0H`4AQHQQ`0000000000l_00000000000000H3Xj00@000000000
;Ad:BBBg6`02UgYf39MjMP000000000000000000000@_5d7000004HF8HL]7@Y90P000:Coool00000
00000000002@0@00000000M002800640K01Y0680LP00000000000000000000000000000000000000
00000000000000000000000000000000I7H02000000U0000300000H0000X0000300000D0001D0000
E00009D2000?0@00/P80044100010000V7j=?m;BS3oU2@00]@@00040001<00000000000000000000
oooooooooom@0000YO142g<0001D0000E00007l3000?0000W0<0044000010000V7j=?m;BS3mV3@00
g0000040001<00000000000000000000oooooooooom@0000YO2ZMW<0001D0000E00004H1000?0000
H`40044000010000V7j=?m;BS3oQ1000g0000040001<00000000000000000000oooooooooom@0000
YO0007<0001B0000L04000D00030o_oo0000000000000000T040000000870P0@D`1i06d0HP1_06`0
00000000000000000000000000000000000000000000000000000000000000000000000000000000
000000000000000000000000NgI8]a/0e4=kMP8000165R676@2`0Me2NgH0000000000000000Z0000
00000>81S0A65R6700000000002Ul26700000000003bl00000000000001P>SX010000000000]7@Y9
000K00:GNWH<UgYf00000000000000000000012lG@L00000AQHQQ`00000002/0B>e42`0000000000
0`0b0000@`0000001d008P00H@3`jd@;0@00043]A0_h]A/0jVHRMdS]A0_Y50[i0000042f6`2Q[JYf
000[0000001TMP08000002D0000<00001@0002P0000<00001P0005@0001D0000=`<00500001T0`00
]0000040002HOXdod]:<?e@<001d0P000@0004`00000000000000000003ooooooooooe00000[l000
/00005@0001D0000TP8006`0002a0P00d0000040002HOXdod]:<?m/9003P0P000@0004`000000000
00000000003ooooooooooe00003gl000N`0005@0001D0000TP800540002a0P00]@000040002HOXdo
d]:<?m/9001f0P000@0004`00000000000000000003ooooooooooe00003gl008N`0005@0001D0000
S`800:40002a0P001@400040002HOXdod]:<?m/9002]0`000@0004`00000000000000000003ooooo
oooooe00003hl4@;N`0005@0001D0000S`8001/0002a0P00O`000040002HOXdod]:<?m/9002Y0@00
0@0004`00000000000000000003ooooooooooe00003fl:UfN`0005@0001D0000fP4006`0003i0@00
d0000040002HOXdod]:<?aX7003P0P000@0004`00000000000000000003ooooooooooe00003Wl000
N`0005@0001D0000fP400540003i0@00]@000040002HOXdod]:<?aX7001f0P000@0004`000000000
00000000003ooooooooooe00003Wl000N`0005@0001D0000fP400:40003m0@001@400040002HOXdo
d]:<?aX7002]0`000@0004`00000000000000000003ooooooooooe00003Xl000N`0005@0001D0000
fP4001/0003m0@00O`000040002HOXdod]:<?aX7002Y0@000@0004`00000000000000000003ooooo
oooooe00003Vl000N`0005@0001D0000l`000=00000P0@00=0400040002HOXdod]:<?j<3001P1000
0@0004`00000000000000000003ooooooooooe00000ml000/00005@0001D0000V@000=D0002h0000
>@400040002HOXdod]:<?dX2001e10000@0004`00000000000000000003ooooooooooe00003gl000
N`0005@0001D0000UP0001`1002h0000P0400040002HOXdod]:<?dX200261@000@0004`000000000
00000000003ooooooooooe00003hl000N`0005@0001D0000UP000:00002h000010400040002HOXdo
d]:<?dX2002X0`000@0004`00000000000000000003ooooooooooe00003fl0XbN`0005@0001D0000
3@000=D0000/0000>@400040002HOXdod]:<?c40001e10000@0004`00000000000000000003ooooo
oooooe00003Wl:YfN`0005@0001D00003@0001`1000`0000P0400040002HOXdod]:<?c4000261@00
0@0004`00000000000000000003ooooooooooe00003Xl27KN`0005@0001D00003@000:00000`0000
10400040002HOXdod]:<?c40002X0`000@0004`00000000000000000003ooooooooooe00003Vl020
N`000580001`0@001P00063oool0000000000000002@0@00000000L2011D06T0K@1U07<0801>06D0
M`0P0580K`1]0640KP00000000000000000000000000000000000000000000000000000000000000
00000000000000000000000000000000001kMTRg6`3D@g]f0P0004HF8HLI0;01000002X000000000
:P><14HF8HL0000000000>K`8HL0000000000:G`8HL0000000000?;`000000000000060j>P000000
000002dM2TT001/00YMjMPbGNWH00000000000000000000000000000:`1XkD@;00000000:`030380
001300000003038000130?3[A0/10000H>e42oRe6`3ZIR9gJ>e42a0R2/L00000@;HK0:6]ZWH002/0
000006Af00P000009@0000`000060000:00000`000050000E00005@0002J0P00g@400:h2000<0P00
0@0009QnSCoBdX`on0T00;d700010000C00000000000000000000?ooooooooooD0000300001A0000
E00005@0001_0`00B`4008<3001j0@000@0009QnSCoBdX`o:`d008/500010000C000000000000000
00000?ooooooooooD0000380001A0000E00005@0001F0`00F04006X300270@000@0009QnSCoBdX`o
bP`00;d500010000C00000000000000000000?ooooooooooD00002X0001A0000E00005@000240`00
g@0009P3000<0@000@0009QnSCoBdX`oN@d00>@300010000C00000000000000000000?oooooooooo
D0000300001A0000E00005@0002n1000B`000=84001j00000@0009QnSCoBdX`o;a800;8100010000
C00000000000000000000?ooooooooooD0000380M01A0000E00005@0002U1000F0000;T400270000
0@0009QnSCoBdX`ocQ400>@100010000C00000000000000000000?ooooooooooD00002X0:P1A0000
E00005@0001N1000F0000784002700000@0009QnSCoBdX`o_A000>@100010000C000000000000000
00000?ooooooooooD0000380OWIA0000E00005@0001;0@00g@0005l1000<0@000@0009QnSCoBdX`o
m0@00>@300010000C00000000000000000000?ooooooooooD0000300:`1A0000E00005@0002g0P00
3`000</2000n00000@0009QnSCoBdX`oJ@X00<X000010000C00000000000000000000?oooooooooo
D0000380001A0000E00005@0001V0P00900007X2001C00000@0009QnSCoBdX`o<0T001`100010000
C00000000000000000000?ooooooooooD00002X0001A0000E00005@0002n0000T`000=8000320000
0@0009QnSCoBdX`of0800<T200010000C00000000000000000000?ooooooooooD0000380001A0000
E00005@0001Z0000g@0007h0000<0@000@0009QnSCoBdX`oV0400>H300010000C000000000000000
00000?ooooooooooD00003@0001A0000DP00070100050000`?koo`00000000000000090100010000
1`8045@0J@1]06D0L`0P04h0I@1g0200DP1_06d0H@1^000000000000000000000000000000000000
000000000000000000000000000000000000000000000000000000000000000000000000000007]f
B;LK0=A3NgH00000i0@00000000b18`4AQHQQ`0000000000=00QQ`0000000000i_0QQ`0000000000
YO0QQ`0000000000l_00000000000000H3Xj000000000000;Ad:B@006`02UgYf39MjMP0000000000
000[04S]A0/00000000[00<0<P0004<0000[00<0<P0004<0l>]42`400010kD@;n;DK0>YV8WM8kD@;
jA@:n`000010]Q/0XJfZMP00:`000000I7H02000000U0000300000D0000X0000300000H0001D0000
F00008X3001I0@00d`<00;H100010000V7j=?m;BS3n@3@00C@H00080001<00000000000000000000
oooooooooom@0000I01j0:40001k0000E00005P0003[0P00F@400503002f0@000@0009QnSCoBdX`o
;`/004d600020000C00000000000000000000?ooooooooooD00004L0M@3V0000X@0005@0001H0000
f@@005T0000R1@00]P000040002HOXdod]:<?i@B001d0P000P0004`00000000000000000003ooooo
oooooe00001T07X0X@0007/0001D0000E00007H4001I0000W`@00;H000010000V7j=?m;BS3lH4@00
M0800040001<00000000000000000000oooooooooom@0000M@2=?j40001D0000F0000=P2001I0000
8@<00;H000010000V7j=?m;BS3oW2P00M0800080001<00000000000000000000oooooooooom@0000
I01j0:40001k0000E00005P0000Q0P00VP0006X2003g00000@0009QnSCoBdX`o:0P006d300020000
C00000000000000000000?ooooooooooD00006@0NP2Q0000N`0005@0001H0000308002D0001O0P00
PP000040002HOXdod]:<?mX7002/0@000P0004`00000000000000000003ooooooooooe00001T07D0
X@000:40001D0000E00009`1001I0000e`400;H000010000V7j=?m;BS3l[1P00M0800040001<0000
0000000000000000oooooooooom@0000A`000>H0001D0000E00004X0000J0@00O@0007L100010000
V7j=?m;BS3lP0@00F@D00040001<00000000000000000000oooooooooom@0000DP0^0<@0001B0000
L04000H0001Poooo0000000000000000T040004000070P0@E01Y06d0I@1c0200CP1U07L0801B06l0
K@1Q06h0000000000000000000000000000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000003T1000000000X5S0A65R67
00000000001B07D000000000000d026700000000003Vl26700000000002Ul26700000000003bl000
00000000001P>SX000000000000]7@Y9000K0000000002/0J>e42`00000002/00`0b0000@`0002/0
0`0b0000@`3`jd@;0@00063]A0_h]A/0jVHRMfS]A0/@8P[90000042f6`2Q[JYf000[0000001TMP08
000002D0000<00001P0002P0000<00001@0005@0001D00002@@00980000H1000`@000040002HOXdo
d]:<?gL?00340P000@0004`00000000000000000003ooooooooooe00001c08`o?@000580001`0@00
1@000<3nool0000000000000002@0@000@000PL2011C07T0K@1R06l0K00000000000000000000000
00000000000000000000000000000000000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000>@400000000HPF<14HF8HL00000
000007<08HL0000000000580M@000000000003@08HL0000000000>K`8HL0000000000:G`8HL00000
00000?;`000000000000060j>P00000000000000:`18kD@;00000000:`03038000130000:`030380
00130?3[A0/10000@>e42oRe6`3ZIR9gB>e42nTD2_d00000@;HK0:6]ZWH002/0000006Af00P00000
9@0000`000050000:00000`000060000E00005@0002o0P00D0400?@2002d0@000@0009QnSCoBdX`o
Q@X004d600010000C00000000000000000000?ooooooooooD000073`002Y0000E00005@0000K1000
D00005<4002d00000@0009QnSCoBdX`o_0l007@200010000C00000000000000000000?oooooooooo
D00007O`00370000E00005P0002Y0`00D00000L4002d00000@0009QnSCoBdX`o1Ph007@200020000
C00000000000000000000?ooooooooooD000073`L_2Y0000[00005@0001D0000L0400500002U0@00
]0000040002HOXdod]:<?h45001d0P000@0004`00000000000000000003ooooooooooe00001`l000
Z@0005@0001D0000>`0009/0001Y0000o`000040002HOXdod]:<?n80002F0`000@0004`000000000
00000000003ooooooooooe00001Wl5T7PP000580001`0@001P0001000007000000000000002l0P00
000000L20R9C07T0L`1d06D0K@000000000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000000000000000000000000000
0000000000000000000002X000000000nPF<14HF8HL00000000006O`L_000000000007<08HL00000
00000580M@000000000003@08HL0000000000>K`8HL0000000000:G`8HL0000000000?;`00000000
00000000:`1XkD@;00000000:`03038000130000:`03038000130?3[A0/10000H>e42oRe6`3ZIR9g
J>e42a0R2//00000@;HK0:6]ZWH002/0000006Af00P000009@0000`000060000:00000`000050000
9@0000`0000700209@0000`000000020<00000`0000?00209@0000`000010000B`00010000000000
1@0002P0000<00000`0002P0000<0000100002P0000<00001P0002D0000<00000000P2D0000<0000
1`00P280000<0000ooooobD0000<00003@00P2P0000<00000P000280000<0000ooooob80000<0000
ooooodH0001`0000I00004E=AR//@000900001P0001DRXDn0000P000083och@noinZ?OooochZ@000
900001P00000080o0000P0000800080o0000P000080Q@0L0300000000004@0003000000000160000
500000P00017A4U30`0002D0000<00003P00P2D0000<00003P00P4H0000d0000:00004E=AR/Z@000
900001P00000080o0000P0000800080o0000P000080Q000020000680000<00000@0004`0001T0000
00000000000o1@004`80000000000000@0D001@2000Y0:X0000000000000080o000000000000080o
0000000000000000000000000000000000000000000R000030000?oooom60000700001000015CDH[
0T0000`0000000003P0001@000000000400001@0
\>"], "Graphics",
           ImageSize->{215, 85},
           ImageMargins->0,
           ImageRegion->{{0., 1.}, {0., 1.}}], ";", "\[IndentingNewLine]", 
          RowBox[{"\[Gamma]1new", "=", 
           SqrtBox[
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{"n1", "+", "n2"}], ")"}], "/", "n3"}], "*", 
             RowBox[{"r0", "^", "2"}]}]]}], ";", "\[IndentingNewLine]", 
          RowBox[{"\[Gamma]4new", "=", 
           SqrtBox[
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{"n1", "+", "n2"}], ")"}], "/", "n4"}], "*", 
             RowBox[{"r0", "^", "2"}]}]]}], ";", "\[IndentingNewLine]", 
          RowBox[{"\[Gamma]1old", "=", "\[Gamma]1"}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"\[Gamma]4old", "=", "\[Gamma]4"}], ";"}]}], 
        "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{"umax", "=", 
        RowBox[{"u", "[", "0", "]"}]}]}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.4978645937622366`*^9, 3.4978645949322367`*^9}, {
  3.609680456333393*^9, 3.609680475603779*^9}, {3.6098234091923466`*^9, 
  3.609823409687396*^9}, {3.609824544348385*^9, 3.609824548408385*^9}, {
  3.609825268779296*^9, 3.6098252696753855`*^9}, {3.609835209336076*^9, 
  3.609835228233786*^9}, {3.612776301752453*^9, 3.612776303742456*^9}, {
  3.612777386830193*^9, 3.612777456007306*^9}, {3.6127780815162077`*^9, 
  3.612778082758211*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"count", "=", "40"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "0"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"(", 
     RowBox[{"count", "+", "20"}], ")"}]}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{
    RowBox[{"f2", "=", 
     RowBox[{"clamped", "[", 
      RowBox[{
       RowBox[{"i", "*", 
        RowBox[{"1", "/", "count"}]}], ",", "20", ",", "1.0"}], "]"}]}], ";", 
    RowBox[{
     RowBox[{"uamp", "[", "i", "]"}], "=", 
     RowBox[{"f2", "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], ";", 
    RowBox[{"Print", "[", 
     RowBox[{"i", ",", "\"\<   \>\"", ",", 
      RowBox[{"uamp", "[", "i", "]"}]}], "]"}]}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.497864510707836*^9, 3.497864533702236*^9}, {
   3.497864564403036*^9, 3.4978645730922365`*^9}, {3.4978646095494366`*^9, 
   3.4978646357730365`*^9}, {3.497864687159436*^9, 3.497864688313836*^9}, {
   3.4978651559238367`*^9, 3.4978651579050364`*^9}, 3.497865231693036*^9, {
   3.4978652631894364`*^9, 3.4978652888670363`*^9}, {3.6096825685280247`*^9, 
   3.6096825702680945`*^9}, {3.6096827298652964`*^9, 
   3.6096827539862614`*^9}, {3.6096828418629637`*^9, 3.609682842482976*^9}, {
   3.609683018942753*^9, 3.6096830195827656`*^9}, {3.6096830692258124`*^9, 
   3.6096830878361845`*^9}, {3.6096836229556627`*^9, 
   3.6096836252257075`*^9}, {3.6096837702388535`*^9, 3.609683801293481*^9}, {
   3.609684034905188*^9, 3.6096840385352607`*^9}, {3.6098202505510087`*^9, 
   3.609820255718009*^9}, {3.609820685290416*^9, 3.609820685609416*^9}, {
   3.609820954715416*^9, 3.609820961194416*^9}, {3.6098212293473787`*^9, 
   3.6098212295543995`*^9}, {3.609821427242815*^9, 3.6098214359778147`*^9}, {
   3.609821577809815*^9, 3.609821578368815*^9}, {3.609821719930437*^9, 
   3.6098217204264865`*^9}, {3.609823249644393*^9, 3.6098232522196503`*^9}, {
   3.60982530091951*^9, 3.609825302231641*^9}, {3.6098257347138844`*^9, 
   3.609825735449958*^9}, {3.6098282708274803`*^9, 3.6098282743638334`*^9}, {
   3.6098320886946335`*^9, 3.6098320892146907`*^9}}],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Abs", "[", 
      RowBox[{"uamp", "[", "i", "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "0", ",", "count"}], "}"}]}], "]"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", " ", 
   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4978649730138364`*^9, 3.4978650841638365`*^9}, {
  3.4978653022986364`*^9, 3.4978653077898364`*^9}, {3.497865509794236*^9, 
  3.4978655905086365`*^9}, {3.609819192590416*^9, 3.609819193404356*^9}, {
  3.609819731508009*^9, 3.609819732148009*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.612776368238552*^9, 3.6127763682485523`*^9}}]
},
WindowSize->{1424, 765},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{1, 1},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{595.1999999999999, 841.8}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 77567, 1430, 2223, "Input"],
Cell[78115, 1452, 2098, 41, 52, "Input"],
Cell[80216, 1495, 639, 14, 31, "Input"],
Cell[80858, 1511, 94, 1, 31, "Input"]
}
]
*)

(* End of internal cache information *)