蒙特卡罗方法(三)：铀235裂变中子谱(瓦特、maxwell)(近似修正、乘减)、方位角、质心系与实验室系散射角、中子弹散射后角度-能量分布

目录内容

1. 铀235裂变中子谱（瓦特谱） 近似修正抽样
2. 铀235裂变中子谱（瓦特谱） 乘减抽样
3. 方位角余弦phi 均匀分布
4. 质心系各向同性散射角余弦theta抽样、转换成实验室系余弦
5. 中子弹性散射后角度、能量
6. 铀235裂变中子谱（maxwell谱）

铀235裂变中子谱（瓦特谱） 近似修正抽样

    '''铀235裂变中子谱（瓦特谱）。P70'''
    # 近似修正抽样P70，比乘减抽样计算量小、效率高。
    a = 0.965
    b = 2.29
    c = 0.453
    Emin = 0
    Emax = inf   # 无穷大
    m = 0.8746
    M = 0.2678
    lamda = 0.5543
    Earr = np.array([])
    for i in range(100000):
        s1 = uniform()
        if s1<=m:
            Efa = -a*log(uniform()*uniform())
            Earr = np.append(Earr,Efa)
        else:
            while True:
                s2 = uniform()
                s3 = uniform()
                Ef1 = -log(s3)*a/(1-a*lamda)
                H1 = a/(1-a