{"nbformat":4,"nbformat_minor":0,"metadata":{"colab":{"name":"Fiche Composition chimique du système dans l’état final.ipynb","provenance":[{"file_id":"1BhdR1vrHIaE53MTxNttNP78U3tOKWTK5","timestamp":1655447008700}],"collapsed_sections":[],"authorship_tag":"ABX9TyMXbu5pAGjX651kHZvz91Kz"},"kernelspec":{"name":"python3","display_name":"Python 3"},"language_info":{"name":"python"}},"cells":[{"cell_type":"markdown","source":["**Déterminer, à l’aide d’un langage de programmation, l’état final d’un système à l’équilibre, siège d’une transformation, modélisée par une ou deux réactions à partir des conditions initiales et valeur(s) de la(es) constante(s) thermodynamique(s) d’équilibre.**"],"metadata":{"id":"rqsA-ngAVU7Z"}},{"cell_type":"markdown","metadata":{"id":"mpgaFvrX0KNZ"},"source":["Déterminer l'avancement volumique à l’équilibre d'une réaction."]},{"cell_type":"code","metadata":{"colab":{"base_uri":"https://localhost:8080/","height":297},"id":"gQgWE6ucjK-J","executionInfo":{"status":"ok","timestamp":1630508219476,"user_tz":-120,"elapsed":890,"user":{"displayName":"Lionnel Malara","photoUrl":"https://lh3.googleusercontent.com/a-/AOh14Gg9uzWMkTid1qlQikWp7IMi29R7n5oSh3ECmY5O=s64","userId":"16342666118884215780"}},"outputId":"6a9108d2-c9f2-467f-c7e7-73bdefb49039"},"source":["import matplotlib.pyplot as plt #Pour tracer des graphiques\n","import numpy as np #Pour faire divers calculs\n","import scipy.optimize as op #Importe la bibliothèque pour faire la résolution d'équation\n","\n","\n","## Constantes d'equilibre réaction\n","\n","K=10**(-1) # valeur de K\n","C_0=10**(-1) # valeur de C0 en mol/L\n","\n","## Quotient reactionnel réaction\n","\n","def Qr(ksi):\n","    # Renvoie la valeur du quotient reactionnel de la reaction pour l'avancements ksi\n","    \n","    return (ksi**2)/(C_0-ksi)\n","\n","## Determination de l'avancements à l'equilibre : resolution de l'equation\n","\n","def f(ksi): # chercher l'avancement a l'equilibre revient a chercher le zero de Qr - K\n","    \n","    return Qr(ksi)-K\n","\n","## Trace de f(ksi)\n","n=100 #nombre de points de l'intervalle 0;C0\n","ksi_i=C_0/n #premier point\n","ksi_f=C_0*(n-1)/n #dernier point\n","ksi = np.linspace(ksi_i,ksi_f,n) # liste des avancements entre 0 et C0\n","\n","y = []# initialisation du np.array des y\n","\n","for x in ksi: #boucle de calcul des y sur l'intervale de ksi\n","\n","    y.append(f(x))# calcul et stockage de la fonction Qr-K dans y\n","\n","plt.figure(1)\n","plt.xlabel('ksi') #Légende de l’axe des abscisses\n","plt.ylabel('f') #Légende de l’axe des ordonnées\n","plt.plot(ksi,y,'b-',label='Qr(ksi)-K°') #Représente y en fonction de x en bleu de manière continue avec une étiquette pour la courbe\n","plt.legend() #Affiche l’étiquette de la courbe\n","plt.grid() #Affiche le quadrillage\n","plt.show() #Affiche le graphique\n","\n","## Determination de l'avancements à l'equilibre : resolution de l'equation\n","\n","def f(ksi):# chercher l'avancement a l'equilibre revient a chercher le zero de Qr - K\n","\n","    return Qr(ksi)-K\n","\n","ksi_eq=op.fsolve(f,C_0/2) #Résolution de la fonction f=0 avec la valeur initiale C0/2\n","\n","print('ksi_eq = ',ksi_eq[0]) #Affiche le résultat\n","   "],"execution_count":null,"outputs":[{"output_type":"display_data","data":{"image/png":"iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAbm0lEQVR4nO3de3hV9Z3v8fc3F5JAAANCWkEM8KAOokATxYPOHBA92uKIcqyntQ7qqDyj5VRbW2daznSm6PTiWB2ttg5jraW1jUcEi2gVR4NajmMlcmkBlYsoUUFBESIkkPA9f6wdyBWy2XuvtXfW5/U869lrr8v+fX8JflxZ67fXMndHRETiIy/qAkREJFwKfhGRmFHwi4jEjIJfRCRmFPwiIjFTEHUB3XHsscd6RUVFt7f/9NNP6dOnT+YKylLqd7yo3/GTbN9ra2u3u/ug9sszFvxm9iBwIfCBu49JLBsAPAJUAJuBy9z94yN9VkVFBcuXL+9220uXLmXSpEnJF53j1O94Ub/jJ9m+m9nbnS3P5Kmeh4AL2i37B+A5dx8FPJd4LyIiIcpY8Lv7i8BH7RZPA36ZmP8lcHGm2hcRkc5ZJr+5a2YVwOJWp3p2uvsxiXkDPm5538m+M4GZAOXl5ZXV1dXdbre+vp7S0tLUis9B6ne8qN/xk2zfJ0+eXOvuVe2XR3Zx193dzLr8v467zwXmAlRVVXn781r79++nrq6OhoaGDvv279+f4uLi9BacA1Lpd3FxMUOHDqWwsDDNVWVeXM/5qt/xk66+hx3828zss+7+vpl9FvjgaD+orq6Ovn37UlFRQfDHwyG7d++mb9++qdaac4623+7Ojh07qKurY/jw4RmoTESySdjj+BcBVybmrwR+d7Qf1NDQwMCBAzuEviTPzBg4cGCnfz2JSM+TseA3s98CLwMnmVmdmV0D/BA4z8zWA+cm3qfSRuqFCqCfpUicZOxUj7t/uYtVUzLVpohIT/GnP8Gjj8KsWTB4cHo/W7dsEBHJQqtXw623ws6d6f9sBX8K6urqmDZtGqNGjWLEiBHMmjWLxsbGTrd9//33ufDCCwF46KGHmDVrVrfaWL58OV/72tcAWLx4Md/97ne73HbMmDFs374dgNraWoYPH86KFSs4cOAAM2bMYOLEiaxZsyaZLopIRFqiJBMDFBX8R8ndmT59OhdffDHr169n/fr17N27l1tuuaXNdk1NTQDceeedXHfddUm3U1VVxT333APA1KlTeeKJJ9izZ89h91m9ejWXXnopjzzyCOPHj2fJkiVMmDCBhQsX8uMf/zjpGkQkfC1jLTIR/Dlxk7YjuekmWLny0Pvm5hLy81P7zHHj4N/+rev1zz//PMXFxVx99dUA5Ofnc9ddd3HCCScwatQolixZQn19Pc3Nzbzwwgs89thj3HbbbR0+58knn+S2227jiSeeoKamhu9973vk5+fTv39/XnzxRZYuXcodd9zB4sWLMTMmTZrE4sWLueyyyzqta926dVx55ZX86le/4owzzkj8PJrJy8sjLy8PPWpTJDco+LPQmjVrqKysbLOsX79+VFRU0NTUxGuvvcbq1asZMGAAb731FmVlZRQVFbXZfuHChdx555089dRTlJWVMWfOHJ555hmGDBnCzi5O7FVVVfHSSy91GfzTpk3j17/+NWefffbBZeeffz5XXHEF8+bNY+7cuSn2XETCoOA/gvZH5rt37438C1znnXceAwYMAILz+4MGtb0z6vPPP8/y5ctZsmQJ/fr1A+Css87iqquu4rLLLmP69Omdfu7gwYN57733umz33HPP5YEHHuD8888nP/FnT0FBAcnc8kJEotcS/Jn4Mr3O8R+l0aNHU1tb22bZrl272Lp1K0VFRW3umV1SUtLhy1EjR45k9+7dvPnmmweX3X///dx2221s2bKFyspKduzY0aHdhoYGSkpKaG5uZty4cYwbN67NBd97770XgBtuuCEt/RSRaDQ0BEf7mfiKjYL/KE2ZMoU9e/Ywb948IDiPfvPNNzNr1ixKSkrabHviiSeyefPmNstOOOEEHnvsMWbMmHFwpM3GjRuZMGECc+bMYdCgQWzZsqVDu2+++SZjxowhPz+flStXsnLlSubMmXNwfV5eHr/5zW94/fXXDzsCSESyW0vwZ4KC/yiZGQsXLmT+/PmMGjWKgQMHkpeXx+zZszts26dPH0aOHMmGDRvaLD/55JN5+OGH+eIXv8jGjRv51re+xamnnsqYMWOYOHEiY8eO7fBZNTU1TJ069bC1FRcXs2jRIhYtWsR9992XWkdFJBKNjZkLftw966fKykpvb+3atR2Wtdi1a1eX6zJl2bJlPmzYMK+tre10/YIFC3z27NkptbF161Y/55xzulyfar8P9zPNZjU1NVGXEAn1u2ebMcO9oqLtsmT7Diz3TjK1R1zczQYTJ07k7bc7fcoZAJdcckmn5+yT8c4772gcvkhMZPJUT04Hv7vn1M3Frr322pT2P/3009NUSUeu8f0iWaWhAdqNAE+bnD3HX1xczI4dOxRYaeCJ+/HH8eE1ItlKR/ydGDp0KHV1dXz44Ycd1jU0NMQyxFLpd8sTuEQkOyj4O1FYWNjl06KWLl3K+PHjQ64oenHtt0hP1NAAZWWZ+eycPdUjItKTZXI4p4JfRCQL6QtcIiIxo+AXEYkZDecUEYkZHfGLiMSMgl9EJEbcNapHRCRW9u0LXhX8IiIxkcnHLoKCX0Qk67QEv0b1iIjEhI74RURiRsEvIhIzCn4RkZhpbAxeFfwiIjGhI34RkZjRqB4RkZjpkUf8ZvZ1M1tjZn82s9+aWfyekygi0oUeF/xmNgT4GlDl7mOAfOBLYdchIpKtelzwJxQAJWZWAPQG3ouoDhGRrJPpUT3m7pn55MM1anYj8C/AXmCJu3+lk21mAjMBysvLK6urq7v9+fX19ZSWlqap2tyhfseL+t1zLVgwhJ/8ZBSPP/4H+vdvOrg82b5Pnjy51t2rOqxw91AnoAx4HhgEFAKPA1ccbp/KykpPRk1NTVLb9xTqd7yo3z3X7be7g3t9fdvlyfYdWO6dZGoUp3rOBd5y9w/dfT+wAJgYQR0iIlmpJw7nfAc408x6m5kBU4B1EdQhIpKVGhogPx8KCjLz+aEHv7u/AswHXgP+lKhhbth1iIhkq0w+dhGC0TWhc/d/Av4pirZFRLJdpoNf39wVEckymXzeLij4RUSyjo74RURipqEhcyN6QMEvIpJ1dMQvIhIzCn4RkZhR8IuIxIyCX0QkZjScU0QkZnTELyISMxrOKSISMzriFxGJGQW/iEjMKPhFRGKkuRmamhT8IiKxkekHrYOCX0Qkq2T6sYug4BcRySotwa8jfhGRmFDwi4jEjIJfRCRmFPwiIjGjUT0iIjGjI34RkZjRcE4RkZjREb+ISMwo+EVEYkbBLyISMwp+EZGY0XBOEZGY0ageEZGYUfCLiMRMQwP06gV5GUznSILfzI4xs/lm9rqZrTOz/xZFHSIi2SbTj10EKMjsx3fpbuBpd7/UzHoBvSOqQ0Qkq/TI4Dez/sBfAVcBuPs+YF/YdYiIZKPGxswHv7l7Zlto36DZOGAusBYYC9QCN7r7p+22mwnMBCgvL6+srq7udhv19fWUlpamreZcoX7Hi/rdM91661+wfn1f5s37Y4d1yfZ98uTJte5e1WGFu4c6AVVAEzAh8f5u4NbD7VNZWenJqKmpSWr7nkL9jhf1u2e6+GL3U0/tfF2yfQeWeyeZGsXF3Tqgzt1fSbyfD3wugjpERLJOGOf4Qw9+d98KbDGzkxKLphCc9hERib0eeXE34X8DDydG9GwCro6oDhGRrNLQAH37ZraNSILf3VcSnOsXEZFWGhpg0KDMtqFv7oqIZJEwhnMq+EVEskhDQ2bv0wMKfhGRrNIjR/WIiEjXFPwiIjGj4BcRiRF3Bb+ISKw0NQXhr+AXEYmJMJ6+BQp+EZGs0RL8OuIXEYkJBb+ISMwo+EVEYkbBLyISMwp+EZGYaWwMXhX8IiIxoeGcIiIxo1M9IiIxo+AXEYmZyIPfzH6VeL0xsyWIiAhkQfADlWZ2HPC3ZlZmZgNaT5ktS0QkfvbsCV4zHfyHe9j6/cBzwAigFrBW6zyxXERE0uTDDyEvDwZk+NC6yyN+d7/H3f8CeNDdR7j78FaTQl9EJM22bYNBgyA/P7PtHPHirrtfn9kSREQEYOtWKC/PfDsa1SMikiW2bVPwi4jEyrZt8JnPZL4dBb+ISBZw16keEZFY2bUruEmbjvhFRGJi69bgVUf8IiIxsW1b8KrgFxGJiZbg16keEZGYiMWpHjPLN7MVZrY4qhpERLLFtm3BN3YHDsx8W1Ee8d8IrIuwfRGRrBHW7RogouA3s6HAVOCBKNoXEck2YY3hBzB3D6el1o2azQd+APQFvunuF3ayzUxgJkB5eXlldXV1tz+/vr6e0tLSNFWbO9TveFG/e5brr/8cffs2cfvtq7vcJtm+T548udbdqzqscPdQJ+BC4KeJ+UnA4iPtU1lZ6cmoqalJavueQv2OF/W7Zxk2zH3GjMNvk2zfgeXeSaZGcarnLOAiM9sMVAPnmNmvI6hDRCQrtNyuIYyhnBDBOX53/7a7D3X3CuBLwPPufkXYdYiIZItPPoF9+8I7x69x/CIiEQtzDD8c/tGLGefuS4GlUdYgIhK1ML+1CzriFxGJXJj36QEFv4hI5FpO9eiIX0QkJlpu1zBgQDjtKfhFRCK2bRsMHgx5ISWygl9EJGJhjuEHBb+ISOS2bQvvwi4o+EVEIrdtm474RURiw11H/CIisbJzZ7i3awAFv4hIpMIeww8KfhGRSIX9rV1Q8IuIRCrs+/SAgl9EJFJh35kTFPwiIpHavBlKSqCsLLw2FfwiIhFatQpOPTW82zWAgl9EJDLuQfCPHRtuuwp+EZGIvPcefPSRgl9EJDZWrQpeTzst3HYV/CIiEVHwi4jEzKpVUFEB/fuH266CX0QkIqtXh3+0Dwp+EZFI7N0Lb7wR/oVdUPCLiERizRo4cEDBLyISG6tXB6861SMiEhOrVkGfPjByZPhtK/hFRCIQxa0aWij4RURC1nKrhihO84CCX0QkdHV1wSMXo7iwCwp+EZHQtXxjV8EvIhITLcF/6qnRtK/gFxEJWU0NjB4N/fpF076CX0QkRDt3wgsvwF//dXQ1hB78Zna8mdWY2VozW2NmN4Zdg4hIVJ5+Gpqa4KKLoquhIII2m4Cb3f01M+sL1JrZs+6+NoJaRERCtWgRDBoEEyZEV0PoR/zu/r67v5aY3w2sA4aEXYeISNj274ff/x6mToX8/OjqMHePrnGzCuBFYIy772q3biYwE6C8vLyyurq6259bX19PaWlp+grNEep3vKjfuWfFimP4xjfGMWfOn/nLv9ye9P7J9n3y5Mm17l7VYYW7RzIBpUAtMP1I21ZWVnoyampqktq+p1C/40X9zj033eReVOS+e/fR7Z9s34Hl3kmmRjKqx8wKgceAh919QRQ1iIiEyT04vz9lCkT9B0sUo3oM+Dmwzt3vDLt9EZEorFsHmzZFO4yzRRRH/GcBfwOcY2YrE9MXIqhDRCQ0ixYFrxdeGG0dEMFwTnf/A2BhtysiEpXmZvj5z+HMM2Ho0Kir0Td3RUQybtEi2LABvv71qCsJKPhFRDLsX/8Vhg+H6dOjriQQxTd3RURiY9kyePll+MlPoCBLEldH/CIiGXTHHTBgAFx9ddSVHKLgFxHJkDffhN/9Dm64IXiwerZQ8IuIZMiPfgS9esGsWVFX0paCX0QkA5YtgwcfDI72y8ujrqYtBb+ISJo1NsJ118GwYTBnTtTVdJQl15hFRHqOH/wguEXDU09Ff1+ezuiIX0Qkjdauhe9/Hy6/HD7/+air6ZyCX0QkTfbuhSuvhL594a67oq6mazrVIyKSBgcOBKFfWwsLF8LgwVFX1DUFv4hIGvzjP8KjjwZf2Jo2LepqDk+nekREUvSLXwTn9WfOhG98I+pqjkzBLyKSgn//d7jmGjjvPLj3XrAcuOm8gl9E5Ci4B2P0/+7v4AtfgMcfh8LCqKvqHp3jFxFJ0r59cNNN8LOfBRd0/+M/cif0QcEvIpKU9evhy18ORu/ccgv88Ie5cXqnNQW/iEg3uMO8ecEN1woLYcECuOSSqKs6OjrHLyJyBGvWwJQpcNVVMH48rFqVu6EPCn4RkS59+GHwnNyxY2HlSvjpT6GmBo4/PurKUqNTPSIi7WzbFjwn92c/C27DcM01wY3Xjj026srSQ8EvIpJQWwv33Qe//W0wcufyy2H2bDj55KgrSy8Fv4jE2o4dwa0WHnoIXnkFeveGGTPg5pvhxBOjri4zFPwiEjvbt8OTTwYjc37/e9i/H0aPhrvvDsbl9+8fdYWZpeAXkR6vuRlWrIBnn4Wnn4Y//CG4m+aQIXDjjfCVrwQXcHNtPP7RUvCLSI+zf38wCuell4LpxRfho4+CdWPHwne+E9xBs7IyPmHfmoJfRHLavn3w+utB0NfWwquvBkf3DQ3B+hEj4KKL4NxzgynbHnweBQW/iOSExkbYsAFeeulYli0LvlS1Zk0Q+vv2Bdv07h0cxX/1qzBhApx1Fhx3XLR1ZyMFv4hkhQMH4IMP4O23g2nzZti0KQj7jRvhnXeCbWAMACecAKecAhdcAOPGBdOoUVCgVDsi/YhEJKPcYefOINTffx+2bg1e33sP3n03mLZsCV5bjtxbHHssjBwJEycGQyxPOgl2717O5ZdX0bdvNP3pCSIJfjO7ALgbyAcecPcfRlGHiCRn//4gxHfuhI8/Di6Ytrxu3x6Mid+xI7jVQcv0wQfBfu0VFQWjaoYMgTPPDG6DcPzxwZF8RUXw2q9fx/2WLq1X6Kco9OA3s3zgPuA8oA541cwWufvasGsR6encg1sOfPpp26m+/tC0e3fHadeuYPrkk0PTzp3BvodzzDEwcCAMGhSE+PjxwcXUwYOD6bOfhc98JpjKyuI5oiYbRHHEfwawwd03AZhZNTANUPBLj+IejB/ft6/j1NjYdupsWUPDodfW0969wevbb4/mzjuD9y3Tnj3B1BL2e/YkV3NpaXCU3bdv8CWmfv2CsD7mmGDq3z8I7LKy4P2AAcFUVhYEvs6v5wZz93AbNLsUuMDdr028/xtggrvParfdTGAmQHl5eWV1dXW326ivr6e0tDR9ReeIVPrt3nqyg68HDgTLIFgWvO+4/tC8JS7AtX3f+rW5ue365uZD6w8csMQUzLesbz0fbHto+d69jRQWlhx83/ozmpvt4BTs0/q17frWU1NTXqv59uuC9S3LW94fmg/e79+f3pvfFhYeoFevAxQWHqCo6AAFBU2UlDhFRcHyoqIDFBU1U1zc8bW4uPngVFJyaOrd+9B8cXEz+flpLTkj4vrfNyTf98mTJ9e6e1X75Vn7/2d3nwvMBaiqqvJJkyZ1e9+lS5cyadIkrr8eXnihO20defnRzreEaWfL269PZr6zqampiby8gsNu0zIdCuwj/3x6kvz84Kg0Pz+YCgsPvW+ZLyjouLy4+ND8kaZevYKpsDA4j90y36tX2/dFRYfet8y3n0pKgrZ79YL8/Dxa30m95d953MS135C+vkcR/O8Cre9mPTSxLO2GDYMxY7q3bVfnGlsvP9p5s66Xt1+fzHz76d13t3L88UO7XN8y5eUdeX3LNu23zc/ver71Pi3zrZcfbr7ls1rm8/IOP9962auv/hdnnXXmwfctU+uALyg41J5I3EUR/K8Co8xsOEHgfwm4PBMNffvbmfjU7LV06QYmTRoadRmh27KlIecfjCESptCD392bzGwW8AzBcM4H3X1N2HWIiMRVJOf43f0p4Kko2hYRiTs9c1dEJGYU/CIiMaPgFxGJGQW/iEjMKPhFRGJGwS8iEjOh36vnaJjZh8DbSexyLLA9Q+VkM/U7XtTv+Em27ye4+6D2C3Mi+JNlZss7uzFRT6d+x4v6HT/p6rtO9YiIxIyCX0QkZnpq8M+NuoCIqN/xon7HT1r63iPP8YuISNd66hG/iIh0QcEvIhIzORf8ZnaBmb1hZhvM7B86WV9kZo8k1r9iZhWt1n07sfwNMzs/zLpTdbT9NrPzzKzWzP6UeD0n7NpTkcrvO7F+mJnVm9k3w6o5HVL8d36amb1sZmsSv/fiMGtPRQr/zgvN7JeJ/q4zs5x6DFM3+v1XZvaamTUlnlveet2VZrY+MV3ZrQbdPWcmgge3bARGAL2AVcDodtvcANyfmP8S8EhifnRi+yJgeOJz8qPuUwj9Hg8cl5gfA7wbdX/C6Her9fOBR4FvRt2fkH7fBcBqYGzi/cCY/Du/HKhOzPcGNgMVUfcpjf2uAE4D5gGXtlo+ANiUeC1LzJcdqc1cO+I/A9jg7pvcfR9QDUxrt8004JeJ+fnAFDOzxPJqd29097eADYnPywVH3W93X+Hu7yWWrwFKzKwolKpTl8rvGzO7GHiLoN+5JJV+/w9gtbuvAnD3He7eHFLdqUql3w70MbMCoATYB+wKp+yUHbHf7r7Z3VcDB9rtez7wrLt/5O4fA88CFxypwVwL/iHAllbv6xLLOt3G3ZuATwiOerqzb7ZKpd+t/U/gNXdvzFCd6XbU/TazUuDvge+FUGe6pfL7PhFwM3smcWrglhDqTZdU+j0f+BR4H3gHuMPdP8p0wWmSSjYd1b6RPHpRwmdmpwA/IjgijIN/Bu5y9/rEHwBxUQCcDZwO7AGeM7Nad38u2rIy7gygGTiO4JTHS2b2n+6+KdqyslOuHfG/Cxzf6v3QxLJOt0n82dcf2NHNfbNVKv3GzIYCC4EZ7r4x49WmTyr9ngDcbmabgZuA75jZrEwXnCap9LsOeNHdt7v7HoJnW38u4xWnRyr9vhx42t33u/sHwDIgV+7nk0o2Hd2+UV/YSPIiSAHBxYvhHLoIckq7bb5K24s//zcxfwptL+5uIncueqXS72MS20+Puh9h9rvdNv9Mbl3cTeX3XQa8RnCBswD4T2Bq1H0Kod9/D/wiMd8HWAucFnWf0tXvVts+RMeLu28lfu9lifkBR2wz6k4fxQ/pC8CbBFfBZyeWzQEuSswXE4zi2AD8ERjRat/Zif3eAD4fdV/C6DfwfwjOfa5sNQ2Ouj9h/L5bfUZOBX+q/QauILig/Wfg9qj7Eka/gdLE8jWJ0P9W1H1Jc79PJ/hr7lOCv3DWtNr3bxM/jw3A1d1pT7dsEBGJmVw7xy8iIilS8IuIxIyCX0QkZhT8IiIxo+AXEYkZBb9IN5hZhZn9uZvb/r9M1yOSCgW/SJq5+8SoaxA5HAW/SJLMbISZrTCz/25mfzSzlWa22sxGJdbXR12jyOHoJm0iSTCzkwhum3sVcC1wt7s/bGa9CO6rLpL1FPwi3TcI+B3BfY/WmtnLwOzETfAWuPv6aMsT6R6d6hHpvk8I7vV+NoC7/wa4CNgLPJVrj7WU+NIRv0j37QMuAZ5JnMf/L2CTu99jZsMIHo33fJQFinSHgl8kCe7+qZldSPCIuyeA/2Vm+4GtwPcjLU6km3R3ThGRmNE5fhGRmFHwi4jEjIJfRCRmFPwiIjGj4BcRiRkFv4hIzCj4RURi5v8D9dB3ad4xiaMAAAAASUVORK5CYII=\n","text/plain":["<Figure size 432x288 with 1 Axes>"]},"metadata":{"needs_background":"light"}},{"output_type":"stream","name":"stdout","text":["ksi_eq =  0.06180339887498969\n"]}]},{"cell_type":"markdown","source":["\n","Déterminer l'avancement volumique à l’équilibre de deux réactions.\n","\n"],"metadata":{"id":"BIN6NKodam5K"}},{"cell_type":"code","source":["import numpy as np #Pour faire divers calculs\n","from scipy.optimize import fsolve #Importe la fonction pour faire la résolution d’équation\n","\n","## Constantes d’équilibre et concentration initiale\n","K1=10**(-4.8) # K° de la réaction 1 HA + H2O = A(-) + H3O(+)\n","K2=10**(-14) # K° de la réaction 2 H2O = HO(-) + H3O(+)\n","C=1e-7 # Concentration initiale de HA en mol/L\n","def f1(x): #définit la première équation x[0]=[A-], [H3O+]= x[0]+x[1]\n","\treturn\tx[0]*(x[0]+x[1])-K1*(C-x[0]) #Qr1-K1\n","def f2(x): #définit la seconde équation x[1]=[HO-], [H3O+]= x[0]+x[1]\n","\treturn\tx[1]*(x[0]+x[1]) - K2 #Qr2-K2\n","def func(x): #définit le système des deux équations\n","\treturn\tf1(x),f2(x) #système\n","root = fsolve(func, [C/2,C/2]) #résolution du système\n","print('les solutions sont :','[A-] = ',format(root[0], \"#.1e\"),'mol/L',' et [HO-] = ',format(root[1],\"#.1e\"), 'mol/L')\n","#affiche la solution\n","pH=-np.log10(root[0]+root[1]) #calcule le pH\n","print('[H3O+] = ',format(root[0]+root[1], \"#.1e\"),'mol/L',' et pH = ',format(pH,\"#.1f\")) #affiche le pH"],"metadata":{"colab":{"base_uri":"https://localhost:8080/"},"id":"A1eIm-gsa5FY","executionInfo":{"status":"ok","timestamp":1655448709770,"user_tz":-120,"elapsed":23,"user":{"displayName":"Lionnel Malara","userId":"16342666118884215780"}},"outputId":"54819f7c-4c39-4197-deea-0a361459217b"},"execution_count":2,"outputs":[{"output_type":"stream","name":"stdout","text":["les solutions sont : [A-] =  9.9e-08 mol/L  et [HO-] =  6.2e-08 mol/L\n","[H3O+] =  1.6e-07 mol/L  et pH =  6.8\n"]}]}]}