calc_epsilon.m

These are matlab and simulink files to model the membrane crystallization system, including the matl

function Fun = Calc_epsilon(x)

global sigma P rhoLac kV kB kG V Wsatcry Wleavecry 

    G = kG*((Wleavecry-Wsatcry)/Wsatcry);
    B = kB*((Wleavecry-Wsatcry)/Wsatcry)^3;
    sigma=2.8375e-4;
  
    Fun = [x(6)-P/(rhoLac*(1-x(7)));
    B*(sigma)^3-x(1)*x(6)/V;
    x(1)*G/sigma-x(2)*x(6)/V;
    2*x(2)*G/sigma-x(3)*x(6)/V;
    3*x(3)*G/sigma-x(4)*x(6)/V;
    4*x(4)*G/sigma-x(5)*x(6)/V;
    1-kV*x(4)-x(7)];