// Die Exponentialfunktion nach Definition durch die Potenzreihe function X = my_exp(A) S= eye(A); // fuer k-ten Summand S_k X= zeros(A); // fuer Summe der S_k for k=0:1000 disp(S); X= X + S; S= S * A/(k+1); if norm(S)<=1E-20 then break; end end // disp(k); endfunction