You should be careful with degrees and radians. Your whole program works with degrees EXCEPT the numpy.sin(x) function which only accepts radians. So I think you just need to replace 360° with 2π everywhere.