1 : C LEVMA1 SOURCE FANDEUR 10/07/07 21:15:16 6708 2 : SUBROUTINE LEVMA1(mle) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : -INC CCOPTIO 6 : -INC SMLENTI 7 : -INC SMLREEL 8 : SEGMENT MLREE4.MLREEL,MLREE5.MLREEL,MLREE6.MLREEL,MLREE8.MLREEL, 9 : + MLREE9.MLREEL,MLRE11.MLREEL,MLRE12.MLREEL,MLRE13.MLREEL 10 : 11 : mlenti = mle 12 : 13 : mlreel = lect(10) 14 : mlree1 = lect(1) 15 : mlree2 = lect(2) 16 : mlree3 = lect(3) 17 : mlree4 = lect(4) 18 : mlree5 = lect(5) 19 : mlree6 = lect(6) 20 : mlree8 = lect(8) 21 : mlree9 = lect(9) 22 : mlre11 = lect(11) 23 : segact mlreel*mod,mlree1,mlree2,mlree3*mod,mlree4,mlree5,mlree6, 24 : + mlree8*mod,mlree9*mod,mlre11*mod 25 : 26 : chi2 = prog(1) 27 : alamda = prog(2) 28 : ndata = mlree1.prog(/1) 29 : ma = mlree3.prog(/1) 30 : jg = ma*ma 31 : segini mlre12 32 : jg = ma 33 : segini mlre13 34 : 34 : call levma2(mlree1.prog,mlree2.prog,mlree6.prog, ndata, 35 : + mlree3.prog,ma,mlree4.prog,mlree5.prog, 36 : + chi2,alamda,mlree8.prog,mlree9.prog,mlre11.prog,mlre12.prog, 37 : + mlre13.prog) 38 : 39 : mlreel = lect(10) 40 : prog(1) = chi2 41 : prog(2) = alamda 42 : 43 : segdes mlree4,mlree5 44 : segsup mlre12,mlre13 45 : 46 : RETURN 47 : END 48 : 49 :
© Cast3M 2003 - All rights reserved.
Disclaimer