Chebyshev polynomial T2 

c chebyshev y2 = 2y1 - y0 + (dx^2) * ( x1*(y1-y0)/dx   -n^2*y1     )
c       /(1-x1^2)     dx<<1
      data n, nstep / 2, 1000/
      kp=int(float(nstep)/50.)
      kount=kp
      y0=-1.
      dx=1./float(nstep)
      dtdx=0.
      y1=dtdx*dx+y0
      print 100 , 0. , y0
      do 10 i=2,nstep
      x=dx*float(i)
      y2= 2.*y1 - y0 +(dx**2)*((x-dx)*(y1-y0)/dx  -float(n)**2*y1)
     $   /(1-(x-dx)**2)
      if(i.eq.kount)then
      print 100 , x , y2
      kount=kount+kp
      endif
      y0=y1
      y1=y2
10    continue
100   format(2x,'x, T =', 2(4x,e12.5))
      stop
      end
      
       x, T =     0.00000E+00    -0.10000E+01
 x, T =     0.20000E-01    -0.99924E+00
 x, T =     0.40000E-01    -0.99689E+00
 x, T =     0.60000E-01    -0.99293E+00
 x, T =     0.80000E-01    -0.98738E+00
 x, T =     0.10000E+00    -0.98023E+00
 x, T =     0.12000E+00    -0.97149E+00
 x, T =     0.14000E+00    -0.96114E+00
 x, T =     0.16000E+00    -0.94920E+00
 x, T =     0.18000E+00    -0.93566E+00
 x, T =     0.20000E+00    -0.92053E+00
 x, T =     0.22000E+00    -0.90380E+00
 x, T =     0.24000E+00    -0.88547E+00
 x, T =     0.26000E+00    -0.86554E+00
 x, T =     0.28000E+00    -0.84401E+00
 x, T =     0.30000E+00    -0.82089E+00
 x, T =     0.32000E+00    -0.79617E+00
 x, T =     0.34000E+00    -0.76985E+00
 x, T =     0.36000E+00    -0.74194E+00
 x, T =     0.38000E+00    -0.71243E+00
 x, T =     0.40000E+00    -0.68132E+00
 x, T =     0.42000E+00    -0.64861E+00
 x, T =     0.44000E+00    -0.61431E+00
 x, T =     0.46000E+00    -0.57841E+00
 x, T =     0.48000E+00    -0.54091E+00
 x, T =     0.50000E+00    -0.50181E+00
 x, T =     0.52000E+00    -0.46112E+00
 x, T =     0.54000E+00    -0.41882E+00
 x, T =     0.56000E+00    -0.37494E+00
 x, T =     0.58000E+00    -0.32945E+00
 x, T =     0.60000E+00    -0.28237E+00
 x, T =     0.62000E+00    -0.23368E+00
 x, T =     0.64000E+00    -0.18340E+00
 x, T =     0.66000E+00    -0.13151E+00
 x, T =     0.68000E+00    -0.78024E-01
 x, T =     0.70000E+00    -0.22930E-01
 x, T =     0.72000E+00     0.33769E-01
 x, T =     0.74000E+00     0.92074E-01
 x, T =     0.76000E+00     0.15199E+00
 x, T =     0.78000E+00     0.21351E+00
 x, T =     0.80000E+00     0.27664E+00
 x, T =     0.82000E+00     0.34138E+00
 x, T =     0.84000E+00     0.40772E+00
 x, T =     0.86000E+00     0.47568E+00
 x, T =     0.88000E+00     0.54526E+00
 x, T =     0.90000E+00     0.61646E+00
 x, T =     0.92000E+00     0.68928E+00
 x, T =     0.94000E+00     0.76374E+00
 x, T =     0.96000E+00     0.83985E+00
 x, T =     0.98000E+00     0.91767E+00
 x, T =     0.10000E+01     0.99756E+00
Hosted by www.Geocities.ws

1