Please show some work.

Consider a single atom, which has the potential energy E(x)=kcx2, where x is atoms coordinate and kc=0.1 is a constant. Develop a Monte Carlo (MC) program by Python code, which computes the average energy of this atom at the temperatures T from 0.1 to 1.0 with the increment of 0.1. For simplicity assume that kB=1. Set maximum position displacement during MC move to dx=10. Also, set at each temperature the initial coordinate value x to 100. If MC code works correctly, then T/2. The deviation of from the theoretical result T/2 should not exceed 10%. Present the results as a plot vs T and as a table. Include in the table the number of MC steps performed to obtain at each T. I also need to see Python code.