function s=mojexp(x,eps) %MOJEXP racuna exp(x) direktno iz Taylorjeve vrste %POZOR!!! Neuporabno za velike x!!! %s=MOJEXP(x) %s je vsota vrste (odrezemo, ko clen v vrsti manjsi od 'eps') %x je argumet eksponentne funkcije, eps pa predpisana toleranca poz=1; %predpostavimo, da je x pozitiven s=1; clen=1; n=1; napaka=eps+1; if x<0 x=abs(x); poz=0; end while napaka>eps*abs(s) clen=clen*x/n; s=s+clen; n=n+1; napaka=abs(clen); end if poz==0 s=1/s; end