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)
ZnJvbSBudW1weSBpbXBvcnQgemVyb3MsIGFycmF5LCBzaW4sIHBpCmRqbiA9IFsxLCAyLCAzLCA0LCA1LCA2LCA3LCA4LCA5LCAxMF0KZGpuID0gYXJyYXkoZGpuKQoKZGVmIHRlc3QoZGpuKToKICAgTSA9IGxlbihkam4pCiAgIGRlbHRhID0gemVyb3MoTSwgZmxvYXQpCiAgIGRlbHRhWzA6TV0gPSAwLjQwOTI3OTcwOSAqIHNpbigoKDIqcGkpLzM2NSkqKDI4NCtkam5bMDpNXSkpCiAgIHJldHVybiBkZWx0YQoKcmVzID0gdGVzdChkam4pCnByaW50KHJlcyk=