Talk:Matlab:Gradual dimming population model

From Puella Magi Wiki
Jump to navigation Jump to search
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.
Note: Please always sign your name when editing talk page by putting four tildes (~~~~) at the end of your comment.

Can I get a few extra simulations with the below values?

First try

C = 10;    % Number of girls contracted by QBe
D = 0.01;    % Proportion of girls (who are fighting) that die
B = 0.1;   % Proportion of girls (who are fighting) to become witches
F = 0.5;   % Proportion of familiars becoming witches
K = 0.2;  % Proportion of witches getting killed by MSes
P = 0.95;   % proportion of MSes fightgins
T = 0.0001; % Number of MSes turning into witches over time

M(1) = 0;   % Number of magical girls at first
W(1) = 0;   % Number of witches at first

Result


Second try

C = 10;    % Number of girls contracted by QBe
D = 0.01;    % Proportion of girls (who are fighting) that die
B = 0.1;   % Proportion of girls (who are fighting) to become witches
F = 0.05;   % Proportion of familiars becoming witches
K = 0.2;  % Proportion of witches getting killed by MSes
P = 0.95;   % proportion of MSes fightgins
T = 0.0001; % Number of MSes turning into witches over time

M(1) = 0;   % Number of magical girls at first
W(1) = 0;   % Number of witches at first

Result


Third try

C = 10;    % Number of girls contracted by QBe
D = 0.01;    % Proportion of girls (who are fighting) that die
B = 0.1;   % Proportion of girls (who are fighting) to become witches
F = 0.5;   % Proportion of familiars becoming witches
K = 0.2;  % Proportion of witches getting killed by MSes
P = 0.95;   % proportion of MSes fightgins
T = 0.01; % Number of MSes turning into witches over time

M(1) = 0;   % Number of magical girls at first
W(1) = 0;   % Number of witches at first

Result

Another interresting result

Here's another result. Third model, in the case K>B+F. The value at which the system seems to become unstable seems to be around t=0.0003595. I have absolutely no clue about such a behavior ...

Also note that reducing T has the effect of reducing the first gap's length, the others don't seem to vary that much ...

Left is with T=0.000359, the system seems to converge to its equilibrium. Right is with T=0.0003596, the system goes crazy as usual.