readme.md

with p = np.linspace(7/8*1j*2*np.pi,1j*2*np.pi,10)

with p = np.linspace(3/8*1j*2*np.pi,5/8*1j*2*np.pi,10)

reference: Book: Problems and Solutions in Introductory and Advanced Matrix Calculus Problem 35

ellipsis_problem.py

import numpy as np
import matplotlib.pyplot as plt


p = np.linspace(0,1j*2*np.pi,50)
θ = np.linspace(0,np.pi/2,50)
P,Θ = np.meshgrid(p,θ)
m = np.exp(P)*np.sin(Θ)*np.cos(Θ) + np.sin(Θ)**2
b = m.flatten()
r,i = np.real(b), np.imag(b)
plt.scatter(r,i)
plt.axis('equal')
plt.show()