import math
import matplotlib.pyplot as plt

Fs = 1000 # samples per second
f = 100   # frequency

phase = 0

time = []
signal = []

for i in range(1000):
    phase += 2 * math.pi * f / Fs
    sample = math.cos(phase)
    
    time.append(i / Fs)
    signal.append(sample)

plt.plot(time, signal)
plt.xlabel("Time (seconds)")
plt.ylabel("Amplitude")
plt.title("Cosine Wave")
plt.grid()
plt.show()