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)