pr5

import numpy as np
import matplotlib.pyplot as plt

alpha = 0.95
Tc = 290
A = 16
B = -230
D = 900

eta = np.linspace(0, 0.4, 200000)
T_values = [286, 289.3, 288.5, 284, 292]
plt.figure(1)
plt.title("Зміна вільної енергії для різних температур")
for T in T_values:
    df = alpha * (T - Tc) * eta**2 + (A / 2) * eta**2 + (B / 3) * eta**3 + (D / 4) * eta**4
    plt.plot(eta, df, label=f'T = {T} K')
plt.grid()
plt.legend()
# plt.xlim([0, 0.4])
plt.ylim([-0.2, 0.2])

plt.figure(2)
plt.title('Залежність eta від температури')
t = np.linspace(0, 450, 200000)
Dis = (B**2 - 4 * D * (2 * alpha * (t-Tc) + A))
eta1=(-B+np.sqrt(Dis))/(2*D)
eta2=(-B-np.sqrt(Dis))/(2*D)
eta3=np.sqrt((-2*alpha*(t-Tc) - A)/D)
eta4=-np.sqrt((-2*alpha*(t-Tc) - A)/D)
plt.plot(t, eta1, 'b')
plt.plot(t, eta2, 'b--')
plt.plot(t, eta3, 'r')
plt.plot(t, eta4, 'r--')
plt.grid()
plt.xlim( [0, 450])
plt.ylim([0, 1.1])
plt.show()