c      main program
       lasl=12
       r=1.0
       g=1.0
       b=1.0
c
       write(10,*) "#VRML V2.0 utf8"
       call linebe(r,g,b)
c
       do 40 ii=1,lasl
       lin=ii
       li1=3*(lin-1)+1
       li2=3*(lin-1)+2
       li3=3*lin
c      X1=po(li1)
c      Y1=po(li2)
c      Z1=po(li3)
       x1=0.1*float(ii-1)
       y1=0.1*float(ii-1)
       z1=0.1*float(ii-1)
c
       if(ii.lt.lasl) then
       write(10,152) x1,z1,y1
  152  format(1h ,f5.1,1x,f5.1,1x,f5.1,1h,)
       else
       write(10,154) x1,z1,y1
  154  format(1h ,f5.1,1x,f5.1,1x,f5.1)
       end if
   40  continue
c
       call lineen(lasl)
c
       stop
       end
c      subroutine linebe(r,g,b)
       subroutine linebe(r,g,b)
c
       write(10,*) " Shape {"
       write(10,*) " geometry IndexedLineSet {"
       write(10,*) " color Color {"
c      write(10,*) " color ["
       write(10,21) r,g,b
   21  format(" color [ ",f3.1,1x,f3.1,1x,f3.1," ]")
       write(10,*) " }"
       write(10,*) " coord Coordinate {"
       write(10,*) " point ["
       return
       end
       subroutine lineen(lasl)
       integer ic(1000)
c
       do 20 i=1,1000
       ic(i)=i-1
   20  continue
       write(10,*) " ]"
       write(10,*) " }"
       write(10,*) " colorIndex [ 0 ]"
       write(10,*) " coordIndex ["
c
       write(10,100) (ic(i),i=1,lasl)
  100  format(1h ,10i2,90i3,900i4)
c
       write(10,*) " ]"
       write(10,*) " colorPerVertex FALSE"
       write(10,*) " }"
       write(10,*) " }"
       return
       end