I am trying to apply a model based on this formula: and so I defined it as a function
def Amp(angular_freq,N,w0,y,c): return ((N)/(((((angular_freq**