1 real FUNCTION amotsa(p,y,psum,mp,np,ndim,pb,yb,funk,
3 INTEGER ihi,mp,ndim,np,NMAX
4 REAL fac,yb,yhi,p(mp,np),pb(np),psum(np),y(mp),funk
9 REAL fac1,fac2,tt,yflu,ytry,ptry(NMAX)
11 COMMON /ambsa/ tt,idum
15 ptry(j)=psum(j)*fac1-p(ihi,j)*fac2
24 yflu=ytry-tt*log(ran1(idum))
29 psum(j)=psum(j)-p(ihi,j)+ptry(j)
36 C (C) Copr. 1986-92 Numerical Recipes Software 0!5,.