Diferencia entre revisiones de «Ecuación del calor CCP»

import numpy as np
import matplotlib.pyplot as plt

# Definimos los tramos
x1 = np.linspace(0, 1/3, 300)
x2 = np.linspace(1/3, 2/3, 300)
x3 = np.linspace(2/3, 1, 300)

# Figura
plt.figure(figsize=(9,5))

# Tramos de la función
plt.plot(x1, 10*np.ones_like(x1), color='red', linewidth=3)
plt.plot(x2, np.zeros_like(x2), color='blue', linewidth=3)
plt.plot(x3, 10*np.ones_like(x3), color='red', linewidth=3)

# Líneas de discontinuidad
plt.axvline(1/3, color='skyblue', linestyle='--', linewidth=1)
plt.axvline(2/3, color='skyblue', linestyle='--', linewidth=1)

# Texto
plt.text(0.15, 10.4, 'Zona caliente', color='red', fontsize=10)
plt.text(0.72, 10.4, 'Zona caliente', color='red', fontsize=10)
plt.text(0.45, 0.6, 'Zona fría', color='blue', fontsize=10)

# Ejes y título
plt.xlabel('$x$')
plt.ylabel('Temperatura')
plt.title('Dato inicial: distribución de temperatura')

# Límites
plt.xlim(0,1)
plt.ylim(-1,11)

# Cuadrícula
plt.grid(alpha=0.3)

plt.tight_layout()
plt.savefig('Dato inicial.png', dpi=300)
plt.show()