****************************************
*** チェビシェフ漸化式のサブルーチン ***
****************************************
      subroutine CHEBYSHEV(n, nmax, che)
      real che(0:n,nmax)
      real xmin, xmax
      real dx
      xmin=-1.0
      xmax=1.0
      dx=(xmax-xmin)/(nmax-1)
*-- 初項の設定 ---
      do 10 i=1,nmax
         che(0,i)=1.0
  10  continue
      if(n.eq.0)then
         go to 100
      end if
      do 20 i=1,nmax
         che(1,i)=xmin+dx*(real(i)-1.0)
  20  continue
      if(n.eq.1)then
         go to 100
      end if
*-- 漸化式の計算 ---
      do 30 j=1,n-1
         do 40 i=1,nmax
            che(j+1,i)=2.0*che(1,i)*che(j,i)-che(j-1,i)
  40  continue
  30  continue
  100 return
      end