求助一个无助的人问题

wnmdyx

计算用参数：ρ=600kg/m3 、c=2830J/(kg·K) 、λ=0.145W/(m2·K) 、T0=300K 、Tf=T∞=500K ；(1)对流换热系数为常数时取：h=100W/(m2·K) 、500100W/(m2·K)、2000100W/(m2·K) 三种、固球半径取rd=100μm ，可构成几个不同的Bi数（Bi=(h·rd)/λ ）。将计算结果转换成无量纲形式画曲线图：
温度：Θ=(T-T∞)/(T0-T∞) ，时间：Fo=(aτ)/rd2 ，a=λ/(ρc) 是热扩散率；径向坐标:r※=r/rd
传热学教科书中（杨世铭、陶文铨第四版）P125，式(3-23a)～(3-23c).
         该方法仅对换热系数为常数的情况适用。
工作内容：将传热学教科书中的公式变为计算程序，其无穷级数至少取前10项，并通过与诺谟图比较，检验程序结果。可选5个点（表面、球心、二分之一半径、四分之一半径、四分之三半径处），检查这5个点在 Fo=(aτ)/rd2 ＝0.1、0.2、0.5、1.0、2.0 五个时刻的无量纲温度与诺模图曲线结果的符合程度，可列表比较说明。