Lista 10 Questão 04

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

# -*- coding: utf-8 -*-
"""
Spyder Editor
Author: Vinícius Hunzicker Ferreira 
"""
import numpy as np
 
def splinelinear(x, f, z):
 i = 0
 assert (z >= x[0]) or (z <= x[-1]), "O valor está fora do intervalo"
 while(x[i+1] < z):
 i = i + 1
 return f[i] + (f[i+1] - f[i])/(x[i+1] - x[i]) * (z - x[i])
 
P, T = np.loadtxt("thermo.dat")
p = splinelinear(T, P, 300.)
 
print "A pressão em que a água entra em ebulição a 300K é", p, "kPa"