1 2 3 4 5 6 7 8 9 10
| function P = plot_densite_probabilite(n1,n2,n3)%,c2,c3)
c2=2.0;
c3=2.0;
L=3.00;
tt=4*gamma((n1+n2)*L)*gamma((n1+n3)*L)*(((n1)^(2*n1*L))*((n2)^(n2*L))*((n3)^(n3*L))/((gamma(n1*L)^2)*gamma(n2*L)*gamma(n3*L)));
t=[0.1 0.2 1];
P =tt*((c2^(2*n2*L)*(1-t)^(2*n1*L-1)/((n1*(1-t)^2+n2*c2^2)^(L*(n1+n2)))+(((1/c2)^(2*n1*L))*(1-t)^(2*n2*L-1))/((n2*(1-t)^2 + n1/(c2^2))^(L*(n1+n2))))*sous_integral(n1,n3,c3,t)+(c3^(2*n3*L)*(1-t)^(2*n1*L-1)/((n1*(1-t)^2+n3*c3^2)^(L*(n1+n3)))+(((1/c3)^(2*n1*L))*(1-t)^(2*n3*L-1))/((n3*(1-t)^2 + n1/(c3^2))^(L*(n1+n3))))*sous_integral(n1,n2,c2,t));
plot(t,P) |
Partager