accidents = [22 22 26 28 27 25 30 29 24];
new_accidents = accidents - 22;
disp(variance(new_accidents), "The variance of the number of accidents is")