Re: Random posts
Posted: October 15th, 2020, 5:44 am
Javascript agm:
e.g. agm(1,2)=1.4567910310469068
PARI lambertw:
e.g. W(-1/4)=-0.35740295618138890306881110405590475332
Code: Select all
function agm(a,b){var c,max_diff=Infinity,diff;for(;;){diff=Math.abs(a-b);c=(a+b)/2;if(diff>=max_diff)return c;max_diff=diff;b=Math.sqrt(a*b);a=c;}}PARI lambertw:
Code: Select all
W(y)=if(y<0,solve(x=-1,0,x*exp(x)-y),lambertw(y))