Code to generate a simple plot (below) runs on my jupyter notebook with no issue, and produces the expected figure 1. But another individual running the same code on their computer (after clearing the kernel) gets a plot that looks completely different 2. Is there some user error (in running the code) that would cause this?
We are both using the anaconda distribution for python.
import numpy as np
import matplotlib.pyplot as plt
sigma = 5.67e-8 #stefan boltzman constant in W/(m^2)/(K^4)
T = np.arange(50,800,) #units in K
F = sigma*T**4 #Flux in units of W/m/m
T1 = 289 #Temperature of the earth in K
T2 = 698 #Temperature of venus in K
T3 = 210 #Temperature of mars in K
T4 = 90 #Temperature of titan in K
T_p = np.array([T1, T2, T3, T4]) #array of planetary surface temperatures
F_p = sigma*T_p**4
print(F_p)
plt.plot(T,F) #plot Stefan-Boltzman law for T = 50 to 800K
plt.plot(T_p[0], F_p[0], 'o', label='Earth', color='red')
plt.plot(T_p[1], F_p[1], 'o', label='Venus', color='blue')
plt.plot(T_p[2], F_p[2], 'o', label='Mars', color='green')
plt.plot(T_p[3], F_p[3], 'o', label='Titan', color='purple')
plt.xlabel('Temperature of planetary surface, T (K)')
plt.ylabel('Flux, F (W/m$^{2})$')
plt.legend()
plt.title('Stefan-Boltzman law')