есть функция вида
proc mini_sun jd {
# Time in julian centuries:
set t [expr ($jd-2451545.0)/36525.0]
set p2 6.283185307; set coseps 0.91748; set sineps 0.39778
set m [expr $p2 * [frac [expr 0.993133 + 99.997361 * $t]]]
set dl [expr 6893.0 * sin($m) + 72.0*sin(2*$m)]
set l [expr $p2 * [frac [expr 0.7859453 + $m/$p2 + (6191.2*$t+$dl)/1296.e3]]]
set sl [expr sin($l)]
set x [expr cos($l)]; set y [expr $coseps*$sl]; set z [expr $sineps*$sl]
set rho [expr sqrt(1.0-$z*$z)]
set sundec [expr (360.0/$p2)*atan($z/$rho)]
set sunra [expr (48./$p2)*atan($y/($x+$rho))]
if {$sunra<0} {set sunra [expr $sunra+24.0]}
return [list $sunra $sundec]
}
как выдать на печать что оно вычисляет??
Спасибо!
Форум —
Development

