CÓDIGOS PYTHON
Presión, clase 11/08/20.
Primero, link para R en diferentes unidades: http://tecno.cruzfierro.com/formularios/R.pdf
Código:
# -*- coding: utf-8 -*-
"""
Spyder Editor
This is a temporary script file.
"""
import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt
# function that returns dP/dt
def model(p,t):
p1=1
v=560
P=0.04
R=0.082
T=313
pm=28
dpdt = P*R*T/(pm*v)*(0.16*0.2-0.00506*0.5*(p*(p-p1))**(0.5))
return dpdt
# initial condition
p0 = 2.72
# time points
t = np.linspace(0,100000000)
# solve ODE
p = odeint(model,p0,t)
# plot results
plt.plot(t,p)
plt.xlabel('time')
plt.ylabel('h(t)')
plt.show()
print(t,p)
Problema 2 (apmonitor) aplicado al ejercicio de P, clase 11/08/20.
import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt
# function that returns dP/dt
def model(p,t):
if t<20000000:
p1=1
