(* Content-type: application/vnd.wolfram.mathematica *)

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

(* CreatedBy='Mathematica 12.1' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     65819,       1269]
NotebookOptionsPosition[     63575,       1226]
NotebookOutlinePosition[     64017,       1243]
CellTagsIndexPosition[     63974,       1240]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Discrete Fourier series", "Section",
 CellChangeTimes->{{3.8250471942961464`*^9, 3.8250472558402195`*^9}, {
  3.825047555742766*^9, 3.825047559409715*^9}, {3.8277096256331263`*^9, 
  3.827709632700748*^9}, {3.827709882245117*^9, 
  3.8277098910601287`*^9}},ExpressionUUID->"86686e14-fc55-42dc-b305-\
e2ed2e4b1a8a"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"n", "=", "200"}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", "list", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"a", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"RandomReal", "[", "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", 
         RowBox[{"n", "-", "1"}]}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.8226211276396384`*^9, 3.82262129486576*^9}, {
   3.822621349911625*^9, 3.8226214833366013`*^9}, {3.822621517147281*^9, 
   3.8226215739258113`*^9}, {3.822621608740857*^9, 3.8226216098937397`*^9}, {
   3.8226216777174053`*^9, 3.8226217894554415`*^9}, {3.8226218227469807`*^9, 
   3.822621875000436*^9}, {3.8226219103329124`*^9, 3.8226219163408504`*^9}, {
   3.8226219542984657`*^9, 3.822622002086941*^9}, {3.822622080352575*^9, 
   3.8226221357673388`*^9}, {3.822645735106282*^9, 3.8226458032017403`*^9}, {
   3.822645834813614*^9, 3.822645837799676*^9}, {3.822645871317651*^9, 
   3.8226459081605883`*^9}, {3.8226459548207536`*^9, 3.822646250191146*^9}, {
   3.8226462896267233`*^9, 3.8226464318411064`*^9}, {3.822646512967785*^9, 
   3.822646516224951*^9}, {3.8226465540576468`*^9, 3.822646564119341*^9}, {
   3.8277098965408745`*^9, 3.8277100592742467`*^9}, {3.8280148398310556`*^9, 
   3.8280148402325363`*^9}, 3.8280149014680333`*^9, {3.8280149529214754`*^9, 
   3.8280150430320406`*^9}, {3.828015120161235*^9, 3.828015158981262*^9}, {
   3.8280151898395896`*^9, 3.8280152089815826`*^9}, {3.8280152398646574`*^9, 
   3.828015354071247*^9}, {3.828015408026306*^9, 3.8280154336339607`*^9}, {
   3.828015472106393*^9, 3.828015472214861*^9}, {3.828015518422146*^9, 
   3.828015669707283*^9}, {3.828015712145762*^9, 3.828015764346391*^9}, {
   3.828015803360965*^9, 3.8280158720273438`*^9}, {3.828015907902726*^9, 
   3.828016019213128*^9}, {3.828148911703225*^9, 3.8281489939384456`*^9}, {
   3.8281493164790792`*^9, 3.8281493205014973`*^9}},
 CellLabel->"In[14]:=",ExpressionUUID->"9174fd4b-30aa-441b-9412-80faf601d72b"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", "modes", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Grid", "[", 
   RowBox[{"Partition", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"ListPlot", "[", 
        RowBox[{
         RowBox[{"Table", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"i", ",", 
             RowBox[{"Sin", "[", 
              RowBox[{"j", " ", "Pi", " ", 
               RowBox[{"i", "/", "n"}]}], "]"}]}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}], ",", 
         RowBox[{"Joined", "\[Rule]", "True"}], ",", 
         RowBox[{"Frame", "\[Rule]", "True"}]}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", 
         RowBox[{"n", "-", "1"}]}], "}"}]}], "]"}], ",", "3"}], "]"}], 
   "]"}]}]], "Input",
 CellLabel->"In[8]:=",ExpressionUUID->"ea83f712-7b09-46d5-9060-e8b8fac5cbd1"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Fourier", " ", "transform"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"\[Alpha]", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"Sum", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Sin", "[", 
            RowBox[{"j", " ", "Pi", " ", 
             RowBox[{"i", "/", "n"}]}], "]"}], 
           RowBox[{"a", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], ",", 
          RowBox[{"{", 
           RowBox[{"i", ",", "1", ",", 
            RowBox[{"n", "-", "1"}]}], "}"}]}], "]"}], "*", 
        RowBox[{"2", "/", "n"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", 
         RowBox[{"n", "-", "1"}]}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.8226211276396384`*^9, 3.82262129486576*^9}, {
   3.822621349911625*^9, 3.8226214833366013`*^9}, {3.822621517147281*^9, 
   3.8226215739258113`*^9}, {3.822621608740857*^9, 3.8226216098937397`*^9}, {
   3.8226216777174053`*^9, 3.8226217894554415`*^9}, {3.8226218227469807`*^9, 
   3.822621875000436*^9}, {3.8226219103329124`*^9, 3.8226219163408504`*^9}, {
   3.8226219542984657`*^9, 3.822622002086941*^9}, {3.822622080352575*^9, 
   3.8226221357673388`*^9}, {3.822645735106282*^9, 3.8226458032017403`*^9}, {
   3.822645834813614*^9, 3.822645837799676*^9}, {3.822645871317651*^9, 
   3.8226459081605883`*^9}, {3.8226459548207536`*^9, 3.822646250191146*^9}, {
   3.8226462896267233`*^9, 3.8226464318411064`*^9}, {3.822646512967785*^9, 
   3.822646516224951*^9}, {3.8226465540576468`*^9, 3.822646564119341*^9}, {
   3.8277098965408745`*^9, 3.8277100592742467`*^9}, {3.8280148398310556`*^9, 
   3.8280148402325363`*^9}, 3.8280149014680333`*^9, {3.8280149529214754`*^9, 
   3.8280150430320406`*^9}, {3.828015120161235*^9, 3.828015158981262*^9}, {
   3.8280151898395896`*^9, 3.8280152089815826`*^9}, {3.8280152398646574`*^9, 
   3.828015354071247*^9}, {3.828015408026306*^9, 3.8280154336339607`*^9}, {
   3.828015472106393*^9, 3.828015472214861*^9}, {3.828015518422146*^9, 
   3.828015669707283*^9}, {3.828015712145762*^9, 3.828015764346391*^9}, {
   3.828015803360965*^9, 3.8280158720273438`*^9}, {3.828015907902726*^9, 
   3.828016019213128*^9}, {3.828148911703225*^9, 3.828148946322839*^9}, {
   3.82814902283848*^9, 3.82814903042307*^9}, 3.828149344617317*^9},
 CellLabel->"In[16]:=",ExpressionUUID->"f6db8324-9aef-4d10-9770-8d2f6a64e224"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"list", " ", "again"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"aa", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Sin", "[", 
           RowBox[{"j", " ", "Pi", " ", 
            RowBox[{"i", "/", "n"}]}], "]"}], 
          RowBox[{"\[Alpha]", "[", 
           RowBox[{"[", "j", "]"}], "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "1", ",", 
           RowBox[{"n", "-", "1"}]}], "}"}]}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", 
         RowBox[{"n", "-", "1"}]}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.8226211276396384`*^9, 3.82262129486576*^9}, {
   3.822621349911625*^9, 3.8226214833366013`*^9}, {3.822621517147281*^9, 
   3.8226215739258113`*^9}, {3.822621608740857*^9, 3.8226216098937397`*^9}, {
   3.8226216777174053`*^9, 3.8226217894554415`*^9}, {3.8226218227469807`*^9, 
   3.822621875000436*^9}, {3.8226219103329124`*^9, 3.8226219163408504`*^9}, {
   3.8226219542984657`*^9, 3.822622002086941*^9}, {3.822622080352575*^9, 
   3.8226221357673388`*^9}, {3.822645735106282*^9, 3.8226458032017403`*^9}, {
   3.822645834813614*^9, 3.822645837799676*^9}, {3.822645871317651*^9, 
   3.8226459081605883`*^9}, {3.8226459548207536`*^9, 3.822646250191146*^9}, {
   3.8226462896267233`*^9, 3.8226464318411064`*^9}, {3.822646512967785*^9, 
   3.822646516224951*^9}, {3.8226465540576468`*^9, 3.822646564119341*^9}, {
   3.8277098965408745`*^9, 3.8277100592742467`*^9}, {3.8280148398310556`*^9, 
   3.8280148402325363`*^9}, 3.8280149014680333`*^9, {3.8280149529214754`*^9, 
   3.8280150430320406`*^9}, {3.828015120161235*^9, 3.828015158981262*^9}, {
   3.8280151898395896`*^9, 3.8280152089815826`*^9}, {3.8280152398646574`*^9, 
   3.828015354071247*^9}, {3.828015408026306*^9, 3.8280154336339607`*^9}, {
   3.828015472106393*^9, 3.828015472214861*^9}, {3.828015518422146*^9, 
   3.828015669707283*^9}, {3.828015712145762*^9, 3.828015764346391*^9}, {
   3.828015803360965*^9, 3.8280158720273438`*^9}, {3.828015907902726*^9, 
   3.828016019213128*^9}, {3.828148911703225*^9, 3.828148946322839*^9}, {
   3.82814902283848*^9, 3.82814903042307*^9}, {3.828149145518894*^9, 
   3.8281491825156827`*^9}, 3.8281493523427954`*^9},
 CellLabel->"In[17]:=",ExpressionUUID->"57ebc5c8-0469-4298-a7a8-c6cc23150f96"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", "comparison", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"gr1", "=", 
     RowBox[{"ListPlot", "[", 
      RowBox[{"a", ",", 
       RowBox[{"Frame", "\[Rule]", "True"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"gr2", "=", 
     RowBox[{"ListPlot", "[", 
      RowBox[{"aa", ",", 
       RowBox[{"Frame", "\[Rule]", "True"}], ",", 
       RowBox[{"PlotStyle", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{"Red", ",", 
          RowBox[{"PointSize", "[", "Large", "]"}]}], "}"}]}]}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Show", "[", 
    RowBox[{"gr2", ",", "gr1"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.8226211276396384`*^9, 3.82262129486576*^9}, {
   3.822621349911625*^9, 3.8226214833366013`*^9}, {3.822621517147281*^9, 
   3.8226215739258113`*^9}, {3.822621608740857*^9, 3.8226216098937397`*^9}, {
   3.8226216777174053`*^9, 3.8226217894554415`*^9}, {3.8226218227469807`*^9, 
   3.822621875000436*^9}, {3.8226219103329124`*^9, 3.8226219163408504`*^9}, {
   3.8226219542984657`*^9, 3.822622002086941*^9}, {3.822622080352575*^9, 
   3.8226221357673388`*^9}, {3.822645735106282*^9, 3.8226458032017403`*^9}, {
   3.822645834813614*^9, 3.822645837799676*^9}, {3.822645871317651*^9, 
   3.8226459081605883`*^9}, {3.8226459548207536`*^9, 3.822646250191146*^9}, {
   3.8226462896267233`*^9, 3.8226464318411064`*^9}, {3.822646512967785*^9, 
   3.822646516224951*^9}, {3.8226465540576468`*^9, 3.822646564119341*^9}, {
   3.8277098965408745`*^9, 3.8277100592742467`*^9}, {3.8280148398310556`*^9, 
   3.8280148402325363`*^9}, 3.8280149014680333`*^9, {3.8280149529214754`*^9, 
   3.8280150430320406`*^9}, {3.828015120161235*^9, 3.828015158981262*^9}, {
   3.8280151898395896`*^9, 3.8280152089815826`*^9}, {3.8280152398646574`*^9, 
   3.828015354071247*^9}, {3.828015408026306*^9, 3.8280154336339607`*^9}, {
   3.828015472106393*^9, 3.828015472214861*^9}, {3.828015518422146*^9, 
   3.828015669707283*^9}, {3.828015712145762*^9, 3.828015764346391*^9}, {
   3.828015803360965*^9, 3.8280158720273438`*^9}, {3.828015907902726*^9, 
   3.828016019213128*^9}, {3.828148911703225*^9, 3.828148946322839*^9}, {
   3.82814902283848*^9, 3.82814903042307*^9}, {3.828149145518894*^9, 
   3.8281491825156827`*^9}},
 CellLabel->"In[18]:=",ExpressionUUID->"97108cef-50e6-4fa2-a79d-60e9f1ea42a0"]
}, Closed]],

Cell[CellGroupData[{

Cell["Modelling assignment", "Section",
 CellChangeTimes->{
  3.4430806364968*^9, {3.4453166713873997`*^9, 3.4453166747414*^9}, {
   3.447462264275*^9, 3.4474622644154*^9}, {3.4695731321116*^9, 
   3.4695731359491997`*^9}, {3.473055192524*^9, 3.473055199133*^9}, {
   3.47305542959*^9, 3.4730554407939997`*^9}, {3.4734696797594*^9, 
   3.4734696967634*^9}, {3.4933064685896*^9, 3.493306471148*^9}, {
   3.5568725291398897`*^9, 3.5568725438967333`*^9}, 3.5568726929512587`*^9, {
   3.5576544837556367`*^9, 3.5576544886852455`*^9}, {3.6852489036101027`*^9, 
   3.68524892333013*^9}, {3.688455478402707*^9, 3.688455486552718*^9}, {
   3.732332061536833*^9, 3.7323320661099243`*^9}, {3.732332251065756*^9, 
   3.732332259663978*^9}, {3.7328755230583534`*^9, 3.7328755307151165`*^9}, {
   3.7648177400838947`*^9, 3.764817746657362*^9}, {3.794895625047804*^9, 
   3.794895629727687*^9}, {3.794895811562347*^9, 3.7948958638747997`*^9}, 
   3.7954055648226833`*^9, {3.7958566176744843`*^9, 3.7958566190993643`*^9}, {
   3.827796552688429*^9, 3.827796552694504*^9}, {3.827799413918418*^9, 
   3.8277994180855923`*^9}},ExpressionUUID->"ddf24a73-a884-42d6-ad85-\
1c9ed46c5997"],

Cell[CellGroupData[{

Cell["Properties [kg,m,s]", "Subsection",
 CellChangeTimes->{{3.473055512584*^9, 3.4730555529630003`*^9}, {
   3.4730575328269997`*^9, 3.473057539222*^9}, {3.4734691424642*^9, 
   3.4734691605758*^9}, {3.4736023142902*^9, 3.4736023273162003`*^9}, {
   3.568715518105007*^9, 3.568715523097016*^9}, 3.5994750893663483`*^9, 
   3.6297916561604834`*^9, {3.63073778491436*^9, 3.630737793931176*^9}, {
   3.6624537703648133`*^9, 3.6624537732749796`*^9}, 3.685248933840145*^9, {
   3.685267058363743*^9, 3.6852671273738394`*^9}, {3.685323461463516*^9, 
   3.6853234663335233`*^9}, {3.6858764211836576`*^9, 3.6858764235236607`*^9}, 
   3.6859482147140565`*^9, {3.686226212775403*^9, 3.6862262151354065`*^9}, {
   3.686466369367293*^9, 3.686466390117322*^9}, {3.732333324403823*^9, 
   3.7323333512571926`*^9}, {3.7330531980857735`*^9, 3.733053200539534*^9}, {
   3.765153863866486*^9, 
   3.76515386489231*^9}},ExpressionUUID->"0f0487a8-da1d-494b-9d29-\
0a7939035da4"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"val", "=", 
    RowBox[{"{", "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{"COLUMN", " ", "DATA"}], "*)"}], "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Es", "\[Rule]", 
       RowBox[{"210", "*", 
        RowBox[{"10", "^", "9"}]}]}], ",", 
      RowBox[{"(*", 
       RowBox[{
        RowBox[{"Young", "'"}], "s", " ", "modulus", " ", "of", " ", 
        "steel"}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"\[Nu]s", "\[Rule]", "0.3"}], ",", "\[IndentingNewLine]", 
      RowBox[{"Gs", "\[Rule]", 
       RowBox[{"Es", "/", 
        RowBox[{"(", 
         RowBox[{"2", "+", 
          RowBox[{"2", "*", "\[Nu]s"}]}], ")"}]}]}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"\[Rho]s", "\[Rule]", "7800"}], ",", " ", 
      RowBox[{"(*", 
       RowBox[{"column", " ", "density"}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"dmaj", "->", "0.006"}], ",", " ", 
      RowBox[{"(*", 
       RowBox[{"nominal", " ", "diamater"}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"P", "\[Rule]", "0.001"}], ",", " ", 
      RowBox[{"(*", "pitch", "*)"}], "\[IndentingNewLine]", 
      RowBox[{"d", "\[Rule]", 
       RowBox[{"dmaj", "-", 
        RowBox[{
         RowBox[{"2", "/", 
          RowBox[{"Sqrt", "[", "3", "]"}]}], " ", "P"}]}]}], ",", " ", 
      RowBox[{"(*", 
       RowBox[{"effective", " ", "diameter"}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"h", "->", "0.156"}], ",", " ", 
      RowBox[{"(*", 
       RowBox[{"floor", "-", "to", "-", 
        RowBox[{"floor", " ", "distance"}]}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"lx", "\[Rule]", "0.410"}], ",", " ", 
      RowBox[{"(*", 
       RowBox[{"column", " ", "spacing"}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"ly", "\[Rule]", "0.243"}], ",", " ", 
      RowBox[{"(*", 
       RowBox[{"column", " ", "spacing"}], "*)"}], "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{"(*", 
       RowBox[{"FLOOR", " ", "DATA"}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"Ea", "\[Rule]", 
       RowBox[{"70", "*", 
        RowBox[{"10", "^", "9"}]}]}], ",", 
      RowBox[{"(*", 
       RowBox[{
        RowBox[{"Young", "'"}], "s", " ", "modulus", " ", "of", " ", 
        "aluminum"}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"\[Nu]a", "\[Rule]", "0.3"}], ",", "\[IndentingNewLine]", 
      RowBox[{"Ga", "\[Rule]", 
       RowBox[{"Ea", "/", 
        RowBox[{"(", 
         RowBox[{"2", "+", 
          RowBox[{"2", "*", "\[Nu]a"}]}], ")"}]}]}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"\[Rho]a", "\[Rule]", "2700"}], ",", " ", 
      RowBox[{"(*", 
       RowBox[{"floor", " ", "density"}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"L", "\[Rule]", "0.440"}], ",", " ", 
      RowBox[{"(*", 
       RowBox[{"floor", " ", "length"}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"W", "\[Rule]", "0.295"}], ",", " ", 
      RowBox[{"(*", 
       RowBox[{"floor", " ", "width"}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"t", "\[Rule]", "0.015"}], ",", " ", 
      RowBox[{"(*", 
       RowBox[{"floor", " ", "thickness"}], "*)"}], "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{"(*", 
       RowBox[{"PLASTIC", " ", "STRIP", " ", "DATA"}], "*)"}], 
      "\[IndentingNewLine]", 
      RowBox[{"Ep", "\[Rule]", 
       RowBox[{"3", "*", 
        RowBox[{"10", "^", "9"}]}]}], ",", "\[IndentingNewLine]", 
      RowBox[{"\[Nu]p", "\[Rule]", "0.3"}], ",", "\[IndentingNewLine]", 
      RowBox[{"Gp", "\[Rule]", 
       RowBox[{"Ep", "/", 
        RowBox[{"(", 
         RowBox[{"2", "+", 
          RowBox[{"2", "*", "\[Nu]p"}]}], ")"}]}]}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"\[Rho]p", "\[Rule]", "1000"}], ",", "\[IndentingNewLine]", 
      RowBox[{"w", "\[Rule]", "0.05"}], ",", " ", 
      RowBox[{"(*", 
       RowBox[{"strip", " ", "width"}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"s", "\[Rule]", "0.002"}]}], " ", 
     RowBox[{"(*", 
      RowBox[{"strip", " ", "thickness"}], "*)"}], "\[IndentingNewLine]", 
     "}"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"Grid", "[", 
  RowBox[{
   RowBox[{"Partition", "[", 
    RowBox[{
     RowBox[{"Flatten", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"val", " ", "/.", 
         RowBox[{"Rule", "\[Rule]", "List"}]}], ")"}], " ", "/.", " ", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"a_", ",", "b_"}], "}"}], "\[Rule]", 
        RowBox[{"{", 
         RowBox[{"a", ",", "b"}], "}"}]}]}], "]"}], ",", "4"}], "]"}], ",", 
   RowBox[{"Alignment", "\[Rule]", "Left"}], ",", 
   RowBox[{"Spacings", "->", 
    RowBox[{"{", 
     RowBox[{"5", ",", "1.2"}], "}"}]}], ",", 
   RowBox[{"Frame", "\[Rule]", "All"}]}], 
  "]"}], "\[IndentingNewLine]"}], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.732963786982458*^9, 3.732963800329712*^9}, {
   3.7329638807173033`*^9, 3.7329639838406887`*^9}, {3.7329640219057136`*^9, 
   3.7329640679494686`*^9}, {3.7329641106095915`*^9, 
   3.7329641238605337`*^9}, {3.7329641586934195`*^9, 
   3.7329641775638523`*^9}, {3.7329642131292753`*^9, 
   3.7329642468139324`*^9}, {3.732964291650645*^9, 3.732964295266472*^9}, {
   3.7329643646012278`*^9, 3.7329645915297546`*^9}, {3.732964625547305*^9, 
   3.732964633879656*^9}, {3.732964682913621*^9, 3.732964776488607*^9}, {
   3.732964812880664*^9, 3.7329648406399565`*^9}, {3.7329649084689856`*^9, 
   3.732965030133172*^9}, {3.732965104873771*^9, 3.7329651443322325`*^9}, {
   3.7329651781396537`*^9, 3.73296520770864*^9}, {3.732965254213273*^9, 
   3.7329652993251047`*^9}, {3.7329653860024896`*^9, 
   3.7329653897065096`*^9}, {3.7329654581300693`*^9, 
   3.7329654674851675`*^9}, {3.7329655079892044`*^9, 3.7329655085784197`*^9}, 
   3.733053175193036*^9, 3.733053206040131*^9, {3.733053281782422*^9, 
   3.7330534732836246`*^9}, {3.733053546072621*^9, 3.7330536355505424`*^9}, {
   3.73305367336732*^9, 3.7330538780699167`*^9}, {3.7330539174964027`*^9, 
   3.7330539761170115`*^9}, {3.7330540174495764`*^9, 
   3.7330540187778425`*^9}, {3.733054221830557*^9, 3.7330542834776363`*^9}, {
   3.7330543376240883`*^9, 3.733054425945921*^9}, {3.7330545145491314`*^9, 
   3.7330545149397078`*^9}, 3.733064445719219*^9, {3.733104688989295*^9, 
   3.7331047545189743`*^9}, {3.7331047887892323`*^9, 
   3.7331049094542694`*^9}, {3.733104944739091*^9, 3.7331049683212028`*^9}, 
   3.7331053896852155`*^9, {3.7331127009548764`*^9, 3.7331127150816717`*^9}, {
   3.733112994611229*^9, 3.7331129958613586`*^9}, {3.733113049569934*^9, 
   3.7331130498667803`*^9}, 3.7331131038881173`*^9, {3.733114899963835*^9, 
   3.733114905386201*^9}, {3.73311625140981*^9, 3.7331162535818663`*^9}, {
   3.73311631318167*^9, 3.733116407363805*^9}, {3.73314636774393*^9, 
   3.7331464169053802`*^9}, 3.7334915575756106`*^9, {3.7334916007179976`*^9, 
   3.7334916015265217`*^9}, 3.733561226688825*^9, {3.7356177526640673`*^9, 
   3.735617761193325*^9}, {3.765112107920904*^9, 3.76511210935957*^9}, {
   3.7651125577057576`*^9, 3.765112628809655*^9}, {3.765152702580862*^9, 
   3.7651528573769875`*^9}, {3.765152921902412*^9, 3.7651529651322827`*^9}, {
   3.7651530067056184`*^9, 3.7651530120276155`*^9}, {3.7651531691641192`*^9, 
   3.7651532270683246`*^9}, {3.7651538681513586`*^9, 
   3.7651539029602003`*^9}, {3.765153944420312*^9, 3.7651539453412943`*^9}, 
   3.765154104220358*^9, {3.765154372827244*^9, 3.765154373048473*^9}, {
   3.765157812723898*^9, 3.765157919729148*^9}, {3.76515796243433*^9, 
   3.765158001596624*^9}, {3.7961782968491864`*^9, 3.796178300896781*^9}, 
   3.7961783600757833`*^9, 3.7961787820454392`*^9, 3.7961794744582567`*^9, {
   3.7962025926308784`*^9, 3.79620259857528*^9}, {3.796346111099641*^9, 
   3.7963461216727734`*^9}, {3.79637298723057*^9, 3.7963729928911247`*^9}, 
   3.79637325914301*^9, {3.7963733017244835`*^9, 3.7963733081905456`*^9}, 
   3.7963733519264555`*^9, {3.7963734957716026`*^9, 3.796373496521605*^9}, {
   3.7964398962391057`*^9, 3.796439896797656*^9}, {3.796783639334712*^9, 
   3.796783652693286*^9}, 3.796783795918022*^9, {3.7969852499549294`*^9, 
   3.7969852616123934`*^9}, {3.7969853925853977`*^9, 3.7969853963543725`*^9}, 
   3.797755086465664*^9, {3.7977575362790594`*^9, 3.797757632579767*^9}, {
   3.7977577228266277`*^9, 3.797757723459545*^9}, {3.7979038267185326`*^9, 
   3.7979038282452393`*^9}, {3.797903962534937*^9, 3.7979040204947233`*^9}, {
   3.797904762968611*^9, 3.797904766733965*^9}, {3.827797382692029*^9, 
   3.8277973833410892`*^9}, 3.827797419388828*^9, 3.82779744998677*^9},
 CellLabel->"In[1]:=",ExpressionUUID->"ce80bae6-fb13-49c7-8dbf-a3af0c00446e"]
}, Closed]],

Cell[CellGroupData[{

Cell["Acceleration signal analysis", "Subsection",
 CellChangeTimes->{{3.473055485005*^9, 3.473055488513*^9}, {3.473057544493*^9,
    3.473057567259*^9}, {3.4734691404674*^9, 3.4734691826186*^9}, {
   3.5674322230551395`*^9, 3.567432229575951*^9}, {3.5674885568810406`*^9, 
   3.5674885586126437`*^9}, 3.5986036084199142`*^9, {3.598603673712649*^9, 
   3.5986036738476562`*^9}, {3.6297775308006935`*^9, 3.629777534217099*^9}, {
   3.6611749978437643`*^9, 3.6611750017959766`*^9}, 3.68524888736008*^9, 
   3.685281845213238*^9, {3.6853234458234944`*^9, 3.685323452933504*^9}, {
   3.6858763883636117`*^9, 3.6858763984036255`*^9}, {3.6862356299073353`*^9, 
   3.6862356300473356`*^9}, {3.6862902438578215`*^9, 3.686290263787849*^9}, {
   3.6862903725280013`*^9, 3.686290373378002*^9}, {3.686370549159279*^9, 
   3.6863705527492843`*^9}, {3.6863950336343946`*^9, 
   3.6863950375944004`*^9}, {3.688458622114133*^9, 3.688458623034134*^9}, {
   3.689506748056718*^9, 3.689506831506835*^9}, {3.690086538897827*^9, 
   3.6900865453778358`*^9}, {3.694943494189322*^9, 3.6949435129193482`*^9}, {
   3.694946432132904*^9, 3.6949465072730093`*^9}, {3.6968528976354847`*^9, 
   3.696852899745488*^9}, {3.696908842156185*^9, 3.6969088815462403`*^9}, {
   3.696908925666302*^9, 3.6969089286863065`*^9}, {3.6969090106584225`*^9, 
   3.6969090154184294`*^9}, {3.7268857155101137`*^9, 
   3.7268857469885635`*^9}, {3.7966133174134192`*^9, 3.796613318932433*^9}, {
   3.796613860802449*^9, 
   3.7966138622196627`*^9}},ExpressionUUID->"0ee45dc8-cf7c-473a-ae94-\
685c39f63702"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"acceleration", " ", "time", " ", "series"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"SetDirectory", "[", 
     RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"series", "=", 
     RowBox[{
      RowBox[{"Import", "[", "\"\<experiment_5.xlsx\>\"", "]"}], "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"g1", "=", 
    RowBox[{"ListPlot", "[", 
     RowBox[{"series", ",", 
      RowBox[{"Joined", "\[Rule]", "True"}], ",", 
      RowBox[{"PlotRange", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"0.67", ",", "1.67"}], "}"}], ",", "All"}], "}"}]}], ",", 
      RowBox[{"Frame", "\[Rule]", "True"}], ",", 
      RowBox[{"FrameLabel", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{
        "\"\<t [s]\>\"", ",", 
         "\"\<a [m/\!\(\*SuperscriptBox[\(s\), \(2\)]\)\>\""}], "}"}]}]}], 
     "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.6895068568568707`*^9, 3.6895068606168756`*^9}, {
   3.689507096377206*^9, 3.689507101947213*^9}, {3.689508864916235*^9, 
   3.6895088669362373`*^9}, 3.689509046516489*^9, 3.6895093142295237`*^9, {
   3.689509404105324*^9, 3.6895094295353594`*^9}, {3.6895120381904573`*^9, 
   3.6895120397004595`*^9}, {3.689512122370575*^9, 3.689512157960625*^9}, {
   3.689512282330799*^9, 3.6895123101408377`*^9}, {3.6895131735720468`*^9, 
   3.689513295222217*^9}, {3.6895133805523367`*^9, 3.689513462322451*^9}, {
   3.689513526522541*^9, 3.689513592872634*^9}, {3.6895147961043186`*^9, 
   3.689514800934325*^9}, 3.689585657689865*^9, {3.689585919908863*^9, 
   3.6895859200938735`*^9}, 3.6896675842177234`*^9, {3.6900856637946014`*^9, 
   3.6900858381548457`*^9}, {3.6900859193649592`*^9, 
   3.6900860779651814`*^9}, {3.6900895242580185`*^9, 
   3.6900895290180254`*^9}, {3.690094654552212*^9, 3.6900947412723336`*^9}, {
   3.690094861452502*^9, 3.690094863612505*^9}, {3.6900949089725685`*^9, 
   3.690095000202696*^9}, {3.690172427710157*^9, 3.6901724615802045`*^9}, {
   3.690172499810258*^9, 3.690172554560334*^9}, {3.6901725949403915`*^9, 
   3.690172619990426*^9}, {3.690172658880481*^9, 3.690172733630585*^9}, {
   3.6901729886409426`*^9, 3.690173006360967*^9}, {3.690173049411027*^9, 
   3.690173089801084*^9}, 3.690173225581274*^9, 3.6901732709113374`*^9, {
   3.6901733763614855`*^9, 3.69017342987156*^9}, 3.6901736413718557`*^9, {
   3.6901736732519007`*^9, 3.690173747172004*^9}, {3.6901738114620943`*^9, 
   3.690173812002095*^9}, {3.690173942672278*^9, 3.6901739433122787`*^9}, {
   3.6947797562079945`*^9, 3.6947797740580196`*^9}, {3.694779859108139*^9, 
   3.6947799021581993`*^9}, {3.6948765488222756`*^9, 
   3.6948767582925687`*^9}, {3.6948771329030933`*^9, 3.694877410933482*^9}, {
   3.694877449223536*^9, 3.694877469123564*^9}, {3.6949175866210814`*^9, 
   3.694917591951089*^9}, {3.6949176842312183`*^9, 3.694917684771219*^9}, {
   3.6949177372612925`*^9, 3.6949178831114964`*^9}, {3.6949179676316147`*^9, 
   3.694917972901622*^9}, {3.694918006341669*^9, 3.694918019701688*^9}, {
   3.694918057981741*^9, 3.6949181203618283`*^9}, {3.6949181744819045`*^9, 
   3.6949181754019055`*^9}, {3.6949182156519623`*^9, 
   3.6949182234419727`*^9}, {3.694918335022129*^9, 3.6949183942422123`*^9}, {
   3.6949184418922787`*^9, 3.6949184896423454`*^9}, 3.6949185895924854`*^9, {
   3.694918662322587*^9, 3.6949186658425922`*^9}, {3.6949188708028793`*^9, 
   3.6949189242029543`*^9}, {3.6949189605230045`*^9, 3.694919119233227*^9}, 
   3.6949191503732705`*^9, {3.69491919269333*^9, 3.694919202403343*^9}, {
   3.694919236243391*^9, 3.694919363013568*^9}, {3.694919403093624*^9, 
   3.6949194079936314`*^9}, {3.6949194522036934`*^9, 3.6949195440338216`*^9}, 
   3.694919577063868*^9, {3.6949196167139235`*^9, 3.694919636613951*^9}, 
   3.694919702364043*^9, {3.694919742224099*^9, 3.6949198399642363`*^9}, {
   3.6949210692072573`*^9, 3.694921117877326*^9}, 3.6949211551773777`*^9, 
   3.694921253767516*^9, {3.694921490017847*^9, 3.69492150647787*^9}, {
   3.696689244189213*^9, 3.696689269359248*^9}, {3.696689311059306*^9, 
   3.6966893218093214`*^9}, 3.6966893534493656`*^9, {3.6966894017894335`*^9, 
   3.6966894024494343`*^9}, 3.6966894403794875`*^9, {3.696689524579605*^9, 
   3.696689545799635*^9}, {3.6966896938898425`*^9, 3.69668974903992*^9}, {
   3.6966899909502583`*^9, 3.6966900131502895`*^9}, {3.696690059300354*^9, 
   3.696690067870366*^9}, {3.6966901023504143`*^9, 3.696690119640438*^9}, 
   3.6966902087005634`*^9, {3.696690269610648*^9, 3.696690285650671*^9}, {
   3.6966904369308825`*^9, 3.6966904480408983`*^9}, 3.6966905310310144`*^9, {
   3.69669056386106*^9, 3.6966905640510607`*^9}, {3.6966906768212185`*^9, 
   3.6966906974212475`*^9}, {3.696848369829146*^9, 3.6968484056391964`*^9}, {
   3.6968484664992814`*^9, 3.6968484989993267`*^9}, {3.6968485344293766`*^9, 
   3.696848561809415*^9}, {3.6968486273195066`*^9, 3.696848628719508*^9}, {
   3.6968486673195624`*^9, 3.6968486728095703`*^9}, {3.6968487103496227`*^9, 
   3.696848745119672*^9}, {3.696848783969726*^9, 3.696848833159795*^9}, {
   3.696849610620883*^9, 3.6968496433509293`*^9}, {3.696852856245427*^9, 
   3.6968528687854443`*^9}, {3.696852904775495*^9, 3.696852947585555*^9}, {
   3.69685299409562*^9, 3.6968530494756975`*^9}, {3.696853309626062*^9, 
   3.6968534102962027`*^9}, {3.6968534477262554`*^9, 
   3.6968534515762606`*^9}, {3.6968535064363375`*^9, 3.696853571396428*^9}, {
   3.696853632906514*^9, 3.696853670556567*^9}, {3.6968537012866096`*^9, 
   3.696853718926635*^9}, {3.6968537656167*^9, 3.696853774756713*^9}, {
   3.696853836246799*^9, 3.6968538409668055`*^9}, {3.6968538924968777`*^9, 
   3.6968540072270384`*^9}, {3.6968569077810993`*^9, 
   3.6968569292611294`*^9}, {3.6968582910830355`*^9, 3.696858359283131*^9}, {
   3.696908932696312*^9, 3.696908949856336*^9}, {3.697719384482682*^9, 
   3.697719385272683*^9}, {3.6977195732529464`*^9, 3.697719628173023*^9}, {
   3.6977197033531284`*^9, 3.6977197182431493`*^9}, 3.6977199911235313`*^9, {
   3.6977204504061747`*^9, 3.6977204574561844`*^9}, {3.6977826595151353`*^9, 
   3.6977826637251415`*^9}, {3.6977827300452337`*^9, 3.697782738075245*^9}, 
   3.7268858695439386`*^9, {3.7268859601130075`*^9, 3.7268860115318284`*^9}, {
   3.796614127709092*^9, 3.796614128038212*^9}, {3.796983648381215*^9, 
   3.7969836524320917`*^9}, {3.7969837026315928`*^9, 
   3.7969837094527044`*^9}, {3.7969837712836246`*^9, 3.7969837798037767`*^9}, 
   3.7969838652120705`*^9, 3.7969839072046843`*^9, 3.796983953777775*^9, 
   3.79698399415638*^9, 3.796984041036425*^9, {3.797755631471511*^9, 
   3.79775578196251*^9}, {3.7977558464918623`*^9, 3.797755859194564*^9}, {
   3.797755909610131*^9, 3.7977559831538444`*^9}, {3.7977560140160637`*^9, 
   3.797756056232184*^9}, {3.797756097527994*^9, 3.7977561454713607`*^9}, 
   3.7977561773991113`*^9, {3.7979049589200907`*^9, 3.7979049602935476`*^9}, 
   3.7979111978974085`*^9, {3.797911413567765*^9, 3.7979114146126385`*^9}, 
   3.8275688637688637`*^9, {3.8277583153794594`*^9, 3.827758316909443*^9}, {
   3.828149669297723*^9, 3.8281496710609474`*^9}, 3.828149886569129*^9},
 CellLabel->"In[3]:=",ExpressionUUID->"88a6ba85-74f9-40a5-aef4-3dbff1896c02"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"frequence", " ", "domain", " ", "representation"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"data", "=", 
     RowBox[{"Transpose", "[", "series", "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"T", "=", 
     RowBox[{
      RowBox[{"Last", "[", 
       RowBox[{"data", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}], "-", 
      RowBox[{"First", "[", 
       RowBox[{"data", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"abs\[Lambda]", "=", 
     RowBox[{"Abs", "[", 
      RowBox[{"Fourier", "[", 
       RowBox[{"data", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}]}], ";"}], "  ", 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "amplitude", " ", "as", " ", "function", " ", "of", " ", "frequency"}], 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f\[Lambda]", "=", 
     RowBox[{"MapIndexed", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"First", "[", "#2", "]"}], "-", "1"}], ")"}], "/", "T"}],
           ",", "#1"}], "}"}], "&"}], ",", "abs\[Lambda]"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"g2", "=", 
    RowBox[{"ListPlot", "[", 
     RowBox[{"f\[Lambda]", ",", 
      RowBox[{"Joined", "\[Rule]", "True"}], ",", 
      RowBox[{"PlotRange", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"0", ",", "100"}], "}"}], ",", "All"}], "}"}]}], ",", 
      RowBox[{"Frame", "\[Rule]", "True"}], ",", 
      RowBox[{"FrameLabel", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{
        "\"\<f [1/s]\>\"", ",", 
         "\"\<\[Alpha] [m/\!\(\*SuperscriptBox[\(s\), \(2\)]\)]\>\""}], 
        "}"}]}]}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.6895068568568707`*^9, 3.6895068606168756`*^9}, {
   3.689507096377206*^9, 3.689507101947213*^9}, {3.689508864916235*^9, 
   3.6895088669362373`*^9}, 3.689509046516489*^9, 3.6895093142295237`*^9, {
   3.689509404105324*^9, 3.6895094295353594`*^9}, {3.6895120381904573`*^9, 
   3.6895120397004595`*^9}, {3.689512122370575*^9, 3.689512157960625*^9}, {
   3.689512282330799*^9, 3.6895123101408377`*^9}, {3.6895131735720468`*^9, 
   3.689513295222217*^9}, {3.6895133805523367`*^9, 3.689513462322451*^9}, {
   3.689513526522541*^9, 3.689513592872634*^9}, {3.6895147961043186`*^9, 
   3.689514800934325*^9}, 3.689585657689865*^9, {3.689585919908863*^9, 
   3.6895859200938735`*^9}, 3.6896675842177234`*^9, {3.6900856637946014`*^9, 
   3.6900858381548457`*^9}, {3.6900859193649592`*^9, 
   3.6900860779651814`*^9}, {3.6900895242580185`*^9, 
   3.6900895290180254`*^9}, {3.690094654552212*^9, 3.6900947412723336`*^9}, {
   3.690094861452502*^9, 3.690094863612505*^9}, {3.6900949089725685`*^9, 
   3.690095000202696*^9}, {3.690172427710157*^9, 3.6901724615802045`*^9}, {
   3.690172499810258*^9, 3.690172554560334*^9}, {3.6901725949403915`*^9, 
   3.690172619990426*^9}, {3.690172658880481*^9, 3.690172733630585*^9}, {
   3.6901729886409426`*^9, 3.690173006360967*^9}, {3.690173049411027*^9, 
   3.690173089801084*^9}, 3.690173225581274*^9, 3.6901732709113374`*^9, {
   3.6901733763614855`*^9, 3.69017342987156*^9}, 3.6901736413718557`*^9, {
   3.6901736732519007`*^9, 3.690173747172004*^9}, {3.6901738114620943`*^9, 
   3.690173812002095*^9}, {3.690173942672278*^9, 3.6901739433122787`*^9}, {
   3.6947797562079945`*^9, 3.6947797740580196`*^9}, {3.694779859108139*^9, 
   3.6947799021581993`*^9}, {3.6948765488222756`*^9, 
   3.6948767582925687`*^9}, {3.6948771329030933`*^9, 3.694877410933482*^9}, {
   3.694877449223536*^9, 3.694877469123564*^9}, {3.6949175866210814`*^9, 
   3.694917591951089*^9}, {3.6949176842312183`*^9, 3.694917684771219*^9}, {
   3.6949177372612925`*^9, 3.6949178831114964`*^9}, {3.6949179676316147`*^9, 
   3.694917972901622*^9}, {3.694918006341669*^9, 3.694918019701688*^9}, {
   3.694918057981741*^9, 3.6949181203618283`*^9}, {3.6949181744819045`*^9, 
   3.6949181754019055`*^9}, {3.6949182156519623`*^9, 
   3.6949182234419727`*^9}, {3.694918335022129*^9, 3.6949183942422123`*^9}, {
   3.6949184418922787`*^9, 3.6949184896423454`*^9}, 3.6949185895924854`*^9, {
   3.694918662322587*^9, 3.6949186658425922`*^9}, {3.6949188708028793`*^9, 
   3.6949189242029543`*^9}, {3.6949189605230045`*^9, 3.694919119233227*^9}, 
   3.6949191503732705`*^9, {3.69491919269333*^9, 3.694919202403343*^9}, {
   3.694919236243391*^9, 3.694919363013568*^9}, {3.694919403093624*^9, 
   3.6949194079936314`*^9}, {3.6949194522036934`*^9, 3.6949195440338216`*^9}, 
   3.694919577063868*^9, {3.6949196167139235`*^9, 3.694919636613951*^9}, 
   3.694919702364043*^9, {3.694919742224099*^9, 3.6949198399642363`*^9}, {
   3.6949210692072573`*^9, 3.694921117877326*^9}, 3.6949211551773777`*^9, 
   3.694921253767516*^9, {3.694921490017847*^9, 3.69492150647787*^9}, {
   3.696689244189213*^9, 3.696689269359248*^9}, {3.696689311059306*^9, 
   3.6966893218093214`*^9}, 3.6966893534493656`*^9, {3.6966894017894335`*^9, 
   3.6966894024494343`*^9}, 3.6966894403794875`*^9, {3.696689524579605*^9, 
   3.696689545799635*^9}, {3.6966896938898425`*^9, 3.69668974903992*^9}, {
   3.6966899909502583`*^9, 3.6966900131502895`*^9}, {3.696690059300354*^9, 
   3.696690067870366*^9}, {3.6966901023504143`*^9, 3.696690119640438*^9}, 
   3.6966902087005634`*^9, {3.696690269610648*^9, 3.696690285650671*^9}, {
   3.6966904369308825`*^9, 3.6966904480408983`*^9}, 3.6966905310310144`*^9, {
   3.69669056386106*^9, 3.6966905640510607`*^9}, {3.6966906768212185`*^9, 
   3.6966906974212475`*^9}, {3.696848369829146*^9, 3.6968484056391964`*^9}, {
   3.6968484664992814`*^9, 3.6968484989993267`*^9}, {3.6968485344293766`*^9, 
   3.696848561809415*^9}, {3.6968486273195066`*^9, 3.696848628719508*^9}, {
   3.6968486673195624`*^9, 3.6968486728095703`*^9}, {3.6968487103496227`*^9, 
   3.696848745119672*^9}, {3.696848783969726*^9, 3.696848833159795*^9}, {
   3.696849610620883*^9, 3.6968496433509293`*^9}, {3.696852856245427*^9, 
   3.6968528687854443`*^9}, {3.696852904775495*^9, 3.696852947585555*^9}, {
   3.69685299409562*^9, 3.6968530494756975`*^9}, {3.696853309626062*^9, 
   3.6968534102962027`*^9}, {3.6968534477262554`*^9, 
   3.6968534515762606`*^9}, {3.6968535064363375`*^9, 3.696853571396428*^9}, {
   3.696853632906514*^9, 3.696853670556567*^9}, {3.6968537012866096`*^9, 
   3.696853718926635*^9}, {3.6968537656167*^9, 3.696853774756713*^9}, {
   3.696853836246799*^9, 3.6968538409668055`*^9}, {3.6968538924968777`*^9, 
   3.6968540072270384`*^9}, {3.6968569077810993`*^9, 
   3.6968569292611294`*^9}, {3.6968582910830355`*^9, 3.696858359283131*^9}, {
   3.696908932696312*^9, 3.696908949856336*^9}, {3.697719384482682*^9, 
   3.697719385272683*^9}, {3.6977195732529464`*^9, 3.697719628173023*^9}, {
   3.6977197033531284`*^9, 3.6977197182431493`*^9}, 3.6977199911235313`*^9, {
   3.6977204504061747`*^9, 3.6977204574561844`*^9}, {3.6977826595151353`*^9, 
   3.6977826637251415`*^9}, {3.6977827300452337`*^9, 3.697782738075245*^9}, 
   3.7268858695439386`*^9, {3.7268859601130075`*^9, 3.7268860115318284`*^9}, {
   3.796614127709092*^9, 3.796614128038212*^9}, {3.796983648381215*^9, 
   3.7969836524320917`*^9}, {3.7969837026315928`*^9, 
   3.7969837094527044`*^9}, {3.7969837712836246`*^9, 3.7969837798037767`*^9}, 
   3.7969838652120705`*^9, 3.7969839072046843`*^9, 3.796983953777775*^9, 
   3.79698399415638*^9, 3.796984041036425*^9, {3.797755631471511*^9, 
   3.79775578196251*^9}, {3.7977558464918623`*^9, 3.797755859194564*^9}, {
   3.797755909610131*^9, 3.7977559831538444`*^9}, {3.7977560140160637`*^9, 
   3.797756056232184*^9}, {3.797756097527994*^9, 3.7977561454713607`*^9}, 
   3.7977561773991113`*^9, {3.7979049589200907`*^9, 3.7979049602935476`*^9}, 
   3.7979111978974085`*^9, {3.7979113115313244`*^9, 3.797911316944172*^9}, {
   3.8277583090377865`*^9, 3.827758310516081*^9}, {3.82775845024944*^9, 
   3.8277584587040834`*^9}},
 CellLabel->"In[6]:=",ExpressionUUID->"fb19b1de-4a78-49bc-8c4f-369353f6c00f"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "find", " ", "the", " ", "peaks", " ", "and", " ", "pick", " ", "the", " ",
     "two", " ", "first"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"peakp", "=", 
     RowBox[{
      RowBox[{"FindPeaks", "[", "abs\[Lambda]", "]"}], "[", 
      RowBox[{"[", 
       RowBox[{
        RowBox[{"1", ";;", "2"}], ",", "1"}], "]"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"f\[Lambda]", "[", 
    RowBox[{"[", 
     RowBox[{"peakp", ",", "1"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.6895068568568707`*^9, 3.6895068606168756`*^9}, {
   3.689507096377206*^9, 3.689507101947213*^9}, {3.689508864916235*^9, 
   3.6895088669362373`*^9}, 3.689509046516489*^9, 3.6895093142295237`*^9, {
   3.689509404105324*^9, 3.6895094295353594`*^9}, {3.6895120381904573`*^9, 
   3.6895120397004595`*^9}, {3.689512122370575*^9, 3.689512157960625*^9}, {
   3.689512282330799*^9, 3.6895123101408377`*^9}, {3.6895131735720468`*^9, 
   3.689513295222217*^9}, {3.6895133805523367`*^9, 3.689513462322451*^9}, {
   3.689513526522541*^9, 3.689513592872634*^9}, {3.6895147961043186`*^9, 
   3.689514800934325*^9}, 3.689585657689865*^9, {3.689585919908863*^9, 
   3.6895859200938735`*^9}, 3.6896675842177234`*^9, {3.6900856637946014`*^9, 
   3.6900858381548457`*^9}, {3.6900859193649592`*^9, 
   3.6900860779651814`*^9}, {3.6900895242580185`*^9, 
   3.6900895290180254`*^9}, {3.690094654552212*^9, 3.6900947412723336`*^9}, {
   3.690094861452502*^9, 3.690094863612505*^9}, {3.6900949089725685`*^9, 
   3.690095000202696*^9}, {3.690172427710157*^9, 3.6901724615802045`*^9}, {
   3.690172499810258*^9, 3.690172554560334*^9}, {3.6901725949403915`*^9, 
   3.690172619990426*^9}, {3.690172658880481*^9, 3.690172733630585*^9}, {
   3.6901729886409426`*^9, 3.690173006360967*^9}, {3.690173049411027*^9, 
   3.690173089801084*^9}, 3.690173225581274*^9, 3.6901732709113374`*^9, {
   3.6901733763614855`*^9, 3.69017342987156*^9}, 3.6901736413718557`*^9, {
   3.6901736732519007`*^9, 3.690173747172004*^9}, {3.6901738114620943`*^9, 
   3.690173812002095*^9}, {3.690173942672278*^9, 3.6901739433122787`*^9}, {
   3.6947797562079945`*^9, 3.6947797740580196`*^9}, {3.694779859108139*^9, 
   3.6947799021581993`*^9}, {3.6948765488222756`*^9, 
   3.6948767582925687`*^9}, {3.6948771329030933`*^9, 3.694877410933482*^9}, {
   3.694877449223536*^9, 3.694877469123564*^9}, {3.6949175866210814`*^9, 
   3.694917591951089*^9}, {3.6949176842312183`*^9, 3.694917684771219*^9}, {
   3.6949177372612925`*^9, 3.6949178831114964`*^9}, {3.6949179676316147`*^9, 
   3.694917972901622*^9}, {3.694918006341669*^9, 3.694918019701688*^9}, {
   3.694918057981741*^9, 3.6949181203618283`*^9}, {3.6949181744819045`*^9, 
   3.6949181754019055`*^9}, {3.6949182156519623`*^9, 
   3.6949182234419727`*^9}, {3.694918335022129*^9, 3.6949183942422123`*^9}, {
   3.6949184418922787`*^9, 3.6949184896423454`*^9}, 3.6949185895924854`*^9, {
   3.694918662322587*^9, 3.6949186658425922`*^9}, {3.6949188708028793`*^9, 
   3.6949189242029543`*^9}, {3.6949189605230045`*^9, 3.694919119233227*^9}, 
   3.6949191503732705`*^9, {3.69491919269333*^9, 3.694919202403343*^9}, {
   3.694919236243391*^9, 3.694919363013568*^9}, {3.694919403093624*^9, 
   3.6949194079936314`*^9}, {3.6949194522036934`*^9, 3.6949195440338216`*^9}, 
   3.694919577063868*^9, {3.6949196167139235`*^9, 3.694919636613951*^9}, 
   3.694919702364043*^9, {3.694919742224099*^9, 3.6949198399642363`*^9}, {
   3.6949210692072573`*^9, 3.694921117877326*^9}, 3.6949211551773777`*^9, 
   3.694921253767516*^9, {3.694921490017847*^9, 3.69492150647787*^9}, {
   3.696689244189213*^9, 3.696689269359248*^9}, {3.696689311059306*^9, 
   3.6966893218093214`*^9}, 3.6966893534493656`*^9, {3.6966894017894335`*^9, 
   3.6966894024494343`*^9}, 3.6966894403794875`*^9, {3.696689524579605*^9, 
   3.696689545799635*^9}, {3.6966896938898425`*^9, 3.69668974903992*^9}, {
   3.6966899909502583`*^9, 3.6966900131502895`*^9}, {3.696690059300354*^9, 
   3.696690067870366*^9}, {3.6966901023504143`*^9, 3.696690119640438*^9}, 
   3.6966902087005634`*^9, {3.696690269610648*^9, 3.696690285650671*^9}, {
   3.6966904369308825`*^9, 3.6966904480408983`*^9}, 3.6966905310310144`*^9, {
   3.69669056386106*^9, 3.6966905640510607`*^9}, {3.6966906768212185`*^9, 
   3.6966906974212475`*^9}, {3.696848369829146*^9, 3.6968484056391964`*^9}, {
   3.6968484664992814`*^9, 3.6968484989993267`*^9}, {3.6968485344293766`*^9, 
   3.696848561809415*^9}, {3.6968486273195066`*^9, 3.696848628719508*^9}, {
   3.6968486673195624`*^9, 3.6968486728095703`*^9}, {3.6968487103496227`*^9, 
   3.696848745119672*^9}, {3.696848783969726*^9, 3.696848833159795*^9}, {
   3.696849610620883*^9, 3.6968496433509293`*^9}, {3.696852856245427*^9, 
   3.6968528687854443`*^9}, {3.696852904775495*^9, 3.696852947585555*^9}, {
   3.69685299409562*^9, 3.6968530494756975`*^9}, {3.696853309626062*^9, 
   3.6968534102962027`*^9}, {3.6968534477262554`*^9, 
   3.6968534515762606`*^9}, {3.6968535064363375`*^9, 3.696853571396428*^9}, {
   3.696853632906514*^9, 3.696853670556567*^9}, {3.6968537012866096`*^9, 
   3.696853718926635*^9}, {3.6968537656167*^9, 3.696853774756713*^9}, {
   3.696853836246799*^9, 3.6968538409668055`*^9}, {3.6968538924968777`*^9, 
   3.6968540072270384`*^9}, {3.6968569077810993`*^9, 
   3.6968569292611294`*^9}, {3.6968582910830355`*^9, 3.696858359283131*^9}, {
   3.696908932696312*^9, 3.696908949856336*^9}, {3.697719384482682*^9, 
   3.697719385272683*^9}, {3.6977195732529464`*^9, 3.697719628173023*^9}, {
   3.6977197033531284`*^9, 3.6977197182431493`*^9}, 3.6977199911235313`*^9, {
   3.6977204504061747`*^9, 3.6977204574561844`*^9}, {3.6977826595151353`*^9, 
   3.6977826637251415`*^9}, {3.6977827300452337`*^9, 3.697782738075245*^9}, 
   3.7268858695439386`*^9, {3.7268859601130075`*^9, 3.7268860115318284`*^9}, {
   3.796614127709092*^9, 3.796614128038212*^9}, {3.796983648381215*^9, 
   3.7969836524320917`*^9}, {3.7969837026315928`*^9, 
   3.7969837094527044`*^9}, {3.7969837712836246`*^9, 3.7969837798037767`*^9}, 
   3.7969838652120705`*^9, 3.7969839072046843`*^9, 3.796983953777775*^9, 
   3.79698399415638*^9, 3.796984041036425*^9, {3.797755631471511*^9, 
   3.79775578196251*^9}, {3.7977558464918623`*^9, 3.797755859194564*^9}, {
   3.797755909610131*^9, 3.7977559831538444`*^9}, {3.7977560140160637`*^9, 
   3.797756056232184*^9}, {3.797756097527994*^9, 3.7977561454713607`*^9}, 
   3.7977561773991113`*^9, {3.7979049589200907`*^9, 3.7979049602935476`*^9}, 
   3.7979111978974085`*^9, {3.7979113115313244`*^9, 3.797911316944172*^9}, 
   3.8281540835524607`*^9},ExpressionUUID->"ccd53af8-ad27-48e8-bd14-\
81fd14600e0f"]
}, Closed]],

Cell[CellGroupData[{

Cell["Simplified models", "Subsection",
 CellChangeTimes->{{3.473055512584*^9, 3.4730555529630003`*^9}, {
   3.4730575328269997`*^9, 3.473057539222*^9}, {3.4734691424642*^9, 
   3.4734691605758*^9}, {3.4736023142902*^9, 3.4736023273162003`*^9}, {
   3.568715518105007*^9, 3.568715523097016*^9}, 3.5994750893663483`*^9, 
   3.6297916561604834`*^9, {3.63073778491436*^9, 3.630737793931176*^9}, {
   3.6624537703648133`*^9, 3.6624537732749796`*^9}, 3.685248933840145*^9, {
   3.685267058363743*^9, 3.6852671273738394`*^9}, {3.685323461463516*^9, 
   3.6853234663335233`*^9}, {3.6858764211836576`*^9, 3.6858764235236607`*^9}, 
   3.6859482147140565`*^9, {3.686226212775403*^9, 3.6862262151354065`*^9}, {
   3.686466369367293*^9, 3.686466390117322*^9}, {3.732333324403823*^9, 
   3.7323333512571926`*^9}, {3.7330531980857735`*^9, 3.733053222557193*^9}, {
   3.7330643803419495`*^9, 3.733064384574917*^9}, {3.8277964716761804`*^9, 
   3.8277964754816017`*^9}, 
   3.8277978043562117`*^9},ExpressionUUID->"1c1abae9-8db2-473b-8232-\
8560ee29afde"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Parameters", " ", "of", " ", "models"}], "*)"}], " ", 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"par", "=", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"mf", "->", 
         RowBox[{"\[Rho]a", " ", "L", " ", "W", " ", "t"}]}], ",", 
        RowBox[{"mb", "\[Rule]", 
         RowBox[{"4", " ", "\[Rho]s", " ", 
          FractionBox[
           RowBox[{
            SuperscriptBox["d", "2"], " ", "\[Pi]"}], "4"], " ", "h"}]}], ",", 
        RowBox[{"kb", "\[Rule]", 
         RowBox[{"4", "*", "12", 
          FractionBox[
           RowBox[{"Es", " ", 
            SuperscriptBox["d", "4"], " ", "\[Pi]"}], 
           RowBox[{"64", " ", 
            SuperscriptBox["h", "3"]}]]}]}]}], " ", "}"}], " ", "//.", 
      "val"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Grid", "[", 
    RowBox[{
     RowBox[{"Partition", "[", 
      RowBox[{
       RowBox[{"Flatten", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"par", "/.", 
           RowBox[{"Rule", "\[Rule]", "List"}]}], ")"}], " ", "/.", " ", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"a_", ",", "b_"}], "}"}], "\[Rule]", 
          RowBox[{"{", 
           RowBox[{"a", ",", "b"}], "}"}]}]}], "]"}], ",", "2"}], "]"}], ",", 
     RowBox[{"Alignment", "\[Rule]", "Left"}], ",", 
     RowBox[{"Spacings", "->", 
      RowBox[{"{", 
       RowBox[{"5", ",", "1.2"}], "}"}]}], ",", 
     RowBox[{"Frame", "\[Rule]", "All"}]}], "]"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.827797788150197*^9, 3.827797797184443*^9}, {
   3.8277978429284835`*^9, 3.8277978695260425`*^9}, {3.8277985001647806`*^9, 
   3.8277985124780045`*^9}, 3.827798544557251*^9, {3.827798577191079*^9, 
   3.8277986598586826`*^9}},
 CellLabel->"In[5]:=",ExpressionUUID->"683dd3b0-b896-4e12-9f2d-79533e986793"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "Mass", " ", "and", " ", "stiffness", " ", "matrices", " ", "by", " ", 
    "PSM"}], "*)"}], " ", "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"M", "=", 
     RowBox[{
      RowBox[{"mf", 
       RowBox[{"(", GridBox[{
          {"1", "0"},
          {"0", "1"}
         }], ")"}]}], "+", 
      RowBox[{
       FractionBox["mb", "2"], 
       RowBox[{"(", GridBox[{
          {"2", "0"},
          {"0", "1"}
         }], ")"}]}]}]}], ";", "   ", 
    RowBox[{"K", "=", 
     RowBox[{"kb", 
      RowBox[{"(", GridBox[{
         {"2", 
          RowBox[{"-", "1"}]},
         {
          RowBox[{"-", "1"}], "1"}
        }], ")"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"{", 
      RowBox[{"vals", ",", "vecs"}], "}"}], "=", 
     RowBox[{
      RowBox[{"Eigensystem", "[", 
       RowBox[{
        RowBox[{"Inverse", "[", "M", "]"}], ".", "K"}], "]"}], " ", "/.", 
      "par"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Sqrt", "[", "vals", "]"}], "/", 
    RowBox[{"(", 
     RowBox[{"2", " ", "Pi"}], ")"}]}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"vecs", "[", 
     RowBox[{"[", "1", "]"}], "]"}], "/", 
    RowBox[{"Norm", "[", 
     RowBox[{"vecs", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"vecs", "[", 
     RowBox[{"[", "2", "]"}], "]"}], "/", 
    RowBox[{"Norm", "[", 
     RowBox[{"vecs", "[", 
      RowBox[{"[", "2", "]"}], "]"}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.732963786982458*^9, 3.732963800329712*^9}, {
   3.7329638807173033`*^9, 3.7329639838406887`*^9}, {3.7329640219057136`*^9, 
   3.7329640679494686`*^9}, {3.7329641106095915`*^9, 
   3.7329641238605337`*^9}, {3.7329641586934195`*^9, 
   3.7329641775638523`*^9}, {3.7329642131292753`*^9, 
   3.7329642468139324`*^9}, {3.732964291650645*^9, 3.732964295266472*^9}, {
   3.7329643646012278`*^9, 3.7329645915297546`*^9}, {3.732964625547305*^9, 
   3.732964633879656*^9}, {3.732964682913621*^9, 3.732964776488607*^9}, {
   3.732964812880664*^9, 3.7329648406399565`*^9}, {3.7329649084689856`*^9, 
   3.732965030133172*^9}, {3.732965104873771*^9, 3.7329651443322325`*^9}, {
   3.7329651781396537`*^9, 3.73296520770864*^9}, {3.732965254213273*^9, 
   3.7329652993251047`*^9}, {3.7329653860024896`*^9, 
   3.7329653897065096`*^9}, {3.7329654581300693`*^9, 
   3.7329654674851675`*^9}, {3.7329655079892044`*^9, 3.7329655085784197`*^9}, 
   3.733053175193036*^9, {3.733054323981842*^9, 3.7330543310918913`*^9}, {
   3.7330544527768817`*^9, 3.733054457152335*^9}, {3.733054573883443*^9, 
   3.7330546234511395`*^9}, {3.7330643435838375`*^9, 
   3.7330643547895327`*^9}, {3.7330643919191055`*^9, 
   3.7330644303402214`*^9}, {3.7331050956954727`*^9, 
   3.7331051189053106`*^9}, {3.7331111086485624`*^9, 
   3.7331111523877835`*^9}, {3.7331116359386663`*^9, 3.733111681037479*^9}, 
   3.7331466985132923`*^9, {3.733146742189515*^9, 3.7331467425179462`*^9}, {
   3.7331882689213552`*^9, 3.7331882705625243`*^9}, {3.733456902620245*^9, 
   3.7334569073552113`*^9}, {3.7961787078657227`*^9, 
   3.7961787392993565`*^9}, {3.7961787880849714`*^9, 
   3.7961788914598694`*^9}, {3.7961789272012496`*^9, 3.796178933116965*^9}, {
   3.796178969908412*^9, 3.7961791130644875`*^9}, 3.7961792235337663`*^9, {
   3.79617937704619*^9, 3.796179391357159*^9}, 3.796179443627824*^9, 
   3.796179486077304*^9, {3.7961796100746794`*^9, 3.796179612132883*^9}, {
   3.7977562054888687`*^9, 3.797756212489262*^9}, {3.7977613590730443`*^9, 
   3.7977614738483243`*^9}, {3.7977616174800205`*^9, 
   3.7977616274995165`*^9}, {3.7977616690287666`*^9, 3.797761757857272*^9}, {
   3.797761790940256*^9, 3.797761803255534*^9}, {3.7977619109048243`*^9, 
   3.797761994990079*^9}, {3.7979038189548445`*^9, 3.7979038355791883`*^9}, 
   3.797904847600319*^9, {3.797909219979147*^9, 3.797909239314906*^9}, {
   3.797909334368704*^9, 3.797909347416195*^9}, 3.8277964392248774`*^9, {
   3.8277964901992283`*^9, 3.8277965171047153`*^9}, {3.8277965644821644`*^9, 
   3.8277966049027786`*^9}, {3.827796646095589*^9, 3.827796914782463*^9}, {
   3.82779700113986*^9, 3.827797056179701*^9}, {3.827797095529213*^9, 
   3.827797336267351*^9}, {3.827797458188809*^9, 3.827797463261814*^9}, {
   3.827797504924575*^9, 3.827797510724432*^9}, {3.8277975601003113`*^9, 
   3.8277978391739416`*^9}, {3.8277986764460735`*^9, 
   3.8277986818887997`*^9}, {3.827798888927939*^9, 3.82779889802024*^9}, {
   3.8280793688420067`*^9, 3.8280794129846125`*^9}},
 CellLabel->"In[7]:=",ExpressionUUID->"1c10060a-870b-4d42-a9fa-219fc56ee982"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "Mass", " ", "and", " ", "stiffness", " ", "matrices", " ", "by", " ", 
    "FDM"}], "*)"}], " ", "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"M", "=", 
     RowBox[{
      RowBox[{"mf", 
       RowBox[{"(", GridBox[{
          {"1", "0"},
          {"0", "1"}
         }], ")"}]}], "+", 
      RowBox[{"mb", 
       RowBox[{"(", GridBox[{
          {"1", "0"},
          {"0", "1"}
         }], ")"}]}]}]}], ";", "   ", 
    RowBox[{"K", "=", 
     RowBox[{"kb", 
      RowBox[{"(", GridBox[{
         {"2", 
          RowBox[{"-", "1"}]},
         {
          RowBox[{"-", "1"}], "1"}
        }], ")"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"vals", ",", "vecs"}], "}"}], "=", 
    RowBox[{
     RowBox[{"Eigensystem", "[", 
      RowBox[{
       RowBox[{"Inverse", "[", "M", "]"}], ".", "K"}], "]"}], " ", "/.", 
     "par"}]}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Sqrt", "[", "vals", "]"}], "/", 
    RowBox[{"(", 
     RowBox[{"2", " ", "Pi"}], ")"}]}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"vecs", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "/", 
     RowBox[{"Norm", "[", 
      RowBox[{"vecs", "[", 
       RowBox[{"[", "1", "]"}], "]"}], "]"}]}], " ", "//", "N"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"vecs", "[", 
      RowBox[{"[", "2", "]"}], "]"}], "/", 
     RowBox[{"Norm", "[", 
      RowBox[{"vecs", "[", 
       RowBox[{"[", "2", "]"}], "]"}], "]"}]}], " ", "//", "N"}]}]}]], "Input",
 CellChangeTimes->{{3.732963786982458*^9, 3.732963800329712*^9}, {
   3.7329638807173033`*^9, 3.7329639838406887`*^9}, {3.7329640219057136`*^9, 
   3.7329640679494686`*^9}, {3.7329641106095915`*^9, 
   3.7329641238605337`*^9}, {3.7329641586934195`*^9, 
   3.7329641775638523`*^9}, {3.7329642131292753`*^9, 
   3.7329642468139324`*^9}, {3.732964291650645*^9, 3.732964295266472*^9}, {
   3.7329643646012278`*^9, 3.7329645915297546`*^9}, {3.732964625547305*^9, 
   3.732964633879656*^9}, {3.732964682913621*^9, 3.732964776488607*^9}, {
   3.732964812880664*^9, 3.7329648406399565`*^9}, {3.7329649084689856`*^9, 
   3.732965030133172*^9}, {3.732965104873771*^9, 3.7329651443322325`*^9}, {
   3.7329651781396537`*^9, 3.73296520770864*^9}, {3.732965254213273*^9, 
   3.7329652993251047`*^9}, {3.7329653860024896`*^9, 
   3.7329653897065096`*^9}, {3.7329654581300693`*^9, 
   3.7329654674851675`*^9}, {3.7329655079892044`*^9, 3.7329655085784197`*^9}, 
   3.733053175193036*^9, {3.733054323981842*^9, 3.7330543310918913`*^9}, {
   3.7330544527768817`*^9, 3.733054457152335*^9}, {3.733054573883443*^9, 
   3.7330546234511395`*^9}, {3.7330643435838375`*^9, 
   3.7330643547895327`*^9}, {3.7330643919191055`*^9, 
   3.7330644303402214`*^9}, {3.7331050956954727`*^9, 
   3.7331051189053106`*^9}, {3.7331111086485624`*^9, 
   3.7331111523877835`*^9}, {3.7331116359386663`*^9, 3.733111681037479*^9}, 
   3.7331466985132923`*^9, {3.733146742189515*^9, 3.7331467425179462`*^9}, {
   3.7331882689213552`*^9, 3.7331882705625243`*^9}, {3.733456902620245*^9, 
   3.7334569073552113`*^9}, {3.7961787078657227`*^9, 
   3.7961787392993565`*^9}, {3.7961787880849714`*^9, 
   3.7961788914598694`*^9}, {3.7961789272012496`*^9, 3.796178933116965*^9}, {
   3.796178969908412*^9, 3.7961791130644875`*^9}, 3.7961792235337663`*^9, {
   3.79617937704619*^9, 3.796179391357159*^9}, 3.796179443627824*^9, 
   3.796179486077304*^9, {3.7961796100746794`*^9, 3.796179612132883*^9}, {
   3.7977562054888687`*^9, 3.797756212489262*^9}, {3.7977613590730443`*^9, 
   3.7977614738483243`*^9}, {3.7977616174800205`*^9, 
   3.7977616274995165`*^9}, {3.7977616690287666`*^9, 3.797761757857272*^9}, {
   3.797761790940256*^9, 3.797761803255534*^9}, {3.7977619109048243`*^9, 
   3.797761994990079*^9}, {3.7979038189548445`*^9, 3.7979038355791883`*^9}, 
   3.797904847600319*^9, {3.797909219979147*^9, 3.797909239314906*^9}, {
   3.797909334368704*^9, 3.797909347416195*^9}, 3.8277964392248774`*^9, {
   3.8277964901992283`*^9, 3.8277965171047153`*^9}, {3.8277965644821644`*^9, 
   3.8277966049027786`*^9}, {3.827796646095589*^9, 3.827796914782463*^9}, {
   3.82779700113986*^9, 3.827797056179701*^9}, {3.827797095529213*^9, 
   3.827797336267351*^9}, {3.827797458188809*^9, 3.827797463261814*^9}, {
   3.827797504924575*^9, 3.827797510724432*^9}, {3.8277975601003113`*^9, 
   3.8277978391739416`*^9}, {3.8277986764460735`*^9, 
   3.8277986818887997`*^9}, {3.8277987396367435`*^9, 
   3.8277987812757683`*^9}, {3.827798841585251*^9, 3.8277989232297945`*^9}, 
   3.8277990090331087`*^9, 3.8280795328646336`*^9, {3.8280795688368926`*^9, 
   3.828079575243995*^9}},
 CellLabel->"In[12]:=",ExpressionUUID->"e8468dda-5343-4252-9f58-23c187ef64e1"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "Mass", " ", "and", " ", "stiffness", " ", "matrices", " ", "by", " ", 
    "FEM"}], "*)"}], " ", "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"M", "=", 
     RowBox[{
      RowBox[{"mf", 
       RowBox[{"(", GridBox[{
          {"1", "0"},
          {"0", "1"}
         }], ")"}]}], "+", 
      RowBox[{
       FractionBox["mb", "6"], 
       RowBox[{"(", GridBox[{
          {"4", "1"},
          {"1", "2"}
         }], ")"}]}]}]}], ";", "   ", 
    RowBox[{"K", "=", 
     RowBox[{"kb", 
      RowBox[{"(", GridBox[{
         {"2", 
          RowBox[{"-", "1"}]},
         {
          RowBox[{"-", "1"}], "1"}
        }], ")"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"{", 
      RowBox[{"vals", ",", "vecs"}], "}"}], "=", 
     RowBox[{
      RowBox[{"Eigensystem", "[", 
       RowBox[{
        RowBox[{"Inverse", "[", "M", "]"}], ".", "K"}], "]"}], " ", "/.", 
      "par"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Sqrt", "[", "vals", "]"}], "/", 
    RowBox[{"(", 
     RowBox[{"2", " ", "Pi"}], ")"}]}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"vecs", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "/", 
     RowBox[{"Norm", "[", 
      RowBox[{"vecs", "[", 
       RowBox[{"[", "1", "]"}], "]"}], "]"}]}], " ", "//", "N"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"vecs", "[", 
      RowBox[{"[", "2", "]"}], "]"}], "/", 
     RowBox[{"Norm", "[", 
      RowBox[{"vecs", "[", 
       RowBox[{"[", "2", "]"}], "]"}], "]"}]}], " ", "//", "N"}]}]}]], "Input",
 CellChangeTimes->{{3.732963786982458*^9, 3.732963800329712*^9}, {
   3.7329638807173033`*^9, 3.7329639838406887`*^9}, {3.7329640219057136`*^9, 
   3.7329640679494686`*^9}, {3.7329641106095915`*^9, 
   3.7329641238605337`*^9}, {3.7329641586934195`*^9, 
   3.7329641775638523`*^9}, {3.7329642131292753`*^9, 
   3.7329642468139324`*^9}, {3.732964291650645*^9, 3.732964295266472*^9}, {
   3.7329643646012278`*^9, 3.7329645915297546`*^9}, {3.732964625547305*^9, 
   3.732964633879656*^9}, {3.732964682913621*^9, 3.732964776488607*^9}, {
   3.732964812880664*^9, 3.7329648406399565`*^9}, {3.7329649084689856`*^9, 
   3.732965030133172*^9}, {3.732965104873771*^9, 3.7329651443322325`*^9}, {
   3.7329651781396537`*^9, 3.73296520770864*^9}, {3.732965254213273*^9, 
   3.7329652993251047`*^9}, {3.7329653860024896`*^9, 
   3.7329653897065096`*^9}, {3.7329654581300693`*^9, 
   3.7329654674851675`*^9}, {3.7329655079892044`*^9, 3.7329655085784197`*^9}, 
   3.733053175193036*^9, {3.733054323981842*^9, 3.7330543310918913`*^9}, {
   3.7330544527768817`*^9, 3.733054457152335*^9}, {3.733054573883443*^9, 
   3.7330546234511395`*^9}, {3.7330643435838375`*^9, 
   3.7330643547895327`*^9}, {3.7330643919191055`*^9, 
   3.7330644303402214`*^9}, {3.7331050956954727`*^9, 
   3.7331051189053106`*^9}, {3.7331111086485624`*^9, 
   3.7331111523877835`*^9}, {3.7331116359386663`*^9, 3.733111681037479*^9}, 
   3.7331466985132923`*^9, {3.733146742189515*^9, 3.7331467425179462`*^9}, {
   3.7331882689213552`*^9, 3.7331882705625243`*^9}, {3.733456902620245*^9, 
   3.7334569073552113`*^9}, {3.7961787078657227`*^9, 
   3.7961787392993565`*^9}, {3.7961787880849714`*^9, 
   3.7961788914598694`*^9}, {3.7961789272012496`*^9, 3.796178933116965*^9}, {
   3.796178969908412*^9, 3.7961791130644875`*^9}, 3.7961792235337663`*^9, {
   3.79617937704619*^9, 3.796179391357159*^9}, 3.796179443627824*^9, 
   3.796179486077304*^9, {3.7961796100746794`*^9, 3.796179612132883*^9}, {
   3.7977562054888687`*^9, 3.797756212489262*^9}, {3.7977613590730443`*^9, 
   3.7977614738483243`*^9}, {3.7977616174800205`*^9, 
   3.7977616274995165`*^9}, {3.7977616690287666`*^9, 3.797761757857272*^9}, {
   3.797761790940256*^9, 3.797761803255534*^9}, {3.7977619109048243`*^9, 
   3.797761994990079*^9}, {3.7979038189548445`*^9, 3.7979038355791883`*^9}, 
   3.797904847600319*^9, {3.797909219979147*^9, 3.797909239314906*^9}, {
   3.797909334368704*^9, 3.797909347416195*^9}, 3.8277964392248774`*^9, {
   3.8277964901992283`*^9, 3.8277965171047153`*^9}, {3.8277965644821644`*^9, 
   3.8277966049027786`*^9}, {3.827796646095589*^9, 3.827796914782463*^9}, {
   3.82779700113986*^9, 3.827797056179701*^9}, {3.827797095529213*^9, 
   3.827797336267351*^9}, {3.827797458188809*^9, 3.827797463261814*^9}, {
   3.827797504924575*^9, 3.827797510724432*^9}, {3.8277975601003113`*^9, 
   3.8277978391739416`*^9}, {3.8277986764460735`*^9, 
   3.8277986818887997`*^9}, {3.8277987396367435`*^9, 
   3.8277987812757683`*^9}, {3.827798841585251*^9, 3.8277989232297945`*^9}, {
   3.8277990194127755`*^9, 3.827799019814051*^9}, {3.8277990559128036`*^9, 
   3.827799085151805*^9}, {3.8280796724105425`*^9, 3.8280796750354853`*^9}},
 CellLabel->"In[17]:=",ExpressionUUID->"b50f9eec-13ff-472c-ba91-4e14a54a0794"]
}, Closed]]
}, Open  ]]
},
WindowSize->{745.875, 373.875},
WindowMargins->{{Automatic, 58.5}, {62.25, Automatic}},
Magnification:>1.25 Inherited,
FrontEndVersion->"12.1 for Microsoft Windows (64-bit) (June 19, 2020)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"07a44198-840f-46f5-9176-46370e18ff27"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 320, 5, 83, "Section",ExpressionUUID->"86686e14-fc55-42dc-b305-e2ed2e4b1a8a"],
Cell[903, 29, 2144, 37, 155, "Input",ExpressionUUID->"9174fd4b-30aa-441b-9412-80faf601d72b"],
Cell[3050, 68, 929, 25, 132, "Input",ExpressionUUID->"ea83f712-7b09-46d5-9060-e8b8fac5cbd1"],
Cell[3982, 95, 2526, 47, 155, "Input",ExpressionUUID->"f6db8324-9aef-4d10-9770-8d2f6a64e224"],
Cell[6511, 144, 2483, 46, 108, "Input",ExpressionUUID->"57ebc5c8-0469-4298-a7a8-c6cc23150f96"],
Cell[8997, 192, 2440, 44, 179, "Input",ExpressionUUID->"97108cef-50e6-4fa2-a79d-60e9f1ea42a0"]
}, Closed]],
Cell[CellGroupData[{
Cell[11474, 241, 1166, 17, 65, "Section",ExpressionUUID->"ddf24a73-a884-42d6-ad85-1c9ed46c5997"],
Cell[CellGroupData[{
Cell[12665, 262, 960, 14, 67, "Subsection",ExpressionUUID->"0f0487a8-da1d-494b-9d29-0a7939035da4"],
Cell[13628, 278, 8657, 171, 845, "Input",ExpressionUUID->"ce80bae6-fb13-49c7-8dbf-a3af0c00446e",
 InitializationCell->True]
}, Closed]],
Cell[CellGroupData[{
Cell[22322, 454, 1558, 22, 46, "Subsection",ExpressionUUID->"0ee45dc8-cf7c-473a-ae94-685c39f63702"],
Cell[23883, 478, 7312, 115, 139, "Input",ExpressionUUID->"88a6ba85-74f9-40a5-aef4-3dbff1896c02"],
Cell[31198, 595, 8174, 144, 234, "Input",ExpressionUUID->"fb19b1de-4a78-49bc-8c4f-369353f6c00f"],
Cell[39375, 741, 6683, 101, 84, "Input",ExpressionUUID->"ccd53af8-ad27-48e8-bd14-81fd14600e0f"]
}, Closed]],
Cell[CellGroupData[{
Cell[46095, 847, 1039, 15, 46, "Subsection",ExpressionUUID->"1c1abae9-8db2-473b-8232-8560ee29afde"],
Cell[47137, 864, 1917, 50, 186, "Input",ExpressionUUID->"683dd3b0-b896-4e12-9f2d-79533e986793"],
Cell[49057, 916, 4691, 98, 229, "Input",ExpressionUUID->"1c10060a-870b-4d42-a9fa-219fc56ee982"],
Cell[53751, 1016, 4845, 101, 249, "Input",ExpressionUUID->"e8468dda-5343-4252-9f58-23c187ef64e1"],
Cell[58599, 1119, 4948, 103, 252, "Input",ExpressionUUID->"b50f9eec-13ff-472c-ba91-4e14a54a0794"]
}, Closed]]
}, Open  ]]
}
]
*)