多元弱酸的分布曲线图

代码

制作多元弱酸的分布曲线，多元弱碱改一下里面的对应关系即可

#2.0增加交点
import numpy as np
import math
import matplotlib.pyplot as plt



n = eval(input("请输入n元弱酸:"))
name = input("请输入酸的全称（如H2CO3）：")
X = input("请输入酸根离子（如CO3）：")
Kn = []
i = 0
while(n>i):
    Kn.append(eval(input("请输入K"+str(i+1)+":")))
    i=i+1

A = []
for i in range(n+1):
    A.append("H"+str(n-i)+X+str(i)+"-")
print(A)

def draw(X,Y,label):
    plt.plot(X,Y,label=label)


def count(n,X,cH,Kn):
    sum = cH**(n-X+1)
    h = X
    while(h>1):
        h -= 1
        sum *= Kn[h-1]
    return sum

step = 0.01
pH = np.arange(0, 14, step)
cH = 10 ** (-pH)
temp = 0
for i in range(1, n + 2):
    temp += count(n, i, cH, Kn)

len  = int(14/0.01)
cross = []
Ylist = []
min = 1.5
def cross1():
    for i in range(1,n+2):
        Ylist.append(count(n, i, cH, Kn) / temp)
    for i in range(n):
        print(i)
        for j in range(len):
            t = []
            if(abs(Ylist[i][j]-Ylist[i+1][j])<0.0055 and Ylist[i][j]>0.1 and Ylist[i+1][j]>0.1):
                t.append(j * 0.01)
                t.append(Ylist[i][j]/2+Ylist[i+1][j]/2)
                cross.append(t)
                break

cross1()
print(cross)

def ffmu(index):
    Y = count(n,index,cH,Kn)/temp
    draw(pH,Y,label=A[index-1])


for i in range(1,n+2):
    ffmu(i)

for i in range(n):
    plt.scatter(cross[i][0],cross[i][1],s=50,label=str(cross[i][0])+",{:.2f}".format(cross[i][1]))
plt.xlabel("PH")
plt.ylabel("δ")
plt.ylim(0, 1)
plt.xlim(0,14)
plt.title(str(name)+"δ-pH")
plt.legend()
plt.show()

效果

H2CO3的分布曲线


HAc

HCN

酒石酸

磷酸