5 * Revision 1.1.1.1 1996/04/01 15:02:24 mclareni
10 SUBROUTINE DLSQP2(N,X,Y,A0,A1,A2,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)
60 DET=(FN*SXXXX-SXX**2)*SXX-FN*SXXX**2
62 A2=(SXX*(FN*SXXY-SXX*SY)-FN*SXXX*SXY)/DET
70 IF(IFAIL .EQ. 0 .AND. N .GT. 3)
71 1 SD=SQRT(MAX(R0,SYY-A0*SY-A1*SXY-A2*SXXY)/(N-3))