function [LSS, VSS, SPSS]=steadystate_helper(bet, tet, chi, tau)
L0=1/2;
V0=1;
SP0=1/3;
x0=[L0; V0; SP0];

F=@(LSS, VSS, SPSS) [(LSS)-(VSS)/(tet-bet*((VSS)*(SPSS))); 1-chi*((SPSS)*(LSS)+(1/bet))-tau*(LSS); (VSS)-chi-(1-chi)*bet*(VSS)*((SPSS)*(LSS)+(1/bet))];
options = optimset ('Display', 'Final', 'Tolx', 1e-10,'Tolfun',1e-10);
fcn = @(b) F(b(1),b(2),b(3));
B=fsolve(fcn, x0, options);
LSS=B(1);
VSS=B(2);
SPSS=B(3);
