        implicit real*8(a-h,o-z)
        dimension x(1000000),y(1000000),z(1000000)
        print'("** 山形綱の載荷節点番号と拘束節点番号")'
         print'("** だしてくれるプログラム")'
         ! nx:x方向分割数(山形綱横方向分割数)
         ! ny:y方向分割数(山形綱縦方向分割数)
         ! nz:z方向分割数(山形綱奥行き方向分割数
         print *,"山形綱の縦の分割数は"
         read *,nx
         print *,"山形綱の横の分割数は"
         read *,ny
         print *,"山形綱の奥行きの分割数は"
         read *,nz
         print *,"山形綱の拘束する節点(載荷点じゃない方は):"
         do i=1,nx+1
         print *,i
         end do
         do j=0,ny-1
         print *,(nx+1)*(nz+1)+j*(nz+1)+1
         end do
         print *,"山形綱の拘束する節点(載荷節点は):"
         do i=1,ny+1
         print *,(nx+1)*nz+i,','
         end do
         do j=0,ny-1
         print *,(nx+1)*(nz+1)+(nz+1)+(nz+1)*j
         end do
         end
