c----------------------------------------------
c      triangb1.f
c    
c      subrtn file : pk1subrtn
c       
c----------------------------------------------
c      main program  triangb1.f
c      triangb1.f    for triangle belt program
c----( nx: x point, ny: y point )
       parameter (nx=20,ny=2)
       parameter (pai=3.1415926,rd=5.0)
c----( p: x,y,z no jun ni set   u: r,g,b no jun ni set )
       dimension p(3*nx*2)
       dimension u(3*nx*2)
c
       call initvrml
c
c################################( test data )
       dd=2.0*pai/float(nx)
c
       do 20  i=1,nx
       do 20  j=1,2
       i1=3*(i-1)+1+3*nx*(j-1)
       i2=i1+1
       i3=i1+2
       p(i1)=i
       p(i2)=ny*(j-1)
       p(i3)=rd*sin(float(i)*dd)
       u(i1)=i*(1/nx)+(j-1)
       u(i2)=i*(1/nx)
       u(i3)=1.0
   20  continue
c#############################################
c
       call triangbk1(nx,ny,u,p)
c
       stop
       end
