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