function [x0,data]=lineariseringspunkt(z,w_Gin,w_Lin,data)

%=========================================================================\
%Funksjon for å finne den stasjonære løsningen x0 det skal lineariseres
%                               rundt.
%=========================================================================/

%=========================================================================\
%Tuningskonstanter for å tilfredstille OLGA-data.
%=========================================================================/
%Disse er funnet fra kap 4.2.1 i paperet "A low-dimensional dynamic model
%of severe slugging for control design and analysis".
data.K1=0.0054
data.K2=6.84
data.K3=0.11
data.n=2.3
data.V_G1=12.64
data.M_G=20.6
%=========================================================================\
%(B.1)-(B.3) settes lik 0 for å finne den (ustabile) stasjonære løsningen
%                       ved gitt ventilåpning z.
%=========================================================================/
w_Lout=w_Lin; %B.1
w_G1=w_Gin;   %B.2  
w_Gout=w_G1;  %B.3
%=========================================================================\
%De tilstandsavhengige interne variablene beregnes
%=========================================================================/
alpha_Lm=w_Lin/(w_Gin+w_Lin); %alpha_Lm er væskefraksjonen (massebasis) oppstrøms choke-valven. Denne relasjonen gjelder KUN i det stasjonære tilfellet.
w_out=w_Gout/(1-alpha_Lm);    %B.23; Nå kan jeg kanskje benytte B.22 for å finne P2 eller rho_T. Men må finne en av disse variablene på en annen måte. HVILKEN/HVORDAN?
                              %Er nok i hovedsak rho_T som går an å finne på en annen måte.  








V_L=V_LR+h1*data.A1; %B.7 må finne h1 og V_LR
%=========================================================================\
%Beregning av x0
%=========================================================================/
mL=V_L*data.rho_L; %B.6 må finne V_L : OK
mG1=P1*data.V_G1*data.M_G/(data.R*data.T); %B.4 her må jeg ha trykket P1
mG2=P2*V_G2*data.M_G/(data.R*data.T);%B.12 her må jeg ha trykket P2 og volumet V_G2

x0=[mL;mG1;mG2];
end