How can this code be improved and made more innovative and commented upon to show every step of the code

import numpy as np

import matplotlib.pyplot as plt

from scipy.signal import bode, TransferFunction

# Step $1$: Define parameters for the ideal design

fc $=2500$# Cutoff frequency in Hz

R $=1000$# Resistance in ohms

L$\_$ideal $=$ R $/(2*$ np$.$pi $*$ fc$)$# Ideal inductance in Henry

# Step $2$: Define the transfer function for the ideal design

num$\_$ideal $=[$L$\_$ideal, $0]$# Numerator: sL $($jL$)$

den$\_$ideal $=[$L$\_$ideal, R$]$# Denominator: sL $+$ R

system$\_$ideal $=$ TransferFunction$($num$\_$ideal, den$\_$ideal$)$

# Step $3$: Adjust component values to match Appendix A

# R $=1$ k$($direct match$)$

# L $=60$H $($series combination of $6$ x $10$H$)+3\mathrm{.}64$H adjustment

L$\_$selected $=63\mathrm{.}64$e$-6$# Selected inductance in Henry

num$\_$selected $=[$L$\_$selected, $0]$# Numerator for selected design

den$\_$selected $=[$L$\_$selected, R$]$# Denominator for selected design

system$\_$selected $=$ TransferFunction$($num$\_$selected, den$\_$selected$)$

# Step $4$: Generate frequency range for the Bode plot

frequencies $=$ np$.$logspace$(2,5,1000)$# Frequency from $100$ Hz to $100$ kHz

w $=2*$ np$.$pi $*$ frequencies $$# Angular frequency

# Step $5$: Calculate Bode response for both designs

$\_,$ mag$\_$ideal, phase$\_$ideal $=$ bode$($system$\_$ideal, w$)$

$\_,$ mag$\_$selected, phase$\_$selected $=$ bode$($system$\_$selected, w$)$

# Step $6$: Plot Bode magnitude and phase for comparison

plt$.$figure$($figsize$=(12,8))$

# Magnitude plot

plt$.$subplot$(2,1,1)$

plt$.$semilogx$($frequencies$,$ mag$\_$ideal, label$=$"Ideal Design", linestyle$="--")$

plt$.$semilogx$($frequencies$,$ mag$\_$selected, label$=$"Selected Design", linestyle$="-")$

plt$.$title$("$Bode Plot for RL High$-$Pass Filter $($Comparison$)")$

plt$.$ylabel$("$Magnitude $($dB$)")$

plt$.$legend$()$

plt$.$grid$($True$,$ which$=$"both", linestyle$="--",$ linewidth$=0\mathrm{.}7)$

# Phase plot

plt$.$subplot$(2,1,2)$

plt$.$semilogx$($frequencies$,$ phase$\_$ideal, label$=$"Ideal Design", linestyle$="--")$

plt$.$semilogx$($frequencies$,$ phase$\_$selected, label$=$"Selected Design", linestyle$="-")$

plt$.$xlabel$("$Frequency $($Hz$)")$

plt$.$ylabel$("$Phase $($degrees$)")$

plt$.$legend$()$

plt$.$grid$($True$,$ which$=$"both", linestyle$="--",$ linewidth$=0\mathrm{.}7)$

plt$.$tight$\_$layout$()$

plt$.$show$()$

# Step $7$: Display calculated component values

print$("$Ideal Inductance $($L$\_$ideal$)$:$",$ L$\_$ideal, $"$H$")$

print$("$Selected Inductance $($L$\_$selected$)$:$",$ L$\_$selected, $"$H$")$