mathematica_ilt.nb

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

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

(* CreatedBy='Mathematica 11.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     57498,       1462]
NotebookOptionsPosition[     56163,       1413]
NotebookOutlinePosition[     56539,       1429]
CellTagsIndexPosition[     56496,       1426]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
Inverse Laplace Transform with Concentrated Matrix-Exponential Functions\
\>", "Title",
 CellChangeTimes->{
  3.802089556196247*^9, {3.802089603785129*^9, 
   3.802089620874042*^9}},ExpressionUUID->"202f27a9-1ca7-4231-b1b7-\
215903ec1cf8"],

Cell[CellGroupData[{

Cell[TextData[{
 "The following function performs the inverse transformation when f(",
 Cell[BoxData[
  FormBox[
   SuperscriptBox["s", "*"], TraditionalForm]],
  FormatType->"TraditionalForm"],
 ")=",
 Cell[BoxData[
  FormBox[
   SuperscriptBox["f", "*"], TraditionalForm]],
  FormatType->"TraditionalForm"],
 "(",
 Cell[BoxData[
  FormBox["s", TraditionalForm]],
  FormatType->"TraditionalForm"],
 ")"
}], "Section",
 CellChangeTimes->{{3.802089576079286*^9, 3.802089598393911*^9}, {
  3.940309007276932*^9, 3.940309028807987*^9}, {3.9403091171697755`*^9, 
  3.9403091495982847`*^9}, {3.9403091857202673`*^9, 3.940309216438936*^9}, {
  3.9403107914107122`*^9, 3.9403108005169706`*^9}},
 ExpressionUUID -> "a8f6ac9c-4991-4d97-98ce-4892bbcd0b84"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"ILT", "[", 
    RowBox[{"fun_", ",", " ", "T_", ",", " ", "maxFnEvals_", ",", " ", 
     RowBox[{"method_:", "\"\<cme\>\""}], ",", 
     RowBox[{"precision_:", "MachinePrecision"}]}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "params", ",", "i", ",", "\[Eta]", ",", "\[Beta]", ",", "\[Xi]", ",", 
       "n", ",", "k", ",", "prec"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"Length", "[", 
           RowBox[{"Position", "[", 
            RowBox[{
             RowBox[{"Names", "[", "\"\<Global`*\>\"", "]"}], ",", 
             "\"\<cmeParams\>\""}], "]"}], "]"}], "\[Equal]", "0"}], " ", "||",
          " ", 
         RowBox[{
          RowBox[{"Length", "[", "cmeParams", "]"}], "\[Equal]", "0"}]}], ",",
         "\[IndentingNewLine]", 
        RowBox[{"(*", 
         RowBox[{
          RowBox[{"Print", "[", "\"\<CME parameters loaded\>\"", "]"}], ";"}],
          "*)"}], "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"cmeParams", " ", "=", " ", 
          RowBox[{"Import", "[", 
           RowBox[{"FileNameJoin", "[", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"NotebookDirectory", "[", "]"}], ",", 
              "\"\<iltcme.json\>\""}], "}"}], "]"}], "]"}]}], ";"}]}], 
       "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"precision", "\[LessEqual]", "0"}], ",", 
        RowBox[{"prec", "=", "MachinePrecision"}], ",", 
        RowBox[{"prec", "=", "precision"}]}], "]"}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"method", "\[Equal]", "\"\<cme\>\""}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"params", " ", "=", " ", 
          RowBox[{"cmeParams", "[", 
           RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
         RowBox[{"For", "[", 
          RowBox[{
           RowBox[{"i", "=", "2"}], ",", 
           RowBox[{"i", "<", 
            RowBox[{"Length", "[", "cmeParams", "]"}]}], ",", 
           RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{
              RowBox[{
               RowBox[{"(", 
                RowBox[{"\"\<cv2\>\"", "/.", 
                 RowBox[{"cmeParams", "[", 
                  RowBox[{"[", "i", "]"}], "]"}]}], ")"}], " ", "<", 
               RowBox[{"(", 
                RowBox[{"\"\<cv2\>\"", "/.", "params"}], ")"}]}], " ", "&&", 
              " ", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"\"\<n\>\"", "+", "1"}], " ", "/.", 
                 RowBox[{"cmeParams", "[", 
                  RowBox[{"[", "i", "]"}], "]"}]}], ")"}], "\[LessEqual]", 
               "maxFnEvals"}]}], ",", 
             RowBox[{"params", "=", 
              RowBox[{"cmeParams", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}]}], " ", "]"}]}], 
          "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
         RowBox[{"\[Eta]", "=", 
          RowBox[{"SetPrecision", "[", 
           RowBox[{
            RowBox[{"Prepend", "[", 
             RowBox[{
              RowBox[{
               RowBox[{
                RowBox[{"\"\<a\>\"", "*", "\"\<mu1\>\""}], "+", 
                RowBox[{
                "\[ImaginaryI]", "*", "\"\<b\>\"", "*", "\"\<mu1\>\""}]}], "/.",
                "params"}], ",", 
              RowBox[{
               RowBox[{"\"\<c\>\"", "*", "\"\<mu1\>\""}], "/.", "params"}]}], 
             "]"}], ",", "prec"}], "]"}]}], ";", "\[IndentingNewLine]", 
         RowBox[{"\[Beta]", "=", 
          RowBox[{"SetPrecision", "[", 
           RowBox[{
            RowBox[{
             RowBox[{
              RowBox[{"Prepend", "[", 
               RowBox[{
                RowBox[{"1", "+", 
                 RowBox[{"\[ImaginaryI]", " ", "*", " ", "\"\<omega\>\"", "*", 
                  RowBox[{"Range", "[", 
                   RowBox[{"\"\<n\>\"", "/.", "params"}], "]"}]}]}], ",", 
                "1"}], "]"}], "*", "\"\<mu1\>\""}], " ", "/.", "params"}], 
            ",", "prec"}], "]"}]}], ";"}], "\[IndentingNewLine]", ",", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{"method", "\[Equal]", "\"\<euler\>\""}], ",", 
           "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"n", "=", 
             RowBox[{"If", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"Mod", "[", 
                 RowBox[{"maxFnEvals", ",", "2"}], "]"}], "\[Equal]", "0"}], 
               ",", 
               RowBox[{"maxFnEvals", "-", "1"}], ",", "maxFnEvals"}], "]"}]}],
             ";", "\[IndentingNewLine]", 
            RowBox[{"\[Beta]", "=", 
             RowBox[{"Table", "[", 
              RowBox[{
               RowBox[{"SetPrecision", "[", 
                RowBox[{
                 RowBox[{
                  FractionBox[
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"n", "-", "1"}], ")"}], 
                    RowBox[{"Log", "[", "10", "]"}]}], "6"], "+", 
                  RowBox[{"\[Pi]", " ", "\[ImaginaryI]", " ", 
                   RowBox[{"(", 
                    RowBox[{"k", "-", "1"}], ")"}]}]}], ",", "prec"}], "]"}], 
               ",", 
               RowBox[{"{", 
                RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";", 
            "\[IndentingNewLine]", 
            RowBox[{"\[Xi]", "=", 
             RowBox[{"ConstantArray", "[", 
              RowBox[{
               RowBox[{"SetPrecision", "[", 
                RowBox[{"1", ",", "prec"}], "]"}], ",", "n"}], "]"}]}], ";", 
            "\[IndentingNewLine]", 
            RowBox[{
             RowBox[{"\[Xi]", "[", 
              RowBox[{"[", "1", "]"}], "]"}], " ", "=", " ", 
             RowBox[{"SetPrecision", "[", 
              RowBox[{
               RowBox[{"1", "/", "2"}], ",", "prec"}], "]"}]}], ";", 
            "\[IndentingNewLine]", 
            RowBox[{
             RowBox[{"\[Xi]", "[", 
              RowBox[{"[", "n", "]"}], "]"}], "=", 
             RowBox[{"SetPrecision", "[", 
              RowBox[{
               FractionBox["1", 
                SuperscriptBox["2", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"n", "-", "1"}], ")"}], "/", "2"}]]], ",", 
               "prec"}], "]"}]}], ";", "\[IndentingNewLine]", 
            RowBox[{"For", "[", 
             RowBox[{
              RowBox[{"k", "=", "1"}], ",", 
              RowBox[{"k", "<", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"n", "-", "1"}], ")"}], "/", "2"}]}], ",", 
              RowBox[{"k", "++"}], ",", 
              RowBox[{
               RowBox[{"\[Xi]", "[", 
                RowBox[{"[", 
                 RowBox[{"n", "-", "k"}], "]"}], "]"}], "=", 
               RowBox[{
                RowBox[{"\[Xi]", "[", 
                 RowBox[{"[", 
                  RowBox[{"n", "-", "k", "+", "1"}], "]"}], "]"}], "+", " ", 
                RowBox[{
                 SuperscriptBox["2", 
                  RowBox[{
                   RowBox[{"-", 
                    RowBox[{"(", 
                    RowBox[{"n", "-", "1"}], ")"}]}], "/", "2"}]], 
                 RowBox[{"Binomial", "[", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"n", "-", "1"}], ")"}], "/", "2"}], ",", "k"}], 
                  "]"}]}]}]}]}], "]"}], ";", "\[IndentingNewLine]", 
            RowBox[{"\[Eta]", "=", 
             RowBox[{"Table", "[", 
              RowBox[{
               RowBox[{"SetPrecision", "[", 
                RowBox[{
                 RowBox[{
                  SuperscriptBox["10", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"n", "-", "1"}], ")"}], "/", "6"}]], 
                  SuperscriptBox[
                   RowBox[{"(", 
                    RowBox[{"-", "1"}], ")"}], 
                   RowBox[{"k", "-", "1"}]], 
                  RowBox[{"\[Xi]", "[", 
                   RowBox[{"[", "k", "]"}], "]"}]}], ",", "prec"}], "]"}], 
               ",", 
               RowBox[{"{", 
                RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
           "\[IndentingNewLine]", ",", "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"If", "[", 
             RowBox[{
              RowBox[{"method", "\[Equal]", "\"\<gaver\>\""}], ",", 
              "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{"n", "=", 
                RowBox[{"If", "[", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"Mod", "[", 
                    RowBox[{"maxFnEvals", ",", "2"}], "]"}], "\[Equal]", 
                   "1"}], ",", 
                  RowBox[{"maxFnEvals", "-", "1"}], ",", "maxFnEvals"}], 
                 "]"}]}], ";", "\[IndentingNewLine]", 
               RowBox[{"\[Beta]", "=", 
                RowBox[{"Table", "[", 
                 RowBox[{
                  RowBox[{"k", " ", 
                   RowBox[{"SetPrecision", "[", 
                    RowBox[{
                    RowBox[{"Log", "[", "2", "]"}], ",", "prec"}], "]"}]}], 
                  ",", 
                  RowBox[{"{", 
                   RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";", 
               "\[IndentingNewLine]", 
               RowBox[{"\[Eta]", "=", 
                RowBox[{"Table", "[", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"SetPrecision", "[", 
                    RowBox[{
                    RowBox[{"Log", "[", "2", "]"}], ",", "prec"}], "]"}], 
                   SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"-", "1"}], ")"}], 
                    RowBox[{
                    RowBox[{"n", "/", "2"}], "+", "k"}]], 
                   RowBox[{"Sum", "[", 
                    RowBox[{
                    RowBox[{
                    FractionBox[
                    SuperscriptBox["j", 
                    RowBox[{
                    RowBox[{"n", "/", "2"}], "+", "1"}]], 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{"n", "/", "2"}], ")"}], "!"}]], 
                    RowBox[{"Binomial", "[", 
                    RowBox[{
                    RowBox[{"n", "/", "2"}], ",", "j"}], "]"}], 
                    RowBox[{"Binomial", "[", 
                    RowBox[{
                    RowBox[{"2", "j"}], ",", "j"}], "]"}], 
                    RowBox[{"Binomial", "[", 
                    RowBox[{"j", ",", 
                    RowBox[{"k", "-", "j"}]}], "]"}]}], ",", 
                    RowBox[{"{", 
                    RowBox[{"j", ",", 
                    RowBox[{"Floor", "[", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{"k", "+", "1"}], ")"}], "/", "2"}], "]"}], ",", 
                    RowBox[{"Min", "[", 
                    RowBox[{"k", ",", 
                    RowBox[{"n", "/", "2"}]}], "]"}]}], "}"}]}], "]"}]}], 
                  "  ", ",", 
                  RowBox[{"{", 
                   RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
               ";"}]}], "\[IndentingNewLine]", "]"}], ";"}]}], 
          "\[IndentingNewLine]", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}],
       ";", "\[IndentingNewLine]", 
      RowBox[{"Return", "[", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{"ListQ", "[", "T", "]"}], ",", 
         RowBox[{"Table", "[", 
          RowBox[{
           RowBox[{"Re", "[", 
            RowBox[{
             RowBox[{"Sum", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"\[Eta]", "[", 
                 RowBox[{"[", "k", "]"}], "]"}], "*", 
                RowBox[{"fun", "[", 
                 RowBox[{
                  RowBox[{"\[Beta]", "[", 
                   RowBox[{"[", "k", "]"}], "]"}], "/", "x"}], "]"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"k", ",", "1", ",", 
                 RowBox[{"Length", "[", "\[Eta]", "]"}]}], "}"}]}], "]"}], 
             "/", "x"}], "]"}], ",", 
           RowBox[{"{", 
            RowBox[{"x", ",", "T"}], "}"}]}], "]"}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"Re", "[", 
           RowBox[{
            RowBox[{"Sum", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"\[Eta]", "[", 
                RowBox[{"[", "k", "]"}], "]"}], "*", 
               RowBox[{"fun", "[", 
                RowBox[{
                 RowBox[{"\[Beta]", "[", 
                  RowBox[{"[", "k", "]"}], "]"}], "/", "x"}], "]"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", "1", ",", 
                RowBox[{"Length", "[", "\[Eta]", "]"}]}], "}"}]}], "]"}], "/",
             "x"}], "]"}], "/.", 
          RowBox[{"x", "\[Rule]", "T"}]}]}], "]"}], "]"}], ";"}]}], 
    "\[IndentingNewLine]", "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.757325469151745*^9, 3.757325502280724*^9}, {
   3.75732570552604*^9, 3.757325767957542*^9}, {3.757325851475429*^9, 
   3.7573259078775578`*^9}, {3.7573259980687437`*^9, 3.75732607799338*^9}, {
   3.75732612019339*^9, 3.757326145684031*^9}, {3.7573261929718733`*^9, 
   3.7573262563262873`*^9}, 3.757396882026514*^9, {3.75739691295605*^9, 
   3.757397014808332*^9}, {3.7573970726919823`*^9, 3.757397221175847*^9}, {
   3.757397601976653*^9, 3.757397666846058*^9}, {3.7573980113072147`*^9, 
   3.757398027672312*^9}, {3.757398063947167*^9, 3.757398074906393*^9}, {
   3.757398144923143*^9, 3.7573981496322947`*^9}, {3.7573993424496107`*^9, 
   3.757399347323094*^9}, {3.757399476063964*^9, 3.75739947780792*^9}, {
   3.757401415453883*^9, 3.757401444829782*^9}, {3.757401556932979*^9, 
   3.757401590183528*^9}, {3.75740192886294*^9, 3.757401929846053*^9}, 
   3.757402167776972*^9, {3.7574022152695227`*^9, 3.7574022283358307`*^9}, {
   3.757402260899828*^9, 3.757402445473905*^9}, {3.757402564187311*^9, 
   3.75740256910693*^9}, {3.75740262917791*^9, 3.757402645140341*^9}, {
   3.7574027762243433`*^9, 3.757402783353739*^9}, {3.757402817032646*^9, 
   3.757402819721209*^9}, {3.757403106078498*^9, 3.7574031183701763`*^9}, {
   3.75740319359807*^9, 3.7574032188856077`*^9}, {3.7574091525588427`*^9, 
   3.757409154870949*^9}, {3.8019831824681797`*^9, 3.8019832067914248`*^9}, {
   3.801983239697904*^9, 3.801983241672761*^9}, {3.801983370480793*^9, 
   3.8019833917103777`*^9}, {3.801983446472941*^9, 3.801983533919042*^9}, {
   3.8020746871383104`*^9, 3.8020746896221056`*^9}, {3.8020747267696257`*^9, 
   3.802074828323705*^9}},
 ExpressionUUID -> "642fd32f-b58f-4afa-be6a-7d7291ce64f9"]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "The following function performs the inverse transformation when f(",
 Cell[BoxData[
  FormBox[
   SuperscriptBox["s", "*"], TraditionalForm]],
  FormatType->"TraditionalForm"],
 ")\[NotEqual]",
 Cell[BoxData[
  FormBox[
   SuperscriptBox["f", "*"], TraditionalForm]],
  FormatType->"TraditionalForm"],
 "(",
 Cell[BoxData[
  FormBox["s", TraditionalForm]],
  FormatType->"TraditionalForm"],
 ")"
}], "Section",
 CellChangeTimes->{{3.802089576079286*^9, 3.802089598393911*^9}, 
   3.940309003891079*^9, {3.940309229415779*^9, 3.9403092526945353`*^9}, {
   3.940309321584119*^9, 3.940309323786025*^9}, 3.940309912386386*^9, {
   3.9403108225636377`*^9, 3.940310830613632*^9}},
 ExpressionUUID -> "a8f6ac9c-4991-4d97-98ce-4892bbcd0b84"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"ILTC", "[", 
    RowBox[{"fun_", ",", " ", "T_", ",", " ", "maxFnEvals_", ",", " ", 
     RowBox[{"method_:", "\"\<cme\>\""}], ",", 
     RowBox[{"precision_:", "MachinePrecision"}]}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "params", ",", "i", ",", "\[Eta]", ",", "\[Beta]", ",", "\[Xi]", ",", 
       "n", ",", "k", ",", "prec"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"Length", "[", 
           RowBox[{"Position", "[", 
            RowBox[{
             RowBox[{"Names", "[", "\"\<Global`*\>\"", "]"}], ",", 
             "\"\<cmeParams\>\""}], "]"}], "]"}], "\[Equal]", "0"}], " ", "||",
          " ", 
         RowBox[{
          RowBox[{"Length", "[", "cmeParams", "]"}], "\[Equal]", "0"}]}], ",",
         "\[IndentingNewLine]", 
        RowBox[{"(*", 
         RowBox[{
          RowBox[{"Print", "[", "\"\<CME parameters loaded\>\"", "]"}], ";"}],
          "*)"}], "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"cmeParams", " ", "=", " ", 
          RowBox[{"Import", "[", 
           RowBox[{"FileNameJoin", "[", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"NotebookDirectory", "[", "]"}], ",", 
              "\"\<iltcme.json\>\""}], "}"}], "]"}], "]"}]}], ";"}]}], 
       "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"precision", "\[LessEqual]", "0"}], ",", 
        RowBox[{"prec", "=", "MachinePrecision"}], ",", 
        RowBox[{"prec", "=", "precision"}]}], "]"}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"method", "\[Equal]", "\"\<cme\>\""}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"params", " ", "=", " ", 
          RowBox[{"cmeParams", "[", 
           RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
         RowBox[{"For", "[", 
          RowBox[{
           RowBox[{"i", "=", "2"}], ",", 
           RowBox[{"i", "<", 
            RowBox[{"Length", "[", "cmeParams", "]"}]}], ",", 
           RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{
              RowBox[{
               RowBox[{"(", 
                RowBox[{"\"\<cv2\>\"", "/.", 
                 RowBox[{"cmeParams", "[", 
                  RowBox[{"[", "i", "]"}], "]"}]}], ")"}], " ", "<", 
               RowBox[{"(", 
                RowBox[{"\"\<cv2\>\"", "/.", "params"}], ")"}]}], " ", "&&", 
              " ", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"\"\<n\>\"", "+", "1"}], " ", "/.", 
                 RowBox[{"cmeParams", "[", 
                  RowBox[{"[", "i", "]"}], "]"}]}], ")"}], "\[LessEqual]", 
               "maxFnEvals"}]}], ",", 
             RowBox[{"params", "=", 
              RowBox[{"cmeParams", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}]}], " ", "]"}]}], 
          "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
         RowBox[{"\[Eta]", "=", 
          RowBox[{"SetPrecision", "[", 
           RowBox[{
            RowBox[{"Prepend", "[", 
             RowBox[{
              RowBox[{
               RowBox[{
                RowBox[{"\"\<a\>\"", "*", "\"\<mu1\>\""}], "+", 
                RowBox[{
                "\[ImaginaryI]", "*", "\"\<b\>\"", "*", "\"\<mu1\>\""}]}], "/.",
                "params"}], ",", 
              RowBox[{
               RowBox[{"\"\<c\>\"", "*", "\"\<mu1\>\""}], "/.", "params"}]}], 
             "]"}], ",", "prec"}], "]"}]}], ";", "\[IndentingNewLine]", 
         RowBox[{"\[Beta]", "=", 
          RowBox[{"SetPrecision", "[", 
           RowBox[{
            RowBox[{
             RowBox[{
              RowBox[{"Prepend", "[", 
               RowBox[{
                RowBox[{"1", "+", 
                 RowBox[{"\[ImaginaryI]", " ", "*", " ", "\"\<omega\>\"", "*", 
                  RowBox[{"Range", "[", 
                   RowBox[{"\"\<n\>\"", "/.", "params"}], "]"}]}]}], ",", 
                "1"}], "]"}], "*", "\"\<mu1\>\""}], " ", "/.", "params"}], 
            ",", "prec"}], "]"}]}], ";"}], "\[IndentingNewLine]", ",", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{"method", "\[Equal]", "\"\<euler\>\""}], ",", 
           "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"n", "=", 
             RowBox[{"If", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"Mod", "[", 
                 RowBox[{"maxFnEvals", ",", "2"}], "]"}], "\[Equal]", "0"}], 
               ",", 
               RowBox[{"maxFnEvals", "-", "1"}], ",", "maxFnEvals"}], "]"}]}],
             ";", "\[IndentingNewLine]", 
            RowBox[{"\[Beta]", "=", 
             RowBox[{"Table", "[", 
              RowBox[{
               RowBox[{"SetPrecision", "[", 
                RowBox[{
                 RowBox[{
                  FractionBox[
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"n", "-", "1"}], ")"}], 
                    RowBox[{"Log", "[", "10", "]"}]}], "6"], "+", 
                  RowBox[{"\[Pi]", " ", "\[ImaginaryI]", " ", 
                   RowBox[{"(", 
                    RowBox[{"k", "-", "1"}], ")"}]}]}], ",", "prec"}], "]"}], 
               ",", 
               RowBox[{"{", 
                RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";", 
            "\[IndentingNewLine]", 
            RowBox[{"\[Xi]", "=", 
             RowBox[{"ConstantArray", "[", 
              RowBox[{
               RowBox[{"SetPrecision", "[", 
                RowBox[{"1", ",", "prec"}], "]"}], ",", "n"}], "]"}]}], ";", 
            "\[IndentingNewLine]", 
            RowBox[{
             RowBox[{"\[Xi]", "[", 
              RowBox[{"[", "1", "]"}], "]"}], " ", "=", " ", 
             RowBox[{"SetPrecision", "[", 
              RowBox[{
               RowBox[{"1", "/", "2"}], ",", "prec"}], "]"}]}], ";", 
            "\[IndentingNewLine]", 
            RowBox[{
             RowBox[{"\[Xi]", "[", 
              RowBox[{"[", "n", "]"}], "]"}], "=", 
             RowBox[{"SetPrecision", "[", 
              RowBox[{
               FractionBox["1", 
                SuperscriptBox["2", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"n", "-", "1"}], ")"}], "/", "2"}]]], ",", 
               "prec"}], "]"}]}], ";", "\[IndentingNewLine]", 
            RowBox[{"For", "[", 
             RowBox[{
              RowBox[{"k", "=", "1"}], ",", 
              RowBox[{"k", "<", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"n", "-", "1"}], ")"}], "/", "2"}]}], ",", 
              RowBox[{"k", "++"}], ",", 
              RowBox[{
               RowBox[{"\[Xi]", "[", 
                RowBox[{"[", 
                 RowBox[{"n", "-", "k"}], "]"}], "]"}], "=", 
               RowBox[{
                RowBox[{"\[Xi]", "[", 
                 RowBox[{"[", 
                  RowBox[{"n", "-", "k", "+", "1"}], "]"}], "]"}], "+", " ", 
                RowBox[{
                 SuperscriptBox["2", 
                  RowBox[{
                   RowBox[{"-", 
                    RowBox[{"(", 
                    RowBox[{"n", "-", "1"}], ")"}]}], "/", "2"}]], 
                 RowBox[{"Binomial", "[", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"n", "-", "1"}], ")"}], "/", "2"}], ",", "k"}], 
                  "]"}]}]}]}]}], "]"}], ";", "\[IndentingNewLine]", 
            RowBox[{"\[Eta]", "=", 
             RowBox[{"Table", "[", 
              RowBox[{
               RowBox[{"SetPrecision", "[", 
                RowBox[{
                 RowBox[{
                  SuperscriptBox["10", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"n", "-", "1"}], ")"}], "/", "6"}]], 
                  SuperscriptBox[
                   RowBox[{"(", 
                    RowBox[{"-", "1"}], ")"}], 
                   RowBox[{"k", "-", "1"}]], 
                  RowBox[{"\[Xi]", "[", 
                   RowBox[{"[", "k", "]"}], "]"}]}], ",", "prec"}], "]"}], 
               ",", 
               RowBox[{"{", 
                RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
           "\[IndentingNewLine]", ",", "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"If", "[", 
             RowBox[{
              RowBox[{"method", "\[Equal]", "\"\<gaver\>\""}], ",", 
              "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{"n", "=", 
                RowBox[{"If", "[", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"Mod", "[", 
                    RowBox[{"maxFnEvals", ",", "2"}], "]"}], "\[Equal]", 
                   "1"}], ",", 
                  RowBox[{"maxFnEvals", "-", "1"}], ",", "maxFnEvals"}], 
                 "]"}]}], ";", "\[IndentingNewLine]", 
               RowBox[{"\[Beta]", "=", 
                RowBox[{"Table", "[", 
                 RowBox[{
                  RowBox[{"k", " ", 
                   RowBox[{"SetPrecision", "[", 
                    RowBox[{
                    RowBox[{"Log", "[", "2", "]"}], ",", "prec"}], "]"}]}], 
                  ",", 
                  RowBox[{"{", 
                   RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";", 
               "\[IndentingNewLine]", 
               RowBox[{"\[Eta]", "=", 
                RowBox[{"Table", "[", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"SetPrecision", "[", 
                    RowBox[{
                    RowBox[{"Log", "[", "2", "]"}], ",", "prec"}], "]"}], 
                   SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"-", "1"}], ")"}], 
                    RowBox[{
                    RowBox[{"n", "/", "2"}], "+", "k"}]], 
                   RowBox[{"Sum", "[", 
                    RowBox[{
                    RowBox[{
                    FractionBox[
                    SuperscriptBox["j", 
                    RowBox[{
                    RowBox[{"n", "/", "2"}], "+", "1"}]], 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{"n", "/", "2"}], ")"}], "!"}]], 
                    RowBox[{"Binomial", "[", 
                    RowBox[{
                    RowBox[{"n", "/", "2"}], ",", "j"}], "]"}], 
                    RowBox[{"Binomial", "[", 
                    RowBox[{
                    RowBox[{"2", "j"}], ",", "j"}], "]"}], 
                    RowBox[{"Binomial", "[", 
                    RowBox[{"j", ",", 
                    RowBox[{"k", "-", "j"}]}], "]"}]}], ",", 
                    RowBox[{"{", 
                    RowBox[{"j", ",", 
                    RowBox[{"Floor", "[", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{"k", "+", "1"}], ")"}], "/", "2"}], "]"}], ",", 
                    RowBox[{"Min", "[", 
                    RowBox[{"k", ",", 
                    RowBox[{"n", "/", "2"}]}], "]"}]}], "}"}]}], "]"}]}], 
                  "  ", ",", 
                  RowBox[{"{", 
                   RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
               ";"}]}], "\[IndentingNewLine]", "]"}], ";"}]}], 
          "\[IndentingNewLine]", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}],
       ";", "\[IndentingNewLine]", 
      RowBox[{"Return", "[", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{"ListQ", "[", "T", "]"}], ",", 
         RowBox[{"Table", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"\[Eta]", "[", 
                RowBox[{"[", "1", "]"}], "]"}], "*", 
               RowBox[{"fun", "[", 
                RowBox[{
                 RowBox[{"\[Beta]", "[", 
                  RowBox[{"[", "1", "]"}], "]"}], "/", "x"}], "]"}]}], "+", 
              RowBox[{"Sum", "[", 
               RowBox[{
                RowBox[{
                 RowBox[{"\[Eta]", "[", 
                  RowBox[{"[", "k", "]"}], "]"}], "*", 
                 RowBox[{
                  RowBox[{"fun", "[", 
                   RowBox[{
                    RowBox[{"\[Beta]", "[", 
                    RowBox[{"[", "k", "]"}], "]"}], "/", "x"}], "]"}], "/", 
                  "2"}]}], ",", 
                RowBox[{"{", 
                 RowBox[{"k", ",", "2", ",", 
                  RowBox[{"Length", "[", "\[Eta]", "]"}]}], "}"}]}], "]"}], 
              "+", 
              RowBox[{"Sum", "[", 
               RowBox[{
                RowBox[{
                 RowBox[{"Conjugate", "[", 
                  RowBox[{"\[Eta]", "[", 
                   RowBox[{"[", "k", "]"}], "]"}], "]"}], "*", 
                 RowBox[{
                  RowBox[{"fun", "[", 
                   RowBox[{
                    RowBox[{"Conjugate", "[", 
                    RowBox[{"\[Beta]", "[", 
                    RowBox[{"[", "k", "]"}], "]"}], "]"}], "/", "x"}], "]"}], 
                  "/", "2"}]}], ",", 
                RowBox[{"{", 
                 RowBox[{"k", ",", "2", ",", 
                  RowBox[{"Length", "[", "\[Eta]", "]"}]}], "}"}]}], "]"}]}], 
             ")"}], "/", "x"}], ",", 
           RowBox[{"{", 
            RowBox[{"x", ",", "T"}], "}"}]}], "]"}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"\[Eta]", "[", 
               RowBox[{"[", "1", "]"}], "]"}], "*", 
              RowBox[{"fun", "[", 
               RowBox[{
                RowBox[{"\[Beta]", "[", 
                 RowBox[{"[", "1", "]"}], "]"}], "/", "x"}], "]"}]}], "+", 
             RowBox[{"Sum", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"\[Eta]", "[", 
                 RowBox[{"[", "k", "]"}], "]"}], "*", 
                RowBox[{
                 RowBox[{"fun", "[", 
                  RowBox[{
                   RowBox[{"\[Beta]", "[", 
                    RowBox[{"[", "k", "]"}], "]"}], "/", "x"}], "]"}], "/", 
                 "2"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"k", ",", "2", ",", 
                 RowBox[{"Length", "[", "\[Eta]", "]"}]}], "}"}]}], "]"}], 
             "+", 
             RowBox[{"Sum", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"Conjugate", "[", 
                 RowBox[{"\[Eta]", "[", 
                  RowBox[{"[", "k", "]"}], "]"}], "]"}], "*", 
                RowBox[{
                 RowBox[{"fun", "[", 
                  RowBox[{
                   RowBox[{"Conjugate", "[", 
                    RowBox[{"\[Beta]", "[", 
                    RowBox[{"[", "k", "]"}], "]"}], "]"}], "/", "x"}], "]"}], 
                 "/", "2"}]}], ",", 
               RowBox[{"{", 
                RowBox[{"k", ",", "2", ",", 
                 RowBox[{"Length", "[", "\[Eta]", "]"}]}], "}"}]}], "]"}]}], 
            ")"}], "/", "x"}], "/.", 
          RowBox[{"x", "\[Rule]", "T"}]}]}], "]"}], "]"}], ";"}]}], 
    "\[IndentingNewLine]", "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.757325469151745*^9, 3.757325502280724*^9}, {
   3.75732570552604*^9, 3.757325767957542*^9}, {3.757325851475429*^9, 
   3.7573259078775578`*^9}, {3.7573259980687437`*^9, 3.75732607799338*^9}, {
   3.75732612019339*^9, 3.757326145684031*^9}, {3.7573261929718733`*^9, 
   3.7573262563262873`*^9}, 3.757396882026514*^9, {3.75739691295605*^9, 
   3.757397014808332*^9}, {3.7573970726919823`*^9, 3.757397221175847*^9}, {
   3.757397601976653*^9, 3.757397666846058*^9}, {3.7573980113072147`*^9, 
   3.757398027672312*^9}, {3.757398063947167*^9, 3.757398074906393*^9}, {
   3.757398144923143*^9, 3.7573981496322947`*^9}, {3.7573993424496107`*^9, 
   3.757399347323094*^9}, {3.757399476063964*^9, 3.75739947780792*^9}, {
   3.757401415453883*^9, 3.757401444829782*^9}, {3.757401556932979*^9, 
   3.757401590183528*^9}, {3.75740192886294*^9, 3.757401929846053*^9}, 
   3.757402167776972*^9, {3.7574022152695227`*^9, 3.7574022283358307`*^9}, {
   3.757402260899828*^9, 3.757402445473905*^9}, {3.757402564187311*^9, 
   3.75740256910693*^9}, {3.75740262917791*^9, 3.757402645140341*^9}, {
   3.7574027762243433`*^9, 3.757402783353739*^9}, {3.757402817032646*^9, 
   3.757402819721209*^9}, {3.757403106078498*^9, 3.7574031183701763`*^9}, {
   3.75740319359807*^9, 3.7574032188856077`*^9}, {3.7574091525588427`*^9, 
   3.757409154870949*^9}, {3.8019831824681797`*^9, 3.8019832067914248`*^9}, {
   3.801983239697904*^9, 3.801983241672761*^9}, {3.801983370480793*^9, 
   3.8019833917103777`*^9}, {3.801983446472941*^9, 3.801983533919042*^9}, {
   3.8020746871383104`*^9, 3.8020746896221056`*^9}, {3.8020747267696257`*^9, 
   3.802074828323705*^9}, {3.940309536119437*^9, 3.940309555345018*^9}, {
   3.9403096084159727`*^9, 3.9403096413859324`*^9}, {3.940310501487999*^9, 
   3.9403105565290623`*^9}, {3.9403106194325895`*^9, 3.940310688230997*^9}},
 ExpressionUUID -> "642fd32f-b58f-4afa-be6a-7d7291ce64f9"]
}, Open  ]],

Cell[CellGroupData[{

Cell["A couple of simple examples", "Section",
 CellChangeTimes->{{3.80208963825071*^9, 3.802089670428808*^9}, {
  3.802089717015321*^9, 
  3.802089721271408*^9}},ExpressionUUID->"5628899b-611e-4983-9b78-\
28bab36e6fdd"],

Cell[CellGroupData[{

Cell["Staircase function at a single point", "Subsection",
 CellChangeTimes->{{3.802089735600309*^9, 
  3.802089745128953*^9}},ExpressionUUID->"c1aed1d9-7428-4139-b73a-\
63783ce34d50"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"h", "[", "t_", "]"}], ":=", " ", 
    RowBox[{"Floor", "[", "t", "]"}]}], " ", ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"time", " ", "domain", " ", "function"}], " ", 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"H", "[", "s_", "]"}], ":=", 
    RowBox[{
     FractionBox["1", "s"], " ", 
     FractionBox["1", 
      RowBox[{
       SuperscriptBox["\[ExponentialE]", "s"], "-", "1"}]]}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"Laplace", " ", "transform", " ", "domain", " ", "function"}], " ",
    "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"n", "=", "50"}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"order", " ", "=", " ", 
    RowBox[{
    "maximum", " ", "number", " ", "of", " ", "function", " ", 
     "evaluations"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"timepoint", "=", "2.5"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"precision", " ", "=", " ", "MachinePrecision"}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"optional", " ", 
    RowBox[{"parameter", ":", " ", 
     RowBox[{
     "number", " ", "of", " ", "precise", " ", "digits", " ", "during", " ", 
      "calculation", " ", "or", " ", "machine", " ", "precision"}]}]}], " ", 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"Perform", " ", 
     RowBox[{"ILT", "!"}], " ", "The", " ", "exact", " ", "result", " ", "is",
      " ", 
     RowBox[{"2.0", ".", " ", "Increasing"}], " ", "\"\<n\>\"", " ", 
     "parameter", " ", "improves", " ", "the", " ", "accuracy"}], ",", " ", 
    RowBox[{
    "at", " ", "least", " ", "in", " ", "case", " ", "of", " ", "the", " ", 
     "cme", " ", "method"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ILT", "[", 
   RowBox[{
   "H", ",", "timepoint", ",", "n", ",", "\"\<euler\>\"", ",", "precision"}], 
   "]"}], " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ILT", "[", 
   RowBox[{
   "H", ",", "timepoint", ",", "n", ",", "\"\<gaver\>\"", ",", "precision"}], 
   "]"}], " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ILT", "[", 
   RowBox[{
   "H", ",", "timepoint", ",", "n", ",", "\"\<cme\>\"", ",", "precision"}], 
   "]"}], " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ILTC", "[", 
   RowBox[{
   "H", ",", "timepoint", ",", "n", ",", "\"\<euler\>\"", ",", "precision"}], 
   "]"}], " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ILTC", "[", 
   RowBox[{
   "H", ",", "timepoint", ",", "n", ",", "\"\<gaver\>\"", ",", "precision"}], 
   "]"}], " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ILTC", "[", 
   RowBox[{
   "H", ",", "timepoint", ",", "n", ",", "\"\<cme\>\"", ",", "precision"}], 
   "]"}], " "}]}], "Input",
 CellChangeTimes->{{3.802074625855833*^9, 3.8020746683458505`*^9}, {
  3.802089673581768*^9, 3.802089674797903*^9}, {3.8020897907097797`*^9, 
  3.802089950789777*^9}, {3.80209016049277*^9, 3.802090170617971*^9}, {
  3.9403099426791377`*^9, 3.9403099498195963`*^9}},
 ExpressionUUID -> "c30545ca-cce9-4ee2-bda3-d42fc1e887de"],

Cell[BoxData["2.0968146963059553`"], "Output",
 CellChangeTimes->{{3.9403099266257296`*^9, 3.940309950560197*^9}, 
   3.940310713696581*^9}],

Cell[BoxData["1.099511627776`*^11"], "Output",
 CellChangeTimes->{{3.9403099266257296`*^9, 3.940309950560197*^9}, 
   3.940310713710579*^9}],

Cell[BoxData["1.9999780120225161`"], "Output",
 CellChangeTimes->{{3.9403099266257296`*^9, 3.940309950560197*^9}, 
   3.940310713731578*^9}],

Cell[BoxData[
 RowBox[{"2.0968146963057737`", "\[VeryThinSpace]", "-", 
  RowBox[{"3.353761712787673`*^-13", " ", "\[ImaginaryI]"}]}]], "Output",
 CellChangeTimes->{{3.9403099266257296`*^9, 3.940309950560197*^9}, 
   3.940310713745577*^9}],

Cell[BoxData["1.099511627776`*^11"], "Output",
 CellChangeTimes->{{3.9403099266257296`*^9, 3.940309950560197*^9}, 
   3.9403107137665777`*^9}],

Cell[BoxData[
 RowBox[{"1.9999780120225166`", "\[VeryThinSpace]", "+", 
  RowBox[{"1.4210854715202005`*^-15", " ", "\[ImaginaryI]"}]}]], "Output",
 CellChangeTimes->{{3.9403099266257296`*^9, 3.940309950560197*^9}, 
   3.9403107137665777`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Staircase function at multiple points", "Subsection",
 CellChangeTimes->{{3.802089754041147*^9, 
  3.802089763473835*^9}},ExpressionUUID->"8d04b0f1-dce2-42fb-b72e-\
e54c0907938d"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"h", "[", "t_", "]"}], ":=", " ", 
    RowBox[{"Floor", "[", "t", "]"}]}], " ", ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"time", " ", "domain", " ", "function"}], " ", 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"H", "[", "s_", "]"}], ":=", 
    RowBox[{
     FractionBox["1", "s"], " ", 
     FractionBox["1", 
      RowBox[{
       SuperscriptBox["\[ExponentialE]", "s"], "-", "1"}]]}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"Laplace", " ", "transform", " ", "domain", " ", "function"}], " ",
    "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"points", " ", "=", " ", 
    RowBox[{"Table", "[", 
     RowBox[{"t", ",", 
      RowBox[{"{", 
       RowBox[{"t", ",", "0.01", ",", "5", ",", "0.01"}], "}"}]}], "]"}]}], 
   ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "list", " ", "of", " ", "time", " ", "points", " ", "where", " ", "the", 
    " ", "ILT", " ", "is", " ", "to", " ", "be", " ", "evaluated"}], " ", 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"n", "=", "30"}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"order", " ", "=", " ", 
    RowBox[{
    "maximum", " ", "number", " ", "of", " ", "function", " ", 
     "evaluations"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"precision", " ", "=", " ", "MachinePrecision"}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"optional", " ", 
    RowBox[{"parameter", ":", " ", 
     RowBox[{
     "number", " ", "of", " ", "precise", " ", "digits", " ", "during", " ", 
      "calculation", " ", "or", " ", "machine", " ", "precision"}]}]}], " ", 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"original", " ", "=", " ", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"h", "[", "t", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0.01", ",", "5", ",", "0.01"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ilteuler", " ", "=", " ", 
   RowBox[{"ILT", "[", 
    RowBox[{
    "H", ",", "points", ",", "n", ",", "\"\<euler\>\"", ",", "precision"}], 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"iltcme", " ", "=", " ", 
   RowBox[{"ILT", "[", 
    RowBox[{
    "H", ",", "points", ",", "n", ",", "\"\<cme\>\"", ",", "precision"}], 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"ListLinePlot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"original", ",", " ", "ilteuler", ",", " ", "iltcme"}], "}"}], 
   ",", 
   RowBox[{"PlotLegends", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<original\>\"", ",", "\"\<Euler\>\"", ",", "\"\<CME\>\""}], 
     "}"}]}], ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<order 30\>\""}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.8020697997642655`*^9, 3.8020699105040226`*^9}, {
   3.8020699534406548`*^9, 3.802070062848937*^9}, {3.8020701255951986`*^9, 
   3.8020701919391794`*^9}, {3.802070226852863*^9, 3.8020703212075124`*^9}, {
   3.802070359731757*^9, 3.8020705177260466`*^9}, {3.8020741369744697`*^9, 
   3.8020741535020375`*^9}, {3.8020742152062206`*^9, 3.802074318693603*^9}, {
   3.80207441499911*^9, 3.80207456596365*^9}, {3.8020746046891236`*^9, 
   3.802074616483035*^9}, 3.8020748602850447`*^9, {3.802089970192419*^9, 
   3.802089971376786*^9}, {3.802090022680962*^9, 3.802090023458226*^9}, {
   3.802090226320466*^9, 3.802090226749135*^9}},
 ExpressionUUID -> "c4e58486-0879-47e3-842c-b44a8ccf623c"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"n", "=", "80"}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "let", " ", "us", " ", "increase", " ", "the", " ", "number", " ", "of", 
    " ", "function", " ", "evaluations"}], " ", 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ilteuler", " ", "=", " ", 
   RowBox[{"ILT", "[", 
    RowBox[{
    "H", ",", "points", ",", "n", ",", "\"\<euler\>\"", ",", "precision"}], 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"iltcme", " ", "=", " ", 
   RowBox[{"ILT", "[", 
    RowBox[{
    "H", ",", "points", ",", "n", ",", "\"\<cme\>\"", ",", "precision"}], 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"ListLinePlot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"original", ",", " ", "ilteuler", ",", " ", "iltcme"}], "}"}], 
   ",", 
   RowBox[{"PlotLegends", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<original\>\"", ",", "\"\<Euler\>\"", ",", "\"\<CME\>\""}], 
     "}"}]}], ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<order 80\>\""}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.802090033658668*^9, 3.802090036674446*^9}, {
  3.8020902168813343`*^9, 3.802090217364505*^9}},
 ExpressionUUID -> "147b2e1c-9a55-4095-b27c-d0ba99914767"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Interactive demo", "Section",
 CellChangeTimes->{{3.8020900642908373`*^9, 
  3.802090066514724*^9}},ExpressionUUID->"68a85931-887d-4cdf-88c1-\
2bb4ed3d40df"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "Function", " ", "definition", " ", "for", " ", "interactive", " ", 
    "plot"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"tFuns", "=", 
    RowBox[{"<|", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"exponential", "\[Rule]", 
       RowBox[{"<|", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"ilt", "\[Rule]", 
          RowBox[{"Function", "[", 
           RowBox[{"x", ",", 
            RowBox[{"Exp", "[", 
             RowBox[{"-", "x"}], "]"}]}], "]"}]}], ",", "\[IndentingNewLine]", 
         RowBox[{"lt", "\[Rule]", 
          RowBox[{"Function", "[", 
           RowBox[{"s", ",", 
            FractionBox["1", 
             RowBox[{"1", "+", "s"}]]}], "]"}]}], ",", "\[IndentingNewLine]", 
         
         RowBox[{"xvals", "\[Rule]", 
          RowBox[{"Range", "[", 
           RowBox[{
            RowBox[{"1", "/", "10"}], ",", "5", ",", 
            RowBox[{"1", "/", "100"}]}], "]"}]}], ",", "\[IndentingNewLine]", 
         
         RowBox[{"yrange", "\[Rule]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "0.5"}], ",", "1.5"}], "}"}]}]}], 
        "\[IndentingNewLine]", "|>"}]}], ",", " ", "\[IndentingNewLine]", 
      RowBox[{"sine", "\[Rule]", 
       RowBox[{"<|", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"ilt", "\[Rule]", 
          RowBox[{"Function", "[", 
           RowBox[{"x", ",", 
            RowBox[{"Sin", "[", "x", "]"}]}], "]"}]}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{"lt", "\[Rule]", 
          RowBox[{"Function", "[", 
           RowBox[{"s", ",", 
            FractionBox["1", 
             RowBox[{"1", "+", 
              SuperscriptBox["s", "2"]}]]}], "]"}]}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{"xvals", "\[Rule]", 
          RowBox[{"Range", "[", 
           RowBox[{
            RowBox[{"1", "/", "10"}], ",", "15", ",", 
            RowBox[{"1", "/", "20"}]}], "]"}]}], ",", "\[IndentingNewLine]", 
         RowBox[{"yrange", "\[Rule]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "1.2"}], ",", "1.2"}], "}"}]}]}], 
        "\[IndentingNewLine]", "|>"}]}], ",", "\[IndentingNewLine]", 
      RowBox[{"heavyside", "\[Rule]", 
       RowBox[{"<|", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"ilt", "\[Rule]", 
          RowBox[{"Function", "[", 
           RowBox[{"x", ",", 
            RowBox[{"If", "[", 
             RowBox[{
              RowBox[{"x", ">", "1"}], ",", "1", ",", "0"}], "]"}]}], "]"}]}],
          ",", "\[IndentingNewLine]", 
         RowBox[{"lt", "\[Rule]", 
          RowBox[{"Function", "[", 
           RowBox[{"s", ",", 
            RowBox[{
             SuperscriptBox["\[ExponentialE]", 
              RowBox[{"-", "s"}]], "/", "s"}]}], "]"}]}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{"xvals", "\[Rule]", 
          RowBox[{"Range", "[", 
           RowBox[{
            RowBox[{"1", "/", "10"}], ",", "3", ",", 
            RowBox[{"1", "/", "100"}]}], "]"}]}], ",", "\[IndentingNewLine]", 
         
         RowBox[{"yrange", "\[Rule]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "0.5"}], ",", "1.5"}], "}"}]}]}], 
        "\[IndentingNewLine]", "|>"}]}], ",", "\[IndentingNewLine]", 
      RowBox[{"expheavyside", "\[Rule]", 
       RowBox[{"<|", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"ilt", "\[Rule]", 
          RowBox[{"Function", "[", 
           RowBox[{"x", ",", 
            RowBox[{
             RowBox[{"If", "[", 
              RowBox[{
               RowBox[{"x", ">", "1"}], ",", "1", ",", "0"}], "]"}], "*", 
             RowBox[{"Exp", "[", 
              RowBox[{"1", "-", "x"}], "]"}]}]}], "]"}]}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{"lt", "\[Rule]", 
          RowBox[{"Function", "[", 
           RowBox[{"s", ",", 
            RowBox[{
             SuperscriptBox["\[ExponentialE]", 
              RowBox[{"-", "s"}]], "/", 
             RowBox[{"(", 
              RowBox[{"1", "+", "s"}], ")"}]}]}], "]"}]}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{"xvals", "\[Rule]", 
          RowBox[{"Range", "[", 
           RowBox[{
            RowBox[{"1", "/", "10"}], ",", "5", ",", 
            RowBox[{"1", "/", "100"}]}], "]"}]}], ",", "\[IndentingNewLine]", 
         
         RowBox[{"yrange", "\[Rule]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "0.5"}], ",", "1.2"}], "}"}]}]}], 
        "\[IndentingNewLine]", "|>"}]}], ",", "\[IndentingNewLine]", 
      RowBox[{"squarewave", "\[Rule]", 
       RowBox[{"<|", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"ilt", "\[Rule]", 
          RowBox[{"Function", "[", 
           RowBox[{"x", ",", 
            RowBox[{"Mod", "[", 
             RowBox[{
              RowBox[{"Floor", "[", "x", "]"}], ",", "2"}], "]"}]}], "]"}]}], 
         ",", "\[IndentingNewLine]", 
         RowBox[{"lt", "\[Rule]", 
          RowBox[{"Function", "[", 
           RowBox[{"s", ",", 
            RowBox[{
             FractionBox["1", "s"], 
             FractionBox["1", 
              RowBox[{
               SuperscriptBox["\[ExponentialE]", "s"], "+", "1"}]]}]}], 
           "]"}]}], ",", "\[IndentingNewLine]", 
         RowBox[{"xvals", "\[Rule]", 
          RowBox[{"Range", "[", 
           RowBox[{
            RowBox[{"1", "/", "10"}], ",", "10", ",", 
            RowBox[{"1", "/", "100"}]}], "]"}]}], ",", "\[IndentingNewLine]", 
         
         RowBox[{"yrange", "\[Rule]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "1"}], ",", "2"}], "}"}]}]}], "\[IndentingNewLine]", 
        "|>"}]}], ",", "\[IndentingNewLine]", 
      RowBox[{"staircase", "\[Rule]", 
       RowBox[{"<|", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"ilt", "\[Rule]", 
          RowBox[{"Function", "[", 
           RowBox[{"x", ",", 
            RowBox[{"Floor", "[", "x", "]"}]}], "]"}]}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{"lt", "\[Rule]", 
          RowBox[{"Function", "[", 
           RowBox[{"s", ",", 
            RowBox[{
             FractionBox["1", "s"], 
             FractionBox["1", 
              RowBox[{
               SuperscriptBox["\[ExponentialE]", "s"], "-", "1"}]]}]}], 
           "]"}]}], ",", "\[IndentingNewLine]", 
         RowBox[{"xvals", "\[Rule]", 
          RowBox[{"Range", "[", 
           RowBox[{
            RowBox[{"1", "/", "10"}], ",", "5", ",", 
            RowBox[{"1", "/", "100"}]}], "]"}]}], ",", "\[IndentingNewLine]", 
         
         RowBox[{"yrange", "\[Rule]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "1"}], ",", "5"}], "}"}]}]}], "\[IndentingNewLine]", 
        "|>"}]}]}], "\[IndentingNewLine]", "|>"}]}], ";"}]}]], "Input",
 CellChangeTimes->{{3.757403719118678*^9, 3.7574037284295692`*^9}, {
  3.757405110766655*^9, 3.7574051141269503`*^9}, {3.757405223006798*^9, 
  3.7574052339533043`*^9}, {3.757405270553302*^9, 3.757405330117991*^9}, {
  3.757405373336069*^9, 3.7574055327190723`*^9}, {3.757405906505804*^9, 
  3.757405935906724*^9}, {3.757405993003598*^9, 3.757406017545661*^9}, {
  3.757406091165361*^9, 3.757406132745515*^9}, {3.757408001760363*^9, 
  3.757408011176996*^9}, {3.757408041620781*^9, 3.757408113191677*^9}, {
  3.757408315515572*^9, 3.757408345522232*^9}, {3.7574085815394087`*^9, 
  3.757408581888969*^9}, {3.757421832722377*^9, 3.757421857121273*^9}, {
  3.757421889275429*^9, 3.757421928219496*^9}, {3.802074098796134*^9, 
  3.8020741046070733`*^9}, {3.802074920708456*^9, 3.80207492381693*^9}},
 ExpressionUUID -> "65ca5cd4-3889-402e-90db-402394ef358c"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"Interactive", " ", "plot"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Manipulate", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"precision", "=", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"highPrecision", "\[Equal]", "0"}], ",", "MachinePrecision", 
        ",", " ", "200"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"euler", " ", "=", " ", 
      RowBox[{"ILT", "[", 
       RowBox[{
        RowBox[{"tFuns", "[", 
         RowBox[{"testFunction", ",", "lt"}], "]"}], ",", 
        RowBox[{"tFuns", "[", 
         RowBox[{"testFunction", ",", "xvals"}], "]"}], ",", " ", "funEvals", 
        ",", " ", "\"\<euler\>\"", ",", "precision"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"cme", " ", "=", " ", 
      RowBox[{"ILT", "[", 
       RowBox[{
        RowBox[{"tFuns", "[", 
         RowBox[{"testFunction", ",", "lt"}], "]"}], ",", 
        RowBox[{"tFuns", "[", 
         RowBox[{"testFunction", ",", "xvals"}], "]"}], ",", " ", "funEvals", 
        ",", " ", "\"\<cme\>\"", ",", "precision"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"gaver", "=", 
      RowBox[{"ILT", "[", 
       RowBox[{
        RowBox[{"tFuns", "[", 
         RowBox[{"testFunction", ",", "lt"}], "]"}], ",", 
        RowBox[{"tFuns", "[", 
         RowBox[{"testFunction", ",", "xvals"}], "]"}], ",", " ", "funEvals", 
        ",", " ", "\"\<gaver\>\"", ",", "precision"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"exact", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"tFuns", "[", 
          RowBox[{"testFunction", ",", "ilt"}], "]"}], "[", "x", "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"x", ",", 
          RowBox[{"tFuns", "[", 
           RowBox[{"testFunction", ",", "xvals"}], "]"}]}], "}"}]}], "]"}]}], 
     ";", "\[IndentingNewLine]", 
     RowBox[{"ListLinePlot", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"exact", ",", "cme", ",", "euler", ",", " ", "gaver"}], "}"}],
        ",", " ", 
       RowBox[{"PlotLegends", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{
         "\"\<Exact\>\"", ",", "\"\<CME\>\"", ",", "\"\<Euler\>\"", ",", 
          "\"\<Gaver\>\""}], "}"}]}], ",", 
       RowBox[{"PlotRange", "->", 
        RowBox[{"tFuns", "[", 
         RowBox[{"testFunction", ",", "yrange"}], "]"}]}], ",", 
       RowBox[{"ImageSize", "\[Rule]", "480"}]}], "]"}]}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"funEvals", ",", "50"}], "}"}], ",", "3", ",", "200", ",", 
      "1"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"highPrecision", ",", "1"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "1"}], "}"}], ",", "Checkbox"}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"testFunction", ",", 
      RowBox[{"{", 
       RowBox[{
       "staircase", ",", "exponential", ",", "sine", ",", " ", "heavyside", 
        ",", "expheavyside", ",", "squarewave"}], "}"}]}], "}"}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->CompressedData["
1:eJxTTMoPSmViYGAQBWIQzZM+PXntv9eOv371NYJo7hVvl4DoIOPfa0H0b/vl
m0D0MstnYNrHaOJOEO1iKrgLRAdG6+4H0YaH1oJpBUPN+nVAWnXy9SYQvd98
dRuInhHY3AGi9+xTXg+iLZ4u3QeiT0mtOgqibwjsOgGiD6x/dwZEexzYfRZE
m/BeugSit0R5XwXRTVHfXoDoLHGltyA6qu7vNxBdxO34HUS/YpF9ux5Ip0Ue
eA+iI07t494A8k/xUT4QreQWJwyiv3+NFgGLH9aQBtEVOTpguk7DXx1EZ/cd
BNPrIngSQDSPUEAGiO7ab71uG8j9FRkbQPSR+ZeePgx+48gg/xBMAwBJG8xx

  "],
 ExpressionUUID -> "689cd0fe-f194-47a6-8a3d-6fd5d464b278"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1680, 972},
WindowMargins->{{0, Automatic}, {0, Automatic}},
Magnification:>1.5 Inherited,
FrontEndVersion->"11.0 for Microsoft Windows (64-bit) (September 21, 2016)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 249, 6, 246, "Title"],
Cell[CellGroupData[{
Cell[854, 32, 746, 21, 106, "Section"],
Cell[1603, 55, 15062, 352, 1255, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16702, 412, 751, 21, 106, "Section"],
Cell[17456, 435, 17429, 408, 1371, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[34922, 848, 220, 4, 106, "Section"],
Cell[CellGroupData[{
Cell[35167, 856, 184, 3, 72, "Subsection"],
Cell[CellGroupData[{
Cell[35376, 863, 3134, 85, 395, "Input"],
Cell[38513, 950, 140, 2, 45, "Output"],
Cell[38656, 954, 140, 2, 45, "Output"],
Cell[38799, 958, 140, 2, 45, "Output"],
Cell[38942, 962, 239, 4, 45, "Output"],
Cell[39184, 968, 142, 2, 45, "Output"],
Cell[39329, 972, 242, 4, 45, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[39620, 982, 185, 3, 72, "Subsection"],
Cell[39808, 987, 3487, 91, 308, "Input"],
Cell[43298, 1080, 1233, 33, 133, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[44580, 1119, 163, 3, 106, "Section"],
Cell[44746, 1124, 7723, 190, 1280, "Input"],
Cell[52472, 1316, 3663, 93, 365, "Input"]
}, Open  ]]
}, Open  ]]
}
]
*)