1 subroutine y12maf(n, z, a, snr, nn, rnr, nn1, pivot, ha, 2 1iha,aflag,iflag,b,ifail) 3 implicit double precision (a-b,g,p,t-y), integer (c,f,h-n,r-s,z) 4 double precision a(nn), pivot(n), aflag(8),b(n) 5 integer snr(nn), rnr(nn1), ha(iha,11), iflag(10) 6 aflag(1)=16.0d0 7 aflag(2)=1.d-12 8 aflag(3)=1.d+16 9 aflag(4)=1.d-12 10 iflag(2)=2 11 iflag(3)=1 12 iflag(4)=0 13 iflag(5)=1 14 call y12mbf(n,z,a,snr,nn,rnr,nn1,ha,iha,aflag,iflag,ifail) 15 if(ifail.ne.0)go to 1 16 call y12mcf(n,z,a,snr,nn,rnr,nn1,pivot,b,ha,iha,aflag,iflag, 17 1 ifail) 18 if(ifail.ne.0)go to 1 19 call y12mdf(n,a,nn,b,pivot,snr,ha,iha,iflag,ifail) 201 return 21 end 22