Lista 03 Questão 6

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

# -*- coding: utf-8 -*-
"""
Created on Sun Oct 08 20:36:55 2017
@author: Vinícius
"""
from raiz import bissecao
from numpy import exp
def PU(t):
 pumax = 80000.
 ku = 0.05
 pumin = 110000.
 return pumax*exp(-ku*t) + pumin
 
def PS(t):
 psmax = 320000.
 p0 = 10000.
 ks = 0.09
 return psmax / (1 + (psmax/(p0 - 1)*exp(-ks*t)))
def g(t):
 return 1.2*PU(t) - PS(t)
tempo = bissecao(g, 1, 100, 1E-2)
print tempo