help resolve the error in my code :import numpy as np

# Given reaction mechanism $($coefficients are already corrected$)$

finput $=["$C$3$H CH$3*+$ C$2$H$5*","$C$3$H$8+$ H H$2+$ C$3$H$7*","$C$3$H C$3$H$6+$ H$*",$

$"$C$2$H C$2$H$4+$ H$*","$C$3$H C$2$H$4+$ CH$3*","$C$2$H$5*+$ C$3$H C$2$H$6+$ C$3$H$7*",$

$"$CH$3+$ C$3$H CH$4+$ C$3$H$7*","$CH$3*+$ C$3$H C$3$H$6+$ CH$4"]$

reactions $=[]$

for line in finput:

stripped$\_$line $=$ line.replace$(\text{'}\text{'},\text{'}\text{'})$ #

if in stripped$\_$line:

left, right $=$ stripped$\_$line.split

elif in stripped$\_$line:

left, right $=$ stripped$\_$line.split

else:

raise ValueError$($f$"$Invalid delimiter in reaction: $\{$line$\}")$

reactions.append$(($left$.$strip$(),$ right.strip$()))$

species$\_$tmp $=[]$

for left, right in reactions:

left$\_$terms $=$ left.split$(\text{'}+\text{'})$

right$\_$terms $=$ right.split$(\text{'}+\text{'})$

terms $=[$t$.$strip$()$ for t in left$\_$terms$]+[$t$.$strip$()$ for t in right$\_$terms$]$

for term in terms:

tmp $=$ term.split$()$

assert len$($tmp$)==1$ or len$($tmp$)==2,$ 'terms $=\%$r$,$ term $=\%$r$,$ tmp $=\%$r$\text{'}\%($terms$,$ term, tmp$)$

if len$($tmp$)==2$:

species$\_$tmp$.$append$($tmp$[1].$strip$())$

else:

species$\_$tmp$.$append$($term$.$strip$())$

species$\_$filtered $=$ set$($species$\_$tmp$)$

species $=$ list$($species$\_$filtered$)$

# Initialize the stoichiometric matrix as zero

s$\_$mtrx $=$ np$.$zeros$(($len$($reactions$),$ len$($species$)))$

for $($i$\_$row, r$)$ in enumerate$($reactions$)$:

left $=$ r$.$splitstrip$()$

right $=$ r$.$splitstrip$()$

left$\_$terms $=$ left.split$(\text{'}+\text{'})$

left$\_$terms $=[$t$.$strip$()$ for t in left$\_$terms$]$ # in$-$place clean up

right$\_$terms $=$ right.split$(\text{'}+\text{'})$

right$\_$terms $=[$t$.$strip$()$ for t in right$\_$terms$]$ # in$-$place clean up

for t in left$\_$terms: # reactants

tmp $=$ t$.$split$(\text{'}\text{'})$ # split stoichiometric coeff and species name

if len$($tmp$)==2$: # stoich coeff and species name

coeff $=$ float$($tmp$[0].$strip$())$

species$\_$member $=$ tmp$[1].$strip$()$

j$\_$col $=$ species.index$($species$\_$member$)$ # find id of species in the species list

assert s$\_$mtrx$[$i$\_$row,j$\_$col$]==0\mathrm{.}0,\backslash$

'duplicates not allowed r$\%$r: $\%$r $\%$r r$\text{'}\%\backslash$

$($i$\_$row,r$,$species$\_$member,s$\_$mtrx$[$i$\_$row,j$\_$col$])$

s$\_$mtrx$[$i$\_$row,j$\_$col$]=-1\mathrm{.}0*$ coeff

else: # only species name

species$\_$member $=$ tmp$[0].$strip$()$

j$\_$col $=$ species.index$($species$\_$member$)$

assert s$\_$mtrx$[$i$\_$row,j$\_$col$]==0\mathrm{.}0,\backslash$

'duplicates not allowed r$\%$r: $\%$r $\%$r r$\text{'}\%\backslash$

$($i$\_$row,r$,$species$\_$member,s$\_$mtrx$[$i$\_$row,j$\_$col$])$

s$\_$mtrx$[$i$\_$row,j$\_$col$]=-1\mathrm{.}0$

for t in right$\_$terms: # products

tmp $=$ t$.$split$(\text{'}\text{'})$

if len$($tmp$)==2$:

coeff $=$ float$($tmp$[0].$strip$())$

species$\_$member $=$ tmp$[1].$strip$()$

j$\_$col $=$ species.index$($species$\_$member$)$

assert s$\_$mtrx$[$i$\_$row,j$\_$col$]==0\mathrm{.}0,\backslash$

'duplicates not allowed r$\%$r: $\%$r $\%$r r$\text{'}\%\backslash$

$($i$\_$row,r$,$species$\_$member,s$\_$mtrx$[$i$\_$row,j$\_$col$])$

s$\_$mtrx$[$i$\_$row,j$\_$col$]=1\mathrm{.}0*$ coeff

else:

species$\_$member $=$ tmp$[0].$strip$()$

j$\_$col $=$ species.index$($species$\_$member$)$

assert s$\_$mtrx$[$i$\_$row,j$\_$col$]==0\mathrm{.}0,\backslash$

'duplicates not allowed r$\%$r: $\%$r $\%$r r$\text{'}\%\backslash$

$($i$\_$row,r$,$species$\_$member,s$\_$mtrx$[$i$\_$row,j$\_$col$])$

s$\_$mtrx$[$i$\_$row,j$\_$col$]=1\mathrm{.}0$

print$(\text{'}$species $=\text{'},$ species$)$

print$(\text{'}$# of species $=\text{'},$ len$($species$))$