function dS=test(x)
R=1.987;
T=1000;
x1=x(1);
x2=x(2);
x3=x(3);
x4=x(4);
x5=x(5);


S1=exp(x1)*(4610/R/T+x1-log(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5)))+exp(x1)*(1-exp(x1)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5)))-exp(x2)*exp(x1)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))-exp(x3)*exp(x1)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))-exp(x4)*exp(x1)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))-exp(x5)*exp(x1)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5));
S2=-exp(x2)*exp(x1)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))+exp(x2)*(-46030/R/T+x2-log(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5)))+exp(x2)*(1-exp(x2)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5)))-exp(x3)*exp(x2)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))-exp(x4)*exp(x2)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))-exp(x5)*exp(x2)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5));
S3=-exp(x3)*exp(x1)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))-exp(x3)*exp(x2)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))+exp(x3)*(-47940/R/T+x3-log(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5)))+exp(x3)*(1-exp(x3)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5)))-exp(x4)*exp(x3)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))-exp(x5)*exp(x3)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5));
S4=-exp(x4)*exp(x1)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))-exp(x4)*exp(x2)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))-exp(x4)*exp(x3)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))+exp(x4)*(-94610/R/T+x4-log(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5)))+exp(x4)*(1-exp(x4)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5)))-exp(x5)*exp(x4)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5));
S5=-exp(x5)*exp(x1)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))-exp(x5)*exp(x2)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))-exp(x5)*exp(x3)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))-exp(x5)*exp(x4)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5))+exp(x5)*(x5-log(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5)))+exp(x5)*(1-exp(x5)/(exp(x1)+exp(x2)+exp(x3)+exp(x4)+exp(x5)));

dS=[S1 S2 S3 S4 S5 ]'

endfunction
