# -*- coding: utf-8 -*-
"""
Created on Mon May  2 10:33:46 2022

@author: neil
"""
# plot of the 1/time behavior of our transient CA
# (based on program )
import matplotlib.pyplot as plt
times = [20,30,40,50,60,80,100,120,140] # 'scale' times
peaks = [100,70,57,45,38,29,23,19,16]   # these are averages of several CA runs
theory= [115,72,57,46,38,29,23,19,17]   # using integers

fig4 = plt.figure()
ax4  = fig4.add_subplot(111)
ax4.plot(times,theory,"r-o")
ax4.plot(times,peaks,"b+")
ax4.set_ylabel("Temperature, red-theory, blue-CA")
ax4.set_xlabel("Time")
ax4.set_title("Temperature decay of a 2D point source: 1/time")