function  v = absorption_geometrie_umgesch (~)


% T in Kelvin
% p in bar
% sigma in cm^2
% k bezogen auf mm

global l

%Beispiel zur Berechnung des gekoppelten Geometrie-Absorptionsfaktors


% siehe Hartinger et al.

    function y=k0(T)
        y=5.61*exp(-3298/T);
    end

    function y=b3(T)
        y=9.9145*exp((T-1091.4)/454.7)^2;
    end

    function y=kO2(x,pO2,T)
        y=(k0(T).*b3(T).^0.5)./(b3(T)+pO2.*x).^0.5*1e-1;
    end

    function y=kCO2(T)
        y=7.59*exp(-3617/T)*1e-1;
    end

    % siehe HVG-Bericht (Monkhouse)

    function y=sigH2O(T)
        y= exp(32.1627-5.7426*(T/1000)+13.9331/(T/1000)+(-16.0154+3.2252*...
            (T/1000)-8.7988/(T/1000))*1.94)*1e-19;
    end

    function y=kH2O(T)
        y= sigH2O(T)/1.381e-23/T*1e-2;
    end

    function y=kCO(T)
        y= 7.59*exp(-3617/T)*1e-3;
    end

    function y=sigH2S(T)
        y= 6.3e-18;
    end

    function y=kH2S(T)
        y= sigH2S(T)/1.381e-23/T*1e-2;
    end

    function  y=kNH3(T)
        y=  230*(1-exp(-1370/T))*(300/T)*9.869e-2;
    end

    function y = f(x,T,pCO2,pH2O,pO2,pCO,pH2S,pNH3)			% Absorptionsfunktion
        y= exp(-(kCO2(T).*pCO2+kH2O(T).*pH2O+kCO(T).*pCO+kH2S(T)*pH2S+kNH3(T)...
            .*pNH3).*x).*exp(-kO2(x,pO2,T).*pO2.*x);
    end


l = 44.4;

function y=g(x)
            
    global T
    global pO2
    global pCO
    global pCO2
    global pH2O
    global pH2S
    global pNH3    
    global a
    global off

    y=f(x,1073.15,23,85,7.8,12.7,42.9,14.8).*8^2./(x+99.2).^2;
    
end

[v] = quad(@g,0,l)


end