Simulated Bioterrorism Attack in SWI-Prolog

This program must be written in SWI-Prolog, and no other language.

Use http://www.sciencedirect.com/science/article/pii/S1201971206000968 as a basis to craft the program after.

You will implement a Time stepping algorithm where each time step divided into phases, each phase representing time spent with various mixing groups

Mixing Groups are collections of individuals and are defined per individual are as follows:

Family Group with size range 1-7

Neighborhood

Two to four random groups (non-family, non-residential)

An initial list of individuals will be created such that

Each individual gets a resistance score with mean of 1.0 and std dev of .15

Initial status of Non-infected, contaminated set to false, disease progression state (incubation, fever, rash, pustules)

Date of initial infection

Generate a number of families for a neighborhood

Sequentially generate the families within the neighborhood

Each person will then assigned two additional random groups

Each group is assigned a level of contact score

Family groups are assigned .05

Other groups are r*.0625*.05 where r is random int 1 to 16

A single non-family mixing group is chosen as the target.

Within the initial target group each individual is then calculated to possibly be infected & contaminated or contaminated or neither

At this point the simulation is ready

In each time period each group is evaluated to see if an individual is either contaminated or infectious

At the end of the time period the contamination score is to 1/(time*time)*initial score

Each individual has a chance of being contaminated by a contaminated individual

Each individual has a chance of being infected by a infectious individual

Governing equation for spread of disease:

The Product term refers to the virus propagation in the first equation and the scores of each mixing in the second

When any individual displays pustules then there is a 52.5% chance the CDC is considered notified and 48 hours later containment protocols are implemented.

After that the simulation stops all data is dumped to file for prolog analysis.

Id number

List of groups for that id

Disease progression status of not-exposed, exposed, infected, symptomatic, diagnosed and dead

In effect you will create a record of the progression of the disease from initial attack until the CDC confirms the outbreak.

For example:

.

Day 15: 4787 Not Exposed, 198 Exposed, 18 symptomatic, 3 dead

.etc

In addition you need to provide functions that will allow the user to examine the status of any individual or any mixing group during runtime using a set of command line options including:

Menu

Do time step

Display results for current day

Display status for a particular person PID

Display status for a particular mixing group MGID

Contamination 1-H(1 Score number of contaminated es stance Infected 1 Score )contaminated +infectious Contamination 1-H(1 Score number of contaminated es stance Infected 1 Score )contaminated +infectious