具体内容请参考： https://www.tabirstrees.top/2022/02/12/Eorupainternalstructure/
parameter.py中更改参数（一般是星球的质量分数）；
ode_model.ipynb是计算的主程序（迭代）以及绘制最后的图表；
eos_no_tem.py是定义计算函数的地方，求解微分方程的位置，具体的计算方式可以看PPT中的描述；
core.txt,curst.txt,low_mantle.txt,mantle.txt分别是定义各个圈层的热力学参数的文件
欢迎讨论！！！(*^_^*)