Benchmark Model 2b

KModel SIRDemography
	attribute: #(status -> S I R);
	parameters: #(beta lambda gamma b mu);
	transitions: #(
		S -- lambda --> I.
		I -- gamma --> R.
		status -- mu --> Empty.
		Empty -- b --> S.
	);
	lambda: #(beta*I/N).
	
KModel AgeStructure
	attribute: #(age -> J A);
	parameters: #(m);
	transitions: #(
		J -- m --> A.
	).
	
Composition Model2
	model: 'SIRDemography';
	model: 'AgeStructure'.
	
Scenario M2Parameters
	 on: 'Model2';
         beta: '0.1*sin(t)';
         gamma: 1/30; 
	 mu_age: {0.01. 0.0025};
	 m: 1/50;
	 b_age: #(0.005 0);
	 lambda: #(beta*I/N sum);
	 S_age: #(0 90);
	 I_age: #(1 9).
	
Simulation M2Sim gillespie
	scenarios: #(M2Parameters);
	from: 0.0; 
	to: 365; 
	step: 0.5.
	
Visualization M2Vis diagram 
	for: 'M2Sim';
	xLabel: 'Time (days)';
        open.