% Fil wvann.m

% Funksjonen beregner Wilsonparametere for sytemet 

% 2-metoxyetanol/vann
% Inn data er damp/vaeske sammensetninger for 2-metoxyetanol/vann

% som er generert i Aspen

% Programmet bruker funksjonen fmins,og et kavdratisk avvik

% Funksjonen bruker funksjonen ycompw.m og wilsonvw.m, dette er funksjoner 

% som tar inn en sammensetningsvektor x og temperaturen T 

% Start av programmet : x=fmins('wiso',W)

% W er en 2x2 matrise som gir start parameter for bergningene 

function f = wvann(W)

global Wglob
Wglob=[0 W(1);W(2) 0]

% trinn 1

T1= 116.7271+273.15;
x1= [0.875; 0.125];
y1= 0.6875;
ycalc1=ycompw(T1,x1)

% trinn 2
T2= 111.1147+273.15;
x2= [0.75; 0.25];
y2= 0.4940;
ycalc2=ycompw(T2,x2)

% trinn 3
T3= 106.9411+273.15;
x3= [0.625; 0.375];
y3= 0.3635;
ycalc3=ycompw(T3,x3)

% trinn 4

T4= 103.8477+273.15;
x4= [0.5; 0.5];
y4= 0.2707;
ycalc4=ycompw(T4,x4)


% trinn 5

T5= 101.6456+273.15;
x5= [0.375; 0.625];
y5= 0.2023;
ycalc5=ycompw(T5,x5)

% trinn 6
T6= 100.4549+273.15;
x6=[0.275;0.725];
y6= 0.1593;
ycalc6=ycompw(T6,x6)

% trinn 7
T7= 99.50169+273.15;
x7=[0.125;0.875];
y7= 0.1015;
ycalc7=ycompw(T7,x7)

% Bergener kvadratavviket


f=(y1-ycalc1(1))^2+(y2-ycalc2(1))^2+(y3-ycalc3(1))^2+(y4-ycalc4(1))^2+(y5-ycalc5(1))^2+(y6-ycalc6(1))^2+(y6-ycalc6(1))^2   ;

