June 24 2002  11:02:19.922 AM
  
 PMETHOD
  
 Solve the two-point boundary value problem
  
   - d/dX (P dU/dX) + Q U  =  F
  
 on the interval [-1,1], with
 U(-1) = U(1)=0.
  
 The P method is used, which represents U as
 a weighted sum of orthogonal polynomials.
  
  
 Highest degree polynomial to use is  2
 Number of points to be used for output =  10
  
 Basis function orthogonality test:
  
    i   j     b(i,j)/a(i)
  
  
   0   0   1.00000    
   0   1  0.223517E-07
   0   2 -0.335276E-07
  
   1   0  0.372529E-07
   1   1   1.00000    
   1   2  0.931323E-08
  
   2   0 -0.122236E-06
   2   1  0.203727E-07
   2   2   1.00000    
  
 Representation of solution:
  
 Basis function coefficients:
  
   0  0.954930    
   1 -0.189175E-08
   2 -0.220787    
  
  
        X     Approximate Solution
  
  -1.00000       0.00000    
 -0.800000      0.308802    
 -0.600000      0.588547    
 -0.400000      0.809559    
 -0.200000      0.950645    
   0.00000      0.999087    
  0.200000      0.950645    
  0.400000      0.809559    
  0.600000      0.588547    
  0.800000      0.308802    
   1.00000       0.00000    
  
 
Comparison of computed and exact solutions:
 
    X        U computed    U exact     Difference
 
 -1.0000   0.00000     -0.437114E-07 -0.437114E-07
 -0.8000  0.308802      0.309017      0.214905E-03
 -0.6000  0.588547      0.587785     -0.761330E-03
 -0.4000  0.809559      0.809017     -0.542462E-03
 -0.2000  0.950645      0.951057      0.411153E-03
  0.0000  0.999087       1.00000      0.912905E-03
  0.2000  0.950645      0.951057      0.411153E-03
  0.4000  0.809559      0.809017     -0.542462E-03
  0.6000  0.588547      0.587785     -0.761330E-03
  0.8000  0.308802      0.309017      0.214905E-03
  1.0000   0.00000     -0.437114E-07 -0.437114E-07
 
Big L2 error =   0.572599E-03
 
PMETHOD
  Normal end of execution.