heart shaped cookie cutter

thingiverse

Here is the code to generate the given plot with all instructions: ``` # This code will draw the above described curve import numpy as np from sympy import sin, cos import matplotlib.pyplot as plt # Number of points on the curves num_points = 15 t = np.linspace(0.00,1.0,num_points) x_ellipses_2d_001 = lambda theta: np.cos(np.radians(4+30*np.tan(np.pi/18+theta/6))) y_ellipses_2d_001 = lambda theta: -np.sin(np.radians(8)) X_ellipses_2d_001, Y_ellipses_2d_001 = [], [] # Points of interest points = [[26.974074740740745,20.4024002402402],[16.174673867335873,-4.8777777777777816],[15.510610061022067,-27.710100110011012]] for i in t: X_ellipses_2d_001.append(x_ellipses_2d_001(i)) Y_ellipses_2d_001.append(y_ellipses_2d_001(i)) plt.figure(figsize=(12,8)) # Ploting points plt.scatter(*zip(*points), color="red", marker = '*') plt.annotate('', xy=points[0], xytext=points[-1], arrowprops=dict(arrowstyle="->",color='black'),zorder=10) # Ploting ellipses x = X_ellipses_2d_001 y = Y_ellipses_2d_001 plt.plot(x,y,color = 'green') # Plotting curves that should mirror m_x = [] m_y = [] def sin_60(n): return n - (7+4*n)/3.464102753114229 for x,y in zip(x[-1:-9:-1], y[-1:-9:-1]): n_x = (y*(np.sin(np.radians(60))/ 2 + np.cos(np.pi/12))) n_y = sin_60(n_x)*2 m_x.append(y-n_y) m_y.append(x-n_x) m_y_ellipses_2d_001, m_X_ellipses_2d_001= [], [] for i in t: n_x = (y_ellipses_2d_001(i)* (np.sin(np.radians(60))/ 2 + np.cos(np.pi/12))) n_y = sin_60(n_x)*2 m_y.append(y_ellipses_2d_001(i)-n_y) m_X_ellipses_2d_001.append(x_ellipses_2d_001(i) - n_x) # Plotting mirror curves plt.plot(m_X_ellipses_2d_001, m_y_ellipses_2d_001, color = 'orange') plt.xlim(-45,25) plt.ylim(-38.21,22.1) for label in ['left','bottom','right','top']: plt.gca().spines[label].set_visible(False) plt.gca().add_patch(plt.Rectangle((20,-37),(7,0), edgecolor='black', facecolor = 'none')) plt.xlabel('X (nm)') plt.ylabel('Y (nm)') plt.show() ```