f[N-1]= -Tb*((rho*vel)/(2*h) - (k)/(h**2))
print('',f)
# La solucion en el arreglo u, incluye las fronteras
U = np.zeros(N+2)
#obtención de solucion del sistema de N x N
U[1:N+1] = np.linalg.solve(A,f)
#Valores de frontera conocidos debido a las condiciones de Dirichlet
U[0]=Ta
U[N+1]=Tb
print('La solución es:',U)
#Se establece el dominio para solución numérica
XNumerica=np.linspace(0,L,N+2)
#Se establece el dominio para solución analítica
XAnalitica = np.linspace(0,L,100)
#Se calcula la solución analítica
SolucionA = func.analyticSol(XAnalitica,rho,vel,k,L,Ta,Tb)
#Salvamos sol
par['solucion']=U
par['x']=XNumerica
hdf5.saveParametros(out_file_name,par)
#Figuras
func.Graf_Sol_A_N(XNumerica,U,XAnalitica,SolucionA)
Ejecucion correcta: python hdf5.py ENTRADA
donde el nombre "ENTRADA" es el nombre del archivo donde
se almacenan los datos del problema.
Los datos que se guardan en el archivo ENTRADA son los que se encuentran en
este archivo hdf5.py
Por ejemplo: python hdf5.py ENTRADA"""
print(mensaje)
sys.exit(1)
#Generamos parametros
Datos = {
'a': 0,
'b': 2.5,
'N': 300,
'Ta': 1,
'Tb': 0,
'k': 0.01,
'q': 0,
'ht': 0.001,
'Tmax': 1,
'vel': 1,
'Tol': 1E-4,
'rho': 1,
'Nt': 500,
}
hdf5.saveParametros(in_file_name, Datos)