一个用于计算协调数的python程序

Question

如何从径向分布函数（RDF）配置文件计算配位数？

我正在尝试从溶质-溶剂相互作用的RDF分布图中获得配位数。 rdf数据有两列，分别对应于r vs g（r）。我使用了一个简单的python代码，将g（r）集成到公式的第一个最小值：

Integral =  r^2 * g(r) * dr
CN = 4 * pi * density * Integral

我在dat文件中具有值g（r）。