5 * Revision 1.1.1.1 1996/04/01 15:02:24 mclareni
10 SUBROUTINE DLSQP1(N,X,Y,A0,A1,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)
54 A0=(SY*SXX-SXY*SX)/DET
61 IF(IFAIL .EQ. 0 .AND. N .GT. 2)
62 1 SD=SQRT(MAX(R0,SYY-A0*SY-A1*SXY)/(N-2))