# Set solver options
m.options.IMODE = 3  # Steady-state optimization
m.options.SOLVER = 1 # APOPT Solver

# Solve the optimization problem
m.solve(disp=True)

# Print optimized maintenance interval
opt = interval.value[0]
print(f"Optimized maintenance interval: {opt:.2f} days")