Написать программу, моделирующую процесс распространения инфекционного заболевания

#include <iostream>
#include <vector>
#include <conio.h>
#include <ctime>
 
typedef std::vector<int> vi;
typedef std::vector<vi> vvi;
 
class Automat
{
    vvi matrix;
    static const int ill_dur;
    static const int immun_dur;
    static const int ill_chance;
    
public:
    
    Automat(){}
    Automat(int n)
    {
        init(n);
    }
    ~Automat(){}
    
    void init(int n)
    {
        matrix.resize(n, vi(n, 0)); 
    }
    
    bool check_range(int x, int y)
    {
        if( x >= matrix[0].size() || y >= matrix.size() || x < 0 || y < 0) return false;
        return true;
    }
    
    void show()
    {
        for(auto i : matrix)
        {
            for(auto j : i) std::cout << (j > 0 ? "#" : (j < 0 ? "@" : "_")) << " ";
            std::cout << "\n";  
        }
    }
    
    void infect(vvi& dest, vvi& src, int x, int y, int duration = ill_dur, int chance = ill_chance )
    {
        if(!check_range(x, y)) return;
        if((chance >= rand()%100 + 1) && !src[y][x]) dest[y][x] = duration;
    }
    
    void primary_infect(int x, int y)
    {
        infect(matrix, matrix, x, y, ill_dur, 100);
    }
    
    void proceed()
    {
        vvi buff(matrix);
        
        for(int i = 0; i < matrix.size(); ++i)
        {
            for(int j = 0; j < matrix[i].size(); ++j)
            {
                if(matrix[i][j] > 0)
                {
                    infect(buff, matrix, j, i-1, (ill_dur + 1));
                    infect(buff, matrix, j, i+1, (ill_dur + 1));
                    infect(buff, matrix, j-1, i, (ill_dur + 1));
                    infect(buff, matrix, j+1, i, (ill_dur + 1));
                }
            }
        }
        
        for(auto& i : buff)
            for(auto& j : i)
            {
                if(j == 1) { j = immun_dur; continue; }
                if(j > 0) { --j; continue; }
                if(j < 0) ++j;
            }
            
        matrix = buff;
    }
    
    void run()
    {
        do
        {
            system("cls");
            show();
            proceed();
        } while(getch() != 27);
    }
    
};
 
const int Automat::ill_dur = 5;
const int Automat::immun_dur = -4;
const int Automat::ill_chance = 50;
 
int main()
{
    srand(time(0));
    
    int n = 30;
    //std::cin >> n;
 
    Automat A(n);
    
    A.primary_infect(rand()%n, rand()%n);
    A.primary_infect(rand()%n, rand()%n);
    A.primary_infect(rand()%n, rand()%n);
 
    A.run();
    
    return 0;
}