function mean_squared_error(predictions, actuals):
   total = 0

   for prediction, actual in zip(predictions, actuals):
      total += (prediction - actual) ** 2
   return total / length(predictions) # squaring punishes large misses more
