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

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

% trinn 1

T1= 383.941;
x1= [0.875; 0.125];
y1= 1-0.4266;
ycalc1=ycompw(T1,x1)

% trinn 2
T2= 371.68756;
x2= [0.625; 0.375];
y2= 0.2840;
ycalc2=ycompw(T2,x2)

% trinn 3
T3= 367.59374;
x3= [0.5; 0.5];
y3= 0.2011;
ycalc3=ycompw(T3,x3)

% trinn 4

T4= 364.05734;
x4= [0.375; 0.625];
y4= 0.1349;
ycalc4=ycompw(T4,x4)


% trinn 5

T5= 358.09155;
x5= [0.125; 0.875];
y5= 0.036;
ycalc5=ycompw(T5,x5)

% trinn 6
T6= 376.7462;
x6=[0.75;0.25];
y6= 0.39547;
ycalc6=ycompw(T6,x6)

% trinn 7
T7= 360.9142;
x7=[0.25;0.75];
y7= 0.08072;
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  ;





