Please fix the following code: import matplotlib.pyplot as plt

def simulate$\_$population$\_$growth$($children$\_$per$\_$female, risk$\_$factor, years$=500,$ initial$\_$population$=4)$:

childbearing$\_$years $=25$

population $=$ initial$\_$population

females $=$ initial$\_$population $/2$

population$\_$over$\_$time $=[]$

for year in range$(1,$ years $+1)$:

safe$\_$children $=$ children$\_$per$\_$female $*(1-$ risk$\_$factor$)$

if year $\%$ childbearing$\_$years $==0$:

females $+=($females $*$ safe$\_$children$)/2$

population $+=$ females $*$ safe$\_$children

population$\_$over$\_$time.append$($population$)$

return population$\_$over$\_$time

# Using the calculated solutions for each risk scenario

children$\_$per$\_$female$\_25=$ solutions$[0\mathrm{.}25]$

children$\_$per$\_$female$\_50=$ solutions$[0\mathrm{.}5]$

children$\_$per$\_$female$\_75=$ solutions$[0\mathrm{.}75]$

# Simulate population growth

population$\_$growth$\_25=$ simulate$\_$population$\_$growth$($children$\_$per$\_$female$\_25,0\mathrm{.}25)$

population$\_$growth$\_50=$ simulate$\_$population$\_$growth$($children$\_$per$\_$female$\_50,0\mathrm{.}5)$

population$\_$growth$\_75=$ simulate$\_$population$\_$growth$($children$\_$per$\_$female$\_75,0\mathrm{.}75)$

# Plotting the results

plt$.$figure$($figsize$=(12,8))$

plt$.$plot$($population$\_$growth$\_25,$ label$=\text{'}25\%$ Risk Scenario'$)$

plt$.$plot$($population$\_$growth$\_50,$ label$=\text{'}50\%$ Risk Scenario'$)$

plt$.$plot$($population$\_$growth$\_75,$ label$=\text{'}75\%$ Risk Scenario'$)$

plt$.$xlabel$(\text{'}$Martian Years'$)$

plt$.$ylabel$(\text{'}$Population$\text{'})$

plt$.$title$(\text{'}$Projected Population Growth on Mars Over $500$ Martian Years'$)$

plt$.$legend$()$

plt$.$grid$($True$)$

plt$.$show$()$

import math

# Initial and final populations

initial$\_$population $=4$

final$\_$population$\_25=$ population$\_$growth$($children$\_$per$\_$female$\_25,0\mathrm{.}25)$

final$\_$population$\_50=$ population$\_$growth$($children$\_$per$\_$female$\_50,0\mathrm{.}5)$

final$\_$population$\_75=$ population$\_$growth$($children$\_$per$\_$female$\_75,0\mathrm{.}75)$

# Time period in Martian years

time$\_$period $=500$

# Calculating the growth coefficients for each scenario

growth$\_$coefficient$\_25=(1/$ time$\_$period$)*$ math.log$($final$\_$population$\_25/$ initial$\_$population$)$

growth$\_$coefficient$\_50=(1/$ time$\_$period$)*$ math.log$($final$\_$population$\_50/$ initial$\_$population$)$

growth$\_$coefficient$\_75=(1/$ time$\_$period$)*$ math.log$($final$\_$population$\_75/$ initial$\_$population$)$

children$\_$per$\_$female$\_25,$ children$\_$per$\_$female$\_50,$ children$\_$per$\_$female$\_75,$ growth$\_$coefficient$\_25,$ growth$\_$coefficient$\_50,$ growth$\_$coefficient$\_75$