function [DFluxRad]=FlammeSolide(Parametres,Entrees)

    %Parametres=[L] et Entrees=[EmissF,TempeF,DiameF,HauteF]

    h=2*Entrees(4)/Entrees(3);
    
    S=2*Parametres(1)/Entrees(3);
    
    A=(h^2+S^2+1)/2/S;
    
    B=(1+S^2)/2/S;
    
    DFluxRad=(B-1/S)/pi/sqrt(B^2-1)*atan(sqrt(((B+1)*(S-1))/((B-1)*(S+1))))...
            -(A-1/S)/pi/sqrt(A^2-1)*atan(sqrt(((A+1)*(S-1))/((A-1)*(S+1))));
    
    DFluxRad=DFluxRad*Entrees(1)*5.67E-8*Entrees(2)^4;
    
end