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
|
void matrice_F(float ***F,int M,int N, int O)
{
float ***Mt=alloue_tableau(vs,vs,s);
matrice_Mt(Mt,vs,vs,s);
float **P=alloue_mtrice(vs,vs);
matrice_Mr(P,vs,vs);
int i,j,k;
for (k = 0; k < O; k++)
{
for (i = 0; i < M; i++)
{
for (j = 0; j < N; j++)
{
if (k<dv)
{
flash[i][j][k]=0;
}
else
{
F[i][j][k]=(2*Q)/(E*sqrt(pow(M_PI,3)*Mt[i][j][k]))*(1/(pow(R,2)+8*D*Mt[i][j][k]))*exp((-2*pow(P[i][j],2))/(pow(R,2)+8*D*Mt[i][j][k]));
}
}
}
}
libre_tab(Mt,vs,vs);
libre_matrice(P);
} |
Partager