Lucas Ferreira Loz

import matplotlib.pyplot as plt import numpy as np import sympy as sym x=sym.symbols('x') #curva de equilibrio a12=2 x=np.linspace(0,1,40) y=(a12*x)/(1+x*(a12-1)) #reta de retificação x1=np.linspace(0.5,0.98,40) b=0.27528 a=0.71910 y1=a*x1+b #reta de alimentação xa=np.ones(40)*0.5 ya=np.linspace(0.5,0.6348,40) #reta de esgotamento xs=np.linspace(xb,xf,40) ys=1.28*xs-0.0056 #pratos teoricos yp=np.ones(40) xp=np.ones(40) i=1 yp[0]=0.9799 xp[0]=0.98 prat=0 while i<22: xp=(yp[i-1])/(a12+yp[i-1]*(1-a12)) yp=yp[i-1] xp[i+1]=xp yp[i+1]=a*xp[i+1]+b i=i+2 i=24 while i<40: xp=(yp[i-1])/(a12+yp[i-1]*(1-a12)) yp=yp[i-1] xp[i+1]=xp yp[i+1]=1.28*xp[i+1]-0.0056 i=i+2 xp usem esse código. O problema é que os valores no segundo while estão dando negativos

Bom, estou tentando criar um programa que crie um vetor e insira valores nesse vetor usando uma função, até determinado ponto depois use outra função import matplotlib.pyplot as plt import numpy as np import sympy as sym x=sym.symbols('x') #curva de equilibrio a12=2 x=np.linspace(0,1,40) y=(a12*x)/(1+x*(a12-1)) #reta de retificação x1=np.linspace(0.5,0.98,40) b=0.27528 a=0.71910 y1=a*x1+b #reta de alimentação xa=np.ones(40)*0.5 ya=np.linspace(0.5,0.6348,40) #reta de esgotamento xs=np.linspace(xb,xf,40) ys=1.28*xs-0.0056 #pratos teoricos yp=np.ones(40)*0.6 xp=np.ones(40)*0.6 i=1 yp[0]=0.9799 xp[0]=0.98 prat=0 while i<25: if xp>0.5: xp=(yp[i-1])/(a12+yp[i-1]*(1-a12)) yp=yp[i-1] xp[i+1]=xp yp[i+1]=a*xp[i+1]+b i=i+2 else: xp=(yp[i-1])/(a12+yp[i-1]*(1-a12)) yp=yp[i-1] xp[i+1]=xp yp[i+1]=1.28*xp[i+1]-0.0056 i=i+2 eu gostaria que em xp<0.5 o meu programa use outra função

isso foi um erro meu. Depois eu mudei a função e esqueci de mudar na legenda

para entradas use 80 e 100 e o valor esperado é 365

Estou fazendo um código para calcular o ponto de bolha de uma mistura binaria, usando o método de bisseção. Estou com alguns problemas, esse código demora pra pedir as entradas e não retorna nenhuma resposta. Alguém poderia me dizer onde estou errando? #calculo da temperatura do ponto de bolha #função objetivo: F=0.55*(exp(14.5463-(2940.46/T-237.22)))+0.45*(exp(14.2724-(2945.47/T-224.0)))-65 import math a=int(input('a:')) b=int(input('b:')) erro=10 while abs(erro)>0.1: c=(a+b)/2 y1=(760 - 0.5*(math.exp(-3848.0840/a + 17.5318) + math.exp(-4328.1150/a + 17.9134))) y2=(760 - 0.5*(math.exp(-3848.0840/c + 17.5318) + math.exp(-4328.1150/c + 17.9134))) if y2==0: print(c) elif y1*y2<0: b=c else: a=c erro=y2 print('O valor da raiz é: {:.2f}'.format(c))