https://camera-optical-zoom.blogspot.com/2023/01/vectorising-equation-using-numpy.html