Mercurial > hg > ishara

annotate general/numerical/scalar/lambertw.m @ 35:f1ce7876346a

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Updated docs.
author samer
date Mon, 21 Jan 2013 11:01:45 +0000
parents db7f4afd27c5
children
rev   line source
samer@4 1 function W = lambertw(x)
samer@4 2 % lambertw - Lambert's W function
samer@4 3 %
samer@4 4 % lambertw :: real -> real.
samer@4 5 %
samer@4 6 % evaluates W(y) for -1 branch of Lambert W function where
samer@4 7 % y = -exp(-x)
samer@4 8 %
samer@4 9 % The W function satisfies W(t)*exp(W(t)) = t
samer@4 10
samer@4 11 if x<1, error('lambertw: argument must be >= 1'); end
samer@4 12 W=fixpoint(@(w)(-x-log(abs(w))),-x,'its',60);