(* Content-type: application/mathematica *)

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

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[      7930,        205]
NotebookOptionsPosition[      6386,        152]
NotebookOutlinePosition[      7572,        192]
CellTagsIndexPosition[      7529,        189]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Ex3: Math. Modelling in Biosciences 1", "Title",
 CellChangeTimes->{{3.415589585204612*^9, 3.415589598615617*^9}, {
  3.415612138399819*^9, 3.415612142124673*^9}, {3.418375570847907*^9, 
  3.418375571000227*^9}, {3.42259016611127*^9, 3.422590166214294*^9}}],

Cell["Summer 2008 ", "Subtitle",
 CellChangeTimes->{{3.41561211409342*^9, 3.415612134422485*^9}}],

Cell["Due Friday 4 July (2.5 weeks)", "Subsubtitle",
 CellChangeTimes->{{3.415605888680184*^9, 3.41560590190364*^9}, {
  3.418375565010216*^9, 3.41837556832417*^9}, {3.418444153818715*^9, 
  3.418444167633437*^9}, {3.418449387119722*^9, 3.41844938826529*^9}, {
  3.422585087353535*^9, 3.422585093446563*^9}},
 FontSize->22],

Cell[TextData[{
 "In this exercise, we use MATCONT to reproduce the bifurcation diagrams \
shown in \n\n",
 StyleBox["Sniffers, buzzers, toggles and blinkers: dynamics of regulatory \
and signaling pathways in the cell.",
  FontSlant->"Italic",
  FontColor->RGBColor[0.6, 0.4, 0.2]],
 StyleBox["  ",
  FontSlant->"Italic"],
 "John J. Tyson, Katherine C. Chen and Bela Novak. ",
 StyleBox["Current Opinion in Cell Biology",
  FontWeight->"Bold",
  FontSlant->"Italic",
  FontColor->RGBColor[0.5, 0, 0.5]],
 StyleBox[" 15:221-231 (2003)",
  FontColor->RGBColor[0.5, 0, 0.5]],
 "\n\nas was also discussed in Lecture 4.\n\nIt would be useful to look at:\n\
\n\t",
 StyleBox["run_BifAnalysis_ML.m\n\tMorrisLecarModel.m",
  FontColor->RGBColor[0.6, 0.4, 0.2]],
 "\n\t\nAlso useful would be information on the usage of MATCONT:\n\t\n\t",
 StyleBox["MATCONT_manual.pdf",
  FontColor->RGBColor[0.6, 0.4, 0.2]]
}], "Text",
 CellChangeTimes->{{3.415605406565242*^9, 3.415605433891287*^9}, {
  3.415605464231621*^9, 3.415605549469606*^9}, {3.415605580453104*^9, 
  3.415605694991905*^9}, {3.415605756477949*^9, 3.415605842835668*^9}, {
  3.41560592575955*^9, 3.415605926779197*^9}, {3.415612161819669*^9, 
  3.415612164172899*^9}, {3.415612206830281*^9, 3.415612226928721*^9}, {
  3.418375652291134*^9, 3.41837568525645*^9}, {3.422585117688557*^9, 
  3.422585195239557*^9}, {3.42258523979635*^9, 3.422585248559625*^9}, {
  3.422591086433815*^9, 3.422591104868098*^9}, {3.422594532355084*^9, 
  3.422594561541268*^9}, {3.422595229232656*^9, 3.422595237729008*^9}}],

Cell[CellGroupData[{

Cell[TextData[{
 "Q1: reproduce the bistability diagram of Figure 1f ",
 StyleBox["(25 points)",
  FontColor->RGBColor[0.6, 0.4, 0.2]]
}], "Item",
 CellChangeTimes->{{3.415605951787572*^9, 3.415605952424048*^9}, {
   3.41560696626531*^9, 3.415607037040669*^9}, {3.415607092775741*^9, 
   3.415607100566137*^9}, {3.415607140563574*^9, 3.415607149990141*^9}, {
   3.415611936111022*^9, 3.41561194400363*^9}, {3.415612033505567*^9, 
   3.415612035785994*^9}, {3.418376177719356*^9, 3.418376196071539*^9}, {
   3.418376255255514*^9, 3.418376277850439*^9}, {3.418376435979911*^9, 
   3.418376446235365*^9}, {3.418376737700236*^9, 3.418376739481401*^9}, {
   3.418377261762019*^9, 3.418377263372706*^9}, {3.418388708011525*^9, 
   3.418388708476753*^9}, 3.42259072538941*^9, {3.422590897224965*^9, 
   3.422590912481908*^9}, {3.422591009091628*^9, 3.422591018071728*^9}, {
   3.422591160919988*^9, 3.422591161055477*^9}, {3.422591192356729*^9, 
   3.422591192447517*^9}},
 FontSize->18],

Cell["\<\
do the continuation of equilibrium to find the saddle-node points\
\>", "Subitem",
 CellChangeTimes->{{3.418376212405069*^9, 3.418376212551148*^9}, {
  3.418376251847973*^9, 3.418376296387643*^9}, {3.418376762144599*^9, 
  3.418376811727863*^9}, {3.42259104218181*^9, 3.422591064204932*^9}}],

Cell[TextData[{
 "Q2: reproduce the bifurcation diagram of Figure 2b ",
 StyleBox["(35 points)",
  FontColor->RGBColor[0.6, 0.4, 0.2]],
 " "
}], "Item",
 CellChangeTimes->{{3.415605951787572*^9, 3.415605952424048*^9}, {
   3.41560696626531*^9, 3.415607037040669*^9}, {3.415607067867762*^9, 
   3.41560707747891*^9}, {3.415608019137917*^9, 3.415608052609592*^9}, {
   3.415608130722612*^9, 3.415608180281544*^9}, {3.415610537357123*^9, 
   3.415610612273495*^9}, {3.415610772280362*^9, 3.415610777369745*^9}, {
   3.418376478958577*^9, 3.418376508286136*^9}, {3.418376654676232*^9, 
   3.418376692528988*^9}, {3.418376890547506*^9, 3.418376894943556*^9}, {
   3.418377044325435*^9, 3.418377045728649*^9}, {3.418377191814391*^9, 
   3.418377193674951*^9}, {3.418377347770495*^9, 3.418377348156679*^9}, {
   3.418377667579614*^9, 3.418377682795572*^9}, {3.422590748417265*^9, 
   3.422590806845214*^9}, {3.422591021931011*^9, 3.422591029055353*^9}, 
   3.422591110280217*^9, {3.422591198163694*^9, 3.422591198254127*^9}},
 FontSize->18],

Cell[CellGroupData[{

Cell["\<\
firstly, do the continuation of equilibrium to find the Hopf points\
\>", "Subitem",
 CellChangeTimes->{{3.418376212405069*^9, 3.418376212551148*^9}, {
  3.418376251847973*^9, 3.418376296387643*^9}, {3.418376762144599*^9, 
  3.418376811727863*^9}, {3.418383181966587*^9, 3.418383228331696*^9}, {
  3.418443664201302*^9, 3.41844377606462*^9}, {3.418443823381574*^9, 
  3.41844385526836*^9}, {3.418443909551971*^9, 3.418443914255221*^9}, {
  3.422590811387737*^9, 3.422590856143824*^9}, {3.422591049745414*^9, 
  3.422591049961051*^9}}],

Cell["\<\
then, get limit cycles off the Hopf point(s) (hint: use init_H_LC)\
\>", "Subitem",
 CellChangeTimes->{{3.422590825519265*^9, 3.422590875159374*^9}, {
  3.422591046579824*^9, 3.422591047437434*^9}, {3.42260263913922*^9, 
  3.422602639154668*^9}}]
}, Open  ]]
}, Open  ]],

Cell[TextData[{
 "\t- to superimpose limit cycle solution on the figure, take a look at:  ",
 StyleBox["SuperimposePlot.m",
  FontColor->RGBColor[0.6, 0.4, 0.2]]
}], "Text",
 CellChangeTimes->{{3.422602640447674*^9, 3.422602667278125*^9}, {
  3.422602799456587*^9, 3.422602840863679*^9}, {3.422603298920088*^9, 
  3.422603303347545*^9}}]
}, Open  ]]
},
WindowSize->{914, 931},
WindowMargins->{{Automatic, 131}, {Automatic, 0}},
PrintingPageRange->{Automatic, Automatic},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{611.25, 789.5625},
"PostScriptOutputFile"->"/home/jameslu/Teaching/MathModelBioSciences_Summer08/\
EX2/EX2_Question.nb.ps"},
ShowSelection->True,
FrontEndVersion->"6.0 for Linux x86 (32-bit) (April 20, 2007)",
StyleDefinitions->Notebook[{
   Cell[
    StyleData[
    StyleDefinitions -> 
     FrontEnd`FileName[{"Creative"}, "PastelColor.nb", CharacterEncoding -> 
       "iso8859-1"]]], 
   Cell[
    StyleData["Input"], FontSize -> 18], 
   Cell[
    StyleData["Output"], FontSize -> 18], 
   Cell[
    StyleData["Text"], FontSize -> 18], 
   Cell[
    StyleData["Item"], FontSize -> 22], 
   Cell[
    StyleData["Subitem"], FontSize -> 18, FontColor -> GrayLevel[0]]}, 
  Visible -> False, FrontEndVersion -> 
  "6.0 for Linux x86 (32-bit) (April 20, 2007)", StyleDefinitions -> 
  "PrivateStylesheetFormatting.nb"]
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[590, 23, 263, 3, 74, "Title"],
Cell[856, 28, 97, 1, 38, "Subtitle"],
Cell[956, 31, 323, 5, 39, "Subsubtitle"],
Cell[1282, 38, 1551, 32, 427, "Text"],
Cell[CellGroupData[{
Cell[2858, 74, 980, 17, 37, "Item"],
Cell[3841, 93, 301, 5, 26, "Subitem"],
Cell[4145, 100, 1033, 18, 37, "Item"],
Cell[CellGroupData[{
Cell[5203, 122, 544, 9, 26, "Subitem"],
Cell[5750, 133, 256, 5, 26, "Subitem"]
}, Open  ]]
}, Open  ]],
Cell[6033, 142, 337, 7, 63, "Text"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)