Voilà un exp
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
format long e
p1=0.07;
TCmin=30000;
for n1=1:5
    for n2=4:7
        for L1=3:5
            for L2=4:5
                n1
                n2
                L1
                L2
                q11=binocdf(L1,n1,p1)
                q12=binocdf(L2,n1,p1)-binocdf(L1,n1,p1)
                q21=binocdf(L1,n2,p1)
                q22=binocdf(L2,n2,p1)-binocdf(L1,n2,p1)
                 Q=[q11 q12;q21 q22];
                 QQ=eye(2)-Q
                 if det(QQ)==0
                     TC=NaN;
                 else
                     A=[0 1]*inv(QQ)*ones(2,1)
                     B=A-2
                     C=(3*A)+B
                     TC=A+C+B
                   end
                 end
             end
         end
     end
si det(QQ)==0 on ne fait pas le calcul de TC.
Lors de l'execution
une erreur s'affiche
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
 
n1 =
 
     5
 
 
n2 =
 
     5
 
 
L1 =
 
     3
 
 
L2 =
 
     5
 
 
q11 =
 
    9.998866727999994e-001
 
 
q12 =
 
    1.133272000005681e-004
 
 
q21 =
 
    9.998866727999994e-001
 
 
q22 =
 
    1.133272000005681e-004
 
 
QQ =
 
    1.133272000005681e-004   -1.133272000005681e-004
   -9.998866727999994e-001    9.998866727999994e-001
 
Warning: Matrix is close to singular or badly scaled.
         Results may be inaccurate. RCOND = 6.776264e-021.
> In C:\MATLAB6p5\work\vssmc\exp55.m at line 22
le determinant dans ce cas =0 et pourtant il fait le calcul
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
A =
 
   -7.379533931401575e+019
 
 
B =
 
   -7.379533931401575e+019
 
 
C =
 
   -2.951813572560630e+020
 
 
TC =
 
   -4.427720358840945e+020
où est le problème?