function X = ssiacov(iniziali,parametri)

BETAH   = parametri(1);
BETAE   = parametri(2); 
BETAB   = parametri(3);
BETAS   = parametri(4);   
JEI     = parametri(5);
CHI     = parametri(6);
PHI      = parametri(7);
ALPHA    = parametri(8);
MH       = parametri(9);
GAMMAB   = parametri(10);
GAMMAE   = parametri(11);
GAMMAS   = parametri(12);
EPSILON  = parametri(13);
Ch      = iniziali(1);
D       = iniziali(2);
W       = iniziali(3);
N       = iniziali(4);
Hh      = iniziali(5);
He      = iniziali(6);
Rh      = iniziali(7);
Y       = iniziali(8);
q       = iniziali(9);
Le      = iniziali(10);
Lb      = iniziali(11);
Ls      = iniziali(12);
S       = iniziali(13);
Re      = iniziali(14);
Rs      = iniziali(15);
X = [Ch+D-Rh*D-W*N;
    Rh-(1/BETAH);
    (q/Ch)-JEI/Hh-BETAH*(q/Ch);
    CHI*N^PHI-W/Ch;
    Y-(He^ALPHA)*(N^(1-ALPHA));
    W-(1-ALPHA)*He^ALPHA*N^(-ALPHA);
    q-MH*q/Re-BETAE*(q+ALPHA*Y/He)+BETAE*MH*q;
    Le-MH*q*He/Re;
    (1-GAMMAB)*(Lb+(1-EPSILON)*S)-Lb-S+D;
    Lb-GAMMAE*D;
    (1-GAMMAS)*Ls-Ls+S;
    1-(1-GAMMAB*GAMMAE)/(EPSILON+GAMMAB-EPSILON*GAMMAB)-GAMMAE-BETAB*Rh+...
    1/(EPSILON+GAMMAB-EPSILON*GAMMAB)*BETAB*Rs*(1-GAMMAB*GAMMAE)+...
    BETAB*GAMMAE*Re;
    (1-GAMMAS)+BETAS*GAMMAS*Rs-BETAS*Re;
    Le-Ls-Lb; 
    He+Hh-1;]