% Fil wilsonvw.m
% Funksjonen brukes for aa generere Wilsonkoeffisienter
% Funksjonen beregner aktivitetskoeffisentene til vaesken ved 
% Wilson metode paa vektorisert form
% Funksjonen er laget av Steinar Hauan

function gamma=wilsonvw(Tk,x)

global Wglob Vglob Xglob Rglob NCglob;

W =Wglob;
V =Vglob;
X=x';
NC=NCglob;
R= Rglob;

% Lager en NCxNC matrise av vektorene V og X

Vmat=V(ones(1,NC),:);
Xmat=X(ones(1,NC),:);

% Vektorisering av en for-loop

Lambda =Vmat./Vmat'.*exp(-W/(R*Tk));
ledd2=sum((Xmat.*Lambda)');
Xs=X./ledd2;
Xsmat=Xs(ones(1,NC),:);
ledd3=sum((Xsmat.*Lambda')');

gamma=(exp(1-log(ledd2)-ledd3))';


