ley de enfriamiento de newtom

# Crear un vector con los tiempos
tiempo <- seq(10, 200, by = 10)

# Crear un vector con las tasas de cambio de temperatura
temperatura_cambio <- c(0.111, 0.111, 0.056, 0.055, 0.056, 0.056, 0.111, 
                        0.055, 0.111, 0.056, 0.000, 0.055, 0.056, 0.056, 
                        0.055, 0.056, 0.055, 0.056, 0.000, 0.000)

# Crear una variable categórica para clasificar si la tasa de cambio es mayor que cero o no
clasificacion <- ifelse(temperatura_cambio > 0, "Positivo", "Cero")

# Crear el gráfico de puntos similar al ejemplo
plot(tiempo, temperatura_cambio, 
     type = "p", # Tipo de gráfico: p para puntos
     main = "Tasa de Cambio de Temperatura con el Tiempo", 
     xlab = "Tiempo (segundos)", 
     ylab = "Cambio de Temperatura (°C/segundo)", 
     pch = ifelse(clasificacion == "Positivo", 15, 0), # Cambiar la forma del punto según la clasificación
     col = ifelse(clasificacion == "Positivo", "green", "black")) # Color de los puntos

# Añadir una leyenda
legend("topright", 
       legend = c("Positivo", "Cero"), 
       col = c("green", "black"), 
       pch = c(15, 0),
       title = "Tasa de Cambio")