Can you help me figure out what's wrong with my matlab code. I am solving a set of differential equations and algebraic equations using matlab. It's running and not ending and telling me there is an issue with fsolve. function $dxds=$ HFcountercurrentflowODE $(s,x)$

$\%$ CROSSFLOW cross flow model

global c$1$ c$2$ z ph pl permN$2$ d do Nf Uf mu$1$ R T di pi

$c1=0\mathrm{.}00000000004$;

$c2=0\mathrm{.}0004696$;

$d=0\mathrm{.}0023876$;

$ph=101325$;

$pl=0\mathrm{.}001$;

permN$2=0\mathrm{.}00000000000000104836$;

$do=0\mathrm{.}0002$;

$z=0\mathrm{.}2$;

$Nf=10000$;

$Uf=0\mathrm{.}1$;

mu$1=0\mathrm{.}00001805$;

$R=8\mathrm{.}314$;

$T=303\mathrm{.}15$;

$di=0\mathrm{.}0001$;

$=3\mathrm{.}14159265$;

permH$2$O ph $\{$:$)$;

pizpermN$2*$ph$)\frac{?}{(d**Uf)}$;

mul $**$ RUfd$)$ pipermN$2$phdodidiNf$)$;

selectivityH$2$O$=$ permH$2$O$/$permN$2$;

fun selectivityH $20**x(1)**y(3)$ selectivityH$20*$x$(4)))$;

$y(2)-\frac{x(2)**y(3)}{1-x(4)+(y(3)**x(4))}$;

selectivityH$20*$y$(2))]$;

$y0=[0\mathrm{.}001$;$0\mathrm{.}001$;$0\mathrm{.}001]$;

$y=$ fsolve $($fun$,$ yo$)$;

selectivityH$20*($x$(1)-(x(4)**y(1))$;

$(\frac{1}{x}(3))**(x(2)-(x(4)**y(2)))**((x\frac{2}{y}(2))-1)$;

$(-\frac{1-x(4)}{y}(3))$;

A$*(1-x(3))$;

endThanks