SIR model with Maxima - Revision history

Mol: /* Problems */

2021-08-12T10:40:48Z

Problems

← Older revision		Revision as of 13:40, 12 August 2021
Line 29:		Line 29:
	=== Problems ===		=== Problems ===

	It was difficult to make <math>\frac{d Ns}{dt} = N \frac {ds}{dt}</math>, but <syntaxhighlight inline>ev(%, nouns)</syntaxhighlight> helps.		It was difficult to make <math>\frac{d Ns}{dt} = N \frac {ds}{dt}</math>, but <syntaxhighlight lang="Bash" inline>ev(%, nouns)</syntaxhighlight> helps.

	Also, Greek letters are not supported for <syntaxhighlight inline>plotdf</syntaxhighlight>, and note that no spaces are allowed in <syntaxhighlight inline>parameters</syntaxhighlight> option.		Also, Greek letters are not supported for <syntaxhighlight lang="Bash" inline>plotdf</syntaxhighlight>, and note that no spaces are allowed in <syntaxhighlight lang="Bash" inline>parameters</syntaxhighlight> option.

Mol: /* Problems */

2021-08-12T10:40:05Z

Problems

← Older revision		Revision as of 13:40, 12 August 2021
Line 29:		Line 29:
	=== Problems ===		=== Problems ===

	It was difficult to make <math>\frac{d Ns}{dt} = N \frac {ds}{dt}</math>, but <~~pre~~ inline>ev(%, nouns)</~~pre~~> helps.		It was difficult to make <math>\frac{d Ns}{dt} = N \frac {ds}{dt}</math>, but <syntaxhighlight inline>ev(%, nouns)</syntaxhighlight> helps.

	Also, Greek letters are not supported for <~~pre~~ inline>plotdf</~~pre~~>, and note that no spaces are allowed in <~~pre~~ inline>parameters</~~pre~~> option.		Also, Greek letters are not supported for <syntaxhighlight inline>plotdf</syntaxhighlight>, and note that no spaces are allowed in <syntaxhighlight inline>parameters</syntaxhighlight> option.

Mol: /* Problems */

2021-08-12T10:39:32Z

Problems

← Older revision		Revision as of 13:39, 12 August 2021
Line 29:		Line 29:
	=== Problems ===		=== Problems ===

	It was difficult to make <math>\frac{d Ns}{dt} = N \frac {ds}{dt}</math>, but <pre>ev(%, nouns)</pre> helps.		It was difficult to make <math>\frac{d Ns}{dt} = N \frac {ds}{dt}</math>, but <pre inline>ev(%, nouns)</pre> helps.

	Also, ~~greek~~ letters are not supported for <pre>plotdf</pre>, and note that no spaces are allowed in <pre>parameters</pre> option.		Also, Greek letters are not supported for <pre inline>plotdf</pre>, and note that no spaces are allowed in <pre inline>parameters</pre> option.

Mol: /* Problems */

2021-08-12T10:38:56Z

Problems

← Older revision		Revision as of 13:38, 12 August 2021
Line 29:		Line 29:
	=== Problems ===		=== Problems ===

	It was difficult to make <math>\frac{d Ns}{dt} = N \frac {ds}{dt}</math>, but <pre>ev, nouns</pre> helps.		It was difficult to make <math>\frac{d Ns}{dt} = N \frac {ds}{dt}</math>, but <pre>ev(%, nouns)</pre> helps.

	Also, greek letters are not supported for <pre>plotdf</~~math~~>, and note that no spaces are allowed in <pre>parameters</pre> option.		Also, greek letters are not supported for <pre>plotdf</pre>, and note that no spaces are allowed in <pre>parameters</pre> option.

Mol: /* Problems */

2021-08-12T10:38:25Z

Problems

← Older revision		Revision as of 13:38, 12 August 2021
Line 29:		Line 29:
	=== Problems ===		=== Problems ===

	It was difficult to make <math>\frac{d Ns}{dt} = N \frac {ds}{dt}, but ev, ~~noun~~ helps.		It was difficult to make <math>\frac{d Ns}{dt} = N \frac {ds}{dt}</math>, but <pre>ev, nouns</pre> helps.

	Also, greek letters are not supported for <pre>plotdf</math>, and note that no spaces are allowed in <pre>parameters</pre> option.		Also, greek letters are not supported for <pre>plotdf</math>, and note that no spaces are allowed in <pre>parameters</pre> option.

Mol: /* Theory */

2021-08-12T10:38:01Z

Theory

← Older revision		Revision as of 13:38, 12 August 2021
Line 26:		Line 26:
	plotdf([subst( [β=b], rhs(eqs)), subst( [α=a, β=b], rhs(eqi))], [s,i], [parameters,"a=10000,b=1.1,N=1000"], [s,0,75], [i,0,60], [sliders,"b=0.5:2.5"] );		plotdf([subst( [β=b], rhs(eqs)), subst( [α=a, β=b], rhs(eqi))], [s,i], [parameters,"a=10000,b=1.1,N=1000"], [s,0,75], [i,0,60], [sliders,"b=0.5:2.5"] );
	</syntaxhighlight>		</syntaxhighlight>

			=== Problems ===

			It was difficult to make <math>\frac{d Ns}{dt} = N \frac {ds}{dt}, but ev, noun helps.

			Also, greek letters are not supported for <pre>plotdf</math>, and note that no spaces are allowed in <pre>parameters</pre> option.

Mol: /* Theory */

2021-08-12T10:36:06Z

Theory

← Older revision		Revision as of 13:36, 12 August 2021
Line 7:		Line 7:
	[[File:Sir.png\|thumb\|SIR model]]		[[File:Sir.png\|thumb\|SIR model]]

	<syntaxhighlight>		<syntaxhighlight lang="powershell">
	depends([S,I,s, i],t);		depends([S,I,s, i],t);

Line 25:		Line 25:

	plotdf([subst( [β=b], rhs(eqs)), subst( [α=a, β=b], rhs(eqi))], [s,i], [parameters,"a=10000,b=1.1,N=1000"], [s,0,75], [i,0,60], [sliders,"b=0.5:2.5"] );		plotdf([subst( [β=b], rhs(eqs)), subst( [α=a, β=b], rhs(eqi))], [s,i], [parameters,"a=10000,b=1.1,N=1000"], [s,0,75], [i,0,60], [sliders,"b=0.5:2.5"] );
	</~~syntaxhiglight~~>		</syntaxhighlight>

Mol: Created page with "== Introduction == Use WxMaxima to calculate SIR model. == Theory == SIR model depends([S,I,s, i],t); eqS : 'diff(S, t) = -βS..."

2021-08-12T10:35:01Z

Created page with "== Introduction == Use WxMaxima to calculate SIR model. == Theory == thumb|SIR model <syntaxhighlight> depends([S,I,s, i],t); eqS : 'diff(S, t) = -β*S*..."

New page

== Introduction ==

Use WxMaxima to calculate SIR model.

== Theory ==

[[File:Sir.png|thumb|SIR model]]

<syntaxhighlight>
depends([S,I,s, i],t);

eqS : 'diff(S, t) = -β*S*I;
eqI : 'diff(I, t) = β*S*I - α*I;

factor( rhs(eqI)>0);

eqsN : subst([S=N*s, I=N*i], eqS );
eqiN : subst([S=N*s, I=N*i], eqI );

eqs1 : ev(eqsN, nouns);
eqi1 : ev(eqiN, nouns);

eqs : eqs1/N;
eqi : ratsimp( eqi1/N );

plotdf([subst( [β=b], rhs(eqs)), subst( [α=a, β=b], rhs(eqi))], [s,i], [parameters,"a=10000,b=1.1,N=1000"], [s,0,75], [i,0,60], [sliders,"b=0.5:2.5"] );
</syntaxhiglight>

SIR model with Maxima - Revision history

Mol: /* Problems */

Mol: /* Problems */

Mol: /* Problems */

Mol: /* Problems */

Mol: /* Problems */

Mol: /* Theory */

Mol: /* Theory */

Mol: Created page with "== Introduction == Use WxMaxima to calculate SIR model. == Theory == SIR model depends([S,I,s, i],t); eqS : 'diff(S, t) = -β*S*..."

Mol: Created page with "== Introduction == Use WxMaxima to calculate SIR model. == Theory == SIR model depends([S,I,s, i],t); eqS : 'diff(S, t) = -βS..."