1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
import numpy as np
import matplotlib.pyplot as plt
def mgc_plot(x, y, sim_name):
"""Plot sim and MGC-plot"""
# simulation
plt.figure(figsize=(8, 8))
ax = plt.gca()
ax.set_title(sim_name + " Simulation", fontsize=20)
ax.scatter(x, y)
ax.set_xlabel('X', fontsize=15)
ax.set_ylabel('Y', fontsize=15)
ax.axis('equal')
ax.tick_params(axis="x", labelsize=15)
ax.tick_params(axis="y", labelsize=15)
plt.show()
rng = np.random.default_rng()
unif = np.array(rng.uniform(0, 5, size=100))
x = unif * np.cos(np.pi * unif)
y = unif * np.sin(np.pi * unif) + 0.4 * rng.random(x.size)
mgc_plot(x, y, "Spiral")
|
