c*********************** problem name: usmap ************************** c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine a1xy(x,y,u,ux,uy,rl,itag,values) c implicit double precision (a-h,o-z) implicit integer (i-n) c double precision + values(*) character*80 + su common /val0/k0,ku,kx,ky,kl,kuu,kux,kuy,kul, + kxu,kxx,kxy,kxl,kyu,kyx,kyy,kyl,klu,klx,kly,kll common /atest2/iu(100),a1x,a1y,a1u,a2x,a2y,a2u, + bux,buy,cu0,cu1,ru(90),su(100) c values(k0)=a1x*ux+a1y*uy+a1u*u values(ku)=a1u values(kx)=a1x values(ky)=a1y return end c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine a2xy(x,y,u,ux,uy,rl,itag,values) c implicit double precision (a-h,o-z) implicit integer (i-n) c double precision + values(*) character*80 + su common /val0/k0,ku,kx,ky,kl,kuu,kux,kuy,kul, + kxu,kxx,kxy,kxl,kyu,kyx,kyy,kyl,klu,klx,kly,kll common /atest2/iu(100),a1x,a1y,a1u,a2x,a2y,a2u, + bux,buy,cu0,cu1,ru(90),su(100) c values(k0)=a2x*ux+a2y*uy+a2u*u values(ku)=a2u values(kx)=a2x values(ky)=a2y return end c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine fxy(x,y,u,ux,uy,rl,itag,values) c implicit double precision (a-h,o-z) implicit integer (i-n) c double precision + values(*) character*80 + su common /val0/k0,ku,kx,ky,kl,kuu,kux,kuy,kul, + kxu,kxx,kxy,kxl,kyu,kyx,kyy,kyl,klu,klx,kly,kll common /atest2/iu(100),a1x,a1y,a1u,a2x,a2y,a2u, + bux,buy,cu0,cu1,ru(90),su(100) c values(k0)= - bux*ux - buy*uy - cu0 - cu1*u values(ku)= - cu1 values(kx) = - bux values(ky) = - buy return end c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine gnxy(x,y,u,rl,itag,values) c implicit double precision (a-h,o-z) implicit integer (i-n) c double precision + values(*) common /val1/k0,ku,kl,kuu,kul,klu,kll c return end c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine gdxy(x,y,rl,itag,values) c implicit double precision (a-h,o-z) implicit integer (i-n) c double precision + values(*) common /val2/k0,kl,kll,klb,kub,kic,kim,kil c return end c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine p1xy(x,y,u,ux,uy,rl,itag,values) c implicit double precision (a-h,o-z) implicit integer (i-n) c double precision + values(*) common /val0/k0,ku,kx,ky,kl,kuu,kux,kuy,kul, + kxu,kxx,kxy,kxl,kyu,kyx,kyy,kyl,klu,klx,kly,kll c values(k0)=1.0d0 return end c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine p2xy(x,y,dx,dy,u,ux,uy,rl,itag,jtag,values) c implicit double precision (a-h,o-z) implicit integer (i-n) c double precision + values(*) common /val0/k0,ku,kx,ky,kl,kuu,kux,kuy,kul, + kxu,kxx,kxy,kxl,kyu,kyx,kyy,kyl,klu,klx,kly,kll c return end c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine qxy(x,y,u,ux,uy,rl,itag,values) c implicit double precision (a-h,o-z) implicit integer (i-n) c double precision + values(*) common /val3/kf,kf1,kf2,ksk,kad c s=x**2+y**2 values(kf)=dfloat(itag) values(ksk)=s values(kad)=s return end c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine usrcmd(vx,vy,xm,ym,itnode,ibndry,ip,rp,sp,iu,ru,su,w) c implicit double precision (a-h,o-z) implicit integer (i-n) integer + itnode(5,*),ibndry(6,*),ip(100),iu(100) double precision + vx(*),vy(*),xm(*),ym(*),rp(100),ru(100),w(*) character*80 + sp(100),su(100),file(80) save len,file c data len/68/ data (file(i),i= 1, 10)/ + 'n i= 1,n=istate,a=is,t=i', 1 'n i= 2,n=icty ,a=ic,t=i', 2 'n i= 3,n=isle ,a=il,t=i', 3 'n i= 1,n=a1x, a=x1,t=r', 4 'n i= 2,n=a1y, a=y1,t=r', 5 'n i= 3,n=a1u, a=u1,t=r', 6 'n i= 4,n=a2x, a=x2,t=r', 7 'n i= 5,n=a2y, a=y2,t=r', 8 'n i= 6,n=a2u, a=u2,t=r', 9 'n i= 7,n=bux, a=bx,t=r'/ data (file(i),i= 11, 20)/ + 'n i= 8,n=buy, a=by,t=r', 1 'n i= 9,n=cu0, a=c0,t=r', 2 'n i=10,n=cu1, a=c1,t=r', 3 's n=istate,v= 1,l="alabama"', 4 's n=istate,v= 2,l="alaska"', 5 's n=istate,v= 3,l="arizona"', 6 's n=istate,v= 4,l="arkansas"', 7 's n=istate,v= 5,l="california"', 8 's n=istate,v= 6,l="colorado"', 9 's n=istate,v= 7,l="connecticut"'/ data (file(i),i= 21, 30)/ + 's n=istate,v= 8,l="deleware"', 1 's n=istate,v= 9,l="florida"', 2 's n=istate,v=10,l="georgia"', 3 's n=istate,v=11,l="hawaii"', 4 's n=istate,v=12,l="idaho"', 5 's n=istate,v=13,l="illinois"', 6 's n=istate,v=14,l="indiana"', 7 's n=istate,v=15,l="iowa"', 8 's n=istate,v=16,l="kansas"', 9 's n=istate,v=17,l="kentucky"'/ data (file(i),i= 31, 40)/ + 's n=istate,v=18,l="louisiana"', 1 's n=istate,v=19,l="maine"', 2 's n=istate,v=20,l="maryland"', 3 's n=istate,v=21,l="massachusetts"', 4 's n=istate,v=22,l="michigan"', 5 's n=istate,v=23,l="minnesota"', 6 's n=istate,v=24,l="mississippi"', 7 's n=istate,v=25,l="missouri"', 8 's n=istate,v=26,l="montana"', 9 's n=istate,v=27,l="nebraska"'/ data (file(i),i= 41, 50)/ + 's n=istate,v=28,l="nevada"', 1 's n=istate,v=29,l="new hampshire"', 2 's n=istate,v=30,l="new jersey"', 3 's n=istate,v=31,l="new mexico"', 4 's n=istate,v=32,l="new york"', 5 's n=istate,v=33,l="north carolina"', 6 's n=istate,v=34,l="north dakota"', 7 's n=istate,v=35,l="ohio"', 8 's n=istate,v=36,l="oklahoma"', 9 's n=istate,v=37,l="oregon"'/ data (file(i),i= 51, 60)/ + 's n=istate,v=38,l="pennsylvania"', 1 's n=istate,v=39,l="rhode island"', 2 's n=istate,v=40,l="south carolina"', 3 's n=istate,v=41,l="south dakota"', 4 's n=istate,v=42,l="tennessee"', 5 's n=istate,v=43,l="texas"', 6 's n=istate,v=44,l="utah"', 7 's n=istate,v=45,l="vermont"', 8 's n=istate,v=46,l="virginia"', 9 's n=istate,v=47,l="washington"'/ data (file(i),i= 61, 68)/ + 's n=istate,v=48,l="west virginia"', 1 's n=istate,v=49,l="wisconsin"', 2 's n=istate,v=50,l="wyoming"', 3 's n=istate,v=51,l="united states"', 4 's n=icty ,v=0,l="no county lines"', 5 's n=icty ,v=1,l="county lines"', 6 's n=isle ,v=0,l="no small islands"', 7 's n=isle ,v=1,l="small islands"'/ c istate=iu(1) icty=iu(2) isle=iu(3) c call usrset(file,len,iu,ru,su) c if(istate.ne.iu(1)) ip(41)=-1 if(icty.ne.iu(2)) ip(41)=-1 if(isle.ne.iu(3)) ip(41)=-1 c ispd=ip(8) a1x=ru(1) a1y=ru(2) a1u=ru(3) a2x=ru(4) a2y=ru(5) a2u=ru(6) bux=ru(7) buy=ru(8) cu0=ru(9) cu1=ru(10) if(ispd.eq.1) then if((a1y.ne.a2x).or.(a1u.ne.0.0d0).or. + (a2u.ne.0.0d0).or.(bux.ne.0.0d0).or. 1 (buy.ne.0.0d0)) ip(8)=0 else if((a1y.eq.a2x).and.(a1u.eq.0.0d0).and. + (a2u.eq.0.0d0).and.(bux.eq.0.0d0).and. 1 (buy.eq.0.0d0)) ip(8)=1 endif if(ispd.ne.ip(8)) ip(6)=max0(ip(26),ip(6),1) c return end c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine gdata(vx,vy,xm,ym,itnode,ibndry,ip,rp,sp,iu,ru,su,w) c implicit double precision (a-h,o-z) implicit integer (i-n) integer + itnode(5,*),ibndry(6,*),ip(100),iu(100), 1 ix(7000),iy(7000),ib1(7000),ib2(7000), 2 ib(5),jb(5),jv(5) double precision + vx(*),vy(*),xm(*),ym(*),rp(100),ru(100),w(*) character*80 + sp(100),su(100),state save hmax,grade c c map of states in the united states c data hmax,grade/0.1d0,1.5d0/ c c if(ip(41).eq.1) then iu(1)=51 iu(2)=1 iu(3)=0 c c initialize as laplacian with dirichlet b.c. c do i=1,10 ru(i)=0.0d0 enddo ru(1)=1.0d0 ru(5)=1.0d0 ru(9)=1.0d0 ip(8)=1 c sp(6)='usmap_mpixxx.rw' sp(7)='usmap.jnl' sp(9)='usmap_mpixxx.out' endif c ip(6)=max0(ip(6),ip(26),1) ip(7)=1 c c istate=iu(1) if(istate.lt.1.or.istate.gt.51) istate=51 icty=iu(2) if(icty.ne.0) icty=1 isle=iu(3) if(isle.ne.1) isle=0 c if (istate.eq.1) then call al(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.2) then call ak(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) c c fixup for alaska c if(isle.eq.0) then isle=1 npart=1 endif else if (istate.eq.3) then call az(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.4) then call ar(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.5) then call ca(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.6) then call co(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.7) then call ct(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.8) then call de(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.9) then call fl(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.10) then call ga(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.11) then call hi(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) c c fixup for hawaii c isle=1 else if (istate.eq.12) then call id(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.13) then call il(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.14) then call in(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.15) then call ia(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.16) then call ks(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.17) then call ky(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.18) then call la(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.19) then call me(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.20) then call ml(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.21) then call ma(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.22) then call mi(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.23) then call mn(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.24) then call ms(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.25) then call mo(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.26) then call mt(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.27) then call nb(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.28) then call nv(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.29) then call nh(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.30) then call nj(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.31) then call nm(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.32) then call ny(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.33) then call nc(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.34) then call nd(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.35) then call oh(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.36) then call ok(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.37) then call og(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.38) then call pa(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.39) then call ri(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.40) then call sc(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.41) then call sd(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.42) then call tn(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.43) then call tx(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.44) then c*** call ut(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) c c fixup for salt lake c call uh(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) if(icty.eq.1) ib(1)=ib(1)+208 else if (istate.eq.45) then call vt(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.46) then call va(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.47) then call wa(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.48) then call wv(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.49) then call wi(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.50) then call wy(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) else if (istate.eq.51) then call us(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) endif c sp(2)=state sp(1)=state sp(3)=state sp(4)=state c nvf=0 nbf=0 ntf=1 ncf=0 do nn=1,npart nvf0=nvf nbf0=nbf nbound=nbf0+ib(nn) if(jb(nn+1)-jb(nn).eq.ib(nn).and.isle.eq.0) go to 10 if(icty.eq.1) then j1=jv(nn+1)-1 j2=jb(nn+1)-1 else j1=jv(nn)+ib(nn)-1 j2=jb(nn)+ib(nn)-1 endif if(irot.eq.0) then do j=jv(nn),j1 nvf=nvf+1 vx(nvf)=dfloat(ix(j)) vy(nvf)=dfloat(iy(j)) enddo else do j=jv(nn),j1 nvf=nvf+1 vx(nvf)=dfloat(iy(j)) vy(nvf)=-dfloat(ix(j)) enddo endif do j=jb(nn),j2 nbf=nbf+1 ibndry(1,nbf)=ib1(j)+nvf0 ibndry(2,nbf)=ib2(j)+nvf0 ibndry(3,nbf)=0 ibndry(4,nbf)=2 if(nbf.gt.nbound) ibndry(4,nbf)=0 ibndry(5,nbf)=0 ibndry(6,nbf)=0 enddo 10 enddo c ip(1)=ntf ip(2)=nvf ip(3)=ncf ip(4)=nbf ip(20)=5 rp(15)=hmax rp(16)=grade c c make itnode c call sklutl(0,vx,vy,xm,ym,itnode,ibndry,ip,w,iflag) c c call box(nvf,vx,vy,ncf,xm,ym) c return end c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine box(nvf,vx,vy,ncf,xm,ym) c implicit double precision (a-h,o-z) implicit integer (i-n) double precision + vx(*),vy(*),xm(*),ym(*) c c rescale domain to lie inside the unit square c xmin=vx(1) xmax=xmin ymin=vy(1) ymax=ymin do i=1,nvf xmin=dmin1(xmin,vx(i)) xmax=dmax1(xmax,vx(i)) ymin=dmin1(ymin,vy(i)) ymax=dmax1(ymax,vy(i)) enddo c c compute scaled coordinates c xmid=(xmin+xmax)/2.0d0 ymid=(ymin+ymax)/2.0d0 scale=1.0d0/dmax1(xmax-xmin,ymax-ymin) do i=1,nvf vx(i)=0.5d0+scale*(vx(i)-xmid) vy(i)=0.5d0+scale*(vy(i)-ymid) enddo if(ncf.le.0) return do i=1,ncf xm(i)=0.5d0+scale*(xm(i)-xmid) ym(i)=0.5d0+scale*(ym(i)-ymid) enddo return end c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine al(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) c implicit double precision (a-h,o-z) implicit integer (i-n) integer + jv(*),jb(*),ib(*),ix(*),iy(*),ib1(*),ib2(*),kv(5), 1 kb(5),iib(5),jx(1960),jy(1960),jb1(2020),jb2(2020) character*80 + state,istate save kv,kb,iib,jx,jy,jb1,jb2,istate,jrot,ipart data istate/'alabama '/ data jrot,ipart/0,2/ data kv/ 1,1879,1952, 0, 0/ data kb/ 1,1945,2018, 0, 0/ data iib/ 309, 73, 0, 0, 0/ data (jx(i),i= 1, 100)/ + 10416,18864,27408,22920, 9984,35088,34344, 9624,36264, 9288, 1 8544,37272,38040, 7848,39024,40416, 6936,39936,39312, 6432, 2 39792, 6672,35904,14904,28944,23976,10800,10152, 9912, 9432, 3 9456,34128,38544,38592,38760,38832,38712,38712,39144,38976, 4 40176,39648,39624,39528,39576,39360,39504,39504,39264,39216, 5 39264,39192,40488,40704,40776,40632,40728,40560,40704,40608, 6 40584,40392,40392,41040,41520,41352,41184,41184,41280,41184, 7 40728,40656,40824,40368,39888,40008,39936,39960,39816,39840, 8 39672,39696,39408,39456,39432,39216,39264,39168,39240,39168, 9 6360,39384,39336,39768,39744,39960,40032,40176,40128,39936/ data (jx(i),i= 101, 200)/ + 39960,39864,39864,39720,39816,39744,39576,39528,39720,39624, 1 39720,39744,39912,40416,40440,40608,40704,23136,14880,16512, 2 16584,16416,16392,16584,16704,16800,16896,16896,16776,16368, 3 16104,16056,15864,15552,15408,15072,14880,14832,14640,14592, 4 14520,14976,14976, 6864, 7080, 7032, 7104, 7176, 7200, 7272, 5 7488, 7680, 7728, 7680, 7800, 7896, 8016, 7992, 8256, 8592, 6 8736, 8976, 9000, 8928, 8856, 8928, 8880, 9288, 9240, 9360, 7 9384, 9480, 9504, 9360, 9456, 9744, 9816, 9768, 9792, 9744, 8 9864, 9816,10032, 9984, 9792, 9720, 9552, 9576, 9672, 9744, 9 9864, 9744, 9840,10008,10272,10272,10320,10320,10536,10512/ data (jx(i),i= 201, 300)/ + 10656,12792,11136,10512,10632,11160,11232,11280,11208,11328, 1 11304,11400,11520,11664,11832,12144,12576,13032,13152,13368, 2 13440,13416,13512,13608,13488,13440,13320,13248,13128,12720, 3 12504,12456,12504,12672,12576,12504,12384,12360,12432,12360, 4 12144,11736,11616,11472,11808,11808,11688,11736,11664,11688, 5 11592,11616,11496,11568,11544,11400,14640,15240,15240,15096, 6 15240,15216,15432,15696,15696,15624,15360,15240,15000,14712, 7 14712,14880,14784,14784,14688,14808,14688,14832,14952,15000, 8 15120,15000,15072,15168,15240,15288,15408,15480,15648,15720, 9 15624,15600,15672,15672,15768,15768,15816,15840,15888,15936/ data (jx(i),i= 301, 400)/ + 16008,16200,16344,16368,16536,16848,16920,16536,16464,16800, 1 16224,16008,15984,15672,15288,15168,14112,13896,13632,13176, 2 12864,12096,11952,11808,11616,10968,10608,18816,18432,18240, 3 18240,18384,17520,17688,17952,18144,15816,16008,15960,16056, 4 16296,16272,16464,16608,16560,16800,16728,16848,16800,16896, 5 16872,16752,18672,19152,19584,19680,19872,19800,27408,27888, 6 27912,27840,27576,27408,27288,27288,27408,27144,27144,27048, 7 27240,27096,27192,27120,27192,26976,26976,27312,27240,27288, 8 23904,23712,23592,23184,22872,20400,20544,21576,21984,15816, 9 22440,25176,25080,24888,24576,24240,24024,29880,29616,29568/ data (jx(i),i= 401, 500)/ + 29496,29496,29376,29328,29280,29208,29040,28992,31968,31944, 1 32280,32280,32640,33552,33624,34008,33984,34080,34008,34128, 2 34200,34272,34248,34368,34344,34464,34440,34536,34464,34632, 3 34656,34776,34872,34992,34992,30408,30792,31008,32448,33312, 4 33528,33648,33888,33936,34176,25320,25656,25944,26112,26352, 5 26496,26640,26712,26688,26784,26904,27024,27240,27288,27240, 6 27456,27480,27264,27264,27552,27480,27504,15792,14784,24864, 7 24864,19824,26112,26136,26064,26040,25872,25872,19416,19800, 8 29376,29448,27552,28656,29448,26208,26184,26568,14760,16056, 9 18600,18792,18936,19032,19200,31992,32064,32280,32424,32472/ data (jx(i),i= 501, 600)/ + 32424,33024,33000,35016,35352,35352,36024,36024,35016,27288, 1 26904,27288,27288,27384,11376,11376,11712,11712,14736,34992, 2 34344,34344,33984,34008,34488,34488,34344,34344,34152,34152, 3 34008,34008,33936,33936,33864,33864,33528,33528,33192,33192, 4 33024,33000,32496,25800,25680,25488,25416,25416,25296,25200, 5 25248,25152,25200,25128,25224,24960,24864,24792,24816,24696, 6 24672,24720,24576,24480,24336,24336,24216,24168,24048,24000, 7 23880,23832,23712,23904,23880,23784,23712,23664,23736,23640, 8 23568,23520,23328,23352,23160,23208,23160,23016,23064,22968, 9 22872,22920,22872,22584,22584,22464,22392,22224,22224,22032/ data (jx(i),i= 601, 700)/ + 21960,22008,21960,21768,21672,21600,21648,21456,21336,21384, 1 21240,21264,21216,19536,19560,19824,19800,19848,20016,20088, 2 19896,19992,20808,20808,27408,27360,27432,27816,27912,28056, 3 28344,28560,28872,28896,29088,29136,29208,29256,29448,29424, 4 29520,29520,29712,29808,29904,29976,29880,32976,32304,31968, 5 31968,31632,31632,31464,31464,31128,31128,30960,30960,30792, 6 30792,30624,30624,30552,30552,30288,30288,30024,30024,21216, 7 21312,21312,27144,27168,26376,26088,25896,25728,25584,25488, 8 24936,24936,21480,21984,23184,24264,24264,29856,30000,29880, 9 29880,30072,30096,29952,29712,29568,29400,29352,29208,29136/ data (jx(i),i= 701, 800)/ + 29112,16776,16776,15768,15744,14736,32520,30576,30576,30336, 1 30264,30120,30048,25272,25272,25368,25368,25440,25464,25368, 2 25368,25296,25296,25200,25200,25296,25296,25536,25536,11712, 3 12744,14760,14760,14424,14424,13680,21312,21144,21144,20808, 4 20808,20640,20640,20568,20568,20472,20472,20328,20304,19968, 5 19968,19464,19464,19128,19128,18960,18888,18960,19080,18936, 6 18816,18576,18432,18360,18312,18240,18216,18288,17784,17784, 7 17280,17280,32496,31968,31968,31800,31800,31632,31632,32640, 8 32640,32808,32808,33984,36888,36864,36792,36792,28992,28896, 9 28848,28776,28584,28560,28848,28944,28896,28608,28608,28752/ data (jx(i),i= 801, 900)/ + 28704,28464,28392,28512,28680,28680,28392,28272,28056,27984, 1 27960,27816,27768,27768,27672,27576,27432,27192,27024,26880, 2 25872,25872,26880,18336,18312,18192,18072,18096,18000,17952, 3 17832,17688,17568,17616,17616,17832,17832,18120,18120,18288, 4 18288,18408,18408,18456,18456,18624,18624,18792,18792,18960, 5 18960,19128,19128,19296,19296,19632,19632,20016,20016,20232, 6 20232,20616,25536,25296,25296,25128,25128,25032,25032,24864, 7 24864,24792,24768,24696,24696,24024,24024,23856,23856,23760, 8 23760,23592,23592,23256,23256,23184,23184,23016,23016,22752, 9 22680,22512,22512,22008,22008,21912,21840,21744,21744,21576/ data (jx(i),i= 901, 1000)/ + 21576,21360,21384,21168,21168,21024,21048,12744,31464,31464, 1 31296,31296,31128,31128,30792,30768,29424,29400,28896,28896, 2 34536,33960,33960,34032,34032,34008,34032,33960,25824,26904, 3 27144,27240,27192,27000,27000,27072,26688,26520,26352,26448, 4 26352,26184,26040,26088,26040,25848,25776,25680,25848,25896, 5 26088,26112,25872,30504,30840,30840,30504,20616,20664,20616, 6 21168,21168,21336,21336,21504,21504,21648,21672,21864,21864, 7 22032,20112,19848,19680,19680,19416,19320,19320,18936,18936, 8 18144,17832,17760,16752,25560,25536,25680,25608,22008,24648, 9 24576,24600,24744,24768,25080,25128,25272,25368,12720,12792/ data (jx(i),i= 1001, 1100)/ + 12768,13944,13896, 9456, 9504, 9744, 9696, 9768,10008,10104, 1 10344,10416,10488,10680,10800,10776,10848,10896,11472,11592, 2 11640,16776,20688,22080,22056,17760,17760,16752,16272,16248, 3 15744,15744,30552,30480,30552,30552,25512,25632,25560,25584, 4 25680,25824,25944,26208,26256,26112,26160,26304,26280,26448, 5 26592,26784,26952,27336,27336,27504,27504,37608,37560,36432, 6 36432,36336,36336,36192,36192,35688,35688,35592,35592,34848, 7 34848,34752,34752,34608,20688,26592,26592,27000,20712,21672, 8 23664,23664,34440,34440,33744,33768,33600,33624,16224,16728, 9 16728,19776,19920,20088,20064,20496,20400,20424,20568,20544/ data (jx(i),i= 1101, 1200)/ + 12528,12432,12312,12336,12096,12120,12288,12456,12408,12480, 1 12528,12648,12480,12288,12024,11880,11784,11280,10560,10320, 2 10320,10560,10584,10536,10704,10584,10512,10200, 9960,10272, 3 10320,10104,10032,10128,10368,10488,10392,10128, 9816, 9744, 4 9648, 9552, 9384, 9336, 9504, 9624, 9360, 9288, 9096, 9192, 5 9384, 9456, 9456, 9360, 9288, 9360, 9504, 9600, 9480,12912, 6 13056,13056,13176,13176,13224,13344,13368,13488,13680,13560, 7 13200,13056,12912,12960,13320,13512,13608,13536,13344,13176, 8 13104,13200,13104,13104,12984,12912,12816,12768,12792,12984, 9 13032,12840,12720,12720,12456,12384,12432,12576,12888,12888/ data (jx(i),i= 1201, 1300)/ + 12936,12888,12912,12864,12912,12744,12792,12888,12768,12744, 1 12888,12864,12960,12912,13008,13032,13224,13296,13392,13464, 2 13320,13320,13416,13464,13560,13512,13656,13752,13752,13584, 3 13680,13608,13560,13728,13800,13968,13992,12888,12816,12816, 4 12888,12792,12600,33456,33456,32616,32640,31776,30888,30888, 5 31056,31080,31800,31800,31656,31704,31536,31488,31608,31728, 6 31680,31776,31704,13752,13728,31920,32088,31872,31800,31848, 7 31728,31344,31104,30888,30840,30912,30816,30744,30696,30744, 8 30720,30576,30432,39840,39864,37176,37224,36192,36144,35592, 9 33960,33960,33624,26640,30288,30216,29904,29808,29520,29496/ data (jx(i),i= 1301, 1400)/ + 29400,29424,29208,29136,29016,28920,28800,28752,28896,28800, 1 28680,28584,28464,28416,28488,28416,28344,28320,28080,28056, 2 28104,27912,27840,27744,27672,27576,27648,27480,27672,27624, 3 27528,27216,26760,26784,26376,26160,26136,26016,25848,25872, 4 25776,21648,21840,21816,22008,22008,22152,22224,22176,22344, 5 22584,22416,22368,22416,22656,22752,22680,22632,22632,22728, 6 22848,22920,23088,22992,23016,23064,23352,23592,23640,23616, 7 23688,24000,24240,24264,24552,24648,24648,24840,24960,25032, 8 25296,25464,12648,12624,12528,12000,11616,11448,11328,11136, 9 10704,10632,10896,11208,11352,11424,11544,11640,11736,15720/ data (jx(i),i= 1401, 1500)/ + 30552,30528,30720,30768,30984,31152,31464,16248,36240,31656, 1 31656,31872,31944,31896,32112,16464,16560,16608,16896,16992, 2 17064,17064,17184,17232,17688,17712,17928,18024,18072,18384, 3 18288,18384,18528,18576,18672,19032,18984,19128,19104,21792, 4 36480,36312,39504,39504,39408,39408,38736,38616,38592,38304, 5 38184,37968,37968,30816,30816,30672,30720,25800,25944,25920, 6 25992,25848,25968,26184,26352,26448,26520,26664,26712,22344, 7 22344,22584,22272,22200,22488,22584,22728,22728,21768,11760, 8 11712,11592,11688,11640,11400,11328,11256,11160,10992,11016, 9 10872,10848,10968,10896,10968,10920,10992,10872,10824,10656/ data (jx(i),i= 1501, 1600)/ + 10536,10608,10368,10128, 9888, 9936, 9888, 9936,10128,10152, 1 10296,15888,15720,15744,14736,14736,14304,14304,30744,28824, 2 27744,27720,26304,26736,22272,22272,36480,36312,36312,36000, 3 35976,35904,35904,35496,35496,35328,35280,34752,34512,34152, 4 34080,32784,32784,31848,31848,31752,31752,31560,31560,31176, 5 31200,31056,31080,30912,30888,30864,10368,10272,10320,10224, 6 10152,10080, 9816, 9864,10008,10128,10344,10488,10440,10272, 7 10152, 9768, 9600, 9480, 9360, 9216, 9432, 9624, 9840, 9960, 8 10080,10128,21768,21792,14304,14592,14736,14760,15936,36504, 9 36504,38448,38592,38736,39024,21792,33240,33336,33480,33504/ data (jx(i),i= 1601, 1700)/ + 33600,33648,33984,34032,33936,33936,33984,32832,29304,28992, 1 29016,28920,28920,28776,28776,29256,29280,29280,29160,29112, 2 28968,28968,26328,26400,26400,26328,26328,26376,26376,26328, 3 26304,25800,25824,23808,22464,22416,22200,21768,32856,28848, 4 28824,29208,29208,29328,26808,26784,27024,27168,27408,27816, 5 27984,28008,28176,28176,36552,15936,15768,15768,15264,15288, 6 14712,14880,15312,15192,14928,14832,14712,14544,14400,14280, 7 14064,13896,13680,13704,13560,13416,13416,13392,13296,13296, 8 13104,13128,13296,21480,21264,21168,20736,20400,20256,20304, 9 20256,20136,19968,19560,19488,19320,19200,19200,18888,18768/ data (jx(i),i= 1701, 1800)/ + 18576,18144,17952,17688,17592,17448,17280,16896,16656,16632, 1 13584,13800,14328,14544,14760,15168,15432,15648,16008,33648, 2 33768,33840,34008,33984,34200,34632,35880,35880,36552,11472, 3 11448,11544,11616,11496,11352,11256,11064,11208,11400,11424, 4 11376,11472,11328,11304,11520,11664,11736,11856,11784,11808, 5 11880,12072,11928,11856,11736,11784,11688,11544,11496,11616, 6 11664,11856,11904,11808,11688,11856,11760,11616,11808,11784, 7 11904,11712,11544,11232,11184,10680,10584,10800,10752,10488, 8 10200,10248,10152,10104,10224,28872,23832,23832,23856,24120, 9 24192,24096,23832,13224,13008,12888,12840,13080,12936,12672/ data (jx(i),i= 1801, 1900)/ + 12600,12408,12552,12480,12576,12528,12384,12288,12432,12336, 1 12240,12192,12168,12072,11952,11880,11808,29640,29640,29544, 2 29544,23208,23112,23208,23136,23184,23136,23232,23184,23280, 3 23448,23520,23496,23664,23640,11520,11208,11040,10728, 7656, 4 7656,35880,14736,23856,23976,11280,11544,11472,11160,11304, 5 11472,11520,11520,11352,11424,11280,11352,11280,11448,11520, 6 11640,11664,11400,11328,11088,11040,11112,11400,11448,11040, 7 10992,10776,10680,10704,10656,10752,10776,10896,10656,10728, 8 10632,10584,15336,15528,15504,15360, 9792, 9648, 9600, 9624, 9 9672, 9888, 9912, 9480, 9408, 9480,15960,16248,16224,15528/ data (jx(i),i= 1901, 1960)/ + 15792, 8568, 8640, 8616, 8544, 8688, 8712, 8736, 8736, 8256, 1 8208, 8232, 8208, 8232, 8328, 7056, 6984, 7056, 6888, 6888, 2 6960, 7896, 7824, 7872, 7080, 7152, 7104, 9816, 9816, 9744, 3 9648,10224,10200,10296,10320, 9840, 9456, 8952, 7824, 7656, 4 9120, 9552, 9552, 9648, 9840, 9984,10056,10056, 7368, 7368, 5 7536, 0, 0, 0, 0, 0, 0, 0, 0, 0/ data (jy(i),i= 1, 100)/ + 60720,61848,61776,61776,57096,56472,60360,54072,50016,51048, 1 45000,44472,40152,38736,36000,33288,30864,28104,24648,23832, 2 19392,17088,15744,15648,15624,15624,12024,60840,61080,61944, 3 62016,61776,37320,37248,37248,37128,36984,36696,36360,36192, 4 33648,34248,34464,34608,34728,34800,34992,35040,35400,35688, 5 35880,35952,32664,32472,32328,32136,31944,31872,31752,31704, 6 31584,31368,31200,30888,30456,30336,30312,30240,30216,29976, 7 29904,29784,29688,29496,28944,28440,28320,27600,27336,27048, 8 26904,26688,26376,26304,26184,26112,25968,25512,25200,24864, 9 26160,24312,23904,23400,23232,23064,22464,22128,21816,21360/ data (jy(i),i= 101, 200)/ + 20832,20760,20424,20040,19608,19224,19104,18984,18456,18024, 1 17952,17712,17664,16968,16704,16440,15864,15672,15696, 9744, 2 9768, 9984,10200,10584,10656,11088,11184,11688,11952,12216, 3 12288,12432,12480,12744,13176,13536,13632,13776,13920,14088, 4 14160,15168,15360, 8688, 8664, 8736, 8808, 8784, 8688, 8856, 5 8856, 8712, 8592, 8472, 8496, 8784, 8688, 8640, 8616, 8472, 6 8352, 8424, 8376, 8160, 8112, 7992, 7848, 7944, 8112, 8160, 7 8040, 8184, 8088, 7896, 7752, 8400, 9072, 9360, 9456, 9504, 8 9600,10008,10440,10680,10776,10944,10896,11088,11088,11280, 9 11280,11040,10872,10752,11136,11592,11664,11760,11904,12120/ data (jy(i),i= 201, 300)/ + 12120, 6792, 6792, 6696, 6864, 6936, 6888, 6936, 7008, 7032, 1 7104, 7152, 7152, 6984, 6912, 6912, 7056, 7104, 7152, 7488, 2 7392, 7176, 7176, 7296, 7464, 7608, 7584, 7512, 7704, 7992, 3 8256, 8424, 8568, 8640, 8904, 8952, 8760, 8688, 8520, 8496, 4 8592, 8904, 9600, 9720,10440,10680,10944,11232,11448,11520, 5 11568,11664,11760,11808,11880,11880, 7080, 7272, 7344, 7320, 6 7392, 7488, 7488, 7608, 7680, 7704, 7560, 7608, 7512, 7560, 7 7632, 7824, 7920, 8016, 8064, 8160, 8232, 8328, 8016, 8160, 8 8160, 7848, 7656, 7776, 7728, 7872, 7680, 7872, 7848, 8064, 9 8040, 8112, 8160, 8304, 8280, 8064, 8064, 8376, 8400, 8136/ data (jy(i),i= 301, 400)/ + 8088, 8280, 8520, 8736, 8928, 9024, 9216, 9504, 9648,59544, 1 59712,59832,59904,59904,59784,59616,59232,58944,58824,58776, 2 58920,59640,60144,60360,60552,60840,60696,59736,59424,59328, 3 59256,59064,59592,59472,59112,58992,56784,57048,57240,57456, 4 57696,57792,57792,57888,58272,58512,58728,58896,58968,59136, 5 59184,59280,59136,58944,58632,58416,58224,53928,61608,61248, 6 61032,60984,61248,61128,61152,61056,60864,60648,60552,60408, 7 59904,59640,59448,59280,59064,58824,57624,57360,57336,57264, 8 57048,57120,57096,56712,56664,57816,57840,57096,56928,56880, 9 56736,56616,56880,57024,57000,56640,56712,55848,55728,55920/ data (jy(i),i= 401, 500)/ + 55992,56352,56352,56448,56400,56496,56400,57264,52680,53184, 1 53136,53376,53376,54360,54528,54696,54792,54864,54984,54960, 2 55152,55200,55272,55296,55368,55488,55824,55800,55992,56064, 3 56208,56184,56328,56376,56496,55896,56136,56424,57576,58968, 4 58992,59136,59616,59976,60240,56112,55896,55824,55920,55848, 5 56040,56064,56136,56328,56352,56256,56352,56208,56256,56304, 6 56544,56832,56856,56952,57024,57048,57192,53832,53880,54624, 7 53832,53760,53328,53496,53496,53664,53664,53808,50232,50208, 8 55272,52656,51480,52464,52488,53304,53208,52656,50352,50280, 9 50256,50232,50136,50352,50376,51648,51648,51816,51864,51816/ data (jy(i),i= 501, 600)/ + 51528,50256,50016,49728,49752,49920,49920,49968,49896,49440, 1 49560,50160,50880,51168,50952,50616,50616,49416,49344,49128, 2 49128,48960,48960,48624,48624,48120,48120,47952,47952,47784, 3 47784,46344,46344,46248,46200,46080,46104,45912,45936,45768, 4 45768,45240,45264,53328,53208,53232,53064,52968,52896,52824, 5 52776,52680,52608,52512,52464,52224,52008,52008,51888,51840, 6 51744,51696,51624,51744,51624,51552,51552,51360,51288,51336, 7 51240,51384,51240,50928,50904,51024,50976,50856,50736,50640, 8 50688,50520,50424,50304,50256,50184,50016,49968,49920,49968, 9 49872,49824,49752,49824,49632,49728,49728,49560,49488,49584/ data (jy(i),i= 601, 700)/ + 49536,49440,49344,49296,49344,49248,48840,48768,48816,48720, 1 48792,48672,48672,50040,49920,49824,49728,49560,49416,49224, 2 49152,49032,49032,48696,49560,49632,49800,50256,50040,49968, 3 50184,50208,49872,49488,49392,49200,49152,49224,49224,49032, 4 48912,48768,48744,48792,48768,48552,48528,49632,49752,49872, 5 49800,49800,49704,49704,49272,49272,49440,49440,49368,49368, 6 49128,49128,48960,48960,48792,48816,48648,48648,48528,48504, 7 48528,48168,49344,49008,48264,48264,48168,48168,48144,48024, 8 48024,47592,48168,48504,48480,47712,47616,48480,48360,48144, 9 48000,47712,47616,47592,47352,47304,47160,47016,47064,46992/ data (jy(i),i= 701, 800)/ + 46632,45720,46728,46776,48792,48864,45600,45576,45744,45744, 1 46272,46368,46584,47616,47256,47256,47088,47112,46344,46344, 2 46104,46104,45936,45936,45264,45264,45168,45168,45072,44856, 3 44856,45816,45480,45480,44808,44784,48000,48024,47832,47856, 4 47688,47688,47352,47352,47184,47184,46680,46680,46344,46344, 5 46008,46008,45840,45840,45672,45672,45432,45264,45168,45096, 6 44952,45000,44688,44736,44880,44880,44760,44688,44688,45552, 7 45552,45720,45096,44592,44496,44496,44328,44328,44232,44256, 8 44424,44424,44592,44568,44616,44688,44664,44568,46560,46608, 9 46776,46872,46824,46632,46440,46104,45936,45864,45768,45648/ data (jy(i),i= 801, 900)/ + 45624,45648,45576,45264,45168,45048,44832,44496,44424,44496, 1 44736,44688,44592,44424,44352,44040,43920,43512,43464,43296, 2 45072,44568,44592,44616,44520,44568,44520,44304,44304,44424, 3 44400,44448,44208,44208,44040,44016,43872,43848,43680,43656, 4 43584,43584,43512,43512,43344,43320,42984,42984,42912,42888, 5 42720,42720,42480,42480,42312,42288,42120,42120,42000,41952, 6 41616,41616,44832,44832,44736,44760,44664,44664,44568,44496, 7 44328,44328,44256,44232,44160,44160,43992,43992,43896,43848, 8 43728,43728,43560,43584,43368,43320,43248,43248,43008,43008, 9 42936,42936,42600,42600,42528,42528,42312,42312,42264,42192/ data (jy(i),i= 901, 1000)/ + 41976,41928,41856,41856,41760,41688,41592,40560,44568,44064, 1 44064,43392,43392,43224,43224,42216,42216,41040,41064,40008, 2 40104,40080,42384,42384,42504,42504,42576,42576,39960,43224, 3 43080,42888,42792,42744,42456,42168,41928,41952,42048,41784, 4 41568,41592,41424,41088,41040,41112,41040,40800,40704,40464, 5 40320,40224,40080,39840,39864,40032,40008,41184,41160,40680, 6 40584,39912,39912,39672,39672,39504,39504,39432,39432,39336, 7 39336,41616,41328,41328,41496,41496,41376,41328,40944,40272, 8 40296,39840,38832,38832,39024,39432,39552,39768,39600,39576, 9 39480,39384,39408,39312,39096,39000,39000,39048,39864,39864/ data (jy(i),i= 1001, 1100)/ + 39000,39000,38856,38784,38976,39144,39240,39288,39312,39240, 1 39240,39408,39408,39240,39288,39504,39432,39552,39744,39840, 2 40032,37320,36888,36864,39336,37320,36816,36816,36816,34800, 3 34800,32784,36768,36720,36648,35976,38808,38616,38448,37968, 4 37776,37704,37296,37200,37056,36840,36624,36624,36528,36384, 5 36336,36000,35952,35952,36096,36096,35928,35952,35784,35760, 6 35856,35832,35928,35928,35760,35760,35832,35832,35760,35736, 7 35808,35808,35736,35712,35640,35400,35904,35904,34872,34872, 8 35376,34872,35712,35544,35520,35352,35352,34176,30768,30792, 9 32808,32856,33240,33504,33840,34176,34416,34632,35016,35160/ data (jy(i),i= 1101, 1200)/ + 33384,33408,33528,33816,34080,34224,34296,34248,34152,34056, 1 34080,34248,34416,34440,34344,34536,34560,34296,34104,34416, 2 34632,34728,34896,35040,35232,35328,35232,35160,35328,35448, 3 35592,35856,36168,36288,36192,36408,36576,36456,36888,36888, 4 36792,36864,36840,36936,36960,37152,37368,37800,37944,38088, 5 38040,38088,38160,38208,38352,38400,38304,38376,38520,33312, 6 33480,33720,33744,33984,34008,33936,33792,33960,34032,34320, 7 34224,34296,34560,34608,34656,34560,34776,34848,34752,34800, 8 34944,35112,35304,35496,35496,35400,35400,35472,35688,35712, 9 35856,36096,36120,36024,36000,36048,36192,36168,36432,36552/ data (jy(i),i= 1201, 1300)/ + 36576,36672,36744,36792,36864,36984,37056,37056,37104,37200, 1 37224,37416,37368,37536,37488,37608,37584,37656,37608,37704, 2 37872,37944,37920,38040,37968,38040,38040,38184,38232,38232, 3 38328,38376,38520,38496,38640,38640,38760,33240,33264,33360, 4 33456,33528,33576,34176,34008,33984,33000,32976,35976,36144, 5 36144,35976,35976,35808,35664,35544,35352,35160,34896,34752, 6 34200,33552,33096,33288,32784,32832,32784,32736,32520,32472, 7 32448,32520,32352,32496,32472,32424,32328,32496,32400,32280, 8 32208,32208,32112,33288,32856,32736,32088,32064,33072,33072, 9 33864,34032,34008,31968,32136,31944,32016,32136,32064,32160/ data (jy(i),i= 1301, 1400)/ + 32136,32256,32160,32280,32328,32328,32232,32328,32400,32520, 1 32448,32712,32640,32664,32784,32904,32856,32952,33024,32952, 2 32880,32856,32712,32736,32880,32856,32640,32352,32088,31992, 3 32088,32232,32160,32016,31848,31896,31584,31440,31440,31272, 4 31200,34824,33840,33336,32952,32784,32496,32496,32256,32280, 5 32112,32040,31944,31800,31536,31320,31152,31080,30840,30768, 6 30792,30840,31152,31488,31752,31776,31728,31896,31848,31512, 7 31392,31296,31560,31800,31920,31848,31536,31392,31560,31560, 8 31488,31128,33312,33240,33192,33288,32928,32640,32544,32520, 9 32184,31608,31536,31224,31248,31104,31056,30888,30840,30768/ data (jy(i),i= 1401, 1500)/ + 31920,31704,31416,31152,31104,30720,30432,30288,30024,30432, 1 30792,30648,30528,30096,29952,30360,30312,30216,29976,30000, 2 29952,29880,29784,29616,29520,29400,29136,29112,28944,28848, 3 28632,28656,28560,28392,28320,28224,28080,27960,27768,27768, 4 29016,29016,28104,28176,28176,28104,28080,28176,28320,28512, 5 28824,28992,29064,30408,30144,30144,27864,31152,31152,31104, 6 31056,30720,30408,30384,30240,30216,30096,30048,27816,31008, 7 30792,30744,30504,30360,29976,29928,29952,29808,29808,30576, 8 30504,30552,30744,30792,30696,30792,30792,30672,30672,30528, 9 30552,30408,30216,30024,29784,29616,29400,28992,28920,28920/ data (jy(i),i= 1501, 1600)/ + 28440,28224,27984,28104,27912,27696,27456,27312,27192,27360, 1 27168,30336,30168,28752,28776,27312,27264,27144,26904,26856, 2 26832,27840,26832,26808,26784,27792,28944,28944,27504,27480, 3 27576,27576,27480,27480,27312,27312,27264,27264,26712,26256, 4 25920,25896,26904,26904,27072,27144,27816,27816,27912,27888, 5 27984,27888,27720,27720,27864,27816,27072,26976,26784,26736, 6 26856,26880,26688,26568,26592,26376,26544,26424,26352,26376, 7 26232,26208,25680,25608,25320,25176,24576,24408,24552,24504, 8 24240,23808,26784,25272,25800,25800,25728,25248,25272,22944, 9 23952,23928,24432,24480,24768,24360,22872,22944,23256,23760/ data (jy(i),i= 1601, 1700)/ + 23880,24264,24600,24816,25080,25488,25776,22872,22848,26448, 1 26016,25680,25368,25032,24816,24840,24528,23352,23328,23112, 2 22920,22824,25512,25416,25320,25296,24624,24624,24456,24456, 3 23256,23256,21744,21744,21744,21840,21984,22968,18000,20784, 4 21816,21816,21912,21912,21768,20904,21072,21312,21360,21816, 5 21816,20952,20952,20784,19416,25080,25080,23736,23736,22080, 6 21624,21456,21192,20688,20424,20400,20424,20592,20616,20592, 7 20400,20352,20016,19632,19536,19584,19824,19896,19896,19704, 8 19488,19272,19128,24288,24120,24120,23904,23952,23760,23496, 9 23424,23400,23208,23088,22560,22344,21840,21672,21120,20976/ data (jy(i),i= 1701, 1800)/ + 20928,20424,20304,20016,19824,19728,19344,19104,18840,18696, 1 19200,19176,18600,18504,18504,18360,18336,18504,18672,18000, 2 18144,18432,18624,18720,18936,18936,18624,19080,19080,17928, 3 18096,18144,18360,18528,18576,18696,18696,18840,18792,18840, 4 18984,19128,19248,19344,19344,19296,19632,19560,19512,19440, 5 19440,19872,20208,20232,20184,20040,20112,20040,20064,20136, 6 20280,20256,20328,20424,20400,20616,20736,20808,20880,21120, 7 21408,21528,21792,21864,21936,22200,22344,22320,22440,22464, 8 22680,22848,23040,23304,23424,17904,17904,19584,19680,19896, 9 20328,20712,20712,19080,19152,19392,19248,18888,18816,18912/ data (jy(i),i= 1801, 1900)/ + 18696,18696,18552,18456,18384,18312,18456,18360,18096,18072, 1 18192,18192,18000,18072,17976,18000,17856,17928,17808,17808, 2 17928,18696,18552,18264,18120,18048,18000,17952,17784,17904, 3 17928,17856,17832,17808,17880,17592,17592,17712,17376,17400, 4 17064,18072,15672,15672,15672,17424,17400,17112,16704,16488, 5 16488,16368,16272,16080,15936,15816,15624,15336,15288,15144, 6 15072,14976,14856,14640,14400,14280,14136,13920,13728,13536, 7 13296,13224,13128,13008,12816,12696,12432,12264,12024,11856, 8 11808,11952, 7272, 7320, 7344, 7344, 7296, 7392, 7320, 7272, 9 7296, 7152, 7176, 7584, 7560, 7536, 7584, 7656, 7608, 7416/ data (jy(i),i= 1901, 1960)/ + 7584, 7920, 7896, 7848, 7872, 7872, 7920, 7896, 7824, 7920, 1 7968, 8160, 8208, 8472, 8088, 8352, 8304, 8400, 8256, 8328, 2 8328, 8496, 8496, 8544, 6720, 6672, 6648, 6672, 6624, 6648, 3 6816, 6360, 6432, 6480, 6408, 6936, 7008, 7008, 6792, 6864, 4 7104, 7080, 7152, 7176, 7080, 7128, 7056, 6984, 8760, 8688, 5 8688, 0, 0, 0, 0, 0, 0, 0, 0, 0/ data (jb1(i),i= 1, 100)/ + 1, 28, 29, 30, 2, 2, 3, 7, 3, 1, 1 6, 5, 8, 6, 10, 9, 12, 13, 33, 34, 2 35, 36, 37, 38, 39, 15, 11, 16, 41, 42, 3 43, 44, 45, 46, 47, 48, 49, 50, 51, 15, 4 16, 53, 54, 55, 56, 57, 58, 59, 60, 61, 5 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 6 72, 73, 74, 75, 76, 18, 14, 18, 78, 79, 7 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 8 19, 17, 20, 19, 92, 93, 94, 95, 96, 97, 9 98, 99, 100, 101, 102, 103, 104, 21, 20, 21/ data (jb1(i),i= 101, 200)/ + 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 1 116, 23, 26, 118, 24, 23, 25, 120, 121, 122, 2 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 3 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 4 24, 22, 144, 145, 146, 147, 148, 149, 150, 151, 5 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 6 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 7 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 8 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 9 192, 193, 194, 195, 196, 197, 198, 199, 200, 27/ data (jb1(i),i= 201, 300)/ + 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 1 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 2 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 3 232, 233, 234, 235, 236, 237, 238, 239, 240, 241, 4 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 5 252, 253, 254, 255, 27, 202, 257, 258, 259, 260, 6 261, 262, 263, 264, 265, 266, 267, 268, 269, 270, 7 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 8 281, 282, 283, 284, 285, 286, 287, 288, 289, 290, 9 291, 292, 293, 294, 295, 296, 297, 298, 299, 300/ data (jb1(i),i= 301, 400)/ + 301, 302, 303, 304, 305, 306, 307, 308, 120, 310, 1 311, 312, 313, 314, 315, 316, 317, 318, 319, 320, 2 321, 322, 323, 324, 325, 326, 327, 2, 328, 329, 3 330, 331, 310, 333, 334, 335, 332, 337, 338, 339, 4 340, 341, 342, 343, 344, 345, 346, 347, 348, 349, 5 350, 351, 310, 332, 353, 354, 355, 356, 357, 3, 6 359, 360, 361, 362, 363, 364, 365, 366, 367, 368, 7 369, 370, 371, 372, 373, 374, 375, 376, 377, 378, 8 379, 381, 382, 383, 384, 357, 386, 387, 388, 385, 9 337, 390, 385, 389, 392, 393, 394, 395, 396, 381/ data (jb1(i),i= 401, 500)/ + 398, 399, 400, 401, 402, 403, 404, 405, 406, 407, 1 380, 409, 410, 411, 412, 413, 414, 415, 416, 417, 2 418, 419, 420, 421, 422, 423, 424, 425, 426, 427, 3 428, 429, 430, 431, 432, 433, 434, 435, 398, 436, 4 437, 438, 439, 440, 441, 442, 443, 444, 445, 392, 5 446, 447, 448, 449, 450, 451, 452, 453, 454, 455, 6 456, 457, 458, 459, 460, 461, 462, 463, 464, 465, 7 466, 380, 337, 469, 468, 392, 470, 358, 468, 358, 8 473, 474, 475, 476, 477, 471, 479, 472, 398, 481, 9 409, 483, 484, 482, 473, 486, 487, 483, 469, 489/ data (jb1(i),i= 501, 600)/ + 490, 491, 492, 493, 494, 479, 409, 496, 497, 498, 1 499, 500, 501, 502, 504, 505, 506, 507, 508, 503, 2 504, 510, 511, 512, 513, 483, 10, 515, 516, 517, 3 489, 518, 504, 520, 521, 522, 523, 524, 525, 526, 4 527, 528, 529, 530, 531, 532, 533, 534, 535, 536, 5 537, 538, 539, 540, 541, 542, 473, 544, 545, 546, 6 547, 548, 549, 550, 551, 552, 553, 554, 555, 556, 7 557, 558, 559, 560, 561, 562, 563, 564, 565, 566, 8 567, 568, 569, 570, 571, 572, 573, 574, 575, 576, 9 577, 578, 579, 580, 581, 582, 583, 584, 585, 586/ data (jb1(i),i= 601, 700)/ + 587, 588, 589, 590, 591, 592, 593, 594, 595, 596, 1 597, 598, 599, 600, 601, 602, 603, 604, 605, 606, 2 607, 608, 609, 610, 611, 612, 479, 614, 615, 616, 3 617, 618, 619, 620, 621, 622, 623, 613, 510, 625, 4 626, 627, 628, 629, 630, 631, 632, 633, 634, 635, 5 636, 637, 638, 639, 640, 641, 642, 643, 644, 645, 6 646, 503, 648, 649, 650, 651, 652, 653, 654, 655, 7 656, 657, 658, 659, 660, 661, 662, 663, 664, 665, 8 666, 667, 668, 647, 613, 670, 671, 510, 673, 674, 9 675, 676, 677, 678, 679, 680, 681, 672, 683, 684/ data (jb1(i),i= 701, 800)/ + 685, 686, 682, 647, 688, 689, 690, 691, 692, 693, 1 694, 695, 696, 697, 698, 699, 700, 702, 703, 704, 2 705, 519, 543, 707, 708, 709, 710, 711, 712, 701, 3 682, 714, 715, 716, 717, 718, 719, 720, 721, 722, 4 723, 724, 725, 726, 727, 728, 730, 518, 730, 702, 5 732, 733, 734, 735, 731, 672, 737, 738, 739, 740, 6 741, 742, 743, 744, 745, 746, 747, 748, 749, 750, 7 751, 752, 753, 754, 755, 756, 757, 758, 759, 760, 8 761, 762, 763, 764, 765, 766, 767, 768, 769, 770, 9 771, 702, 543, 773, 774, 775, 776, 777, 778, 779/ data (jb1(i),i= 801, 900)/ + 780, 781, 782, 783, 12, 785, 786, 787, 784, 701, 1 789, 790, 791, 792, 793, 794, 795, 796, 797, 798, 2 799, 800, 801, 802, 803, 804, 805, 806, 807, 808, 3 809, 810, 811, 812, 813, 814, 815, 816, 817, 818, 4 819, 729, 821, 822, 820, 768, 824, 825, 826, 827, 5 828, 829, 830, 831, 832, 833, 834, 835, 836, 837, 6 838, 839, 840, 841, 842, 843, 844, 845, 846, 847, 7 848, 849, 850, 851, 852, 853, 854, 855, 856, 857, 8 858, 859, 860, 861, 729, 863, 864, 865, 866, 867, 9 868, 869, 870, 871, 872, 873, 874, 875, 876, 877/ data (jb1(i),i= 901, 1000)/ + 878, 879, 880, 881, 882, 883, 884, 885, 886, 887, 1 888, 889, 890, 891, 892, 893, 894, 895, 896, 897, 2 898, 899, 900, 901, 902, 903, 904, 905, 906, 862, 3 731, 774, 909, 910, 911, 912, 913, 914, 915, 916, 4 917, 918, 919, 921, 921, 922, 923, 924, 925, 926, 5 927, 784, 920, 820, 930, 931, 932, 933, 934, 935, 6 936, 937, 938, 939, 940, 941, 942, 943, 944, 945, 7 946, 947, 948, 949, 950, 951, 952, 929, 954, 955, 8 922, 920, 954, 862, 958, 959, 960, 961, 962, 963, 9 964, 965, 966, 967, 968, 969, 970, 861, 972, 973/ data (jb1(i),i= 1001, 1100)/ + 974, 975, 976, 977, 978, 979, 980, 981, 982, 983, 1 985, 986, 987, 929, 971, 989, 990, 991, 992, 993, 2 994, 995, 996, 997, 985, 908, 999, 1000, 1001, 1002, 3 984, 1004, 1005, 1006, 1007, 1008, 1009, 1010, 1011, 1012, 4 1013, 1014, 1015, 1016, 1017, 1018, 1019, 1020, 908, 1004, 5 984, 1023, 1024, 971, 1022, 1026, 1023, 1022, 1028, 1029, 6 1030, 1031, 954, 1033, 1034, 1035, 985, 1037, 1038, 1039, 7 1040, 1041, 1042, 1043, 1044, 1045, 1046, 1047, 1048, 1049, 8 1050, 1051, 1052, 1053, 1054, 1055, 1056, 1036, 15, 1058, 9 1059, 1060, 1061, 1062, 1063, 1064, 1065, 1066, 1067, 1068/ data (jb1(i),i= 1101, 1200)/ + 1069, 1070, 1071, 1072, 1073, 921, 1023, 1076, 1077, 1053, 1 1075, 1079, 1076, 1081, 1080, 1074, 1083, 1084, 1085, 1086, 2 1087, 1089, 1090, 1091, 1092, 1093, 1094, 1095, 1096, 1097, 3 1098, 1099, 1075, 1101, 1102, 1103, 1104, 1105, 1106, 1107, 4 1108, 1109, 1110, 1111, 1112, 1113, 1114, 1115, 1116, 1117, 5 1118, 1119, 1120, 1121, 1122, 1123, 1124, 1125, 1126, 1127, 6 1128, 1129, 1130, 1131, 1132, 1133, 1134, 1135, 1136, 1137, 7 1138, 1139, 1140, 1141, 1142, 1143, 1144, 1145, 1146, 1147, 8 1148, 1149, 1150, 1151, 1152, 1153, 1154, 1155, 1156, 1157, 9 1158, 1004, 1160, 1161, 1162, 1163, 1164, 1165, 1166, 1167/ data (jb1(i),i= 1201, 1300)/ + 1168, 1169, 1170, 1171, 1172, 1173, 1174, 1175, 1176, 1177, 1 1178, 1179, 1180, 1181, 1182, 1183, 1184, 1185, 1186, 1187, 2 1188, 1189, 1190, 1191, 1192, 1193, 1194, 1195, 1196, 1197, 3 1198, 1199, 1200, 1201, 1202, 1203, 1204, 1205, 1206, 1207, 4 1208, 1209, 1210, 1211, 1212, 1213, 1214, 1215, 1216, 1217, 5 1218, 1219, 1220, 1221, 1222, 1223, 1224, 1225, 1226, 1227, 6 1228, 1229, 1230, 1231, 1232, 1233, 1234, 1235, 1236, 1003, 7 1160, 1238, 1239, 1240, 1241, 1242, 1101, 1088, 1244, 1245, 8 1246, 1247, 1036, 1249, 1250, 1251, 1252, 1253, 1254, 1255, 9 1256, 1257, 1258, 1259, 1260, 1261, 1262, 1248, 1160, 1264/ data (jb1(i),i= 1301, 1400)/ + 1032, 1248, 1266, 1267, 1268, 1269, 1270, 1271, 1272, 1273, 1 1274, 1275, 1276, 1277, 1278, 1279, 1280, 1281, 1282, 16, 2 1284, 1285, 1286, 1287, 1288, 1289, 1290, 1291, 1292, 1088, 3 1076, 1283, 1295, 1296, 1297, 1298, 1299, 1300, 1301, 1302, 4 1303, 1304, 1305, 1306, 1307, 1308, 1309, 1310, 1311, 1312, 5 1313, 1314, 1315, 1316, 1317, 1318, 1319, 1320, 1321, 1322, 6 1323, 1324, 1325, 1326, 1327, 1328, 1329, 1330, 1331, 1332, 7 1333, 1294, 1294, 1335, 1336, 1337, 1338, 1339, 1340, 1080, 8 1342, 1343, 1344, 1345, 1346, 1347, 1348, 1349, 1350, 1351, 9 1352, 1353, 1354, 1355, 1356, 1357, 1358, 1359, 1360, 1361/ data (jb1(i),i= 1401, 1500)/ + 1362, 1363, 1364, 1365, 1366, 1367, 1368, 1369, 1370, 1371, 1 1372, 1373, 1374, 1375, 1376, 1377, 1378, 1379, 1380, 1381, 2 1341, 1101, 1383, 1384, 1385, 1386, 1387, 1388, 1389, 1390, 3 1391, 1392, 1393, 1394, 1395, 1396, 1397, 1398, 1399, 1089, 4 1032, 1283, 1401, 1402, 1403, 1404, 1405, 1406, 1089, 1288, 5 1407, 1410, 1411, 1412, 1413, 1414, 1409, 1408, 1416, 1417, 6 1418, 1419, 1420, 1421, 1422, 1423, 1424, 1425, 1426, 1427, 7 1428, 1429, 1430, 1431, 1432, 1433, 1434, 1435, 1436, 1437, 8 1438, 1439, 1441, 1409, 18, 1443, 1444, 1445, 1446, 1447, 9 1448, 1449, 1450, 1451, 1452, 1441, 1407, 1454, 1455, 1456/ data (jb1(i),i= 1501, 1600)/ + 1341, 1458, 1459, 1460, 1461, 1462, 1463, 1464, 1465, 1466, 1 1467, 1468, 1357, 1470, 1471, 1472, 1473, 1474, 1475, 1476, 2 1477, 1478, 1440, 1399, 1480, 1481, 1482, 1483, 1484, 1485, 3 1486, 1487, 1488, 1489, 1490, 1491, 1492, 1493, 1494, 1495, 4 1496, 1497, 1498, 1499, 1500, 1501, 1502, 1503, 1504, 1505, 5 1506, 1507, 1508, 1509, 1510, 1408, 1512, 1513, 1514, 1515, 6 1516, 1517, 1511, 1457, 1519, 1520, 1521, 1469, 1523, 1469, 7 1525, 1440, 1523, 1441, 1527, 1528, 1529, 1530, 1531, 1532, 8 1533, 1534, 1535, 1536, 1537, 1538, 1539, 1540, 1541, 1542, 9 1543, 1544, 1545, 1546, 1547, 1548, 1549, 1550, 1551, 1552/ data (jb1(i),i= 1601, 1700)/ + 1553, 1554, 1555, 1457, 1511, 1557, 1558, 1559, 1560, 1561, 1 1562, 1563, 1564, 1565, 1566, 1567, 1568, 1569, 1570, 1571, 2 1572, 1573, 1574, 1575, 1576, 1577, 1578, 1579, 1580, 1581, 3 1525, 1583, 1518, 1585, 1586, 1587, 1588, 1584, 1590, 1591, 4 1592, 1593, 1594, 1595, 1584, 1582, 1590, 1597, 1598, 1599, 5 1600, 1601, 1602, 1603, 1604, 1605, 1606, 1541, 1597, 1608, 6 1520, 1610, 1611, 1612, 1613, 1614, 1615, 1616, 1617, 1618, 7 1619, 1620, 1621, 1609, 1523, 1623, 1624, 1625, 1626, 1627, 8 1628, 1629, 1630, 1631, 1632, 1633, 1634, 1635, 1636, 1637, 9 1596, 1608, 1640, 1641, 1642, 1643, 1609, 1633, 1645, 1646/ data (jb1(i),i= 1701, 1800)/ + 1647, 1648, 1649, 1650, 1651, 1652, 1653, 1640, 1590, 1655, 1 1589, 1656, 1657, 1658, 1659, 1660, 1661, 1662, 1663, 1664, 2 1665, 1666, 1667, 1668, 1669, 1670, 1671, 1672, 1673, 1674, 3 1675, 1676, 1677, 1678, 1679, 1680, 1681, 1682, 1596, 1684, 4 1685, 1686, 1687, 1688, 1689, 1690, 1691, 1692, 1693, 1694, 5 1695, 1696, 1697, 1698, 1699, 1700, 1701, 1702, 1703, 1704, 6 1705, 1706, 1707, 1708, 1709, 1683, 1711, 1712, 1713, 1714, 7 1715, 1716, 1717, 1718, 1710, 1720, 1721, 1722, 1723, 1724, 8 1725, 1726, 1727, 1728, 1655, 1639, 1730, 1731, 1732, 1733, 9 1734, 1735, 1736, 1737, 1738, 1739, 1740, 1741, 1742, 1743/ data (jb1(i),i= 1801, 1900)/ + 1744, 1745, 1746, 1747, 1748, 1749, 1750, 1751, 1752, 1753, 1 1754, 1755, 1756, 1757, 1758, 1759, 1760, 1761, 1762, 1763, 2 1764, 1765, 1766, 1767, 1768, 1769, 1770, 1771, 1772, 1773, 3 1774, 1775, 1776, 1777, 1778, 1779, 1780, 1781, 1782, 1783, 4 1784, 1582, 1640, 1787, 1788, 1789, 1790, 1791, 1792, 1634, 5 1683, 1794, 1795, 1796, 1797, 1798, 1799, 1800, 1801, 1802, 6 1803, 1804, 1805, 1806, 1807, 1808, 1809, 1810, 1811, 1812, 7 1813, 1814, 1815, 1816, 1730, 1639, 1818, 1819, 1820, 1786, 8 1710, 1822, 1823, 1824, 1825, 1826, 1827, 1828, 1829, 1830, 9 1831, 1832, 1833, 1834, 1787, 1730, 1836, 1837, 1838, 1839/ data (jb1(i),i= 1901, 2000)/ + 1840, 1841, 23, 1720, 1715, 1843, 1786, 1787, 1844, 1845, 1 1837, 1846, 1847, 1848, 1849, 1850, 1851, 1852, 1853, 1854, 2 1855, 1856, 1857, 1858, 1859, 1860, 1861, 1862, 1863, 1864, 3 1865, 1866, 1867, 1868, 1869, 1870, 1871, 1872, 1873, 1874, 4 1875, 1876, 1877, 1878, 1, 2, 3, 1, 5, 6, 5 7, 5, 9, 10, 11, 12, 13, 14, 9, 16, 6 17, 16, 19, 20, 21, 22, 19, 24, 25, 26, 7 24, 28, 29, 30, 28, 32, 33, 34, 35, 36, 8 32, 38, 39, 38, 41, 42, 41, 44, 45, 44, 9 47, 48, 47, 50, 51, 52, 50, 54, 55, 56/ data (jb1(i),i= 2001, 2020)/ + 54, 58, 59, 60, 61, 62, 63, 64, 65, 66, 1 67, 68, 69, 58, 71, 72, 71, 0, 0, 0/ data (jb2(i),i= 1, 100)/ + 28, 29, 30, 31, 31, 4, 4, 32, 32, 5, 1 7, 8, 10, 9, 11, 12, 13, 33, 34, 35, 2 36, 37, 38, 39, 40, 40, 14, 41, 42, 43, 3 44, 45, 46, 47, 48, 49, 50, 51, 52, 52, 4 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 5 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 6 73, 74, 75, 76, 77, 77, 17, 78, 79, 80, 7 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 8 90, 91, 91, 92, 93, 94, 95, 96, 97, 98, 9 99, 100, 101, 102, 103, 104, 105, 105, 22, 106/ data (jb2(i),i= 101, 200)/ + 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 1 117, 117, 118, 119, 119, 25, 26, 121, 122, 123, 2 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 3 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 4 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 5 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 6 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 7 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 8 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 9 193, 194, 195, 196, 197, 198, 199, 200, 201, 201/ data (jb2(i),i= 201, 300)/ + 203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 1 213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 2 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 3 233, 234, 235, 236, 237, 238, 239, 240, 241, 242, 4 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 5 253, 254, 255, 256, 256, 257, 258, 259, 260, 261, 6 262, 263, 264, 265, 266, 267, 268, 269, 270, 271, 7 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 8 282, 283, 284, 285, 286, 287, 288, 289, 290, 291, 9 292, 293, 294, 295, 296, 297, 298, 299, 300, 301/ data (jb2(i),i= 301, 400)/ + 302, 303, 304, 305, 306, 307, 308, 309, 309, 311, 1 312, 313, 314, 315, 316, 317, 318, 319, 320, 321, 2 322, 323, 324, 325, 326, 327, 1, 328, 329, 330, 3 331, 332, 333, 334, 335, 336, 336, 338, 339, 340, 4 341, 342, 343, 344, 345, 346, 347, 348, 349, 350, 5 351, 352, 352, 353, 354, 355, 356, 357, 358, 359, 6 360, 361, 362, 363, 364, 365, 366, 367, 368, 369, 7 370, 371, 372, 373, 374, 375, 376, 377, 378, 379, 8 380, 382, 383, 384, 385, 386, 387, 388, 389, 4, 9 390, 5, 391, 391, 393, 394, 395, 396, 397, 397/ data (jb2(i),i= 401, 500)/ + 399, 400, 401, 402, 403, 404, 405, 406, 407, 408, 1 408, 410, 411, 412, 413, 414, 415, 416, 417, 418, 2 419, 420, 421, 422, 423, 424, 425, 426, 427, 428, 3 429, 430, 431, 432, 433, 434, 435, 6, 436, 437, 4 438, 439, 440, 441, 442, 443, 444, 445, 7, 446, 5 447, 448, 449, 450, 451, 452, 453, 454, 455, 456, 6 457, 458, 459, 460, 461, 462, 463, 464, 465, 466, 7 467, 467, 468, 8, 469, 470, 471, 471, 472, 472, 8 474, 475, 476, 477, 478, 478, 480, 480, 481, 482, 9 482, 484, 485, 485, 486, 487, 488, 488, 489, 490/ data (jb2(i),i= 501, 600)/ + 491, 492, 493, 494, 495, 495, 496, 497, 498, 499, 1 500, 501, 502, 503, 505, 506, 507, 508, 9, 509, 2 509, 511, 512, 513, 514, 514, 515, 516, 517, 518, 3 519, 519, 520, 521, 522, 523, 524, 525, 526, 527, 4 528, 529, 530, 531, 532, 533, 534, 535, 536, 537, 5 538, 539, 540, 541, 542, 543, 544, 545, 546, 547, 6 548, 549, 550, 551, 552, 553, 554, 555, 556, 557, 7 558, 559, 560, 561, 562, 563, 564, 565, 566, 567, 8 568, 569, 570, 571, 572, 573, 574, 575, 576, 577, 9 578, 579, 580, 581, 582, 583, 584, 585, 586, 587/ data (jb2(i),i= 601, 700)/ + 588, 589, 590, 591, 592, 593, 594, 595, 596, 597, 1 598, 599, 600, 601, 602, 603, 604, 605, 606, 607, 2 608, 609, 610, 611, 612, 613, 614, 615, 616, 617, 3 618, 619, 620, 621, 622, 623, 624, 624, 625, 626, 4 627, 628, 629, 630, 631, 632, 633, 634, 635, 636, 5 637, 638, 639, 640, 641, 642, 643, 644, 645, 646, 6 647, 648, 649, 650, 651, 652, 653, 654, 655, 656, 7 657, 658, 659, 660, 661, 662, 663, 664, 665, 666, 8 667, 668, 669, 669, 670, 671, 672, 673, 674, 675, 9 676, 677, 678, 679, 680, 681, 682, 683, 684, 685/ data (jb2(i),i= 701, 800)/ + 686, 687, 687, 688, 689, 690, 691, 692, 693, 694, 1 695, 696, 697, 698, 699, 700, 701, 703, 704, 705, 2 706, 706, 707, 708, 709, 710, 711, 712, 713, 713, 3 714, 715, 716, 717, 718, 719, 720, 721, 722, 723, 4 724, 725, 726, 727, 728, 729, 731, 730, 11, 732, 5 733, 734, 735, 736, 736, 737, 738, 739, 740, 741, 6 742, 743, 744, 745, 746, 747, 748, 749, 750, 751, 7 752, 753, 754, 755, 756, 757, 758, 759, 760, 761, 8 762, 763, 764, 765, 766, 767, 768, 769, 770, 771, 9 772, 772, 773, 774, 775, 776, 777, 778, 779, 780/ data (jb2(i),i= 801, 900)/ + 781, 782, 783, 784, 785, 786, 787, 788, 788, 789, 1 790, 791, 792, 793, 794, 795, 796, 797, 798, 799, 2 800, 801, 802, 803, 804, 805, 806, 807, 808, 809, 3 810, 811, 812, 813, 814, 815, 816, 817, 818, 819, 4 820, 821, 822, 823, 823, 824, 825, 826, 827, 828, 5 829, 830, 831, 832, 833, 834, 835, 836, 837, 838, 6 839, 840, 841, 842, 843, 844, 845, 846, 847, 848, 7 849, 850, 851, 852, 853, 854, 855, 856, 857, 858, 8 859, 860, 861, 862, 863, 864, 865, 866, 867, 868, 9 869, 870, 871, 872, 873, 874, 875, 876, 877, 878/ data (jb2(i),i= 901, 1000)/ + 879, 880, 881, 882, 883, 884, 885, 886, 887, 888, 1 889, 890, 891, 892, 893, 894, 895, 896, 897, 898, 2 899, 900, 901, 902, 903, 904, 905, 906, 907, 907, 3 908, 909, 910, 911, 912, 913, 914, 915, 916, 917, 4 918, 919, 920, 13, 922, 923, 924, 925, 926, 927, 5 928, 928, 929, 930, 931, 932, 933, 934, 935, 936, 6 937, 938, 939, 940, 941, 942, 943, 944, 945, 946, 7 947, 948, 949, 950, 951, 952, 953, 953, 955, 956, 8 956, 957, 957, 958, 959, 960, 961, 962, 963, 964, 9 965, 966, 967, 968, 969, 970, 971, 972, 973, 974/ data (jb2(i),i= 1001, 1100)/ + 975, 976, 977, 978, 979, 980, 981, 982, 983, 984, 1 986, 987, 988, 988, 989, 990, 991, 992, 993, 994, 2 995, 996, 997, 998, 998, 999, 1000, 1001, 1002, 1003, 3 1003, 1005, 1006, 1007, 1008, 1009, 1010, 1011, 1012, 1013, 4 1014, 1015, 1016, 1017, 1018, 1019, 1020, 1021, 1021, 14, 5 1022, 1024, 1025, 1025, 1026, 1027, 1027, 1028, 1029, 1030, 6 1031, 1032, 1033, 1034, 1035, 1036, 1037, 1038, 1039, 1040, 7 1041, 1042, 1043, 1044, 1045, 1046, 1047, 1048, 1049, 1050, 8 1051, 1052, 1053, 1054, 1055, 1056, 1057, 1057, 1058, 1059, 9 1060, 1061, 1062, 1063, 1064, 1065, 1066, 1067, 1068, 1069/ data (jb2(i),i= 1101, 1200)/ + 1070, 1071, 1072, 1073, 1074, 1074, 1075, 1077, 1078, 1078, 1 1079, 1080, 1081, 1082, 1082, 1083, 1084, 1085, 1086, 1087, 2 1088, 1090, 1091, 1092, 1093, 1094, 1095, 1096, 1097, 1098, 3 1099, 1100, 1100, 1102, 1103, 1104, 1105, 1106, 1107, 1108, 4 1109, 1110, 1111, 1112, 1113, 1114, 1115, 1116, 1117, 1118, 5 1119, 1120, 1121, 1122, 1123, 1124, 1125, 1126, 1127, 1128, 6 1129, 1130, 1131, 1132, 1133, 1134, 1135, 1136, 1137, 1138, 7 1139, 1140, 1141, 1142, 1143, 1144, 1145, 1146, 1147, 1148, 8 1149, 1150, 1151, 1152, 1153, 1154, 1155, 1156, 1157, 1158, 9 1159, 1159, 1161, 1162, 1163, 1164, 1165, 1166, 1167, 1168/ data (jb2(i),i= 1201, 1300)/ + 1169, 1170, 1171, 1172, 1173, 1174, 1175, 1176, 1177, 1178, 1 1179, 1180, 1181, 1182, 1183, 1184, 1185, 1186, 1187, 1188, 2 1189, 1190, 1191, 1192, 1193, 1194, 1195, 1196, 1197, 1198, 3 1199, 1200, 1201, 1202, 1203, 1204, 1205, 1206, 1207, 1208, 4 1209, 1210, 1211, 1212, 1213, 1214, 1215, 1216, 1217, 1218, 5 1219, 1220, 1221, 1222, 1223, 1224, 1225, 1226, 1227, 1228, 6 1229, 1230, 1231, 1232, 1233, 1234, 1235, 1236, 1237, 1237, 7 1238, 1239, 1240, 1241, 1242, 1243, 1243, 1244, 1245, 1246, 8 1247, 1248, 1249, 1250, 1251, 1252, 1253, 1254, 1255, 1256, 9 1257, 1258, 1259, 1260, 1261, 1262, 1263, 1263, 1264, 1265/ data (jb2(i),i= 1301, 1400)/ + 1265, 1266, 1267, 1268, 1269, 1270, 1271, 1272, 1273, 1274, 1 1275, 1276, 1277, 1278, 1279, 1280, 1281, 1282, 1283, 1284, 2 1285, 1286, 1287, 1288, 1289, 1290, 1291, 1292, 1293, 1293, 3 1294, 1295, 1296, 1297, 1298, 1299, 1300, 1301, 1302, 1303, 4 1304, 1305, 1306, 1307, 1308, 1309, 1310, 1311, 1312, 1313, 5 1314, 1315, 1316, 1317, 1318, 1319, 1320, 1321, 1322, 1323, 6 1324, 1325, 1326, 1327, 1328, 1329, 1330, 1331, 1332, 1333, 7 1334, 1334, 1335, 1336, 1337, 1338, 1339, 1340, 1341, 1342, 8 1343, 1344, 1345, 1346, 1347, 1348, 1349, 1350, 1351, 1352, 9 1353, 1354, 1355, 1356, 1357, 1358, 1359, 1360, 1361, 1362/ data (jb2(i),i= 1401, 1500)/ + 1363, 1364, 1365, 1366, 1367, 1368, 1369, 1370, 1371, 1372, 1 1373, 1374, 1375, 1376, 1377, 1378, 1379, 1380, 1381, 1382, 2 1382, 1383, 1384, 1385, 1386, 1387, 1388, 1389, 1390, 1391, 3 1392, 1393, 1394, 1395, 1396, 1397, 1398, 1399, 17, 1400, 4 1400, 1401, 1402, 1403, 1404, 1405, 1406, 1407, 1408, 1409, 5 1410, 1411, 1412, 1413, 1414, 1415, 1415, 1416, 1417, 1418, 6 1419, 1420, 1421, 1422, 1423, 1424, 1425, 1426, 1427, 1428, 7 1429, 1430, 1431, 1432, 1433, 1434, 1435, 1436, 1437, 1438, 8 1439, 1440, 1442, 1442, 1443, 1444, 1445, 1446, 1447, 1448, 9 1449, 1450, 1451, 1452, 1453, 1453, 1454, 1455, 1456, 1457/ data (jb2(i),i= 1501, 1600)/ + 1458, 1459, 1460, 1461, 1462, 1463, 1464, 1465, 1466, 1467, 1 1468, 1469, 1470, 1471, 1472, 1473, 1474, 1475, 1476, 1477, 2 1478, 1479, 1479, 1480, 1481, 1482, 1483, 1484, 1485, 1486, 3 1487, 1488, 1489, 1490, 1491, 1492, 1493, 1494, 1495, 1496, 4 1497, 1498, 1499, 1500, 1501, 1502, 1503, 1504, 1505, 1506, 5 1507, 1508, 1509, 1510, 1511, 1512, 1513, 1514, 1515, 1516, 6 1517, 1518, 1518, 1519, 1520, 1521, 1522, 1522, 1524, 1524, 7 1526, 1526, 1525, 1527, 1528, 1529, 1530, 1531, 1532, 1533, 8 1534, 1535, 1536, 1537, 1538, 1539, 1540, 1541, 1542, 1543, 9 1544, 1545, 1546, 1547, 1548, 1549, 1550, 1551, 1552, 1553/ data (jb2(i),i= 1601, 1700)/ + 1554, 1555, 1556, 1556, 1557, 1558, 1559, 1560, 1561, 1562, 1 1563, 1564, 1565, 1566, 1567, 1568, 1569, 1570, 1571, 1572, 2 1573, 1574, 1575, 1576, 1577, 1578, 1579, 1580, 1581, 1582, 3 1583, 1584, 1585, 1586, 1587, 1588, 1589, 1589, 1591, 1592, 4 1593, 1594, 1595, 19, 1596, 20, 1597, 1598, 1599, 1600, 5 1601, 1602, 1603, 1604, 1605, 1606, 1607, 1607, 1608, 1609, 6 1610, 1611, 1612, 1613, 1614, 1615, 1616, 1617, 1618, 1619, 7 1620, 1621, 1622, 1622, 1623, 1624, 1625, 1626, 1627, 1628, 8 1629, 1630, 1631, 1632, 1633, 1634, 1635, 1636, 1637, 1638, 9 1638, 1639, 1641, 1642, 1643, 1644, 1644, 1645, 1646, 1647/ data (jb2(i),i= 1701, 1800)/ + 1648, 1649, 1650, 1651, 1652, 1653, 1654, 1654, 1655, 21, 1 1656, 1657, 1658, 1659, 1660, 1661, 1662, 1663, 1664, 1665, 2 1666, 1667, 1668, 1669, 1670, 1671, 1672, 1673, 1674, 1675, 3 1676, 1677, 1678, 1679, 1680, 1681, 1682, 1683, 1684, 1685, 4 1686, 1687, 1688, 1689, 1690, 1691, 1692, 1693, 1694, 1695, 5 1696, 1697, 1698, 1699, 1700, 1701, 1702, 1703, 1704, 1705, 6 1706, 1707, 1708, 1709, 1710, 1711, 1712, 1713, 1714, 1715, 7 1716, 1717, 1718, 1719, 1719, 1721, 1722, 1723, 1724, 1725, 8 1726, 1727, 1728, 1729, 1729, 1720, 1731, 1732, 1733, 1734, 9 1735, 1736, 1737, 1738, 1739, 1740, 1741, 1742, 1743, 1744/ data (jb2(i),i= 1801, 1900)/ + 1745, 1746, 1747, 1748, 1749, 1750, 1751, 1752, 1753, 1754, 1 1755, 1756, 1757, 1758, 1759, 1760, 1761, 1762, 1763, 1764, 2 1765, 1766, 1767, 1768, 1769, 1770, 1771, 1772, 1773, 1774, 3 1775, 1776, 1777, 1778, 1779, 1780, 1781, 1782, 1783, 1784, 4 1785, 1785, 1786, 1788, 1789, 1790, 1791, 1792, 1793, 1793, 5 1794, 1795, 1796, 1797, 1798, 1799, 1800, 1801, 1802, 1803, 6 1804, 1805, 1806, 1807, 1808, 1809, 1810, 1811, 1812, 1813, 7 1814, 1815, 1816, 1817, 1817, 1818, 1819, 1820, 1821, 1821, 8 1822, 1823, 1824, 1825, 1826, 1827, 1828, 1829, 1830, 1831, 9 1832, 1833, 1834, 1835, 1835, 1836, 1837, 1838, 1839, 1840/ data (jb2(i),i= 1901, 2000)/ + 1841, 22, 1842, 1842, 1843, 24, 25, 1844, 1845, 26, 1 1846, 1847, 1848, 1849, 1850, 1851, 1852, 1853, 1854, 1855, 2 1856, 1857, 1858, 1859, 1860, 1861, 1862, 1863, 1864, 1865, 3 1866, 1867, 1868, 1869, 1870, 1871, 1872, 1873, 1874, 1875, 4 1876, 1877, 1878, 27, 2, 3, 4, 4, 6, 7, 5 8, 8, 10, 11, 12, 13, 14, 15, 15, 17, 6 18, 18, 20, 21, 22, 23, 23, 25, 26, 27, 7 27, 29, 30, 31, 31, 33, 34, 35, 36, 37, 8 37, 39, 40, 40, 42, 43, 43, 45, 46, 46, 9 48, 49, 49, 51, 52, 53, 53, 55, 56, 57/ data (jb2(i),i= 2001, 2020)/ + 57, 59, 60, 61, 62, 63, 64, 65, 66, 67, 1 68, 69, 70, 70, 72, 73, 73, 0, 0, 0/ c state=istate irot=jrot npart=ipart do i=1,5 jv(i)=kv(i) jb(i)=kb(i) ib(i)=iib(i) enddo do i=1,jv(npart+1)-1 ix(i)=jx(i) iy(i)=jy(i) enddo do i=1,jb(npart+1)-1 ib1(i)=jb1(i) ib2(i)=jb2(i) enddo return end c----------------------------------------------------------------------- c c piecewise lagrange triangle multi grid package c c edition 10.0 - - - september, 2007 c c----------------------------------------------------------------------- subroutine ak(state,irot,npart,jv,jb,ib,ix,iy,ib1,ib2) c implicit double precision (a-h,o-z) implicit integer (i-n) integer + jv(*),jb(*),ib(*),ix(*),iy(*),ib1(*),ib2(*),kv(5), 1 kb(5),iib(5),jx(1900),jy(1900),jb1(1900),jb2(1900) character*80 + state,istate save kv,kb,iib,jx,jy,jb1,jb2,istate,jrot,ipart data istate/'alaska '/ data jrot,ipart/0,2/ data kv/ 1,1172,1901, 0, 0/ data kb/ 1,1172,1901, 0, 0/ data iib/1171, 729, 0, 0, 0/ data (jx(i),i= 1, 100)/ + 11304,11304,11184,11280,11352,11304,11328,11304,11376,11304, 1 11304,11280,11256,11280,11280,11304,11376,11256,11280,11280, 2 11304,11256,11256,11232,11250,11232,11304,11304,11400,11352, 3 11328,11328,11304,11232,11232,11184,11208,11160,11232,11184, 4 11136,11112,10968,10752,10752,10776,10704,10656,10656,10632, 5 10560,10536,10584,10608,10536,10488,10608,10608,10656,10656, 6 10704,10704,10704,10800,10800,10728,10872,10920,10896,10992, 7 10992,11016,10896,10992,10944,10968,10872,10728,10728,10656, 8 10656,10632,10608,10584,10704,10656,10632,10608,10560,10608, 9 10608,10632,10632,10704,10752,10584,10588,10550,10550,10536/ data (jx(i),i= 101, 200)/ + 10488,10488,10512,10440,10440,10488,10488,10464,10392,10392, 1 10416,10416,10512,10440,10440,10560,10536,10560,10752,10790, 2 10790,10848,10848,10848,10800,10824,10800,10920,10824,10824, 3 10848,10824,10872,10968,10920,10968,11040,11112,11136,11160, 4 11088,11160,11208,11280,11376,11376,11472,11472,11544,11496, 5 11736,11832,11856,11832,11832,11784,11784,11712,11736,11832, 6 11808,11880,11904,11832,11808,11760,11736,11688,11640,11592, 7 11472,11472,11424,11400,11424,11376,11352,11304,11376,11328, 8 11256,11088,10992,11016,10968,10968,10944,10982,10920,10680, 9 10632,10656,10608,10584,10584,10608,10594,10632,10656,10680/ data (jx(i),i= 201, 300)/ + 10680,10704,10776,10656,10488,10368,10368,10440,10392,10392, 1 10440,10512,10512,10560,10632,10656,10656,10680,10704,10824, 2 10896,10968,10896,10944,11064,11088,11184,11232,11256,11352, 3 11400,11328,11376,11328,11328,11256,11304,11328,11352,11544, 4 11568,11688,11736,11736,11760,11832,11904,11904,11928,11904, 5 11784,11756,11760,11736,11736,11760,11760,11712,11640,11640, 6 11664,11688,11760,11736,11832,11904,11976,12096,12120,12096, 7 12120,12048,11976,11928,11904,11832,11760,11808,11856,11808, 8 11712,11664,11712,11712,11640,11424,11424,11400,11280,11184, 9 11064,11016,10920,11016,11016,11064,11064,11088,11448,11616/ data (jx(i),i= 301, 400)/ + 11688,11664,11712,11736,11712,11760,11760,11760,11976,12024, 1 12000,12000,11976,11952,12024,12048,12168,12168,12312,12360, 2 12360,12312,12360,12312,12336,12360,12360,12384,12408,12480, 3 12384,12380,12336,12408,12480,12504,12504,12456,12480,12480, 4 12528,12528,12600,12648,12624,12744,12840,12936,12984,13008, 5 12970,12984,13056,13056,13056,13080,13128,13152,13104,13032, 6 13056,13056,13080,13080,13128,13152,13176,13200,13200,13224, 7 13224,13248,13320,13296,13272,13392,13440,13560,13632,13704, 8 13728,13680,13680,13752,13656,13824,13800,13824,13776,13920, 9 13920,13992,14064,14064,14160,14232,14376,14400,14448,14520/ data (jx(i),i= 401, 500)/ + 14568,14904,14952,15072,15216,15264,15360,15504,15576,15672, 1 15744,15768,15768,15768,15816,15936,16200,16608,16776,16896, 2 16920,17040,17112,17256,17256,17400,17448,17784,17856,17976, 3 18048,18048,18072,18024,18096,18048,18192,18240,18384,18408, 4 18384,18432,18432,18456,18528,18576,18600,18648,18696,18816, 5 18984,19008,19080,19104,19440,19488,19488,19584,19584,19656, 6 19704,19776,19944,19968,20040,20088,20160,20184,20280,20328, 7 20328,20376,20448,20424,20424,20376,20376,20352,20304,20328, 8 20304,20280,20256,20208,20256,20280,20256,20232,20184,20160, 9 20184,20280,20256,20208,20136,20088,20112,20088,19992,20016/ data (jx(i),i= 501, 600)/ + 19968,19920,19824,19800,19824,19824,19776,19800,19824,19728, 1 19680,19704,19728,19704,19728,19714,19728,19734,19824,19704, 2 19656,19608,19584,19536,19536,19488,19488,19464,19488,19464, 3 19368,19296,19344,19320,19296,19224,19224,19200,19176,19176, 4 19104,19104,19176,19200,19176,19104,19104,19080,19056,19104, 5 19032,18990,18984,19008,18984,18984,18960,18960,18936,18864, 6 18840,18840,18840,18850,18864,18816,18816,18696,18624,18624, 7 18528,18552,18528,18504,18456,18432,18336,18336,18288,18340, 8 18384,18336,18312,18360,18384,18480,18480,18552,18576,18504, 9 18432,18456,18408,18336,18336,18312,18336,18312,18336,18240/ data (jx(i),i= 601, 700)/ + 18288,18240,18240,18216,18168,18192,18216,18192,18144,18168, 1 18120,18096,18072,18072,17948,18000,17952,17952,17976,18072, 2 18120,18144,18120,18168,18144,18192,18192,18264,18240,18264, 3 18312,18288,18240,18216,18168,18168,18192,18216,18144,18168, 4 18168,18144,18144,18096,18072,18072,18048,18024,17856,17880, 5 17832,17736,17712,17616,17544,17568,17520,17352,17376,17328, 6 17328,17160,17208,17232,17184,17184,17256,17256,17294,17280, 7 17304,17304,17266,17376,17280,17164,17164,17112,16992,16824, 8 16704,16728,16728,16680,16704,16632,16344,16056,15984,16032, 9 15960,15960,15864,15816,15768,15792,15792,15768,15792,15768/ data (jx(i),i= 701, 800)/ + 15696,15648,15576,15504,15552,15528,15504,15504,15480,15480, 1 15408,15456,15456,15408,15384,15336,15288,15288,15432,15384, 2 15360,15288,15288,15264,15312,15264,15288,15312,15288,15384, 3 15360,15288,15216,15192,15192,15168,15144,15120,15144,15120, 4 15096,15096,15096,15072,15072,15048,15048,15024,15024,15000, 5 14976,14928,15000,15000,14928,14904,14832,14856,14832,14904, 6 14880,14856,14880,14880,14856,14784,14832,14784,14784,14784, 7 14856,14856,14904,14904,14880,14832,14808,14856,14904,14904, 8 14904,14928,14928,14976,14952,14928,14928,14904,14928,14904, 9 14904,14880,14952,14928,14952,14952,14928,14904,14904,14880/ data (jx(i),i= 801, 900)/ + 14880,14856,14736,14712,14736,14712,14688,14664,14640,14640, 1 14602,14616,14592,14592,14568,14568,14544,14496,14544,14544, 2 14472,14496,14496,14472,14448,14448,14424,14424,14376,14424, 3 14352,14328,14352,14304,14328,14304,14328,14232,14256,14232, 4 14208,14160,14184,14112,14040,14016,13968,13968,13992,13968, 5 14040,14112,14088,14160,14160,14208,14088,14112,13992,14016, 6 14088,14088,14112,14112,14088,14160,14328,14424,14472,14688, 7 14640,14592,14424,14496,14616,14496,14448,14424,14304,14280, 8 14280,14184,14184,14136,13992,13968,14016,13968,13920,13848, 9 13848,13872,13824,13776,13704,13704,13800,13752,13632,13680/ data (jx(i),i= 901, 1000)/ + 13680,13632,13570,13608,13560,13560,13488,13536,13488,13488, 1 13368,13440,13392,13392,13344,13368,13368,13416,13488,13584, 2 13608,13584,13608,13584,13560,13464,13440,13464,13416,13440, 3 13344,13416,13368,13392,13392,13368,13368,13344,13344,13296, 4 13320,13272,13296,13272,13248,13224,13152,13128,13128,13080, 5 13080,12984,12984,12984,12936,12960,12864,12864,12840,12864, 6 12816,12720,12720,12768,12720,12744,12744,12720,12696,12672, 7 12696,12624,12600,12576,12576,12528,12504,12504,12480,12432, 8 12432,12384,12408,12408,12360,12336,12240,12216,12216,12144, 9 12072,12096,12096,12120,12192,12144,12144,12216,12144,12120/ data (jx(i),i= 1001, 1100)/ + 12072,12072,12048,12048,12024,12048,12048,12000,11976,11832, 1 11808,11808,11784,11760,11736,11712,11736,11760,11736,11760, 2 11688,11688,11640,11616,11592,11592,11520,11520,11472,11472, 3 11400,11448,11424,11400,11376,11256,11184,11184,11232,11304, 4 11184,11136,11112,11040,11016,11016,10992,10992,10944,10944, 5 10872,10848,10896,10848,10848,10800,10848,10824,10776,10752, 6 10704,10728,10656,10656,10680,10584,10584,10632,10608,10608, 7 10608,10704,10848,10848,10800,11088,11256,11208,11232,11352, 8 11352,11304,11328,11328,11352,11376,11376,11448,11472,11544, 9 11712,11904,11880,11880,12000,11976,12024,12048,12048,12072/ data (jx(i),i= 1101, 1200)/ + 12168,12240,12288,12360,12360,12384,12432,12384,12360,12384, 1 12456,12528,12480,12432,12408,12408,12504,12600,12552,12600, 2 12624,12648,12624,12604,12552,12528,12240,12264,12240,12240, 3 12168,12192,12216,12216,12288,12336,12312,12216,12192,12216, 4 12144,12130,12144,12120,12120,12120,12096,12110,12096,12072, 5 12120,12096,12024,11952,11880,11904,11856,11856,11808,11784, 6 11736,11712,11736,11712,11568,11544,11448,11448,11424,11400, 7 11424,19008,18840,18864,18840,18840,18816,18816,18792,18816, 8 18768,18768,18744,18696,18672,18696,18648,18624,18648,18624, 9 18624,18576,18538,18552,18558,18600,18600,18576,18648,18648/ data (jx(i),i= 1201, 1300)/ + 18696,18720,18696,18720,18768,18912,18552,18528,18552,18528, 1 18504,18552,18600,18576,18600,18600,19176,19128,19128,19104, 2 19104,19080,19080,19080,19104,19104,19128,19128,19108,19056, 3 19104,19056,19056,19032,19008,18984,19032,19008,19080,19008, 4 19104,19128,19135,19128,19152,19152,19200,19176,19224,19224, 5 19152,19176,19176,19224,19200,19224,19176,19200,19224,19224, 6 19224,19200,19176,19176,19152,19176,13752,13632,13632,13728, 7 13656,13680,13752,13752,13728,13776,13776,13824,13848,13872, 8 13944,13920,13968,13968,13944,13920,13896,13872,13896,13872, 9 13896,13824,13824,13800,13752, 8640, 8712, 8784, 8784, 8832/ data (jx(i),i= 1301, 1400)/ + 8928, 8904, 8952, 9072, 9120, 9096, 8952, 8856, 6024, 6024, 1 6096, 6072, 6120, 6120, 6144, 6144, 6216, 6192, 6168, 6216, 2 6216, 6192, 6072, 6072, 6792, 6600, 6696, 6744, 6840, 6864, 3 6864, 6960, 6912, 7032, 7056, 7008, 6984, 6984, 6888, 6912, 4 6888, 6864, 6864, 9144, 9096, 9144, 9192, 9216, 9264, 9264, 5 9312, 9312, 9384, 9336, 9384, 9360, 9384, 9408, 9432, 9432, 6 9480, 9384, 9360, 9504, 9528, 9528, 9528, 9576, 9600, 9600, 7 9648, 9648, 9624, 9528, 9552, 9528, 9552, 9576, 9576, 9600, 8 9600, 9528, 9504, 9480, 9480, 9456, 9456, 9432, 9432, 9408, 9 9312, 9192, 9634, 9610, 9682, 5736, 5688, 5736, 5688, 5784/ data (jx(i),i= 1401, 1500)/ + 5784, 5856, 5808,10176,10104,10080,10176,10248,10344,10368, 1 10488,10560,10584,10656,10680,10584,10560,10488,10320,10248, 2 19800,19752,19704,19704,19608,19584,19608,19560,19584,19656, 3 19608,19752,19776,19752,13416,13368,13488,13416,13488,13464, 4 13392,13392,13416,13296,13296,13296,13344,13392,13368,13344, 5 13320,13320,13272,13272,13248,13200,13200,13224,13200,13272, 6 13272,13296,13320,13368,13416,13416,13392,13440,13440,13512, 7 13464,13464,13512,13440,13536,13440,13440,13512,13560,13560, 8 13584,13560,13608,13560,13632,13656,13632,13680,13608,13752, 9 13728,13728,13738,13800,13800,13896,13824,13848,13848,13872/ data (jx(i),i= 1501, 1600)/ + 13848,13944,13896,13704,13728,13680,13752,13824,13800,13752, 1 13728,13656,13656,13704,13584,13584,13560,13536,13560,13488, 2 13512,13536,13536,13512,13488,13512,19992,19896,19920,19872, 3 19896,19848,19824,19824,19800,19800,19776,19776,19800,19752, 4 19752,19704,19728,19728,19680,19656,19608,19632,19536,19512, 5 19536,19632,19584,19560,19584,19536,19536,19488,19440,19416, 6 19464,19464,19440,19392,19272,19272,19272,19272,19344,19368, 7 19392,19440,19464,19440,19464,19560,19656,19704,19752,19824, 8 19824,19680,19680,19704,19680,19728,19800,19728,19800,19824, 9 19872,19824,19848,19872,19896,19896,19920,19896,19920,19872/ data (jx(i),i= 1601, 1700)/ + 19944,19944,19968,19944,19992,19968,19992,19944,19920,19920, 1 19992,19968,20064,20040,20064,20016,20040,20088,20112,20088, 2 20064,20064,20136,20098,20112,20088,20088,20040,20064,20016, 3 20050,20040,19992,19992,20016,20016,19992,19920,19944,19944, 4 19920,19920,19944,19896,19944,19872,19944,19896,19982,20136, 5 20160,20136,20160,19320,19272,19296,19344,19296,19392,19392, 6 19344,19680,19656,19656,19584,19560,19536,19560,19560,19608, 7 19608,19608,19632,19688,19664,19702,19678,19698,19632,19632, 8 19584,19656,19704,19752,19488,19416,19440,19488,19536,19526, 9 19248,19254,19230,19224,19216,19152,19152,19104,19152,19056/ data (jx(i),i= 1701, 1800)/ + 19080,19248,19272,19248,19280,19348,19368,19392,19392,19248, 1 19296,19320,19320,19344,19368,19368,19296,19296,19464,19416, 2 19398,19416,19388,19416,19512,19464,19488,18888,18864,18840, 3 18888,18840,18888,18840,18912,18912,18864,18816,18864,18816, 4 18768,18672,18672,18624,18624,18600,18576,18624,18696,18816, 5 18840,18960,18984,19008,18960,18864,18840,18792,18816,18816, 6 18840,18864,18936,18936,18960,19032,18960,18960,18984,18984, 7 19008,19032,19008,18936,18984,18936,18792,18696,18672,18696, 8 15024,15072,15048,15072,15192,15168,15240,15216,15240,15192, 9 15144,15168,15120,15144,10248,10008, 9960, 9936,10104,10104/ data (jx(i),i= 1801, 1900)/ + 10200,10224,10296,10296,10392,10368,10392,10368,10416,10368, 1 10200,10248,15336,15312,15360,15288,15336,15456, 8736, 8640, 2 8568, 8568, 8616, 8616, 8712, 9720, 9696, 9672, 9624, 9504, 3 9456, 9360, 9288, 9288, 9336, 9360, 9360, 9384, 9360, 9384, 4 9408, 9504, 9576, 9648, 9696, 9696, 9672, 9864, 9984, 9936, 5 9816,18384,18552,18576,18576,18696,18768,18720,18456,18480, 6 18456,18528,18528,18552,18552,18528,18480,18432,18504,18408, 7 18384,18384,18336,18264,18264,18264,18384,18600,18648,18672, 8 18648,18672,18648,18504,18456,18408,18456,18408,18456,18480, 9 18384,18360,18312,18288,18336,18288,18264,18240,18264,18240/ data (jy(i),i= 1, 100)/ + 14304,14328,14352,14352,14424,14424,14520,14544,14568,14592, 1 14568,14616,14688,14688,14664,14736,14784,15096,15144,15168, 2 15192,15192,15144,15144,15168,15240,15336,15360,15432,15408, 3 15456,15408,15384,15288,15264,15240,15216,15216,15168,15072, 4 15096,15048,15000,15048,15120,15144,15192,15288,15336,15312, 5 15384,15408,15408,15456,15432,15480,15528,15576,15600,15624, 6 15624,15600,15576,15552,15504,15432,15528,15504,15432,15456, 7 15504,15504,15576,15576,15600,15624,15576,15600,15648,15648, 8 15672,15648,15648,15672,15744,15744,15720,15768,15816,15840, 9 15864,15888,15912,15936,16008,15888,15840,15816,15792,15744/ data (jy(i),i= 101, 200)/ + 15768,15840,15864,15888,15936,15960,15984,16008,16008,16056, 1 16104,16080,16080,16104,16152,16152,16248,16296,16488,16488, 2 16440,16416,16440,16464,16464,16488,16488,16584,16512,16608, 3 16608,16632,16728,16704,16776,16824,16824,16800,16752,16752, 4 16680,16728,16680,16704,16776,16800,16896,16920,16872,16848, 5 16872,16992,17040,17160,17256,17328,17352,17352,17424,17400, 6 17424,17448,17496,17616,17616,17544,17520,17568,17520,17520, 7 17448,17376,17352,17424,17424,17472,17520,17496,17448,17400, 8 17496,17520,17496,17520,17496,17520,17496,17486,17472,17616, 9 17688,17760,17832,17904,17952,18000,17928,17904,17904,17928/ data (jy(i),i= 201, 300)/ + 17952,17952,17952,18024,18096,18192,18240,18288,18240,18216, 1 18264,18240,18288,18312,18312,18336,18360,18336,18384,18432, 2 18360,18384,18432,18480,18528,18504,18552,18576,18576,18576, 3 18528,18552,18528,18456,18384,18360,18360,18312,18288,18264, 4 18216,18240,18168,18192,18240,18336,18312,18240,18312,18360, 5 18408,18384,18336,18336,18360,18432,18480,18576,18600,18672, 6 18696,18696,18600,18552,18408,18456,18360,18360,18432,18480, 7 18504,18480,18504,18504,18476,18456,18576,18672,18696,18744, 8 18720,18720,18792,18816,18744,18816,18960,19056,19224,19344, 9 19416,19512,19560,19560,19584,19656,19728,19824,19776,19872/ data (jy(i),i= 301, 400)/ + 19992,19944,19968,20064,20088,20136,20160,20184,20424,20400, 1 20400,20376,20376,20352,20352,20400,20400,20424,20544,20544, 2 20496,20472,20472,20400,20424,20328,20424,20462,20472,20496, 3 20472,20554,20554,20640,20664,20664,20640,20616,20568,20616, 4 20616,20640,20616,20616,20640,20640,20736,20856,20904,20904, 5 20880,20832,20832,20808,20784,20808,20784,20760,20712,20664, 6 20640,20568,20592,20616,20616,20664,20712,20688,20736,20736, 7 20760,20760,20712,20664,20640,20592,20640,20664,20640,20640, 8 20592,20544,20520,20496,20472,20472,20448,20448,20424,20376, 9 20424,20448,20448,20400,20472,20472,20424,20376,20400,20400/ data (jy(i),i= 401, 500)/ + 20352,20352,20304,20280,20352,20376,20400,20352,20304,20280, 1 20208,20208,20232,20256,20232,19608,18288,16248,15432,15408, 2 15456,15408,15504,15552,15408,15336,15288,15072,14904,15096, 3 15096,15168,15168,15264,15312,15336,15456,15504,15432,15360, 4 15336,15336,15312,15288,15288,15240,15240,15168,15120,15096, 5 14976,14904,14856,14856,14520,14472,14400,14400,14304,14304, 6 14208,14232,14184,14184,14184,14160,14160,14112,14136,14040, 7 13944,13848,13728,13704,13512,13440,13464,13416,13512,13416, 8 13392,13392,13488,13488,13560,13560,13584,13536,13536,13584, 9 13632,13680,13704,13656,13848,13872,13896,13896,13944,14040/ data (jy(i),i= 501, 600)/ + 13944,13944,13824,13752,13752,13704,13728,13680,13632,13656, 1 13728,13728,13728,13776,13776,13848,13872,13968,14040,14016, 2 14064,14064,14088,14112,14184,14160,14184,14184,14256,14208, 3 14232,14280,14304,14352,14304,14328,14352,14352,14352,14304, 4 14328,14400,14400,14472,14424,14424,14472,14472,14520,14568, 5 14520,14544,14592,14640,14712,14736,14760,14808,14712,14736, 6 14832,14880,14904,14976,15000,14952,14784,14856,14832,14904, 7 14976,14976,15048,15000,15048,15120,15240,15264,15264,15216, 8 15144,15192,15216,15144,15096,14928,14904,14832,14736,14736, 9 14832,14784,14784,14760,14784,14784,14808,14856,14880,14952/ data (jy(i),i= 601, 700)/ + 15024,15024,14976,15024,15072,15024,14976,14928,14928,14952, 1 14952,15000,14976,14952,15024,14952,14952,14880,14928,14928, 2 14880,14856,14904,14880,14832,14832,14856,14832,14808,14808, 3 14736,14712,14688,14736,14736,14712,14712,14688,14688,14640, 4 14616,14640,14664,14664,14688,14664,14688,14664,14736,14784, 5 14760,14832,14880,14952,14952,15000,14952,15000,15024,15048, 6 15024,15096,15144,15192,15288,15336,15288,15216,15158,15144, 7 15144,15168,15264,15288,15288,15360,15312,15240,15144,15168, 8 15192,15240,15288,15312,15264,15216,15240,15120,15144,15144, 9 15192,15240,15216,15240,15264,15336,15360,15336,15456,15408/ data (jy(i),i= 701, 800)/ + 15288,15264,15312,15312,15432,15408,15408,15432,15384,15432, 1 15384,15456,15480,15432,15456,15408,15408,15432,15504,15528, 2 15480,15480,15504,15528,15528,15552,15576,15576,15600,15624, 3 15648,15624,15528,15528,15552,15552,15576,15576,15528,15504, 4 15504,15552,15576,15528,15600,15648,15552,15576,15600,15552, 5 15552,15600,15744,15768,15624,15672,15648,15648,15600,15648, 6 15600,15600,15576,15528,15528,15504,15504,15480,15456,15432, 7 15504,15408,15504,15432,15384,15384,15336,15370,15336,15360, 8 15384,15360,15408,15336,15312,15312,15288,15288,15264,15216, 9 15192,15168,15216,15192,15168,15144,15168,15072,15048,15072/ data (jy(i),i= 801, 900)/ + 15048,15024,15024,15048,15072,15072,14976,15048,15120,15048, 1 14952,14904,14904,15000,15024,14952,14880,14880,14856,14784, 2 14856,14808,14784,14760,14808,14736,14736,14712,14688,14808, 3 14664,14712,14736,14712,14688,14664,14664,14592,14568,14520, 4 14568,14592,14544,14544,14496,14544,14568,14592,14592,14640, 5 14664,14640,14688,14736,14760,14832,14760,14736,14808,14952, 6 15048,15144,15144,15240,15336,15360,15504,15432,15480,15432, 7 15480,15466,15576,15696,15768,15744,15696,15624,15624,15672, 8 15624,15600,15576,15504,15432,15384,15336,15336,15288,15192, 9 15168,15144,15096,15048,15072,15048,14976,14904,14856,14856/ data (jy(i),i= 901, 1000)/ + 14784,14736,14760,14784,14832,14736,14736,14688,14688,14640, 1 14616,14592,14568,14520,14472,14472,14424,14448,14448,14400, 2 14352,14328,14328,14256,14208,14208,14160,14136,14136,14088, 3 14040,14016,14016,13992,13968,13968,13944,13968,13920,13992, 4 13920,13944,13920,13896,13920,13872,13896,13848,13824,13800, 5 13752,13776,13728,13704,13704,13656,13656,13632,13632,13584, 6 13608,13536,13512,13512,13488,13464,13440,13440,13392,13392, 7 13344,13392,13320,13368,13320,13272,13296,13272,13248,13296, 8 13248,13224,13200,13176,13176,13224,13152,13128,13104,13104, 9 13008,13008,13032,13056,13008,12984,12960,12984,12936,12960/ data (jy(i),i= 1001, 1100)/ + 12984,12960,12936,12912,12888,12888,12840,12888,12840,12816, 1 12768,12840,12840,12792,12744,12768,12792,12816,12840,12864, 2 12888,12864,12864,12816,12840,12792,12792,12768,12720,12792, 3 12768,12744,12720,12720,12744,12672,12696,12744,12792,12768, 4 12816,12768,12720,12624,12648,12576,12552,12528,12504,12456, 5 12480,12528,12504,12600,12624,12576,12552,12504,12456,12432, 6 12480,12480,12552,12504,12456,12408,12432,12456,12456,12480, 7 12552,12600,12644,12672,12672,12888,12912,12888,12888,12912, 8 12864,12864,12840,12816,12744,12792,12840,12816,12840,13008, 9 13128,13200,13200,13224,13296,13272,13248,13224,13344,13368/ data (jy(i),i= 1101, 1200)/ + 13512,13536,13584,13622,13584,13536,13536,13608,13632,13848, 1 13910,13896,13920,13920,13968,14016,14136,14184,14208,14280, 2 14328,14400,14376,14328,14256,14256,14136,14232,14256,14304, 3 14352,14424,14424,14448,14448,14376,14448,14496,14544,14472, 4 14376,14424,14448,14448,14400,14328,14304,14352,14352,14328, 5 14184,14136,14160,14352,14376,14448,14448,14400,14352,14424, 6 14424,14472,14472,14520,14424,14448,14424,14400,14376,14352, 7 14328,13728,13872,13968,13944,13992,13892,13920,13944,13992, 8 13944,14016,13992,14016,14064,14112,14112,14136,14160,14146, 9 14208,14232,14256,14280,14304,14352,14328,14280,14328,14304/ data (jy(i),i= 1201, 1300)/ + 14304,14280,14280,14256,14232,13992,14040,14112,14160,14136, 1 14208,14232,14184,14184,14136,14064,13704,13728,13776,13800, 2 13848,13824,13848,13896,13872,13896,13872,13944,13920,13944, 3 13992,13968,14016,13992,14016,14064,14064,14112,14088,14136, 4 14112,14088,14040,13992,14064,14088,14088,14040,14016,13992, 5 13992,13982,13944,13920,13896,13824,13872,13824,13800,13776, 6 13728,13728,13824,13752,13752,13728,13876,13968,13992,13968, 7 14016,14040,14016,14064,14088,14112,14136,14136,14064,14136, 8 14088,14040,14088,14064,13992,13992,14040,14040,14016,14016, 9 13968,13944,14016,13944,13944,11712,11756,11808,11856,11904/ data (jy(i),i= 1301, 1400)/ + 11880,11928,11976,12000,11952,11904,11856,11784,11808,11832, 1 11880,11904,11880,11904,11904,11952,11952,11928,11904,11880, 2 11808,11832,11784,11856,11784,11832,11818,11856,11808,11808, 3 11832,11856,11904,11928,11856,11808,11832,11760,11808,11784, 4 11760,11798,11760,11808,11856,11904,11880,11904,11880,11904, 5 11904,11928,11928,11976,11976,12000,12024,12000,11952,12000, 6 12000,12048,12096,12168,12120,12096,12048,12072,12144,12096, 7 12120,12096,12072,12024,12024,12000,11976,12000,11976,12024, 8 11952,11952,11904,11928,11904,11928,11880,11928,11880,11904, 9 11880,11808,11976,12024,12024,11904,11976,12000,12096,12072/ data (jy(i),i= 1401, 1500)/ + 12024,12000,12000,12240,12264,12360,12384,12504,12480,12504, 1 12552,12528,12384,12360,12312,12360,12336,12312,12336,12264, 2 13152,13152,13200,13224,13248,13320,13344,13368,13416,13320, 3 13320,13224,13200,13176,13224,13224,13344,13344,13416,13440, 4 13344,13368,13416,13416,13392,13368,13406,13406,13320,13296, 5 13272,13320,13344,13440,13488,13512,13536,13536,13584,13632, 6 13656,13656,13680,13704,13680,13656,13632,13632,13560,13488, 7 13560,13656,13632,13680,13680,13728,13776,13824,13680,13728, 8 13728,13752,13728,13800,13776,13728,13824,13800,13872,13848, 9 13800,13728,13800,13800,13848,13776,13728,13728,13704,13728/ data (jy(i),i= 1501, 1600)/ + 13680,13680,13584,13632,13608,13584,13608,13560,13512,13488, 1 13536,13536,13512,13488,13464,13440,13416,13440,13392,13368, 2 13368,13344,13320,13344,13296,13272,13224,13224,13272,13248, 3 13296,13320,13368,13296,13344,13296,13320,13344,13368,13392, 4 13416,13464,13416,13392,13416,13440,13392,13440,13416,13440, 5 13464,13488,13488,13512,13512,13536,13584,13560,13560,13632, 6 13656,13728,13728,13848,13848,13872,13896,13920,13920,13896, 7 13944,13920,13872,13848,13824,13824,13752,13656,13632,13608, 8 13584,13656,13608,13608,13560,13584,13584,13536,13560,13560, 9 13512,13488,13464,13488,13536,13512,13440,13440,13392,13368/ data (jy(i),i= 1601, 1700)/ + 13386,13368,13344,13320,13296,13272,13488,13512,13584,13608, 1 13560,13536,13416,13488,13512,13560,13560,13560,13464,13440, 2 13488,13464,13584,13584,13608,13656,13584,13584,13680,13776, 3 13656,13608,13704,13680,13608,13584,13584,13632,13704,13728, 4 13728,13752,13800,13776,13824,13824,13872,13872,13944,13824, 5 13680,13680,13608,13704,13680,13776,13800,13800,13824,13752, 6 13752,13800,13824,13872,13848,13872,13944,13944,14016,14016, 7 13968,13944,13958,13934,13920,13848,13848,13944,13968,13992, 8 14078,13990,14004,13992,13920,13968,14016,14040,14016,13944, 9 13968,14016,14016,14040,14094,14160,14136,14160,14184,14208/ data (jy(i),i= 1701, 1800)/ + 14256,14256,14232,14280,14280,14232,14208,14112,14088,14184, 1 14126,14126,14088,14088,14040,14016,14016,13992,14064,14064, 2 14088,14112,14184,14208,14136,14088,14064,14160,14232,14256, 3 14280,14304,14328,14352,14352,14376,14376,14400,14448,14424, 4 14448,14664,14688,14736,14760,14784,14832,14808,14736,14780, 5 14736,14616,14544,14520,14544,14688,14712,14736,14712,14688, 6 14592,14616,14544,14520,14520,14448,14448,14424,14424,14400, 7 14424,14400,14352,14376,14328,14208,14784,14784,14808,14832, 8 14904,14952,14976,15000,15144,15216,15216,15168,15192,15096, 9 15024,15000,15000,14952,15096,15264,15312,15384,15360,15408/ data (jy(i),i= 1801, 1900)/ + 15432,15408,15432,15384,15360,15288,15240,15240,15192,15144, 1 15144,15120,15168,15192,15240,15240,15288,15264,15744,15768, 2 15888,15912,15936,15864,15792,16920,17016,17040,17136,17256, 3 17256,17208,17256,17352,17424,17472,17424,17400,17376,17376, 4 17376,17328,17352,17352,17280,17256,17232,17088,17064,16992, 5 17040,14544,14496,14472,14448,14520,14376,14352,14448,14400, 6 14376,14376,14352,14352,14304,14256,14256,14304,14280,14352, 7 14352,14376,14424,14424,14448,14520,14562,14496,14520,14544, 8 14568,14592,14640,14664,14554,14568,14592,14592,14616,14664, 9 14712,14640,14664,14640,14568,14568,14640,14640,14592,14558/ data (jb1(i),i= 1, 100)/ + 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 2 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 3 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 4 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 5 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 6 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 7 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 8 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 9 91, 92, 93, 94, 95, 96, 97, 98, 99, 100/ data (jb1(i),i= 101, 200)/ + 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 1 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 2 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 3 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 4 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 5 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 6 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 7 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 8 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 9 191, 192, 193, 194, 195, 196, 197, 198, 199, 200/ data (jb1(i),i= 201, 300)/ + 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, 1 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 2 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 3 231, 232, 233, 234, 235, 236, 237, 238, 239, 240, 4 241, 242, 243, 244, 245, 246, 247, 248, 249, 250, 5 251, 252, 253, 254, 255, 256, 257, 258, 259, 260, 6 261, 262, 263, 264, 265, 266, 267, 268, 269, 270, 7 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 8 281, 282, 283, 284, 285, 286, 287, 288, 289, 290, 9 291, 292, 293, 294, 295, 296, 297, 298, 299, 300/ data (jb1(i),i= 301, 400)/ + 301, 302, 303, 304, 305, 306, 307, 308, 309, 310, 1 311, 312, 313, 314, 315, 316, 317, 318, 319, 320, 2 321, 322, 323, 324, 325, 326, 327, 328, 329, 330, 3 331, 332, 333, 334, 335, 336, 337, 338, 339, 340, 4 341, 342, 343, 344, 345, 346, 347, 348, 349, 350, 5 351, 352, 353, 354, 355, 356, 357, 358, 359, 360, 6 361, 362, 363, 364, 365, 366, 367, 368, 369, 370, 7 371, 372, 373, 374, 375, 376, 377, 378, 379, 380, 8 381, 382, 383, 384, 385, 386, 387, 388, 389, 390, 9 391, 392, 393, 394, 395, 396, 397, 398, 399, 400/ data (jb1(i),i= 401, 500)/ + 401, 402, 403, 404, 405, 406, 407, 408, 409, 410, 1 411, 412, 413, 414, 415, 416, 417, 418, 419, 420, 2 421, 422, 423, 424, 425, 426, 427, 428, 429, 430, 3 431, 432, 433, 434, 435, 436, 437, 438, 439, 440, 4 441, 442, 443, 444, 445, 446, 447, 448, 449, 450, 5 451, 452, 453, 454, 455, 456, 457, 458, 459, 460, 6 461, 462, 463, 464, 465, 466, 467, 468, 469, 470, 7 471, 472, 473, 474, 475, 476, 477, 478, 479, 480, 8 481, 482, 483, 484, 485, 486, 487, 488, 489, 490, 9 491, 492, 493, 494, 495, 496, 497, 498, 499, 500/ data (jb1(i),i= 501, 600)/ + 501, 502, 503, 504, 505, 506, 507, 508, 509, 510, 1 511, 512, 513, 514, 515, 516, 517, 518, 519, 520, 2 521, 522, 523, 524, 525, 526, 527, 528, 529, 530, 3 531, 532, 533, 534, 535, 536, 537, 538, 539, 540, 4 541, 542, 543, 544, 545, 546, 547, 548, 549, 550, 5 551, 552, 553, 554, 555, 556, 557, 558, 559, 560, 6 561, 562, 563, 564, 565, 566, 567, 568, 569, 570, 7 571, 572, 573, 574, 575, 576, 577, 578, 579, 580, 8 581, 582, 583, 584, 585, 586, 587, 588, 589, 590, 9 591, 592, 593, 594, 595, 596, 597, 598, 599, 600/ data (jb1(i),i= 601, 700)/ + 601, 602, 603, 604, 605, 606, 607, 608, 609, 610, 1 611, 612, 613, 614, 615, 616, 617, 618, 619, 620, 2 621, 622, 623, 624, 625, 626, 627, 628, 629, 630, 3 631, 632, 633, 634, 635, 636, 637, 638, 639, 640, 4 641, 642, 643, 644, 645, 646, 647, 648, 649, 650, 5 651, 652, 653, 654, 655, 656, 657, 658, 659, 660, 6 661, 662, 663, 664, 665, 666, 667, 668, 669, 670, 7 671, 672, 673, 674, 675, 676, 677, 678, 679, 680, 8 681, 682, 683, 684, 685, 686, 687, 688, 689, 690, 9 691, 692, 693, 694, 695, 696, 697, 698, 699, 700/ data (jb1(i),i= 701, 800)/ + 701, 702, 703, 704, 705, 706, 707, 708, 709, 710, 1 711, 712, 713, 714, 715, 716, 717, 718, 719, 720, 2 721, 722, 723, 724, 725, 726, 727, 728, 729, 730, 3 731, 732, 733, 734, 735, 736, 737, 738, 739, 740, 4 741, 742, 743, 744, 745, 746, 747, 748, 749, 750, 5 751, 752, 753, 754, 755, 756, 757, 758, 759, 760, 6 761, 762, 763, 764, 765, 766, 767, 768, 769, 770, 7 771, 772, 773, 774, 775, 776, 777, 778, 779, 780, 8 781, 782, 783, 784, 785, 786, 787, 788, 789, 790, 9 791, 792, 793, 794, 795, 796, 797, 798, 799, 800/ data (jb1(i),i= 801, 900)/ + 801, 802, 803, 804, 805, 806, 807, 808, 809, 810, 1 811, 812, 813, 814, 815, 816, 817, 818, 819, 820, 2 821, 822, 823, 824, 825, 826, 827, 828, 829, 830, 3 831, 832, 833, 834, 835, 836, 837, 838, 839, 840, 4 841, 842, 843, 844, 845, 846, 847, 848, 849, 850, 5 851, 852, 853, 854, 855, 856, 857, 858, 859, 860, 6 861, 862, 863, 864, 865, 866, 867, 868, 869, 870, 7 871, 872, 873, 874, 875, 876, 877, 878, 879, 880, 8 881, 882, 883, 884, 885, 886, 887, 888, 889, 890, 9 891, 892, 893, 894, 895, 896, 897, 898, 899, 900/ data (jb1(i),i= 901, 1000)/ + 901, 902, 903, 904, 905, 906, 907, 908, 909, 910, 1 911, 912, 913, 914, 915, 916, 917, 918, 919, 920, 2 921, 922, 923, 924, 925, 926, 927, 928, 929, 930, 3 931, 932, 933, 934, 935, 936, 937, 938, 939, 940, 4 941, 942, 943, 944, 945, 946, 947, 948, 949, 950, 5 951, 952, 953, 954, 955, 956, 957, 958, 959, 960, 6 961, 962, 963, 964, 965, 966, 967, 968, 969, 970, 7 971, 972, 973, 974, 975, 976, 977, 978, 979, 980, 8 981, 982, 983, 984, 985, 986, 987, 988, 989, 990, 9 991, 992, 993, 994, 995, 996, 997, 998, 999, 1000/ data (jb1(i),i= 1001, 1100)/ + 1001, 1002, 1003, 1004, 1005, 1006, 1007, 1008, 1009, 1010, 1 1011, 1012, 1013, 1014, 1015, 1016, 1017, 1018, 1019, 1020, 2 1021, 1022, 1023, 1024, 1025, 1026, 1027, 1028, 1029, 1030, 3 1031, 1032, 1033, 1034, 1035, 1036, 1037, 1038, 1039, 1040, 4 1041, 1042, 1043, 1044, 1045, 1046, 1047, 1048, 1049, 1050, 5 1051, 1052, 1053, 1054, 1055, 1056, 1057, 1058, 1059, 1060, 6 1061, 1062, 1063, 1064, 1065, 1066, 1067, 1068, 1069, 1070, 7 1071, 1072, 1073, 1074, 1075, 1076, 1077, 1078, 1079, 1080, 8 1081, 1082, 1083, 1084, 1085, 1086, 1087, 1088, 1089, 1090, 9 1091, 1092, 1093, 1094, 1095, 1096, 1097, 1098, 1099, 1100/ data (jb1(i),i= 1101, 1200)/ + 1101, 1102, 1103, 1104, 1105, 1106, 1107, 1108, 1109, 1110, 1 1111, 1112, 1113, 1114, 1115, 1116, 1117, 1118, 1119, 1120, 2 1121, 1122, 1123, 1124, 1125, 1126, 1127, 1128, 1129, 1130, 3 1131, 1132, 1133, 1134, 1135, 1136, 1137, 1138, 1139, 1140, 4 1141, 1142, 1143, 1144, 1145, 1146, 1147, 1148, 1149, 1150, 5 1151, 1152, 1153, 1154, 1155, 1156, 1157, 1158, 1159, 1160, 6 1161, 1162, 1163, 1164, 1165, 1166, 1167, 1168, 1169, 1170, 7 1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 8 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 9 20, 21, 22, 23, 24, 25, 26, 27, 28, 29/ data (jb1(i),i= 1201, 1300)/ + 30, 31, 32, 33, 34, 1, 36, 37, 38, 39, 1 40, 41, 42, 43, 44, 36, 46, 47, 48, 49, 2 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 3 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 4 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 5 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 6 90, 91, 92, 93, 94, 46, 96, 97, 98, 99, 7 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 8 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 9 120, 121, 122, 123, 96, 125, 126, 127, 128, 129/ data (jb1(i),i= 1301, 1400)/ + 130, 131, 132, 133, 134, 135, 136, 125, 138, 139, 1 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 2 150, 151, 152, 138, 154, 155, 156, 157, 158, 159, 3 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 4 170, 171, 154, 173, 174, 175, 176, 177, 178, 179, 5 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 6 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 7 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 8 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 9 220, 173, 222, 223, 222, 225, 226, 227, 228, 229/ data (jb1(i),i= 1401, 1500)/ + 230, 231, 225, 233, 234, 235, 236, 237, 238, 239, 1 240, 241, 242, 243, 244, 245, 246, 247, 248, 233, 2 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 3 260, 261, 262, 250, 264, 265, 266, 267, 268, 269, 4 270, 271, 272, 273, 274, 275, 276, 277, 278, 279, 5 280, 281, 282, 283, 284, 285, 286, 287, 288, 289, 6 290, 291, 292, 293, 294, 295, 296, 297, 298, 299, 7 300, 301, 302, 303, 304, 305, 306, 307, 308, 309, 8 310, 311, 312, 313, 314, 315, 316, 317, 318, 319, 9 320, 321, 322, 323, 324, 325, 326, 327, 328, 329/ data (jb1(i),i= 1501, 1600)/ + 330, 331, 332, 333, 334, 335, 336, 337, 338, 339, 1 340, 341, 342, 343, 344, 345, 346, 347, 348, 349, 2 350, 351, 352, 353, 354, 264, 356, 357, 358, 359, 3 360, 361, 362, 363, 364, 365, 366, 367, 368, 369, 4 370, 371, 372, 373, 374, 375, 376, 377, 378, 379, 5 380, 381, 382, 383, 384, 385, 386, 387, 388, 389, 6 390, 391, 392, 393, 394, 395, 396, 397, 398, 399, 7 400, 401, 402, 403, 404, 405, 406, 407, 408, 409, 8 410, 411, 412, 413, 414, 415, 416, 417, 418, 419, 9 420, 421, 422, 423, 424, 425, 426, 427, 428, 429/ data (jb1(i),i= 1601, 1700)/ + 430, 431, 432, 433, 434, 356, 436, 437, 438, 439, 1 440, 436, 442, 443, 444, 445, 446, 447, 448, 449, 2 450, 442, 452, 453, 454, 455, 456, 457, 458, 459, 3 460, 461, 462, 463, 464, 465, 466, 467, 468, 469, 4 470, 471, 472, 473, 474, 475, 476, 477, 478, 479, 5 480, 481, 452, 483, 484, 485, 486, 487, 488, 489, 6 483, 491, 492, 493, 494, 495, 496, 497, 498, 499, 7 500, 501, 502, 503, 504, 505, 491, 507, 508, 509, 8 510, 511, 512, 507, 514, 515, 516, 517, 518, 514, 9 520, 521, 522, 523, 524, 525, 526, 527, 528, 529/ data (jb1(i),i= 1701, 1800)/ + 530, 531, 532, 533, 534, 535, 536, 537, 538, 539, 1 540, 541, 542, 543, 544, 545, 546, 520, 548, 549, 2 550, 551, 552, 553, 554, 555, 548, 557, 558, 559, 3 560, 561, 562, 563, 564, 565, 566, 567, 568, 569, 4 570, 571, 572, 573, 574, 575, 576, 577, 578, 579, 5 580, 581, 582, 583, 584, 585, 586, 587, 588, 589, 6 590, 591, 592, 593, 594, 595, 596, 597, 598, 599, 7 600, 601, 602, 603, 604, 557, 606, 607, 608, 606, 8 610, 611, 612, 613, 614, 615, 616, 617, 618, 619, 9 620, 621, 622, 610, 624, 625, 626, 627, 628, 629/ data (jb1(i),i= 1801, 1900)/ + 630, 631, 632, 633, 634, 635, 636, 637, 638, 639, 1 640, 624, 642, 643, 644, 645, 646, 642, 648, 649, 2 650, 651, 652, 653, 648, 655, 656, 657, 658, 659, 3 660, 661, 662, 663, 664, 665, 666, 667, 668, 669, 4 670, 671, 672, 673, 674, 675, 676, 677, 678, 679, 5 655, 681, 682, 683, 684, 685, 686, 687, 688, 689, 6 690, 691, 692, 693, 694, 695, 696, 697, 698, 699, 7 700, 701, 702, 703, 704, 681, 706, 707, 708, 709, 8 710, 711, 712, 713, 714, 715, 716, 717, 718, 719, 9 720, 721, 722, 723, 724, 725, 726, 727, 728, 706/ data (jb2(i),i= 1, 100)/ + 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 1 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 2 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 3 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 4 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 5 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 6 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 7 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 8 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 9 92, 93, 94, 95, 96, 97, 98, 99, 100, 101/ data (jb2(i),i= 101, 200)/ + 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 1 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 2 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 3 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 4 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 5 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 6 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 7 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 8 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 9 192, 193, 194, 195, 196, 197, 198, 199, 200, 201/ data (jb2(i),i= 201, 300)/ + 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 1 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 2 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 3 232, 233, 234, 235, 236, 237, 238, 239, 240, 241, 4 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 5 252, 253, 254, 255, 256, 257, 258, 259, 260, 261, 6 262, 263, 264, 265, 266, 267, 268, 269, 270, 271, 7 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 8 282, 283, 284, 285, 286, 287, 288, 289, 290, 291, 9 292, 293, 294, 295, 296, 297, 298, 299, 300, 301/ data (jb2(i),i= 301, 400)/ + 302, 303, 304, 305, 306, 307, 308, 309, 310, 311, 1 312, 313, 314, 315, 316, 317, 318, 319, 320, 321, 2 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 3 332, 333, 334, 335, 336, 337, 338, 339, 340, 341, 4 342, 343, 344, 345, 346, 347, 348, 349, 350, 351, 5 352, 353, 354, 355, 356, 357, 358, 359, 360, 361, 6 362, 363, 364, 365, 366, 367, 368, 369, 370, 371, 7 372, 373, 374, 375, 376, 377, 378, 379, 380, 381, 8 382, 383, 384, 385, 386, 387, 388, 389, 390, 391, 9 392, 393, 394, 395, 396, 397, 398, 399, 400, 401/ data (jb2(i),i= 401, 500)/ + 402, 403, 404, 405, 406, 407, 408, 409, 410, 411, 1 412, 413, 414, 415, 416, 417, 418, 419, 420, 421, 2 422, 423, 424, 425, 426, 427, 428, 429, 430, 431, 3 432, 433, 434, 435, 436, 437, 438, 439, 440, 441, 4 442, 443, 444, 445, 446, 447, 448, 449, 450, 451, 5 452, 453, 454, 455, 456, 457, 458, 459, 460, 461, 6 462, 463, 464, 465, 466, 467, 468, 469, 470, 471, 7 472, 473, 474, 475, 476, 477, 478, 479, 480, 481, 8 482, 483, 484, 485, 486, 487, 488, 489, 490, 491, 9 492, 493, 494, 495, 496, 497, 498, 499, 500, 501/ data (jb2(i),i= 501, 600)/ + 502, 503, 504, 505, 506, 507, 508, 509, 510, 511, 1 512, 513, 514, 515, 516, 517, 518, 519, 520, 521, 2 522, 523, 524, 525, 526, 527, 528, 529, 530, 531, 3 532, 533, 534, 535, 536, 537, 538, 539, 540, 541, 4 542, 543, 544, 545, 546, 547, 548, 549, 550, 551, 5 552, 553, 554, 555, 556, 557, 558, 559, 560, 561, 6 562, 563, 564, 565, 566, 567, 568, 569, 570, 571, 7 572, 573, 574, 575, 576, 577, 578, 579, 580, 581, 8 582, 583, 584, 585, 586, 587, 588, 589, 590, 591, 9 592, 593, 594, 595, 596, 597, 598, 599, 600, 601/ data (jb2(i),i= 601, 700)/ + 602, 603, 604, 605, 606, 607, 608, 609, 610, 611, 1 612, 613, 614, 615, 616, 617, 618, 619, 620, 621, 2 622, 623, 624, 625, 626, 627, 628, 629, 630, 631, 3 632, 633, 634, 635, 636, 637, 638, 639, 640, 641, 4 642, 643, 644, 645, 646, 647, 648, 649, 650, 651, 5 652, 653, 654, 655, 656, 657, 658, 659, 660, 661, 6 662, 663, 664, 665, 666, 667, 668, 669, 670, 671, 7 672, 673, 674, 675, 676, 677, 678, 679, 680, 681, 8 682, 683, 684, 685, 686, 687, 688, 689, 690, 691, 9 692, 693, 694, 695, 696, 697, 698, 699, 700, 701/ data (jb2(i),i= 701, 800)/ + 702, 703, 704, 705, 706, 707, 708, 709, 710, 711, 1 712, 713, 714, 715, 716, 717, 718, 719, 720, 721, 2 722, 723, 724, 725, 726, 727, 728, 729, 730, 731, 3 732, 733, 734, 735, 736, 737, 738, 739, 740, 741, 4 742, 743, 744, 745, 746, 747, 748, 749, 750, 751, 5 752, 753, 754, 755, 756, 757, 758, 759, 760, 761, 6 762, 763, 764, 765, 766, 767, 768, 769, 770, 771, 7 772, 773, 774, 775, 776, 777, 778, 779, 780, 781, 8 782, 783, 784, 785, 786, 787, 788, 789, 790, 791, 9 792, 793, 794, 795, 796, 797, 798, 799, 800, 801/ data (jb2(i),i= 801, 900)/ + 802, 803, 804, 805, 806, 807, 808, 809, 810, 811, 1 812, 813, 814, 815, 816, 817, 818, 819, 820, 821, 2 822, 823, 824, 825, 826, 827, 828, 829, 830, 831, 3 832, 833, 834, 835, 836, 837, 838, 839, 840, 841, 4 842, 843, 844, 845, 846, 847, 848, 849, 850, 851, 5 852, 853, 854, 855, 856, 857, 858, 859, 860, 861, 6 862, 863, 864, 865, 866, 867, 868, 869, 870, 871, 7 872, 873, 874, 875, 876, 877, 878, 879, 880, 881, 8 882, 883, 884, 885, 886, 887, 888, 889, 890, 891, 9 892, 893, 894, 895, 896, 897, 898, 899, 900, 901/ data (jb2(i),i= 901, 1000)/ + 902, 903, 904, 905, 906, 907, 908, 909, 910, 911, 1 912, 913, 914, 915, 916, 917, 918, 919, 920, 921, 2 922, 923, 924, 925, 926, 927, 928, 929, 930, 931, 3 932, 933, 934, 935, 936, 937, 938, 939, 940, 941, 4 942, 943, 944, 945, 946, 947, 948, 949, 950, 951, 5 952, 953, 954, 955, 956, 957, 958, 959, 960, 961, 6 962, 963, 964, 965, 966, 967, 968, 969, 970, 971, 7 972, 973, 974, 975, 976, 977, 978, 979, 980, 981, 8 982, 983, 984, 985, 986, 987, 988, 989, 990, 991, 9 992, 993, 994, 995, 996, 997, 998, 999, 1000, 1001/ data (jb2(i),i= 1001, 1100)/ + 1002, 1003, 1004, 1005, 1006, 1007, 1008, 1009, 1010, 1011, 1 1012, 1013, 1014, 1015, 1016, 1017, 1018, 1019, 1020, 1021, 2 1022, 1023, 1024, 1025, 1026, 1027, 1028, 1029, 1030, 1031, 3 1032, 1033, 1034, 1035, 1036, 1037, 1038, 1039, 1040, 1041, 4 1042, 1043, 1044, 1045, 1046, 1047, 1048, 1049, 1050, 1051, 5 1052, 1053, 1054, 1055, 1056, 1057, 1058, 1059, 1060, 1061, 6 1062, 1063, 1064, 1065, 1066, 1067, 1068, 1069, 1070, 1071, 7 1072, 1073, 1074, 1075, 1076, 1077, 1078, 1079, 1080, 1081, 8 1082, 1083, 1084, 1085, 1086, 1087, 1088, 1089, 1090, 1091, 9 1092, 1093, 1094, 1095, 1096, 1097, 1098, 1099, 1100, 1101/ data (jb2(i),i= 1101, 1200)/ + 1102, 1103, 1104, 1105, 1106, 1107, 1108, 1109, 1110, 1111, 1 1112, 1113, 1114, 1115, 1116, 1117, 1118, 1119, 1120, 1121, 2 1122, 1123, 1124, 1125, 1126, 1127, 1128, 1129, 1130, 1131, 3 1132, 1133, 1134, 1135, 1136, 1137, 1138, 1139, 1140, 1141, 4 1142, 1143, 1144, 1145, 1146, 1147, 1148, 1149, 1150, 1151, 5 1152, 1153, 1154, 1155, 1156, 1157, 1158, 1159, 1160, 1161, 6 1162, 1163, 1164, 1165, 1166, 1167, 1168, 1169, 1170, 1171, 7 1171, 2, 3, 4, 5, 6, 7, 8, 9, 10, 8 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 9 21, 22, 23, 24, 25, 26, 27, 28, 29, 30/ data (jb2(i),i= 1201, 1300)/ + 31, 32,