from numpy import zeros, array, sin, pi
djn = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
djn = array(djn)

def test(djn):
   M = len(djn)
   delta = zeros(M, float)
   delta[0:M] = 0.409279709 * sin(((2*pi)/365)*(284+djn[0:M]))
   return delta

res = test(djn)
print(res)