(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 10.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 26343, 1048] NotebookOptionsPosition[ 22344, 905] NotebookOutlinePosition[ 22905, 927] CellTagsIndexPosition[ 22862, 924] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ StyleBox[ RowBox[{"Lotka", "-", RowBox[{"Volterra", " ", "mudel", " ", RowBox[{"(", RowBox[{"kahe", " ", "liigi", " ", "kooslus"}], ")"}]}]}], "Subsubsection", FontSize->24]], "Input", CellChangeTimes->{{3.7818315058479276`*^9, 3.781831505849928*^9}}], Cell[CellGroupData[{ Cell[TextData[{ StyleBox["Probleemi p\[UDoubleDot]stitus: ", FontWeight->"Bold"], "Olgu ", StyleBox["x", FontSlant->"Italic"], "(", StyleBox["t", FontSlant->"Italic"], ") j\[ADoubleDot]neste ja ", StyleBox["y", FontSlant->"Italic"], " (", StyleBox["t", FontSlant->"Italic"], ") huntide arv ajamomendil ", StyleBox["t.", FontSlant->"Italic"], " J\[ADoubleDot]neste toidulaual pole praktilisi piiranguid, sest hundid ei \ lase neil liialt paljuneda. Huntide toidulaud s\[OTilde]ltub \ j\[ADoubleDot]neste arvukusest. Teeme lihtsustava eelduse, et \ j\[ADoubleDot]nestel pole teisi looduslikke vaenlasi peale huntide, ning j\ \[ADoubleDot]nesed on huntide ainsaks toiduks. Kui poleks hunte, alluks j\ \[ADoubleDot]neste arvukus Malthuse seadusele: " }], "Text", CellChangeTimes->{{3.78183844438999*^9, 3.781838448635233*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"x", "'"}], RowBox[{"(", "t", ")"}]}], "=", RowBox[{"a", " ", "x", RowBox[{"(", "t", ")"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"kus", " ", "a"}], ">", "0"}], ";"}]}]], "DisplayFormula"], Cell[CellGroupData[{ Cell["Kui poleks j\[ADoubleDot]neseid,alluks huntide arvukus seadusele", \ "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "'"}], RowBox[{"(", "t", ")"}]}], "=", RowBox[{ RowBox[{"-", "b"}], " ", "y", RowBox[{"(", "t", ")"}]}]}], ",", RowBox[{ RowBox[{"kus", " ", "b"}], ">", "0."}]}]], "DisplayFormula"], Cell[CellGroupData[{ Cell["\<\ J\[ADoubleDot]neste ja huntide kooseksisteerimist v\[OTilde]ib kirjeldada \ lihtsa diferentsiaalv\[OTilde]rrandite s\[UDoubleDot]steemiga\ \>", "Text"], Cell[BoxData[ RowBox[{ RowBox[{ StyleBox[ FractionBox["dx", "dt"], FontFamily->"Times New Roman", FontSlant->"Italic"], "=", RowBox[{ StyleBox["ax", FontSlant->"Italic"], StyleBox[" ", FontSlant->"Italic"], StyleBox["-", FontSlant->"Italic"], StyleBox[" ", FontSlant->"Italic"], RowBox[{ FormBox[ SubscriptBox["a", "1"], TraditionalForm], StyleBox["xy", FontSlant->"Italic"]}]}]}], StyleBox[",", FontSlant->"Italic"]}]], "DisplayFormula"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ StyleBox[ FractionBox["dy", "dt"], FontFamily->"Times New Roman", FontSlant->"Italic"], "=", RowBox[{ RowBox[{"-", StyleBox["by", FontSlant->"Italic"]}], StyleBox[" ", FontSlant->"Italic"], StyleBox["+", FontSlant->"Italic"], StyleBox[" ", FontSlant->"Italic"], RowBox[{ FormBox[ SubscriptBox["b", "1"], TraditionalForm], StyleBox["xy", FontSlant->"Italic"]}]}]}], StyleBox[",", FontSlant->"Italic"]}]], "DisplayFormula", CellChangeTimes->{3.7818395199485083`*^9}], Cell[CellGroupData[{ Cell[TextData[{ "kus ", StyleBox["a, b, ", FontSlant->"Italic"], Cell[BoxData[ FormBox[ RowBox[{ SubscriptBox["a", RowBox[{"1", " "}]], "ja", " ", StyleBox[ SubscriptBox["b", StyleBox["1", FontSlant->"Plain"]], FontSlant->"Italic"]}], TraditionalForm]]], " on positiivsed konstandid. Siin liige ", StyleBox["x", FontSlant->"Italic"], "(", StyleBox["t", FontSlant->"Italic"], ")", StyleBox["y", FontSlant->"Italic"], "(", StyleBox["t", FontSlant->"Italic"], ") on proportsionaalne j\[ADoubleDot]nese ja hundi kohtumise t\[OTilde]en\ \[ADoubleDot]osusele ning sellised kohtumised m\[OTilde]juvad negatiivselt j\ \[ADoubleDot]neste arvukusele ning positiivselt huntide arvukusele." }], "Text"], Cell[CellGroupData[{ Cell[TextData[{ "Siin ", StyleBox["x", FontSlant->"Italic"], "(", StyleBox["t", FontSlant->"Italic"], ") ja ", StyleBox["y", FontSlant->"Italic"], "(", StyleBox["t", FontSlant->"Italic"], ") on liikide arvukused, seega probleemi lahendamisel vaatleme neid v\ \[OTilde]rrandeid mitte kogu ", StyleBox["xy ", FontSlant->"Italic"], "tasandil vaid selle esimeses veerandis, s.t.", StyleBox[" x ", FontSlant->"Italic"], "\[GreaterEqual]", StyleBox[" ", FontSlant->"Italic"], "0 ja ", StyleBox["y ", FontSlant->"Italic"], "\[GreaterEqual]", StyleBox[" ", FontSlant->"Italic"], "0. " }], "Text"], Cell[CellGroupData[{ Cell["\<\ Jagades teise v\[OTilde]rrandi esimesega, saame \[UDoubleDot]he hariliku \ diferentsiaalv\[OTilde]rrandi\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ StyleBox[ FractionBox["dy", "dx"], FontSlant->"Italic"], StyleBox["=", FontSlant->"Italic"], FractionBox[ RowBox[{ RowBox[{"-", "by"}], "+", RowBox[{ SubscriptBox["b", "1"], "xy"}]}], RowBox[{"ax", "-", RowBox[{ SubscriptBox["a", "1"], "xy"}]}]]}]], "DisplayFormula"], Cell[CellGroupData[{ Cell["ehk ", "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"y", RowBox[{"(", RowBox[{ RowBox[{"-", StyleBox["b", FontSlant->"Italic"]}], StyleBox[" ", FontSlant->"Italic"], StyleBox["+", FontSlant->"Italic"], StyleBox[" ", FontSlant->"Italic"], RowBox[{ FormBox[ SubscriptBox["b", "1"], TraditionalForm], StyleBox["x", FontSlant->"Italic"]}]}], StyleBox[")", FontSlant->"Italic"]}], StyleBox["dx", FontSlant->"Italic"]}], StyleBox["=", FontSlant->"Italic"], RowBox[{"x", RowBox[{"(", RowBox[{"a", StyleBox[" ", FontSlant->"Italic"], StyleBox["-", FontSlant->"Italic"], StyleBox[" ", FontSlant->"Italic"], RowBox[{ FormBox[ SubscriptBox["a", "1"], TraditionalForm], StyleBox["y", FontSlant->"Italic"]}]}], StyleBox[")", FontSlant->"Italic"]}], StyleBox[ RowBox[{"dy", "."}], FontSlant->"Italic"]}]}]], "DisplayFormula"], Cell[CellGroupData[{ Cell["See on eralduvate muutujatega v\[OTilde]rrand. Eraldame muutujad", \ "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"xy", "[", RowBox[{ RowBox[{ FractionBox[ RowBox[{ RowBox[{"-", "b"}], "+", RowBox[{ SubscriptBox["b", "1"], "x"}]}], "x"], "dx"}], "-", RowBox[{ FractionBox[ RowBox[{"a", "-", RowBox[{ SubscriptBox["a", "1"], "y"}]}], "y"], "dy"}]}], "]"}], "=", "0."}]], "DisplayFormula"], Cell[CellGroupData[{ Cell[TextData[{ "Lahendid ", StyleBox["x", FontSlant->"Italic"], " = 0 ja ", StyleBox["y", FontSlant->"Italic"], " = 0 j\[ADoubleDot]tame k\[OTilde]rvale kui \[UDoubleDot]lesande sisust l\ \[ADoubleDot]htudes ebahuvitavad", ". Arvutame" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Integral]", RowBox[{ FractionBox[ RowBox[{ RowBox[{"-", "b"}], "+", RowBox[{ SubscriptBox["b", "1"], "x"}]}], "x"], "dx"}]}], "+", RowBox[{"\[Integral]", RowBox[{ FractionBox[ RowBox[{ RowBox[{"-", "a"}], "+", RowBox[{ SubscriptBox["a", "1"], "y"}]}], "y"], "dy"}]}]}], " ", "=", " ", StyleBox["C", FontSlant->"Italic"]}], StyleBox[",", FontSlant->"Italic"]}]], "DisplayFormula"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", StyleBox["b", FontSlant->"Italic"]}], StyleBox[" ", FontSlant->"Italic"], StyleBox["ln", FontSlant->"Plain"], StyleBox[" ", FontSlant->"Italic"], StyleBox["x", FontSlant->"Italic"]}], StyleBox[" ", FontSlant->"Italic"], StyleBox["+", FontSlant->"Plain"], StyleBox[" ", FontSlant->"Plain"], RowBox[{ SubscriptBox["b", "1"], "x"}], "-", RowBox[{ StyleBox["a", FontSlant->"Italic"], StyleBox[" ", FontSlant->"Italic"], StyleBox[ RowBox[{"ln", StyleBox["y", FontSlant->"Italic"]}]]}], StyleBox[" ", FontSlant->"Italic"], StyleBox["+", FontSlant->"Italic"], StyleBox[" ", FontSlant->"Italic"], RowBox[{ SubscriptBox["a", "1"], "y"}]}], " ", "=", " ", StyleBox["C", FontSlant->"Italic"]}]], "DisplayFormula", CellChangeTimes->{3.8468164542957644`*^9, 3.8468165088958607`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ StyleBox["ln", FontSlant->"Plain"], StyleBox[" ", FontSlant->"Italic"], RowBox[{ StyleBox["(", FontSlant->"Italic"], RowBox[{ SuperscriptBox["x", RowBox[{"-", "b"}]], SuperscriptBox[ SuperscriptBox["e", SubscriptBox["b", "1"]], "x"], SuperscriptBox["y", RowBox[{"-", "a"}]], SuperscriptBox[ SuperscriptBox["e", SubscriptBox["a", "1"]], "y"]}], ")"}]}], "=", RowBox[{"ln", " ", SubscriptBox["C", "1"]}]}], ",", " ", RowBox[{ RowBox[{"kus", " ", "ln", " ", SubscriptBox["C", "1"]}], " ", "=", " ", StyleBox[ RowBox[{"C", "."}], FontSlant->"Italic"]}]}]], "DisplayFormula", CellChangeTimes->{3.846816433344928*^9, 3.8468165088958607`*^9}], Cell[CellGroupData[{ Cell["Siit ", "Text", CellChangeTimes->{{3.7818315217888393`*^9, 3.781831525959078*^9}, { 3.846816407948083*^9, 3.846816420178505*^9}, 3.8468165088958607`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ FractionBox[ SuperscriptBox["e", RowBox[{ SubscriptBox["b", "1"], "x"}]], SuperscriptBox["x", "b"]], ".", FractionBox[ SuperscriptBox["e", RowBox[{ SubscriptBox["a", "1"], "y"}]], SuperscriptBox["y", "a"]]}], "=", RowBox[{ SubscriptBox["C", "1"], "."}]}]], "DisplayFormula", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{{3.846816407948083*^9, 3.846816420178505*^9}, { 3.8468165088958607`*^9, 3.846816532623502*^9}}], Cell[CellGroupData[{ Cell["\<\ Viimane seos esitabki lahendeid ilmutamata kujul. Selle l\[ADoubleDot]hemal \ uurimisel ilmneb, et lahendid kujutavad endast kinnisi k\[OTilde]veraid \ \[UDoubleDot]mber punkti, mille koordinaatideks on\ \>", "Text", CellGroupingRules->{GroupTogetherGrouping, 10001.}, CellChangeTimes->{ 3.846816380585635*^9, 3.8468164873054223`*^9, {3.846816532623502*^9, 3.8468165501735325`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"x", "=", FractionBox["b", SubscriptBox["b", "1"]]}], ",", " ", RowBox[{ StyleBox["y", FontSlant->"Italic"], StyleBox[" ", FontSlant->"Plain"], StyleBox["=", FontSlant->"Plain"], StyleBox[" ", FontSlant->"Plain"], RowBox[{ FractionBox["a", SubscriptBox["a", "1"]], "."}]}]}]], "DisplayFormula", CellGroupingRules->{GroupTogetherGrouping, 10001.}, CellChangeTimes->{ 3.846816380585635*^9, 3.8468164873054223`*^9, {3.846816532623502*^9, 3.8468165501735325`*^9}}] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Konkreetne mudel", "Subsection"], Cell[CellGroupData[{ Cell[TextData[{ "Oletame, et ilma huntideta suureneks j\[ADoubleDot]neste arvukus \ \[UDoubleDot]hes ajavahemikus neljakordseks. J\[ADoubleDot]neste niisugust \ kasvu saab modelleerida diferentsiaalv\[OTilde]rrandi ", Cell[BoxData[ RowBox[{ RowBox[{ StyleBox[ FractionBox["dx", "dt"], "InlineFormula", FontFamily->"Times New Roman", FontSize->16, FontSlant->"Italic"], StyleBox["=", "InlineFormula"], RowBox[{ StyleBox["1.4", FontSlant->"Plain"], StyleBox["x", FontFamily->"Times New Roman", FontSlant->"Italic"]}]}], StyleBox[" ", FontSlant->"Italic"], StyleBox[",", FontSlant->"Italic"]}]], "DisplayFormula"], " kui ", StyleBox["y", FontSlant->"Italic"], " = 0 abil. " }], "Text", CellGroupingRules->{GroupTogetherGrouping, 10001.}, CellChangeTimes->{3.84681635506399*^9}], Cell[CellGroupData[{ Cell[TextData[{ "Selle v\[OTilde]rrandi t\[ADoubleDot]pne lahend on ", StyleBox["x", FontSlant->"Italic"], "(", StyleBox["t", FontSlant->"Italic"], ") = ", StyleBox["x", FontSlant->"Italic"], "(0) ", Cell[BoxData[ FormBox[ RowBox[{ StyleBox[ SuperscriptBox["e", RowBox[{"1", ",", RowBox[{"4", "t"}]}]], FontSize->16], ",", " "}], TraditionalForm]]], "mis t\[ADoubleDot]hendab eksponentsiaalset kasvu teguriga ", Cell[BoxData[ FormBox[ StyleBox[ SuperscriptBox["e", RowBox[{"1", ",", "4"}]], FontSize->16], TraditionalForm]]], " mis on umbes 4,0552. Seega ", StyleBox["a", FontSlant->"Italic"], " = 1,4." }], "Text", CellGroupingRules->{GroupTogetherGrouping, 10001.}, CellChangeTimes->{3.84681635506399*^9}], Cell[CellGroupData[{ Cell[TextData[{ "Oletame,et ilma j\[ADoubleDot]nesteta v\[ADoubleDot]heneks huntide arvukus \ \[UDoubleDot]hes ajavahemikus umbes 10 % v\[OTilde]rra. \nHuntide arvu \ niisugust kahanemist saab modelleerida diferentsiaalv\[OTilde]rrandi \n ", Cell[BoxData[ RowBox[{ RowBox[{ StyleBox[ FractionBox["dy", "dt"], "InlineFormula", FontFamily->"Times New Roman", FontSize->16, FontSlant->"Italic"], StyleBox["=", "InlineFormula"], RowBox[{ RowBox[{"-", StyleBox["0.1", FontSlant->"Plain"]}], StyleBox["y", FontFamily->"Times New Roman", FontSlant->"Italic"]}]}], StyleBox[",", FontSlant->"Italic"]}]], "DisplayFormula"], " kui ", StyleBox["x", FontSlant->"Italic"], " = 0 abil. " }], "Text", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{3.846816341975567*^9}], Cell[CellGroupData[{ Cell[BoxData[{ StyleBox[ RowBox[{ RowBox[{ "Selle", " ", "v\[OTilde]rrandi", " ", "t\[ADoubleDot]pne", " ", "lahend", " ", "on", " ", StyleBox["y", FontSlant->"Italic"], RowBox[{"(", StyleBox["t", FontSlant->"Italic"], ")"}]}], "=", RowBox[{ StyleBox["y", FontSlant->"Italic"], RowBox[{"(", "0", ")"}], FormBox[ RowBox[{ StyleBox[ SuperscriptBox["e", RowBox[{ RowBox[{"-", "0"}], ",", RowBox[{"1", "t"}]}]], FontFamily->"Times New Roman", FontSize->16], ",", " "}], TraditionalForm], "mis", " ", "t\[ADoubleDot]hendab"}]}], FontFamily->"Times New Roman"], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{ RowBox[{"eksponentsiaalset", " ", "kahanemist", " ", "teguriga", FormBox[ StyleBox[ SuperscriptBox["e", RowBox[{ RowBox[{"-", "0"}], ",", "1", " "}]], FontSize->16], TraditionalForm], "mis", " ", "on", " ", "umbes", " ", RowBox[{"0.904837", ".", " ", "Seega"}], " ", StyleBox["b", FontSlant->"Italic"]}], " ", "=", "0"}], ",", "1."}], FontFamily->"Times New Roman"]}], "Text", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{3.8468162953002853`*^9}], Cell[CellGroupData[{ Cell[TextData[{ "Olgu ", Cell[BoxData[ FormBox[ SubscriptBox["a", "1"], TraditionalForm]]], " = 0.02 katsetest v\[OTilde]etud konstant, mis n\[ADoubleDot]itab, kui suur \ osa j\[ADoubleDot]nestest tapetakse \[UDoubleDot]he hundi poolt \ \[UDoubleDot]he ajaperioodi jooksul." }], "Text", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{3.8468162808234596`*^9}], Cell[CellGroupData[{ Cell[TextData[{ "Olgu ", Cell[BoxData[ FormBox[ SubscriptBox["b", "1"], TraditionalForm]]], " = 0.001 katsetest v\[OTilde]etud konstant, mis n\[ADoubleDot]itab, kui \ palju s\[UDoubleDot]nnib hunte juurde \[UDoubleDot]he ajaperioodi jooksul iga \ \[ADoubleDot]ras\[ODoubleDot]\[ODoubleDot]dud j\[ADoubleDot]nese kohta." }], "Text", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{3.8468162808234596`*^9}], Cell[CellGroupData[{ Cell["\<\ Seega j\[ADoubleDot]neste ja huntide kooseksisteerimist v\[OTilde]ib \ kirjeldagu diferentsiaalv\[OTilde]rrandite s\[UDoubleDot]steem\ \>", "Text", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{3.8468162294213696`*^9}], Cell[BoxData[ RowBox[{ RowBox[{ StyleBox[ FractionBox["dx", "dt"], FontFamily->"Times New Roman", FontSlant->"Italic"], "=", StyleBox["1", FontFamily->"Times New Roman"]}], StyleBox[",", FontFamily->"Times New Roman"], StyleBox[ RowBox[{ RowBox[{"4", StyleBox["x", FontSlant->"Italic"]}], StyleBox[" ", FontSlant->"Italic"], StyleBox["-", FontSlant->"Italic"], StyleBox[" ", FontSlant->"Italic"], RowBox[{ FormBox[ RowBox[{"0", ",", "02"}], TraditionalForm], StyleBox["xy", FontSlant->"Italic"]}]}], FontFamily->"Times New Roman"], StyleBox[",", FontFamily->"Times New Roman", FontSlant->"Italic"]}]], "DisplayFormula", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{3.8468162294213696`*^9}], Cell[BoxData[ RowBox[{ RowBox[{ StyleBox[ FractionBox["dy", "dt"], FontFamily->"Times New Roman", FontSlant->"Italic"], "=", StyleBox[ RowBox[{"-", "0"}], FontFamily->"Times New Roman"]}], StyleBox[",", FontFamily->"Times New Roman"], StyleBox[ RowBox[{ RowBox[{"1", StyleBox["y", FontSlant->"Italic"]}], StyleBox[" ", FontSlant->"Italic"], "+", StyleBox[" ", FontSlant->"Italic"], "0"}], FontFamily->"Times New Roman"], StyleBox[",", FontFamily->"Times New Roman"], StyleBox[ RowBox[{"001", StyleBox["xy", FontSlant->"Italic"]}], FontFamily->"Times New Roman"], StyleBox[",", FontFamily->"Times New Roman", FontSlant->"Italic"]}]], "DisplayFormula", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{3.7818395120780582`*^9, 3.8468162294213696`*^9, 3.846859750300012*^9}] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Konkreetse mudeli realisatsioon programmina", "Subsection"], Cell[CellGroupData[{ Cell["\<\ a = 1.4; a1 = 0.02; b = 0.1; b1 = 0.001; sol = NDSolve[ {x'[t] == a x[t] - a1 x[t] y[t], y'[t] == - b y[t] + b1 x[t] y[t], x[0] == 100, y[0] == 60 }, {y, x}, {t, 0, 80} ]; Plot[Evaluate[{x[t], y[t]} /. sol], {t, 0, 80}, PlotStyle -> \ {RGBColor[1,0,0], RGBColor[0,0,1]}, PlotRange->All] \ \>", "Input", CellChangeTimes->{{3.781831557500882*^9, 3.7818315595569997`*^9}, 3.8468161581292443`*^9, 3.8468162122769394`*^9, 3.8468598213260746`*^9, { 3.8468599663723707`*^9, 3.8468599671214137`*^9}, {3.8468601263545213`*^9, 3.8468601444185543`*^9}, {3.846860196374526*^9, 3.846860201962846*^9}, { 3.8468602578760443`*^9, 3.8468602593571286`*^9}, {3.8468603887440352`*^9, 3.846860400537656*^9}}], Cell[CellGroupData[{ Cell["\<\ gr1 = ParametricPlot[Evaluate[{x[t], y[t]} /. sol], {t, 0, 20}, Epilog -> \ {PointSize[0.02],Point[{b/b1, a/a1}], \tDashing[{0.03,0.03}], Line[{{b/b1, 0},{b/b1, 120}}], Line[{{0, a/a1},{420, \ a/a1}}] \t}]\ \>", "Input", CellChangeTimes->{{3.7818315405849147`*^9, 3.781831552310585*^9}, 3.8468161361800056`*^9, 3.8468162122769394`*^9}], Cell[BoxData[ RowBox[{"TableForm", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", RowBox[{ RowBox[{"(", RowBox[{"x", "/.", RowBox[{"sol", "[", RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "[", "t", "]"}], ",", RowBox[{ RowBox[{"(", RowBox[{"y", "/.", RowBox[{"sol", "[", RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "80"}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"TableHeadings", "\[Rule]", RowBox[{"{", RowBox[{"None", ",", RowBox[{"{", RowBox[{ "\"\\"", ",", " ", "\"\\"", ",", " ", "\"\\""}], "}"}]}], "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.781832338470551*^9, 3.7818323399956384`*^9}, 3.8468161361800056`*^9, 3.8468162122769394`*^9, 3.8468600530293274`*^9}] }, Closed]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Tasakaalupunkt konkreetse mudeli jaoks", "Subsection"], Cell[CellGroupData[{ Cell["\<\ Tasakaalupunkti antud algparameetrite puhul saame varemvaadeldud \ v\[OTilde]rdustest:\ \>", "Text", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{ 3.846815803618622*^9, {3.8468596135661917`*^9, 3.846859621396639*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"x", "=", RowBox[{ FractionBox["b", SubscriptBox["b", "1"]], "=", RowBox[{ FractionBox["0.1", "0.001"], "=", "100"}]}]}], ",", " ", RowBox[{ StyleBox["y", FontSlant->"Italic"], StyleBox[" ", FontSlant->"Plain"], StyleBox["=", FontSlant->"Plain"], StyleBox[" ", FontSlant->"Plain"], RowBox[{ FractionBox["a", SubscriptBox["a", "1"]], "=", RowBox[{ FractionBox["1.4", "0.02"], "=", "70."}]}]}]}]], "Input", CellChangeTimes->{{3.8468597317399507`*^9, 3.8468597881821785`*^9}}] }, Closed]] }, Closed]] }, CellGrouping->Manual, WindowSize->{1358, 697}, WindowMargins->{{-4, Automatic}, {Automatic, 0}}, PrintingCopies->1, PrintingPageRange->{32000, 32000}, PrintingOptions->{"Magnification"->1., "PaperOrientation"->"Portrait", "PaperSize"->{595.3199999999999, 841.92}}, Magnification:>2.2 Inherited, FrontEndVersion->"10.0 for Microsoft Windows (32-bit) (July 1, 2014)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[580, 22, 289, 8, 107, "Input"], Cell[CellGroupData[{ Cell[894, 34, 847, 25, 233, "Text"], Cell[CellGroupData[{ Cell[1766, 63, 269, 10, 50, "DisplayFormula"], Cell[CellGroupData[{ Cell[2060, 77, 82, 1, 65, "Text"], Cell[CellGroupData[{ Cell[2167, 82, 269, 10, 50, "DisplayFormula"], Cell[CellGroupData[{ Cell[2461, 96, 161, 3, 107, "Text"], Cell[2625, 101, 539, 23, 88, "DisplayFormula"], Cell[CellGroupData[{ Cell[3189, 128, 605, 25, 88, "DisplayFormula"], Cell[CellGroupData[{ Cell[3819, 157, 758, 29, 150, "Text"], Cell[CellGroupData[{ Cell[4602, 190, 625, 30, 151, "Text"], Cell[CellGroupData[{ Cell[5252, 224, 128, 3, 65, "Text"], Cell[CellGroupData[{ Cell[5405, 231, 339, 14, 81, "DisplayFormula"], Cell[CellGroupData[{ Cell[5769, 249, 20, 0, 65, "Text"], Cell[CellGroupData[{ Cell[5814, 253, 1031, 45, 50, "DisplayFormula"], Cell[CellGroupData[{ Cell[6870, 302, 82, 1, 65, "Text"], Cell[CellGroupData[{ Cell[6977, 307, 384, 15, 81, "DisplayFormula"], Cell[CellGroupData[{ Cell[7386, 326, 258, 10, 108, "Text"], Cell[CellGroupData[{ Cell[7669, 340, 550, 21, 87, "DisplayFormula"], Cell[CellGroupData[{ Cell[8244, 365, 1003, 42, 50, "DisplayFormula"], Cell[CellGroupData[{ Cell[9272, 411, 832, 30, 58, "DisplayFormula"], Cell[CellGroupData[{ Cell[10129, 445, 165, 2, 65, "Text"], Cell[CellGroupData[{ Cell[10319, 451, 534, 17, 87, "DisplayFormula", CellGroupingRules->{GroupTogetherGrouping, 10000.}], Cell[CellGroupData[{ Cell[10878, 472, 400, 8, 149, "Text", CellGroupingRules->{GroupTogetherGrouping, 10001.}], Cell[11281, 482, 561, 20, 80, "DisplayFormula", CellGroupingRules->{GroupTogetherGrouping, 10001.}] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[12131, 528, 38, 0, 82, "Subsection"], Cell[CellGroupData[{ Cell[12194, 532, 869, 29, 163, "Text", CellGroupingRules->{GroupTogetherGrouping, 10001.}], Cell[CellGroupData[{ Cell[13088, 565, 786, 32, 121, "Text", CellGroupingRules->{GroupTogetherGrouping, 10001.}], Cell[CellGroupData[{ Cell[13899, 601, 873, 28, 207, "Text", CellGroupingRules->{GroupTogetherGrouping, 10000.}], Cell[CellGroupData[{ Cell[14797, 633, 1294, 42, 122, "Text", CellGroupingRules->{GroupTogetherGrouping, 10000.}], Cell[CellGroupData[{ Cell[16116, 679, 392, 10, 107, "Text", CellGroupingRules->{GroupTogetherGrouping, 10000.}], Cell[CellGroupData[{ Cell[16533, 693, 437, 10, 107, "Text", CellGroupingRules->{GroupTogetherGrouping, 10000.}], Cell[CellGroupData[{ Cell[16995, 707, 254, 5, 107, "Text", CellGroupingRules->{GroupTogetherGrouping, 10000.}], Cell[17252, 714, 846, 33, 88, "DisplayFormula", CellGroupingRules->{GroupTogetherGrouping, 10000.}], Cell[18101, 749, 911, 34, 88, "DisplayFormula", CellGroupingRules->{GroupTogetherGrouping, 10000.}] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[19133, 795, 65, 0, 82, "Subsection"], Cell[CellGroupData[{ Cell[19223, 799, 728, 19, 446, "Input"], Cell[CellGroupData[{ Cell[19976, 822, 350, 7, 130, "Input"], Cell[20329, 831, 990, 29, 108, "Input"] }, Closed]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[21380, 867, 60, 0, 82, "Subsection"], Cell[CellGroupData[{ Cell[21465, 871, 256, 6, 66, "Text", CellGroupingRules->{GroupTogetherGrouping, 10000.}], Cell[21724, 879, 592, 22, 102, "Input"] }, Closed]] }, Closed]] } ] *) (* End of internal cache information *)