Given stoichiometric$\_$matrix $=$ np$.$array$([[-1,1,0,0,0,0],)[1,-2,-1,-1,0,0],)[0,2,1,1,0,0],[0,-1,-0\mathrm{.}5,0,0,0],[0,0,0,-1,0,-1],[0,0,0,2,-2,1]$

And the species production rate density vector $[$mol$/($s cc$)]$ below,

g $=[1\mathrm{.}70838,1\mathrm{.}38896,1\mathrm{.}44727,0\mathrm{.}80652,2\mathrm{.}27862,2\mathrm{.}03963,2\mathrm{.}81438,0\mathrm{.}30002]$

in the order: C$($g$),$ C$($s$),$ CH$4,$ CO$,$ CO$2,$ H$2,$ H$2$O$,$ O$2$

compute the reaction rates vector r on Python and make a plot Given that S$^$T$*$r $=$ g$.$ Use your own algorithms and explain your work, that is$,$ what is the meaning of your computed reaction rate vector? What problem are you solving to obtain the reaction rate vector from the given production rate vector? Do the reaction rate equations have a solution? If so$,$ is it unique? Why or