Python 代码片段

计算收敛阶

import numpy as np

def order(h1, h2, E1, E2):
    return np.log(E1/E2) / np.log(h1/h2)

p = order(h1=1/32, h2=1/64, E1=0.01, E2=0.0025)
print(p)