function  fnk = g(Z, T)
%Einheiten Umrechnen
atm=1.01325e5; %pa
l=10e-3; %Volumen
%Konstanten
a=-9.2211e-3*atm*l^3;   %pa*m^9/mol^3
A=9.7319e-2*atm*l^2;    %pa*m^6/mol^2
b=1.7976e-4*l^2;        %m^6/mol^2
B=1.9041e-2*l;          %m^3/mol
c=-2.4613e2*atm*l^3;    %pa*K^2*m^9/mol^3
C=3.8914e2*atm*l^2;     %pa*K*m^6/mol^2
alpha=-3.4215e-6*l^3;   %m^9/mol^3
gamma=1.89e-3*l^2;      %m^6/mol^2

p= 1.02925e5;             % p vor der 1. Verdichterstufe


fnk= -Z+1 + (B -(A/(R*T))-(C/T^3))*(p/(Z*R*T)) + (b-(a/(R*T)))*(p/(Z*R*T))^2+((a*alpha)/(R*T))*(p/(Z*R*T))^5+((c*(p/(Z*R*T))^2)/(R*T^3))*(1+gamma*(p/(Z*R*T))^2)*exp(-gamma*(p/(Z*R*T))^2);
 end