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

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

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[      9775,        252]
NotebookOptionsPosition[      9080,        229]
NotebookOutlinePosition[      9434,        245]
CellTagsIndexPosition[      9391,        242]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
Seasonal migration : u_t = -A*c*Cos[c, t] u_x, u[x, 0] = u0 (x) >= 0

Instructions: Evaluate notebook to start the module! \
\>", "Section",
 Evaluatable->False,
 CellChangeTimes->{{3.7603674368785467`*^9, 3.7603674698813343`*^9}, {
  3.76129844411535*^9, 3.761298451634179*^9}, {3.761298517443467*^9, 
  3.761298550224895*^9}},ExpressionUUID->"f774301f-fd8a-424c-adac-\
1746364c2bef"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.7612985296444273`*^9, 3.761298534688897*^9}},
 CellLabel->"In[31]:=",ExpressionUUID->"b661f96d-5af1-430c-8100-4b7e990bc0e9"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"DynamicModule", "[", 
  RowBox[{
   RowBox[{"{", "U", "}"}], ",", 
   RowBox[{
    RowBox[{"USeasonal", "=", 
     RowBox[{
      RowBox[{"U", "[", 
       RowBox[{"x", ",", "t"}], "]"}], "/.", 
      RowBox[{"First", "[", 
       RowBox[{"DSolve", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"D", "[", 
             RowBox[{
              RowBox[{"U", "[", 
               RowBox[{"x", ",", "t"}], "]"}], ",", "t"}], "]"}], "\[Equal]", 
            
            RowBox[{
             RowBox[{"-", "A"}], "*", "c", "*", 
             RowBox[{"Cos", "[", 
              RowBox[{"c", "*", "t"}], "]"}], "*", 
             RowBox[{"D", "[", 
              RowBox[{
               RowBox[{"U", "[", 
                RowBox[{"x", ",", "t"}], "]"}], ",", "x"}], "]"}]}]}], ",", 
           RowBox[{
            RowBox[{"U", "[", 
             RowBox[{"x", ",", "0"}], "]"}], "\[Equal]", 
            RowBox[{"PDF", "[", 
             RowBox[{
              RowBox[{"NormalDistribution", "[", 
               RowBox[{"50", ",", "6"}], "]"}], ",", "x"}], "]"}]}]}], "}"}], 
         ",", 
         RowBox[{"U", "[", 
          RowBox[{"x", ",", "t"}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"x", ",", "t"}], "}"}]}], "]"}], "]"}]}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Manipulate", "[", 
     RowBox[{
      RowBox[{"Plot", "[", 
       RowBox[{
        RowBox[{"USeasonal", "/.", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"A", "\[Rule]", "amplitude"}], ",", 
           RowBox[{"c", "\[Rule]", "speed"}], ",", 
           RowBox[{"t", "\[Rule]", "time"}]}], "}"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"x", ",", "0", ",", "100"}], "}"}], ",", 
        RowBox[{"AxesLabel", "\[Rule]", 
         RowBox[{"{", 
          RowBox[{"x", ",", "u"}], "}"}]}], ",", 
        RowBox[{"PlotLabel", "\[Rule]", 
         RowBox[{"Row", "[", 
          RowBox[{"{", 
           RowBox[{"\"\<Seasonal Migration\\n\>\"", ",", 
            RowBox[{"HoldForm", "@", 
             RowBox[{"D", "[", 
              RowBox[{"u", ",", "t"}], "]"}]}], ",", "\"\<=\>\"", ",", 
            RowBox[{"-", 
             RowBox[{"AcCos", "[", "ct", "]"}]}], ",", 
            RowBox[{"HoldForm", "@", 
             RowBox[{"D", "[", 
              RowBox[{"u", ",", "x"}], "]"}]}]}], "}"}], "]"}]}], ",", 
        RowBox[{"Ticks", "->", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"25", ",", "\"\< \>\""}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{"50", ",", " ", "\"\< \>\""}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{"75", ",", " ", "\"\< \>\""}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{"100", ",", " ", "\"\< \>\""}], "}"}]}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{".02", ",", " ", "\"\< \>\""}], "}"}], ",", " ", 
             RowBox[{"{", 
              RowBox[{".04", ",", "\"\< \>\""}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{".06", ",", " ", "\"\< \>\""}], "}"}]}], "}"}]}], 
          "}"}]}]}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"amplitude", ",", "10", ",", "\"\<A\>\""}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"10", ",", "30"}], "}"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"speed", ",", " ", "10", ",", " ", "\"\<speed c\>\""}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"10", ",", "20"}], "}"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"time", ",", "0", ",", "\"\<time t\>\""}], "}"}], ",", " ", 
        "0", ",", 
        RowBox[{"2", " ", "*", 
         RowBox[{"Pi", "/", "10"}]}]}], "}"}], ",", 
      RowBox[{"ControlType", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{"SetterBar", ",", "SetterBar", ",", "Trigger"}], "}"}]}], ",", 
      RowBox[{"SaveDefinitions", "\[Rule]", "True"}]}], "]"}]}]}], 
  "]"}]], "Input",
 CellOpen->False,
 CellChangeTimes->{{3.760367485683138*^9, 3.760367845646244*^9}},
 CellLabel->"In[32]:=",ExpressionUUID->"6624747d-c676-4a2d-b674-edb90d6b189e"],

Cell[BoxData[
 DynamicModuleBox[{$CellContext`U$$}, 
  TagBox[
   StyleBox[
    DynamicModuleBox[{$CellContext`amplitude$$ = 10, $CellContext`speed$$ = 
     10, $CellContext`time$$ = 0, Typeset`show$$ = True, 
     Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", 
     Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = 
     "\"untitled\"", Typeset`specs$$ = {{{
        Hold[$CellContext`amplitude$$], 10, "A"}, {10, 30}}, {{
        Hold[$CellContext`speed$$], 10, "speed c"}, {10, 20}}, {{
        Hold[$CellContext`time$$], 0, "time t"}, 0, Rational[1, 5] Pi}}, 
     Typeset`size$$ = {360., {143., 148.}}, Typeset`update$$ = 0, 
     Typeset`initDone$$, Typeset`skipInitDone$$ = 
     False, $CellContext`amplitude$33198$$ = 
     False, $CellContext`speed$33199$$ = False, $CellContext`time$33200$$ = 
     0}, 
     DynamicBox[Manipulate`ManipulateBoxes[
      1, StandardForm, 
       "Variables" :> {$CellContext`amplitude$$ = 10, $CellContext`speed$$ = 
         10, $CellContext`time$$ = 0}, "ControllerVariables" :> {
         Hold[$CellContext`amplitude$$, $CellContext`amplitude$33198$$, 
          False], 
         Hold[$CellContext`speed$$, $CellContext`speed$33199$$, False], 
         Hold[$CellContext`time$$, $CellContext`time$33200$$, 0]}, 
       "OtherVariables" :> {
        Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
         Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
         Typeset`specs$$, Typeset`size$$, Typeset`update$$, 
         Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> Plot[
         ReplaceAll[$CellContext`USeasonal, {$CellContext`A -> \
$CellContext`amplitude$$, $CellContext`c -> $CellContext`speed$$, \
$CellContext`t -> $CellContext`time$$}], {$CellContext`x, 0, 100}, 
         AxesLabel -> {$CellContext`x, $CellContext`u}, PlotLabel -> 
         Row[{"Seasonal Migration\n", 
            HoldForm[
             D[$CellContext`u, $CellContext`t]], 
            "=", -$CellContext`AcCos[$CellContext`ct], 
            HoldForm[
             D[$CellContext`u, $CellContext`x]]}], 
         Ticks -> {{{25, " "}, {50, " "}, {75, " "}, {100, " "}}, {{
            0.02, " "}, {0.04, " "}, {0.06, " "}}}], 
       "Specifications" :> {{{$CellContext`amplitude$$, 10, "A"}, {10, 
          30}}, {{$CellContext`speed$$, 10, "speed c"}, {10, 
          20}}, {{$CellContext`time$$, 0, "time t"}, 0, Rational[1, 5] Pi}}, 
       "Options" :> {ControlType -> {SetterBar, SetterBar, Trigger}}, 
       "DefaultOptions" :> {}],
      ImageSizeCache->{405., {214., 220.}},
      SingleEvaluation->True],
     Deinitialization:>None,
     DynamicModuleValues:>{},
     Initialization:>({$CellContext`USeasonal = 
        1/(6 E^((-50 + $CellContext`x - $CellContext`A 
             Sin[$CellContext`c $CellContext`t])^2/72) 
         Sqrt[2 Pi]), $CellContext`u[$CellContext`x, 
          0] = $CellContext`u0 $CellContext`x >= 0}; 
      Typeset`initDone$$ = True),
     SynchronousInitialization->True,
     UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$},
     UnsavedVariables:>{Typeset`initDone$$},
     UntrackedVariables:>{Typeset`size$$}], "Manipulate",
    Deployed->True,
    StripOnInput->False],
   Manipulate`InterpretManipulate[1]],
  DynamicModuleValues:>{}]], "Output",
 CellChangeTimes->{3.760367846510125*^9, 3.760433170683612*^9, 
  3.761298476461625*^9, 3.761298522830193*^9, 3.761298560812953*^9},
 CellLabel->"Out[32]=",ExpressionUUID->"ae096931-d2c4-4126-9e1f-411898cdb3f6"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{808, 911},
WindowMargins->{{Automatic, 390}, {Automatic, 0}},
FrontEndVersion->"11.3 for Mac OS X x86 (32-bit, 64-bit Kernel) (March 5, \
2018)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 395, 9, 143, "Section",ExpressionUUID->"f774301f-fd8a-424c-adac-1746364c2bef",
 Evaluatable->False],
Cell[978, 33, 173, 2, 30, "Input",ExpressionUUID->"b661f96d-5af1-430c-8100-4b7e990bc0e9"],
Cell[CellGroupData[{
Cell[1176, 39, 4355, 117, 19, "Input",ExpressionUUID->"6624747d-c676-4a2d-b674-edb90d6b189e",
 CellOpen->False],
Cell[5534, 158, 3518, 67, 487, "Output",ExpressionUUID->"ae096931-d2c4-4126-9e1f-411898cdb3f6"]
}, Open  ]]
}, Open  ]]
}
]
*)