#
# PYTHON • MyCompiler
# Zhukovsky
# https://[Log in to view URL]
# Transformação de Zhukovsky
# 
# 2026 Maio Dom, 17 — Seg, 18
# by 
# j
# 

import matplotlib.pyplot as plt
import numpy as np

def graf():

    plt.style.use('dark_background')

    x = []
    y = []
    w = []
    z = []

    x0 = -0.15
    y0 = +0.35
    r  = 1.28
    a  = 1.00

    for t in range(0,3600):

        trad = t*np.pi/1800

        xx = x0 + r*np.sin(trad)
        yy = y0 + r*np.cos(trad)

        zz = xx + yy*1j

        zkw = zz + (a**2)/zz

        x.append(zkw.real)
        y.append(zkw.imag)

        w.append(xx)
        z.append(yy)

    plt.plot(x,y)
    plt.plot(w,z)

    sx0 = str(x0)
    sy0 = str(y0)
    sr  = str(r)
    sa  = str(a)

    ss  = 'x0 = ' + sx0 + ' y0 = ' + sy0 + ' r = ' + sr + ' a = ' + sa
    tit = 'Zhukovsky | ' + ss

    plt.axis('equal') 
    plt.title(tit)
    # plt.xticks([])  
    # plt.yticks([])
    plt.grid(False)
    plt.show()

graf()