https://keroseneheaterreview.blogspot.com/2023/10/vectorising-equation-using-numpy.html