(* 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[     40928,        882]
NotebookOptionsPosition[     40559,        867]
NotebookOutlinePosition[     41014,        884]
CellTagsIndexPosition[     40971,        881]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", 
    RowBox[{"-", "5"}]}], ",", 
   RowBox[{"i", "<", "2"}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"r0", "=", "1"}], ";", "\[IndentingNewLine]", 
    RowBox[{"m", "=", "2500"}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Mu]0", "=", "0"}], ";", "\[IndentingNewLine]", 
    RowBox[{"l", "=", 
     RowBox[{"20", "*", "r0"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"\[Mu]", "[", "z_", "]"}], "=", 
     RowBox[{"\[Mu]0", "+", 
      RowBox[{"m", "*", "z"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Nu]", "=", "0.49"}], ";", "\[IndentingNewLine]", 
    RowBox[{"P0", "=", "100"}], ";", "\[IndentingNewLine]", 
    RowBox[{"M0", "=", "0"}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Rho]s", "=", "1"}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Rho]p", "=", "10"}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Lambda]", "=", 
     RowBox[{"2", "*", "\[Nu]", "*", 
      RowBox[{
       RowBox[{"\[Mu]", "[", "z", "]"}], "/", 
       RowBox[{"(", 
        RowBox[{"1", "-", 
         RowBox[{"2", "*", "\[Nu]"}]}], ")"}]}]}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Ip", "=", 
     RowBox[{
      RowBox[{"\[Pi]", "/", "4"}], "*", 
      RowBox[{"r0", "^", "4"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Omega]", "=", "0"}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Gamma]1new", "=", "1"}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Gamma]1old", "=", "0.001"}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Gamma]4new", "=", "1"}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Gamma]4old", "=", "0.001"}], ";", "\[IndentingNewLine]", 
    RowBox[{"rm", "=", 
     RowBox[{"15", "*", "r0"}]}], ";", "\[IndentingNewLine]", 
    "\[IndentingNewLine]", 
    RowBox[{"\[Lambda]\[Mu]", "=", 
     RowBox[{"2", "*", 
      RowBox[{"\[Nu]", "/", 
       RowBox[{"(", 
        RowBox[{"1", "-", 
         RowBox[{"2", "*", "\[Nu]"}]}], ")"}]}]}]}], ";", 
    "\[IndentingNewLine]", "\[IndentingNewLine]", 
    RowBox[{"kr", "=", 
     RowBox[{"10", "^", 
      RowBox[{"(", "i", ")"}]}]}], ";", "\[IndentingNewLine]", 
    "\[IndentingNewLine]", 
    RowBox[{"Ep", "=", 
     RowBox[{"kr", "/", 
      RowBox[{"(", 
       RowBox[{"Ip", "/", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"\[Mu]", "[", "l", "]"}], "*", "2", "*", 
          RowBox[{"(", 
           RowBox[{"1", "+", "\[Nu]"}], ")"}], "*", 
          RowBox[{"l", "^", "4"}]}], ")"}]}], ")"}]}]}], ";", 
    "\[IndentingNewLine]", "\[IndentingNewLine]", "\[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]", "\[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]", "\[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[{"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", "]"}], "-", 
               "\[IndentingNewLine]", 
               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"}]}], "}"}],
            ",", "\[IndentingNewLine]", 
           RowBox[{"{", 
            RowBox[{"\[Phi]1", ",", "\[Phi]2"}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{"r", ",", "r0", ",", "rm"}], "}"}], ",", 
           RowBox[{"MaxSteps", "\[Rule]", "\[Infinity]"}]}], "]"}], "]"}]}], 
       ";", "\[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"}], "}"}], ",", 
             "\[IndentingNewLine]", 
             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"}], "}"}], ",", 
             "\[IndentingNewLine]", 
             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]", ")"}], "*", "\[IndentingNewLine]", 
           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]", "*", "\[IndentingNewLine]", 
           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]", 
       RowBox[{"h", "=", 
        RowBox[{"cs", "*", "m"}]}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{"u", "=", 
        RowBox[{"First", "[", 
         RowBox[{"y", "/.", 
          RowBox[{"NDSolve", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{
               RowBox[{
                RowBox[{"Ep", "*", "Ip", "*", 
                 RowBox[{
                  RowBox[{"y", "''''"}], "[", "z", "]"}]}], "-", 
                RowBox[{"c", "*", 
                 RowBox[{"\[Mu]", "[", "z", "]"}], "*", 
                 RowBox[{
                  RowBox[{"y", "''"}], "[", "z", "]"}]}], "+", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"k", "*", 
                    RowBox[{"\[Mu]", "[", "z", "]"}]}], "-", 
                   RowBox[{"m1", "*", 
                    RowBox[{"\[Omega]", "^", "2"}]}]}], ")"}], "*", 
                 RowBox[{"y", "[", "z", "]"}]}]}], "\[Equal]", "0"}], ",", 
              "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{
                RowBox[{"Ep", "*", "Ip", "*", 
                 RowBox[{
                  RowBox[{"y", "'''"}], "[", "0", "]"}]}], "-", 
                RowBox[{"c", "*", 
                 RowBox[{"\[Mu]", "[", "0", "]"}], "*", 
                 RowBox[{
                  RowBox[{"y", "'"}], "[", "0", "]"}]}], "-", "P0"}], 
               "\[Equal]", "0"}], ",", 
              RowBox[{
               RowBox[{
                RowBox[{"y", "''"}], "[", "l", "]"}], "\[Equal]", "0"}], ",", 
              RowBox[{
               RowBox[{
                RowBox[{"Ep", "*", "Ip", "*", 
                 RowBox[{
                  RowBox[{"y", "''"}], "[", "0", "]"}]}], "+", "M0"}], 
               "\[Equal]", "0"}], ",", 
              RowBox[{
               RowBox[{
                RowBox[{"Ep", "*", "Ip", "*", 
                 RowBox[{
                  RowBox[{"y", "'''"}], "[", "l", "]"}]}], "-", 
                RowBox[{"c", "*", 
                 RowBox[{"\[Mu]", "[", "l", "]"}], "*", 
                 RowBox[{
                  RowBox[{"y", "'"}], "[", "l", "]"}]}], "+", 
                "\[IndentingNewLine]", 
                RowBox[{
                 FractionBox["1", 
                  SqrtBox["cs"]], 
                 SqrtBox["k"], " ", "*", 
                 RowBox[{"\[Mu]", "[", "l", "]"}], " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "cs"}], "-", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{"h", "-", 
                    FractionBox[
                    RowBox[{
                    SqrtBox["cs"], " ", "m2", " ", 
                    SuperscriptBox["\[Omega]", "2"]}], 
                    SqrtBox["k"]]}], ")"}], " ", 
                    RowBox[{"HypergeometricU", "[", 
                    RowBox[{
                    RowBox[{"1", "+", 
                    FractionBox[
                    RowBox[{"h", "-", 
                    FractionBox[
                    RowBox[{
                    SqrtBox["cs"], " ", "m2", " ", 
                    SuperscriptBox["\[Omega]", "2"]}], 
                    SqrtBox["k"]]}], 
                    RowBox[{"2", " ", "cs", " ", "m"}]]}], ",", "2", ",", 
                    FractionBox[
                    RowBox[{"2", " ", 
                    SqrtBox["k"], " ", "*", 
                    RowBox[{"\[Mu]", "[", "l", "]"}]}], 
                    RowBox[{
                    SqrtBox["cs"], " ", "m"}]]}], "]"}]}], ")"}], "/", 
                    "\[IndentingNewLine]", 
                    RowBox[{"(", 
                    RowBox[{"m", " ", 
                    RowBox[{"HypergeometricU", "[", 
                    RowBox[{
                    FractionBox[
                    RowBox[{"h", "-", 
                    FractionBox[
                    RowBox[{
                    SqrtBox["cs"], " ", "m2", " ", 
                    SuperscriptBox["\[Omega]", "2"]}], 
                    SqrtBox["k"]]}], 
                    RowBox[{"2", " ", "cs", " ", "m"}]], ",", "1", ",", 
                    FractionBox[
                    RowBox[{"2", " ", 
                    SqrtBox["k"], "*", 
                    RowBox[{"\[Mu]", "[", "l", "]"}]}], 
                    RowBox[{
                    SqrtBox["cs"], " ", "m"}]]}], "]"}]}], ")"}]}]}], ")"}], 
                 " ", 
                 RowBox[{"y", "[", "l", "]"}]}]}], "\[Equal]", "0"}]}], "}"}],
             ",", "y", ",", 
            RowBox[{"{", 
             RowBox[{"z", ",", "0", ",", "l"}], "}"}], ",", 
            RowBox[{"MaxSteps", "\[Rule]", "\[Infinity]"}]}], "]"}]}], 
         "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", Cell[
       GraphicsData["Metafile", "\<\
CF5dJ6E]HGAYHf4PEfU^I6mgLb15CDHPAVmbKF5d0@00022<0@0006`00000000000000<H4003S0000
00000000000g5000`0<00215CDH00040S20007H0000600000000000000000000Wa80090J003:0000
8040000000000000000003dD0`1kI0@0AP0002`0000P0000ADe6:`500@0L00004000008@`=/00000
F08005P200160000G00005000015CDH[8T0400`0000000007T0900`000000000940100`000000000
<4020100000400000020?b501`0<0000000000A0000<0000000001H0000<0000600000X0000@0000
000000000009000040000<H4003R00009@0000`0000>00209@0000`0000>00204P0000`000010000
DP00070100010000Y?ooo`0000000000000009010000000014008T<0H@1/06T0HP1b06T000000000
00000000000000000000000000000000000000000000000000000000000002`0?>4/0100002Pi2`0
8>8/0>A@dVBPi2`0V>4/01000008hb`0Q>@/08Y?dVBPi2`0V>4/0200001iMOIOV>4/0:3T;00P0000
oooooa`Wb07dMOIOooooooooooooo`60;@81P?ooool0000000P00008001300000@000000001H0P00
9@0003L^T040000000000000003_0P2PNb00@00000000000W`00000000130640K01Y0680LP000000
0`00063R;02o8OIO:7kGH<3U;03<hB`0Y3K_G`H0000100002>8/00SR;02`NNeO1P00033R;00L9lP1
I7H02000000U000030000040000U000030000040000U000030000040000H000030000000009D0000
E000000000000000=@0006l000010000DYV7@8R7Qd000000E`000040001<00001000000000000000
aP@00><0001@0000800003H000160000:00001`00017A4U30P000?ooooooooooa`@00><000000000
AP000405000d1@00ADe6:bY0000T000060000000P3l000200000P000P3l000200000P2501@0<0000
000000Q000FL1000T0@0008@`=/200000P0006H4003GcLJJ00000000H19P0`490000011?000100T0
00<c0P000`0L0000000500002@8000001@0000810@0000D000010_ooo`050000;P4H00001@0000/2
000000D0000<0V03H18B00009PH?01X0ooooo`0040000<3oooneoooo818001D3000;00009PH?00`0
CF5dJ5AiL6D00:0020000?X2000@0000000000@0000]0@001@0001@2d04<0`D0000C0]01/@D50000
50;@0J/81@0001<2d06g2PD0000D0]0190h500004`;@0C0@70000?/2`?h00000002@0@0000040P0@
E6U]IG<PCVEg859_KF5^0>S`4P0<f4Mg`61:Ml<]I/<40000;@4100P0000b2R02MQ410000<7TL0000
n`9Po`0000000901000000@2011DJFeULb1>IGLPDVm]HFh0j?0B00cHAgO0H4Yg`beV``@0000]0@80
10000?010@080000<PZ;0/D40@0003=i2000038:a@001040000cNA`0003k0/3n00000000T0400002
108045=iKF9_K01gPal:J01k7`3Xl1803=Q7Ml1PBWO3;FK3100002d10@040000l04200P0000b2R02
Pa010000?GT80000<PXP0PX;0@0002]i2000038:80841P40000]NA`0003k0V3o00000000T0400002
108045=iKF9_K01g7b8:>>1j7`3Xl1803=Q7Ml1PBWO3;FK3100002d10P040000l04100P0000b2/D0
X@l10000:WT80000<P[502P:0@0002Yi2000038:a@0R1@40000ZNA`0003k0/3n00000000T0410000
108045AYKFEc84iUMb1BKfeQKP3Xl1803=Q7Ml1PBWO3;FK3100002d10@040000l04200P0000b2Q/3
V@h20000I7X80000<PYE0Cl>0P0006Ae2000038:80;l30400017M@P0000b2R020P`10000HgD80000
<PXK0b090P0006Aj2000038:E@762080001TM@P0000b2R02m`H20000HdL80000<PXK0h<30P0006Aj
2000038:E@5P1040001eNPP0000b2UD19`<10000I7X80000<PXP0Zd10@0004Uj2000038:80950040
0015NQ`0003k0V3o00000000T0410000108045AYKFEc84iUMb1BKfeQKP3Xl1803=Q7Ml1PBWO3;FK3
100002d10P040000l04100P0000b2W02X``10000LgX80000<PY`0U420@00071j2000038:L08_0@40
001`NPX0000V1Pl02P3ooooo0@00000070000?/240070000002l0P0000010P8RDgUcM6E]0<?3;FK3
000:04D0RP4000000@0000Cc4P040000;@4100@0003`0@800`0000000008@048900001P000024<3K
0@0000<00000000000000000000K@000@00003@0000100000P00000000000000002CA@00F4@30000
o_oo?_nWZSgngiQ4ojNZ?OkoociBEF938@0000P0001R000030000040000Q000020000680000<0000
0@000140000<0000200000X0000@000000000000000900004000060B001P0`003000010000000000
000000/0000@0000a`@00>80000Q0000200002D0000<00001`00P2D0000<00000000P1H0000<0000
000001P0000<0000000001T0000<0000oooo01@0000<00003@000580001`0@000P000:Coool00000
00000000002@0@00000000M002930640K01Y0680LP1Y000000000000000000000000000000000000
00000000000000000000000000000000002`0Me2AgL20000000003HDS@41`XL1M@P000<0000`_2`0
8m[^GoKo0018_2`020000000001<_2`0JFG`GbYVl5mT_2`0OfC`Gg`T602;0hSZT<T/08g>kel001P0
1@000;@m000kWZEf]lk_GkSAQ0G8ImUP2R`QIA00T0400?l600000000003ooooo0000000000000000
00000000000000000@0000000000000000000000002?0P0000000000000000000000000000000000
U@0000000008_B`03Ca7M`Rm;014?4Mg504T2]RSbPOmJ2]g9`bIM6Af00P000009@0000`000020000
600000`0000000004P0000`0000100006@0000`0003oool05P0000`0000H00002P00010000000000
000000T0000@0000H1800603000V0000700000<00000000040000000000000009@0000`000030000
6`000100000<0`00d04003H0000@0000/@D00=01000K000040000:/8003@0@00=P000100002g2P00
d04001/0000@000090h00=01000f00004000030@003@0@00DP00070100040000`?koo`0000000000
00000901000000001`8045@0J@1]06D0L`0P04h0I@1g0200DP1_06d0H@1^00000000000000000000
000000000000000000000000000004MgM;d/0=A3AgL200002R`QIAT0/07M@TMg@lk_G`0060050000
;<d/0<m^l5l0D?l600000:Rk;030H4Yg/II6MhbFAWNX^b`0I04000000000000046/@2P@000000000
WQD:hE2m;002UdIg39M6M`000000000000000000002hdH@5000000X/8FFN5@[Q0P000:Coool00000
00000000002@0@00000000M002930640K01Y0680LP00000000000000000000000000000000000000
00000000000000000000000000000000002`0Me2AgL20000I7H02000000U0000300000@0001D0000
E00008X400130000/`@00:4000010000[0B=?lg=SCmf4@0080800040001<00000000000000000000
oooooooooom@0000<0000:80001B0000L04000D0001Poooo0000000000000000T040000000070P0@
E01Y06d0I@1c0200CP1U07L0801B06l0K@1Q06h00000000000000000000000000000000000000000
000000000000000000000000AgMd_B`0e4=7M`80000:;25U6@2`0Me2AgM3c^mO000H00D0000/cB`0
cfk`G`00003T1000000002X1m@D:;25U00000000000`000000000000000@Ja0:10000000002N5@[Q
D;d/00:GAWL<UdIg000000000000000000000;SAQ0D000002R`QIIhE2^420000Y?ooo`0000000000
00000901000000001d008P00H@1/06T0HP1b00000000000000000000000000000000000000000000
000000000000000000000000001TMP08000002D0000<00001@0002P0000<0000100005@0001D0000
?@4008@0001A0@00/`000040002/18docLf=?lD4002;0P000@0004`00000000000000000003ooooo
oooooe00000c06EUD@0005@0001D00002P4000h0000N0@00?@000040002/18docLf=?`0400350000
0@0004`00000000000000000003ooooooooooe00000c0008D@000580001`0@0010000<3nool00000
00000000002@0@0000000PL2011C07T0K@1R06l0K000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000000000000017MgBm;03D@dMg
0P0000X/8FDI0;01gD97Md?>kel001P000000>@400000000TP7e1@X/8FD00000000003<08FD00000
00000300000000000000011[40X40000000009hE2^4002`00YM6M`bGAWL000000000000000000000
^=641@00000:;25U00000000A010M4X00000000000000?oo003e00000007@00R001Q02QcBP010000
>7A:02Bl;03ZIR]g@7A:0=lB2RH00000K;`/0:6]YGH004@0000006Af00P000009@0000`000040000
:00000`000050000E00005@0001;1000>P0007P4002O00000@000:`4SCo=cHdoPa00020200010000
C00000000000000000000?ooooooooooD00003g`002a0000E00005@0003O0P00>P0000`3002O0000
0@000:`4SCo=cHdo2P/0020200010000C00000000000000000000?ooooooooooD00002_`002a0000
E00005@0002@0@00>P000;d1002O00000@000:`4SCo=cHdo10H0020200010000C000000000000000
00000?ooooooooooD00002g`002a0000DP00070100050000H?ooo`00000000000000090100000002
1`8045<0N@1]0680K`1/000000000000000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000000004MgM;d/0=A3AgL20000
2R`QIAT0/0400000:P000000000:0_D52R`QI@0000000000;O0QI@0000000000<`0QI@0000000000
<00000000000000046/@2P0000000000WQD:h@00;002UdIg39M6M`00000000000000000000000000
00140>1cBP00000000140000ool00?D000000000ool00?D0:7=:0040003HLdX09;`/0>YV:gOPLdX0
@b/:W000001/_2`0XJfUMP00A0000000I7H02000000U0000300000D0000X0000300000@0001D0000
E0000104000:000090@003`000010000[0B=?lg=SCnQ3`00a@000040001<00000000000000000000
oooooooooom@0000:_2<jU40001D0000E0000:@2000:0000^08003`000010000[0B=?lg=SClX2P00
a@000040001<00000000000000000000oooooooooom@0000:_000540001D0000E00005H1000:0000
JP4003`000010000[0B=?lg=SClR1@00a@000040001<00000000000000000000oooooooooom@0000
:_000540001B0000L04000@00030o_oo0000000000000000T040004000070P0@E01Y06d0I@1c0200
CP1U07L0801B06l0K@1Q06h000000000000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000AgMd_B`0e4=7M`00000Z0000
00000882m@D:;25U00000000000Zl25U00000000000]l25U00000000000c025U00000000000`0000
00000000000@Ja0:00000000002N5@[Q000/00:GAWL<UdIg00000000000004@0@7A:0000000004@0
003oo`00m@0004@0003oo`00m@0XLdX00@0003QdBP0T_2`0jVH[Md1dBP3O4PXX000006bl;02Q[JEf
00140000001TMP08000002D0000<0000100002P0000<00001@0005@0001H0000c0<008D0000F1000
h`000040002/18docLf=?iT>000K0`000P0004`00000000000000000003ooooooooooe00001T07X0
XP0007h0001D0000F0000;@3000>00001`@006`000010000[0B=?lg=SClo3P00E@400080001<0000
0000000000000000oooooooooom@0000I01e0:80002Q0000E00005@0001P0`00@`0009`3002Q0000
0@000:`4SCo=cHdoo0`0020200010000C00000000000000000000?ooooooooooD00004L0B@3[0000
E00005@0000O0`00@`0004<3002Q00000@000:`4SCo=cHdo0P`0020200010000C000000000000000
00000?ooooooooooD00006<0002>0000E00005P0001O0P00Q@000:T2003S00000@000:`4SCo=cHdo
80T001/300020000C00000000000000000000?ooooooooooD00006@0NP2R0000OP0005@0001H0000
B08000h0002K0P00K0000040002/18docLf=?lH8001E0@000P0004`00000000000000000003ooooo
oooooe00001T07D0XP000:40001D0000F0000=0100130000<@800:4000010000[0B=?lg=SCog1P00
80800080001<00000000000000000000oooooooooom@0000H`1708h0003[0000E00005P0003Z0000
Q@0003@1003S00000@000:`4SCo=cHdoP`<001/300020000C00000000000000000000?oooooooooo
D00006@0NP2R0000OP0005@0001D00008`4000h0001<0@00K0000040002/18docLf=?f04001E0@00
0@0004`00000000000000000003ooooooooooe00001e0?ooXP0005@0001D0000dP0000h000010@00
K0000040002/18docLf=?bL3001E0@000@0004`00000000000000000003ooooooooooe00001T0?oo
XP0005@0001D0000K`0004<0002C0000X@000040002/18docLf=?jd1000P0P000@0004`000000000
00000000003ooooooooooe0000190000K00005@0001D00004P0004<000160000X@000040002/18do
cLf=?dD0000P0P000@0004`00000000000000000003ooooooooooe0000150000a0000580001`0@00
1@00063oool0000000000000002@0@000@0000L2011D06T0K@1U07<0801>06D0M`0P0580K`1]0640
KP000000000000000000000000000000000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000>@400000000RP?e1@X/8FD00000
000004D0NP000000000002[`8FD00000000002g`8FD00000000003<08FD000000000030000000000
0000011[40X00000000009hE2^4002`000000000A03PLdX000000000A0000?oo003e0000A0000?oo
003e02QcBP010000f7=:02Bl;03ZIR]gh7=:04<[2Yh00000K;`/0:6]YGH004@0000006Af00P00000
9@0000`000050000:00000`000040000E00005@000190`00O@0005P3002/00000@000:`4SCo=cHdo
X``0070200010000C00000000000000000000?ooooooooooD00007<02TXn0000E00005@0002F0000
O@000:h0002/00000@000:`4SCo=cHdoD@80070200010000C00000000000000000000?oooooooooo
D000070000QA0000E00005@0001;0000O@0006<0002/00000@000:`4SCo=cHdo;`40070200010000
C00000000000000000000?ooooooooooD0000700001A0000DP00070100040000400000L000000000
00000;`2000000001`828U<0N@1c07@0I@1]00000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000000000000000000000000000
00000000000000000000000000000000i0@0000000021?D52R`QI@0000000000L00QI@0000000000
A@1j000000000000:_0QI@0000000000;O0QI@0000000000<`0QI@0000000000<000000000000000
46/@2P00000000000014041dBP00000000140000ool00?D000140000ool00?D0:7=:0040000hM4X0
9;`/0>YV:gM0M4X0ga8::P00001/_2`0XJfUMP00A0000000I7H02000000U0000300000@0000X0000
300000D0000U0000300000L0080U000030000000080`0000300000l0080U000030000040001;0000
4000000000050000:00000`000030000:00000`0000400009@0000`0000000209@0000`000070020
8P0000`0003ooooo9@0000`0000=0020:00000`0000200008P0000`0003ooooo8P0000`0003ooooo
AP000700001T0000ADe6:ba0000T000060000:/?QCh000200000P?DBQSknoolnojNZ?BY0000T0000
60000000P3l000200000P000P3l000200000P2501`0<0000000000A0000<0000000004H0000D0000
200004M4BD<300009@0000`0000>00209@0000`0000>0020AP0003@0000X0000ADe6:bY0000T0000
60000000P3l000200000P000P3l000200000P24000080000HP0000`000010000C00006@000000000
00000<H4003Q00000000000000371000hP0002T0ZP00000000000000P3l0000000000000P3l00000
0000000000000000000000000000000000000280000<0000ooooodH0000L0000400004E=AR/2@000
30000000000>000050000000000@00005000
\>"], "Graphics",
        ImageSize->{195, 36},
        ImageMargins->0,
        ImageRegion->{{0., 1.}, {0., 1.}}], ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"u2", "[", "z_", "]"}], "=", 
        FractionBox[
         RowBox[{
          SuperscriptBox["\[ExponentialE]", 
           FractionBox[
            RowBox[{
             SqrtBox["k"], " ", 
             RowBox[{"(", 
              RowBox[{"l", "-", "z"}], ")"}]}], 
            SqrtBox["cs"]]], " ", 
          RowBox[{"HypergeometricU", "[", 
           RowBox[{
            FractionBox[
             RowBox[{"h", "-", 
              FractionBox[
               RowBox[{
                SqrtBox["cs"], " ", "m2", " ", 
                SuperscriptBox["\[Omega]", "2"]}], 
               SqrtBox["k"]]}], 
             RowBox[{"2", " ", "cs", " ", "m"}]], ",", "1", ",", 
            FractionBox[
             RowBox[{"2", " ", 
              SqrtBox["k"], " ", 
              RowBox[{"\[Mu]", "[", "z", "]"}]}], 
             RowBox[{
              SqrtBox["cs"], " ", "m"}]]}], "]"}], " ", 
          RowBox[{"u", "[", "l", "]"}]}], 
         RowBox[{"HypergeometricU", "[", 
          RowBox[{
           FractionBox[
            RowBox[{"h", "-", 
             FractionBox[
              RowBox[{
               SqrtBox["cs"], " ", "m2", " ", 
               SuperscriptBox["\[Omega]", "2"]}], 
              SqrtBox["k"]]}], 
            RowBox[{"2", " ", "cs", " ", "m"}]], ",", "1", ",", 
           FractionBox[
            RowBox[{"2", " ", 
             SqrtBox["k"], " ", 
             RowBox[{"\[Mu]", "[", "l", "]"}]}], 
            RowBox[{
             SqrtBox["cs"], " ", "m"}]]}], "]"}]]}], ";", 
       "\[IndentingNewLine]", "\[IndentingNewLine]", 
       RowBox[{"n1", "=", 
        RowBox[{
         RowBox[{
          RowBox[{"Integrate", "[", 
           RowBox[{
            RowBox[{"2", "*", 
             RowBox[{"\[Mu]", "[", "z", "]"}], "*", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                SubscriptBox["\[PartialD]", "z"], " ", 
                RowBox[{"u", "[", "z", "]"}]}], ")"}], "^", "2"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"z", ",", "0", ",", "l"}], "}"}]}], "]"}], "+", 
          RowBox[{"Integrate", "[", 
           RowBox[{
            RowBox[{"2", "*", 
             RowBox[{"\[Mu]", "[", "z", "]"}], "*", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                SubscriptBox["\[PartialD]", "z"], " ", 
                RowBox[{"u2", "[", "z", "]"}]}], ")"}], "^", "2"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"z", ",", "l", ",", "\[Infinity]"}], "}"}]}], "]"}]}], "//",
          "N"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"n2", "=", 
        RowBox[{
         RowBox[{
          RowBox[{"Integrate", "[", 
           RowBox[{
            RowBox[{"2", "*", "\[Rho]s", "*", 
             RowBox[{"\[Omega]", "^", "2"}], "*", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"u", "[", "z", "]"}], ")"}], "^", "2"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"z", ",", "0", ",", "l"}], "}"}]}], "]"}], "+", 
          RowBox[{"Integrate", "[", 
           RowBox[{
            RowBox[{"2", "*", "\[Rho]s", "*", 
             RowBox[{"\[Omega]", "^", "2"}], "*", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"u2", "[", "z", "]"}], ")"}], "^", "2"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"z", ",", "l", ",", "\[Infinity]"}], "}"}]}], "]"}]}], "//",
          "N"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
       RowBox[{"n3", "=", 
        RowBox[{
         RowBox[{
          RowBox[{"Integrate", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"\[Lambda]", "+", 
               RowBox[{"2", "*", 
                RowBox[{"\[Mu]", "[", "z", "]"}]}]}], ")"}], "*", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"u", "[", "z", "]"}], ")"}], "^", "2"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"z", ",", "0", ",", "l"}], "}"}]}], "]"}], "+", 
          RowBox[{"Integrate", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"\[Lambda]", "+", 
               RowBox[{"2", "*", 
                RowBox[{"\[Mu]", "[", "z", "]"}]}]}], ")"}], "*", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"u2", "[", "z", "]"}], ")"}], "^", "2"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"z", ",", "l", ",", "\[Infinity]"}], "}"}]}], "]"}]}], "//",
          "N"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"n4", "=", 
        RowBox[{"n3", "/", 
         RowBox[{"(", 
          RowBox[{"(", 
           RowBox[{"\[Lambda]\[Mu]", "+", "2"}], ")"}], ")"}]}]}], ";", 
       "\[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[{
     RowBox[{"Umax", "[", "i", "]"}], "=", 
     RowBox[{"u", "[", "0", "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Ih", "[", "i", "]"}], "=", 
     RowBox[{
      RowBox[{"u", "[", "0", "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"P0", "/", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"\[Mu]", "[", "l", "]"}], "*", "2", "*", 
          RowBox[{"(", 
           RowBox[{"1", "+", "\[Nu]"}], ")"}], "*", "l"}], ")"}]}], ")"}]}]}],
     ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"kri", "[", "i", "]"}], "=", "kr"}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Mmax", "[", "i", "]"}], "=", 
     RowBox[{"First", "[", 
      RowBox[{"FindMaximum", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Ep", "*", "Ip", "*", 
           RowBox[{
            RowBox[{"u", "''"}], "[", "z", "]"}]}], ",", 
          RowBox[{"0", "\[LessEqual]", "z", "\[LessEqual]", "l"}]}], "}"}], 
        ",", " ", 
        RowBox[{"{", 
         RowBox[{"z", ",", 
          RowBox[{"l", "/", "2"}]}], "}"}]}], "]"}], "]"}]}]}]}], 
  "\[IndentingNewLine]", "]"}]], "Input",
 CellChangeTimes->{{3.4810538072042*^9, 3.4810538450186*^9}, {
   3.4810539294146*^9, 3.4810539417074003`*^9}, 3.4810954480842*^9, {
   3.4810957540782003`*^9, 3.4810957828914003`*^9}, {3.4811123328988*^9, 
   3.4811124917358*^9}, {3.4811125217918*^9, 3.4811125427208*^9}, {
   3.4811125738608*^9, 3.4811126166808*^9}, {3.4811126675418*^9, 
   3.4811127012218*^9}, {3.4811128237668*^9, 3.4811128497438*^9}, {
   3.4811130397808*^9, 3.4811130484958*^9}, {3.4811131090358*^9, 
   3.4811131164358*^9}, 3.4811132051558*^9, {3.4811132822658*^9, 
   3.4811132856558*^9}, {3.4811133824018*^9, 3.4811133828898*^9}, 
   3.4811134408458*^9, 3.4811134720538*^9, {3.4811135466328*^9, 
   3.4811135554748*^9}, {3.4811136439077997`*^9, 3.4811136492918*^9}, {
   3.4811137201228*^9, 3.4811137486768*^9}, {3.4811138131948*^9, 
   3.4811138141818*^9}, {3.4811140508708*^9, 3.4811140637268*^9}, {
   3.4811141306998*^9, 3.4811141320428*^9}, {3.4811141897368*^9, 
   3.4811142342108*^9}, {3.4811145485568*^9, 3.4811145559498*^9}, 
   3.4811147149538*^9, {3.481138140817*^9, 3.4811381556682*^9}, {
   3.481138346508*^9, 3.481138348586*^9}, {3.481138517092*^9, 
   3.481138519428*^9}, {3.481138728112*^9, 3.481138730278*^9}, {
   3.481138786105*^9, 3.481138796848*^9}, {3.481138915725*^9, 
   3.481138928488*^9}, {3.4811389792139997`*^9, 3.481138979699*^9}, 
   3.481139062439*^9, 3.481139104743*^9, {3.481139138426*^9, 
   3.481139154164*^9}, {3.481139376148*^9, 3.481139429175*^9}, {
   3.4811395052130003`*^9, 3.481139507035*^9}, {3.4811396348780003`*^9, 
   3.4811396444969997`*^9}, 3.481139696078*^9, {3.48114008244*^9, 
   3.481140109597*^9}, {3.481140250592*^9, 3.48114025429*^9}, {
   3.481140322001*^9, 3.481140325299*^9}, {3.4811857600109997`*^9, 
   3.481185798824*^9}, {3.612776475006717*^9, 3.612776478276722*^9}, {
   3.612847946337908*^9, 3.612847950767914*^9}, {3.612848066346081*^9, 
   3.6128481241961737`*^9}}]
},
WindowSize->{1424, 765},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> True}},
Magnification->0.8999999761581421,
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, 40010, 845, 1744, "Input"]
}
]
*)

(* End of internal cache information *)