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

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

(* CreatedBy='Mathematica 10.4' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[      9181,        230]
NotebookOptionsPosition[      8645,        206]
NotebookOutlinePosition[      9052,        224]
CellTagsIndexPosition[      9009,        221]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Exercise 1.82", "Subsubsection",
 CellChangeTimes->{{3.791440344522746*^9, 3.7914403637039137`*^9}, {
  3.791440528281019*^9, 3.791440539308681*^9}, {3.791440593421216*^9, 
  3.7914406129403152`*^9}, {3.791709135382835*^9, 3.791709145103059*^9}, {
  3.79170953911546*^9, 3.791709551345149*^9}, {3.792399815781332*^9, 
  3.79239982712624*^9}, {3.792418754873685*^9, 3.792418764520961*^9}, {
  3.825435357666072*^9, 3.8254353598893642`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"tsim", "=", "10"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Lambda]", "[", "t_", "]"}], "=", 
   RowBox[{"1", "/", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"t", "^", "2"}], "+", "1"}], ")"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"path", "=", 
   RowBox[{"RandomFunction", "[", 
    RowBox[{
     RowBox[{"InhomogeneousPoissonProcess", "[", 
      RowBox[{
       RowBox[{"\[Lambda]", "[", "t", "]"}], ",", "t"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "tsim"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"f", "=", 
   RowBox[{"path", "[", "\"\<PathFunction\>\"", "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"jumps", "=", 
    RowBox[{"path", "[", "\"\<Times\>\"", "]"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{"f", "[", "t", "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"t", ",", "0", ",", "tsim"}], "}"}], ",", 
    RowBox[{"Exclusions", "\[Rule]", "jumps"}], ",", 
    RowBox[{"ExclusionsStyle", "\[Rule]", "Dashed"}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"\"\<t\>\"", ",", "\"\<N(t)\>\""}], "}"}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", "Black"}]}], "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[ScriptCapitalP]", "=", 
   RowBox[{"InhomogeneousPoissonProcess", "[", 
    RowBox[{
     RowBox[{"\[Lambda]", "[", "t", "]"}], ",", "t"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Mean", "[", 
  RowBox[{"\[ScriptCapitalP]", "[", "t", "]"}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"PDF", "[", 
    RowBox[{
     RowBox[{"PoissonDistribution", "[", 
      RowBox[{"Mean", "[", 
       RowBox[{"\[ScriptCapitalP]", "[", 
        RowBox[{
         SqrtBox["2"], "-", "1"}], "]"}], "]"}], "]"}], ",", "2"}], "]"}], 
   "*", 
   RowBox[{"CDF", "[", 
    RowBox[{
     RowBox[{"PoissonDistribution", "[", 
      RowBox[{
       RowBox[{"Mean", "[", 
        RowBox[{"\[ScriptCapitalP]", "[", "1", "]"}], "]"}], "-", 
       RowBox[{"Mean", "[", 
        RowBox[{"\[ScriptCapitalP]", "[", 
         RowBox[{
          SqrtBox["2"], "-", "1"}], "]"}], "]"}]}], "]"}], ",", "1"}], 
    "]"}]}], "//", "N"}]}], "Input",
 CellChangeTimes->{{3.792399880589603*^9, 3.79239990929224*^9}, {
   3.7924187781792603`*^9, 3.792418817393754*^9}, {3.7924188627234592`*^9, 
   3.792418904363263*^9}, {3.792418943050221*^9, 3.792418954848938*^9}, {
   3.792418994057446*^9, 3.792418994610649*^9}, {3.792419054929047*^9, 
   3.7924190864956408`*^9}, {3.792419131317869*^9, 3.792419181676249*^9}, {
   3.7924194165892878`*^9, 3.792419417421748*^9}, {3.82566394106176*^9, 
   3.825664037120438*^9}, {3.825665998546997*^9, 3.8256660048506308`*^9}, {
   3.825666035575835*^9, 3.825666080940946*^9}, {3.825666121259947*^9, 
   3.825666124185342*^9}, {3.8256662502028217`*^9, 3.825666281770584*^9}, {
   3.825666314140428*^9, 3.8256663957756844`*^9}, 3.825666517801463*^9, {
   3.825666578665436*^9, 3.825666624372158*^9}}],

Cell[BoxData[
 GraphicsBox[{{{{}, {}, 
     {GrayLevel[0], AbsoluteThickness[1.6], Opacity[1.], 
      LineBox[{{2.0408163265306121`*^-7, 0.}, {2.0408163265306121`*^-7, 0.}}],
       LineBox[{{0.00041436299956834155`, 0.}, {0.0018552058964517784`, 0.}, {
       0.0035063301685803582`, 0.}, {0.006808578712837517, 0.}, {
       0.013413075801351836`, 0.}, {0.026622069978380472`, 0.}, {
       0.053040058332437745`, 0.}, {0.054691182604566325`, 0.}, {
       0.056342306876694906`, 0.}, {0.05964455542095207, 0.}, {
       0.06624905250946639, 0.}, {0.07945804668649502, 0.}, {
       0.10481389926529329`, 0.}}], 
      LineBox[{{0.10734633406445769`, 1.}, {0.30497945502612234`, 1.}, {
       0.36750516656754845`, 1.}, {0.42704618422309365`, 1.}}], 
      LineBox[{{0.4334237352435018, 2.}, {0.641747236125647, 2.}, {
       0.8389537374769044, 2.}, {1.0526707890598728`, 2.}, {
       1.2624902284422326`, 2.}, {1.4582079744078238`, 2.}, {
       1.6704362706051261`, 2.}, {1.8685628733856596`, 2.}, {
       2.0627918639655847`, 2.}, {2.2735314047772204`, 2.}, {
       2.3353454854589053`, 2.}, {2.39417487225471, 2.}}], 
      LineBox[CompressedData["
1:eJxTTMoPSmViYGBQB2IQfSzN6EqICbMDAxhwOGxPZFXe/4gFzl/4ivWCVz8b
nP9X6pjyVzkOOD9wosCB/AuccL74eilmzlRuOP+q0Pq+ImleOL/97+yFD07y
IeTdl6zdpisA5z/evF8vlV0Qzl/12Yp58SUE/3VZkvTjmUJwvobs+XDuImE4
X2767OeHrEXgfAGTtbOTxETh/BqR7MksbxH8hSc67pluFYPzD81/7evfJw7n
z/FaJ/AjQgLO/2AhK75YXhLOd/C66Gj2HcH31H+6P+KAFJzvZ26XlDxLGs7n
Ef4e/jdVBs53MYqZOFNXFs6vKGF7L8smB+endgksMLyA4Psdji/VXS0P59+Y
r3tc1FEBzn/qHl0bsxTBL2hc7SoupIhQv2T19pnlCP7fl43esrcR/AUHqj6K
uCvB+fKNyo5XViP4MpffLWKXVIbzr7lNULNpQfC95t7mlzuE4O/7rPbz/38E
HwBFsopd
       "]]}, 
     {Dashing[{Small, Small}], 
      LineBox[{{2.0408163265306121`*^-7, 0.}, {0.00041436299956834155`, 0.}}],
       LineBox[{{0.10481389926529329`, 0.}, {0.10734633406445769`, 1.}}], 
      LineBox[{{0.42704618422309365`, 1.}, {0.4334237352435018, 2.}}], 
      LineBox[{{2.39417487225471, 2.}, {2.4005524232751183`, 
       3.}}]}}, {{}, {}, {}}}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{
    FormBox["\"t\"", TraditionalForm], 
    FormBox["\"N(t)\"", TraditionalForm]},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None},
  PlotRange->{{0, 10}, {0., 3.}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.792418916263425*^9, 3.792418955364818*^9}, {
   3.792419070463434*^9, 3.792419133408869*^9}, 3.792419185858501*^9, 
   3.7924194183806543`*^9, {3.825435395596754*^9, 3.825435443395748*^9}, {
   3.825663913453311*^9, 3.8256639174404716`*^9}, {3.825666060231928*^9, 
   3.825666081645631*^9}, 3.8256661251332912`*^9, 3.825666266157131*^9, 
   3.8256663249524384`*^9, {3.825666381488139*^9, 3.825666397139192*^9}, {
   3.8256665798101397`*^9, 3.825666634680894*^9}}],

Cell[BoxData[
 RowBox[{"ArcTan", "[", "t", "]"}]], "Output",
 CellChangeTimes->{{3.792418916263425*^9, 3.792418955364818*^9}, {
   3.792419070463434*^9, 3.792419133408869*^9}, 3.792419185858501*^9, 
   3.7924194183806543`*^9, {3.825435395596754*^9, 3.825435443395748*^9}, {
   3.825663913453311*^9, 3.8256639174404716`*^9}, {3.825666060231928*^9, 
   3.825666081645631*^9}, 3.8256661251332912`*^9, 3.825666266157131*^9, 
   3.8256663249524384`*^9, {3.825666381488139*^9, 3.825666397139192*^9}, {
   3.8256665798101397`*^9, 3.82566663519592*^9}}],

Cell[BoxData["0.048961304060546995`"], "Output",
 CellChangeTimes->{{3.792418916263425*^9, 3.792418955364818*^9}, {
   3.792419070463434*^9, 3.792419133408869*^9}, 3.792419185858501*^9, 
   3.7924194183806543`*^9, {3.825435395596754*^9, 3.825435443395748*^9}, {
   3.825663913453311*^9, 3.8256639174404716`*^9}, {3.825666060231928*^9, 
   3.825666081645631*^9}, 3.8256661251332912`*^9, 3.825666266157131*^9, 
   3.8256663249524384`*^9, {3.825666381488139*^9, 3.825666397139192*^9}, {
   3.8256665798101397`*^9, 3.825666635753625*^9}}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{663, 595},
WindowMargins->{{4, Automatic}, {Automatic, 4}},
PrintingCopies->1,
PrintingPageRange->{1, Automatic},
FrontEndVersion->"10.4 for Mac OS X x86 (32-bit, 64-bit Kernel) (April 11, \
2016)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 445, 6, 35, "Subsubsection"],
Cell[CellGroupData[{
Cell[1050, 32, 3137, 82, 289, "Input"],
Cell[4190, 116, 3342, 67, 249, "Output"],
Cell[7535, 185, 545, 8, 28, "Output"],
Cell[8083, 195, 534, 7, 28, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)