S = 412
r = 0.09
M = 4
c_k = 2/M

d_0_1 = 1/(1+r/M)

storage_costs = c_k / d_0_1^3 + c_k / d_0_1^2 + c_k / d_0_1

F = S / d_0_1^3 + storage_costs

print(sprintf('storage_costs= %f; F= %f', storage_costs, F))