h2.append(H2)
# t ning vaqt bo'yicha o'zgarish qiymatlari
T = []
T.append(T0)
# P10 ning vaqt bo'yicha o'zgarish qiymatlari
p1010 = []
# P11 ning vaqt bo'yicha o'zgarish qiymatlari
p1111 = []
h11 = H1
h22 = H2
t = 0
p88 = []
p99 = []
v = []
while i<=n:
p10 = (po*H1G)/(H1G-h11)
p8 = p10 + p*g*h11
p88.append(p8)
v1 = k1*pow(abs((p1-p8)), 0.5)*sign(p1, p8)
v2 = k2*pow(abs((p2-p8)), 0.5)*sign(p2, p8)
v3 = k3*pow(abs((p3-p8)), 0.5)*sign(p3, p8)
v6 = k6*pow(abs((p8-p6)), 0.5)*sign(p8, p6)
v7 = k7*pow(abs((p8-p7)), 0.5)*sign(p8, p7)
p11 = (po*H2G)/(H2G-h22)
p9 = p*g*h22+p8
p99.append(p9)
v4 = k4*pow(abs((p9-p4)), 0.5)*sign(p9, p4)
v5 = k5*pow(abs((p9-p5)), 0.5)*sign(p9, p5)
v8 = k8*pow(abs((p9-p8)), 0.5)*sign(p9, p8)
h11 = h11+delT*(v1+v2+v3-v6-v7-v8)/S1
h22 = h22+delT*(v8-v4-v5)/S2
h1.append(round(h11, 3))
h2.append(round(h22, 3))
i = i+1
t = t+delT
T.append(t)
v.append(t)
# h1 va t qiymatlarini grafikga joylanishi
plt.plot(T, h2)
# t ning qiymatlari
plt.xlabel('t - ning qiymatlari')
# h1 ning qiymatlari
plt.ylabel('h1 - ning qiymatlari')
# Grafik haqida malumot
plt.title('h1 ning t ga bog\'liqlik grafigi!')
# Natijani ekranga chop etish
plt.show()
# h2 va t qiymatlarini grafikga joylanishi
plt.plot(T, h2)
# t ning qiymatlari
plt.xlabel('t - ning qiymatlari')
# h2 ning qiymatlari
Do'stlaringiz bilan baham: |