make simulation code pyhon of a system in contact with a thermal reservoir at temperature T, the system is composed of N particles that can be found in two states with zero energy or greater than zero energy. calculate using Monte Carlo method with metropole algorithms to calculate average energy and specific

USE MATLAB PLEASE!

MATLAB CODE as if not reiterate..

How do I input the code into MATLAB to solve this problem? Thank you!

Using python or matlab to compute this.

Write a MATLAB function that calculates the approximate value of six(x) using the Maclaurin series approximation from x-0 to xi+1 six(x) = x -+-+. 3 5 The function should accept 3 parameters: value of x, number of significant figures accuracy i.e. n, and the maximum number of iterations. In the function, use ɛ = (0.5 × 10²-" )% in order to continue until the ɛa falls below this criteria. The function should return 3 values: the approximate value of arctan(x) at the end of the program, final ɛg and the number of iterations it took. Submit short snapshots showing the program running.

In MATLAB. Please do not change varaible name. will give like

In matlab code, please

2. calculates the trajectory r(t) and stores the coordinates for time steps At as a nested list trajectory that contains [[xe, ye, ze], [x1, y1, z1], [x2, y2, z2], ...]. Start from time t = 0 and use a time step At = 0.01; the last data point in the trajectory should be the time when the oscillator "hits the ground", i.e., when z(t) ≤ 0; 3. stores the time for hitting the ground (i.e., the first time t when z(t) ≤ 0) in the variable t_contact and the corresponding positions in the variables x_contact, y_contact, and z_contact. Print t_contact = 1.430 X_contact = 0.755 y contact = -0.380 z_contact = (Output floating point numbers with 3 decimals using format (), e.g., "t_contact = {:.3f}" .format(t_contact).) The partial example output above is for ze = 10. 4. calculates the average x- and y-coordinates 1 y = Yi N where the x, y, are the x(t), y(t) in the trajectory and N is the number of data points that you calculated. Store the result as a list in the variable center = [x_avg, y_avg]…