#! /usr/bin/env python
import sys
from typing import Final
import numpy as np
from numpy import pi
import matplotlib.pyplot as plt


amplitude: tuple[float, float] = (1, 1)
# frequent, some values (-3, 5) (-3, 8) (3, 11)
omega: tuple[float, float] = (3, 11)

# phase
phi: tuple[float, float] = (pi/2, 3*pi/4)

N: Final[int] = 360*2
T = np.linspace(0, 2*pi, num=N)
x = amplitude[0] * np.cos(omega[0]*T + phi[0])
y = amplitude[1] * np.cos(omega[1]*T + phi[1])

# plot
fig, ax = plt.subplots()
ax.plot(x, y, linewidth=2.0)

write_to_file = False
if write_to_file:
    output_dir = "../../2024/python-output"
    image_filename = f"{output_dir}/{sys.argv[0].replace('.py','.pdf')}"
    plt.savefig(image_filename)

plt.show()