编辑代码

# coding:utf-8
#JSRUN引擎2.0，支持多达30种语言在线运行，全仿真在线交互输入输出。 
print("Hello world!   -  python.jsrun.net .")

TR0 = 0.39
CR0 = 2.669283

def FT(m):
    a = 1.5
    b = -1.329
    c = 0.829
    return a * m**2 + b * m + c

def FC(m):
    a = -0.732
    b = 0.845
    c = 0.886
    return a * m**2 + b * m + c

dt = 0.001
Pcin = 7492.649
norm_old = 1.0
norm = 1.0
Pt_out0 = 7637.135
print("Pci    Pco    CR     Pti     Pto     TR     nor_m")
for i in range(100):
    CR = FC(norm)*CR0
    Pco = Pcin*CR

    Ptin = Pco*0.99*0.99
    TR = FT(norm)*TR0
    Pto = Ptin*TR

    dp = Pto - Pt_out0
    norm = norm_old + dt*dp*1000*3.1415926*0.01*0.01/4/5
    norm_old = norm
    # Pcin = Pto*0.99*0.99

    print("{:10.3f} {:10.3f} {:10.3f} {:10.3f} {:10.3f} {:10.3f} {:10.6f}".format(Pcin,Pco ,CR,Ptin,Pto,TR,norm))