1
subroutine trbize(nmaxi,maxdeg,nzm,nzero,sm,zm,zzr,zzi)
3
c bilinear transformation of the zeros and the locations of the extrema
6
implicit double precision (a-h,o-z)
8
dimension sm(maxdeg,*), nzero(*)
9
double precision zm(maxdeg,*), zzr(*), zzi(*)
10
external slamch, dlamch
13
flma=2.0d+0**(int(slamch('l'))-2)
14
flmi=2.0d+0*dlamch('p')
20
zm(i,j) = 2.0d+0*atan(qi)
22
if (qi.ge.flma) go to 10
23
if (qi.lt.flmi) go to 20
26
zzr(i) = (1.0d+0-qqi)/q