Révision 210 TrouNoir/TrouNoir.py
| TrouNoir.py (revision 210) | ||
|---|---|---|
| 158 | 158 |
// Initial version |
| 159 | 159 |
flux_int=2.*planck/c2*pow(nu_em,3)/(exp(planck*nu_em/(k*temp_em))-1.); |
| 160 | 160 |
// Version with log used |
| 161 |
//flux_int=2.*exp(lplanck-lc2+3.*log(nu_em))/(exp(exp(lplanck-lk+log(nu_em/temp_em)))-1.);
|
|
| 161 |
flux_int=2.*exp(lplanck-lc2+3.*log(nu_em))/(exp(exp(lplanck-lk+log(nu_em/temp_em)))-1.); |
|
| 162 | 162 |
flux_int*=pow(rf,3)*b*db*h; |
| 163 | 163 |
flx+=flux_int; |
| 164 | 164 |
} |
Formats disponibles : Unified diff