c      main program  point1.fi
c      and symblv test program
c      parameter (nx=10,ny=10,nz=10)
       parameter (nx=20,ny=20,nz=20)
       parameter (n1=nx,n2=n1*ny,n3=n2*nz)
       parameter (x0=0.0,y0=0.0,z0=0.0)
       parameter (xl=9.0,yl=9.0,zl=9.0)
c      parameter (ico=5,rd=1.0)
       parameter (ico=7,rd=1.0)
       dimension u(n3)
c
       hx=xl/float(nx-1)
       hy=yl/float(ny-1)
       hz=zl/float(nz-1)
c
       i1=101
       write(6,*) i1
       do 10 k=1,nz
       do 10 j=1,ny
       do 10 i=1,nx
c      x=hx*float(i-1)
c      y=hy*float(j-1)
c      z=hz*float(k-1)
       x=0.5*hx*float(2*i-1-nx)
       y=0.5*hy*float(2*j-1-ny)
       z=0.5*hz*float(2*k-1-nz)
       r2=x*x+y*y+z*z
       r1=sqrt(r2)
       i1=i+n1*(j-1)+n2*(k-1)
       u(i1)=exp(-(r1/rd)**2)
       u(i1)=1.0/(1.0+r2/rd)
c      u(i1)=x/xl
   10  continue
c
       i1=101
       write(6,*) i1
c
       call initvrml
c
       call point3d(nx,ny,nz,x0,y0,z0,xl,yl,zl,ico,u)
c
       stop
       end