def seir(seir, t, b, a, g, N):
    """
    seir: seir[0]: S, seir[1]: E, seir[2]: I, seir[3]: R
    t: 時間
    b: 感染率
    a: 発症率
    g: 回復率
    N: 全人口
    """
    return [
        -b * seir[0] * seir[2] / N,
        b * seir[0] * seir[2] / N - a * seir[1],
        a * seir[1] - g * seir[2],
        g * seir[2]
    ]