subroutine llx(p,m,el,lel)                         
      implicit real*8 (a-h,o-z)                          
      dimension p(m)                                     
      common /i721/ n
      common /r721/ x(11)
      data pi /3.14159265d0/
      lel=0
      if(p(2) .gt. 0.d0) then
         sum=0.d0
         do 1 i=1,n
            sum=sum+(x(i)-p(1))**2
    1    continue
         el= -n*0.5*( dlog(pi*2)+dlog(p(2)) ) - 0.5*sum/p(2)
         lel=1 
      end if
      return                                             
      end                                                

1.3. 数式