function [erg] = intExp_ax(a,h)
s = size(a);
if s ~= 2
    error('a muss 2 werte haben');
end
for i = 1:2
    if a(i) < 0 
    error('a muss ueberall groesser gleich als wie 0 sein!!');
    end
end
if h <= 0
    error('h muss echt groesser als wie die null sein');
elseif a(1) == 0 & a(2) == 0  % Fall a=[0 0]
    erg = (h^2)/2;
elseif a(1) == a(2) & a(1) <= 0.0000001
    erg = (h^2)/2;
elseif a(1) == 0 & a(2) ~= 0  % Fall a1=0 und a2!=0 und h egal
    if a(2) <= 0.0001
        erg = 1/2*h^2+1/6*a(2)*h^3+1/24*a(2)^2*h^4+1/120*a(2)^3*h^5;
    else
        erg = -(-exp(a(2)*h)+a(2)*h+1)/(a(2)^2);
    end
elseif a(2) == 0 & a(1) ~= 0  % Fall a2=0 und a1!=0 und h egal
    if a(1) <= 0.0001
        erg = 1/2*h^2+1/6*a(1)*h^3+1/24*a(1)^2*h^4+1/120*a(1)^3*h^5;
    else
        erg = -(-exp(a(1)*h)+a(1)*h+1)/(a(1)^2);
    end 
%Folgende vier Zeilen bestehen aus reinem Nonsens
%elseif a(1) ~= 0 & a(2) == 0 & h >= 0.000001 %Fall a1!=0 und a2=0 und h gross
%    erg = -(a(1)*h+1-exp(a(1)*h))/(a(1)^2);
%elseif a(1) ~= 0 & a(2) == 0 & h < 0.000001 % Fall a1!=0 und a2=0 und h klein also Taylor
%    erg = 1/2*h^2+1/6*a(1)*h^3+1/24*a(1)^2*h^4+1/120*a(1)^3*h^5;
elseif a(1) == a(2) & h >= 0.000001 % Fall a1 = a2 und h gross
    erg = (1-exp(a(1)*h)+exp(a(1)*h)*a(1)*h)/a(1)^2;
elseif a(1) == a(2) & h < 0.000001 % Fall a1 = a2 und h klein also Taylor
    erg = 1/2*h^2+1/3*a(1)*h^3+1/8*a(1)^2*h^4+1/30*a(1)^3*h^5;
elseif abs((a(1)-a(2))) <= 0.00001
    erg = -(1/2*a(2)^2*a(1)-1/2*a(1)^2*a(2))/(a(2)*a(1)*(a(1)-a(2)))*h^2-(-1/6*a(1)^3*a(2)+1/6*a(2)^3*a(1))/(a(2)*a(1)*(a(1)-a(2)))*h^3-(1/24*a(2)^4*a(1)-1/24*a(1)^4*a(2))/(a(2)*a(1)*(a(1)-a(2)))*h^4-(1/120*a(2)^5*a(1)-1/120*a(1)^5*a(2))/(a(2)*a(1)*(a(1)-a(2)))*h^5;
elseif h < 0.000001 % Hier wird die 'Schneider'-Formel benutzt, weil die e-Fkt aufgrund von Auslöschung Scheiße baut
    erg = -(1/2*a(2)^2*a(1)-1/2*a(1)^2*a(2))/(a(2)*a(1)*(a(1)-a(2)))*h^2-(-1/6*a(1)^3*a(2)+1/6*a(2)^3*a(1))/(a(2)*a(1)*(a(1)-a(2)))*h^3-(1/24*a(2)^4*a(1)-1/24*a(1)^4*a(2))/(a(2)*a(1)*(a(1)-a(2)))*h^4-(1/120*a(2)^5*a(1)-1/120*a(1)^5*a(2))/(a(2)*a(1)*(a(1)-a(2)))*h^5;
else % Hier wird die e-Fkt benutzt, weil der 'Schneider' Scheiße baut
    erg = -(-a(1)+a(2)+exp(a(2)*h)*a(1)-exp(a(1)*h)*a(2))/(a(2)*a(1)*(a(1)-a(2)));
end