5 * Revision 1.1.1.1 1996/04/01 15:02:24 mclareni
10 SUBROUTINE DLSQPM(N,X,Y,M,A,SD,IFAIL)
11 #if !defined(CERNLIB_DOUBLE)
14 CALL MTLPRT(NAME,'E201',
15 +'not available on this machine - see documentation')
19 #if defined(CERNLIB_DOUBLE)
20 IMPLICIT DOUBLE PRECISION (A-H,O-Z)
22 PARAMETER (IDIM = 21, R0 = 0)
24 DIMENSION X(*),Y(*),A(0:*),B(IDIM,IDIM),XY(0:IDIM)
28 IF(N .LE. 1 .OR. M .LT. 0 .OR. M1 .GT. IDIM .OR. M1 .GT. N) THEN
31 XY(0)=DVSUM(N,Y(1),Y(2))
33 SYY=DVMPY(N,Y(1),Y(2),Y(1),Y(2))
65 CALL DSEQN(M1,B,IDIM,IFAIL,1,A)
74 SD=SQRT(MAX(R0,SD)/(N-M1))
77 CALL DVSET(M1,R0,A(0),A(1))