这种计算A_N系数和SINF需要的N项SIR逼近
给定为在巴洛温斯坦，物理学d 408，1（2020）的公式（12a）的（https://doi.org/10.1016/j.physd.2020.132540）
输入对应于SIR参数和初始条件如
由等式（1）中指定。可在Python代码https://github.com/nsbsma/SIR-approximant

%%%例如，在巴洛温斯坦，物理学d 408，1再现图2b（2020）
阿尔法= 2.73;R = 0.0178;I0 = 7;S0 = 254;％的输入参数
t = 0时：0.01：6％的时间间隔为0〜6的0.01的增量
N = 25;在逼近方面的％次数，增加直到答案停止变化
[A，SINF，A0，A] = ApproximantCoefficientsSIR（N，α，R，S0，I0）;使用提供获得所需下面的东西代码％
卡帕= R * SINF-α;
SSUM = 1;
对于j = 1：N
SSUM = SSUM + A（j）的* EXP（卡帕* T * j）的;％近似值（在纸方程4.7A）
结束
S = SINF / SSUM。
情节（T，S， '显示名称'， 'S'）;坚持，稍等
I =阿尔法/ R *日志（S / S 0）+ S0-S + I0;
图（T，I， '显示名称'， 'I'）;
R = S0 + I0-S-I;％B /​​ C S + I + R =常数
情节（T，R， '显示名称'， 'R'）;xlabel（ 'T'）;传奇表演