# -*- coding: utf-8 -*-
"""
Created on Tue Dec  6 10:40:12 2022

@author: CK
"""
#1- importation des bibliothèques
import matplotlib.pyplot as plt
import numpy as np
import scipy.optimize as op

#2- données
DrH0 = -150E3 #kJ/mol
CpmE = 2.1E3 #J/K/kg
Te=273+200 #K
T0=273+20 #K
Dv = 3E-3/3600 # m3/s soit 3 L/h
A=1.0E15 #s-1
Ea=157E3 #J/mol
R=8.314 #J/K/mol
V=500E-6 #m3 (volume du réacteur)
tau = V/Dv #s
ro=0.9E3 #kg/m3 masse volumique du liquide
h=80 #W/(m²K) coeff conducto-convectif
S=300E-4 #m²
M=146E-3 #kg/mol
eps=0.001

#3- définitions des fonctions
def k(T):
    


def Xa_BM(T):


    
def Xa_BT(T):
        


def f(T) : #fonction à annuler pour trouver les points de fonctionnement


def deriv(f,x,epsy): #dérivée de la fonction f au point x avec un pas = epsy 


def stabil(Ts):



#4-tracé des courbes pour déterminer visuellement les encadrements des zéros




#4- calcul des points de fonctionnement