new_contents: # Write your mean_fits function here:
  from astropy.io import fits
  import numpy as np
  
  def mean_fits(fitsfiles):
    count = 0
    
    for fitsfile in fitsfiles:
      hdulist = fits.open(fitsfile)
      newdata = hdulist[0].data
      
      if count == 0:
        data = newdata
      else:
       data = data + newdata #adds matrixes
      count = count + 1
  
    data = data/count #calculate the average for each value
  
    return data
  
  if __name__ == '__main__':
      # Test your function with examples from the question
    data  = mean_fits(['image0.fits', 'image1.fits', 'image2.fits'])
    print(data[100, 100])
  
    # You can also plot the result:
    #import matplotlib.pyplot as plt
    #plt.imshow(data.T, cmap=plt.cm.viridis)
    #plt.colorbar()
    #plt.show()
  
subject: Read a set of FITS files and calculate the mean (average). Based on code from two previous exercizes.