update.jl

using DataFrames

s = DataFrame([[true], [true], [false], [false], [false], [false], [false], [false], [false]],
[:MBF, :SBF, :VOX1, :HCM1, :YHP1, :SFF, :ACE2, :SWI5, :CLN3])

for step in 1:30
    last = s[end, :]
    new = proceed(last)
    push!(s, new)
end

function proceed(row)
    MBF = row.CLN3
    SBF = (row.CLN3 || row.MBF) && !(row.VOX1 || row.YHP1)
    VOX1 = row.MBF && row.SBF
    HCM1 = row.MBF && row.SBF
    YHP1 = row.MBF && row.SBF
    SFF = row.SBF && row.HCM1
    ACE2 = row.SFF
    SWI5 = row.SFF
    CLN3 = (row.SWI5 && row.ACE2) && !(row.VOX1 || row.YHP1)

    [MBF, SBF, VOX1, HCM1, YHP1, SFF, ACE2, SWI5, CLN3]
end