% Fil wmet.m
% Funksjonen beregner Wilsonparametere for sytemet 
% 2-metoxyetanol/metanol
% Inn data er damp/vaeske sammensetninger for 2-metoxyetanol/metanol
% 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 = wmet(W)

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

% trinn 1

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

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

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

% trinn 4

T4= 80.63078+273.15;
x4= [0.5; 0.5];
y4= 0.1063;
ycalc4=ycompw(T4,x4)


% trinn 5

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

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

% trinn 7
T7= 67.58888+273.15;
x7=[0.125;0.875];
y7= 0.01491;
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   ;





