<!DOCTYPE html>
<html lang="en">
<head>
<meta charset="UTF-8">
<meta name="generator" content="TempleOS V5.03">
<meta name="viewport" content="width=device-width">
<link rel="stylesheet" href="/style/templeos.css">
<script src="/script/templeos.js"></script>
<style type="text/css">
.cF0{color:#000000;background-color:#ffffff;}
.cF1{color:#0000aa;background-color:#ffffff;}
.cF2{color:#00aa00;background-color:#ffffff;}
.cF3{color:#00aaaa;background-color:#ffffff;}
.cF4{color:#aa0000;background-color:#ffffff;}
.cF5{color:#aa00aa;background-color:#ffffff;}
.cF6{color:#aa5500;background-color:#ffffff;}
.cF7{color:#aaaaaa;background-color:#ffffff;}
.cF8{color:#555555;background-color:#ffffff;}
.cF9{color:#5555ff;background-color:#ffffff;}
.cFA{color:#55ff55;background-color:#ffffff;}
.cFB{color:#55ffff;background-color:#ffffff;}
.cFC{color:#ff5555;background-color:#ffffff;}
.cFD{color:#ff55ff;background-color:#ffffff;}
.cFE{color:#ffff55;background-color:#ffffff;}
.cFF{color:#ffffff;background-color:#ffffff;}
</style>
</head>
<body>
<pre id="content">
<a name="l1"></a><span class=cF2>//Classic 2nd order nonlinear ordinary differential equation</span><span class=cF0>
<a name="l2"></a>
<a name="l3"></a></span><span class=cF9>I64</span><span class=cF0> mp_not_done_flags;
<a name="l4"></a>
<a name="l5"></a>#</span><span class=cF1>define</span><span class=cF0> OX 10
<a name="l6"></a>#</span><span class=cF1>define</span><span class=cF0> OY </span><span class=cF3>GR_HEIGHT</span><span class=cF0>/2-</span><span class=cF3>FONT_HEIGHT</span><span class=cF0>
<a name="l7"></a>
<a name="l8"></a>#</span><span class=cF1>define</span><span class=cF0> MICRO_STEP 0.000001
<a name="l9"></a>
<a name="l10"></a></span><span class=cF1>class</span><span class=cF0> State
<a name="l11"></a>{
<a name="l12"></a> </span><span class=cF1>F64</span><span class=cF0> wabbits,hawks,
<a name="l13"></a> d_wabbits,d_hawks; </span><span class=cF2>//Technically, these are not part of the state.</span><span class=cF0>
<a name="l14"></a>};
<a name="l15"></a>
<a name="l16"></a></span><span class=cF1>U0</span><span class=cF0> Steps(State *s)
<a name="l17"></a>{
<a name="l18"></a> </span><span class=cF9>I64</span><span class=cF0> i;
<a name="l19"></a> </span><span class=cF1>for</span><span class=cF0> (i=0;i<10000;i++) </span><span class=cF7>{</span><span class=cF0>
<a name="l20"></a> s->d_wabbits=s->wabbits*(1.0-0.10*s->hawks);
<a name="l21"></a> s->d_hawks =-s->hawks *(1.0-0.01*s->wabbits);
<a name="l22"></a> s->hawks +=s->d_hawks *MICRO_STEP;
<a name="l23"></a> s->wabbits+=s->d_wabbits*MICRO_STEP;
<a name="l24"></a> </span><span class=cF7>}</span><span class=cF0>
<a name="l25"></a>}
<a name="l26"></a>
<a name="l27"></a></span><span class=cF1>U0</span><span class=cF0> PlotTrajectory(</span><span class=cF9>CTask</span><span class=cF0> *parent)
<a name="l28"></a>{
<a name="l29"></a> </span><span class=cF9>I64</span><span class=cF0> i;
<a name="l30"></a> </span><span class=cF9>CDC</span><span class=cF0> *dc=</span><span class=cF5>DCAlias</span><span class=cF0>(,parent);
<a name="l31"></a> State s;
<a name="l32"></a>
<a name="l33"></a> </span><span class=cF5>MemSet</span><span class=cF0>(&s,0,</span><span class=cF1>sizeof</span><span class=cF7>(</span><span class=cF0>State</span><span class=cF7>)</span><span class=cF0>);
<a name="l34"></a> s.wabbits=</span><span class=cF5>RandU16</span><span class=cF0>%100 +10;
<a name="l35"></a> s.hawks =</span><span class=cF5>RandU16</span><span class=cF0>%10 +1;
<a name="l36"></a>
<a name="l37"></a> </span><span class=cF1>for</span><span class=cF0> (i=0;i<1000;i++) </span><span class=cF7>{</span><span class=cF0>
<a name="l38"></a> dc->color=</span><span class=cF3>LTGREEN</span><span class=cF0>;
<a name="l39"></a> </span><span class=cF5>GrPlot</span><span class=cF0>(dc,s.wabbits+OX,OY-s.d_wabbits);
<a name="l40"></a> dc->color=</span><span class=cF3>LTRED</span><span class=cF0>;
<a name="l41"></a> </span><span class=cF5>GrPlot</span><span class=cF0>(dc,s.hawks +OX,OY-s.d_hawks);
<a name="l42"></a> Steps(&s);
<a name="l43"></a> </span><span class=cF7>}</span><span class=cF0>
<a name="l44"></a> </span><span class=cF5>DCDel</span><span class=cF0>(dc);
<a name="l45"></a> </span><span class=cF5>LBtr</span><span class=cF0>(&mp_not_done_flags,</span><span class=cF5>Gs</span><span class=cF0>->num);
<a name="l46"></a>}
<a name="l47"></a>
<a name="l48"></a></span><span class=cF1>U0</span><span class=cF0> PredatorPrey()
<a name="l49"></a>{
<a name="l50"></a> </span><span class=cF9>I64</span><span class=cF0> i;
<a name="l51"></a> </span><span class=cF9>CDC</span><span class=cF0> *dc=</span><span class=cF5>DCAlias</span><span class=cF0>;
<a name="l52"></a>
<a name="l53"></a> </span><span class=cF5>PopUpOk</span><span class=cF0>(</span><span class=cF6>"This will plot multiple predator-prey\n"</span><span class=cF0>
<a name="l54"></a> </span><span class=cF6>"trajectories. It restarts many times\n"</span><span class=cF0>
<a name="l55"></a> </span><span class=cF6>"with different, random, initial populations.\n"</span><span class=cF0>);
<a name="l56"></a> </span><span class=cF5>SettingsPush</span><span class=cF0>; </span><span class=cF2>//See </span><a href="/Wb/Adam/TaskSettings.HC.HTML#l3"><span class=cF4>SettingsPush</span></a><span class=cF0>
<a name="l57"></a> </span><span class=cF1>try</span><span class=cF0> </span><span class=cF7>{</span><span class=cF0>
<a name="l58"></a> </span><span class=cF5>AutoComplete</span><span class=cF0>;
<a name="l59"></a> </span><span class=cF5>WinBorder</span><span class=cF0>;
<a name="l60"></a> </span><span class=cF5>WinMax</span><span class=cF0>;
<a name="l61"></a> </span><span class=cF5>DocClear</span><span class=cF0>;
<a name="l62"></a> </span><span class=cF5>Refresh</span><span class=cF0>;
<a name="l63"></a>
<a name="l64"></a> dc->color=</span><span class=cF3>BLACK</span><span class=cF0>;
<a name="l65"></a> </span><span class=cF5>GrLine</span><span class=cF0>(dc,OX,0,OX,</span><span class=cF3>GR_HEIGHT</span><span class=cF0>-</span><span class=cF3>FONT_HEIGHT</span><span class=cF0>-1);
<a name="l66"></a> </span><span class=cF5>GrLine</span><span class=cF0>(dc,0,OY,</span><span class=cF3>GR_WIDTH</span><span class=cF0>-1,OY);
<a name="l67"></a> </span><span class=cF1>while</span><span class=cF0> (!</span><span class=cF5>ScanChar</span><span class=cF0>) {
<a name="l68"></a> mp_not_done_flags=1<<</span><span class=cFB>mp_cnt</span><span class=cF0>-1;
<a name="l69"></a> </span><span class=cF1>for</span><span class=cF0> (i=0;i<</span><span class=cFB>mp_cnt</span><span class=cF0>;i++)
<a name="l70"></a> </span><span class=cF5>JobQue</span><span class=cF0>(&PlotTrajectory,</span><span class=cF5>Fs</span><span class=cF0>,i);
<a name="l71"></a> </span><span class=cF1>do</span><span class=cF0> </span><span class=cF5>Yield</span><span class=cF0>;
<a name="l72"></a> </span><span class=cF1>while</span><span class=cF0> (mp_not_done_flags);
<a name="l73"></a> }
<a name="l74"></a> </span><span class=cF7>}</span><span class=cF0> </span><span class=cF1>catch</span><span class=cF0>
<a name="l75"></a> </span><span class=cF5>PutExcept</span><span class=cF0>;
<a name="l76"></a> </span><span class=cF5>SettingsPop</span><span class=cF0>;
<a name="l77"></a> </span><span class=cF5>DCFill</span><span class=cF0>(dc);
<a name="l78"></a> </span><span class=cF5>DCDel</span><span class=cF0>(dc);
<a name="l79"></a>}
<a name="l80"></a>
<a name="l81"></a>PredatorPrey;
</span></pre></body>
</html>