#!/usr/bin/gnuplot # teplota_rovnovazna.plt # Teplota pri termodynamicke rovnovaze. # Miroslav Broz (miroslav.broz@email.cz), Jan 28th 2006 AU = 149.6e9; # m, astronomicka jednotka sigma = 5.76e-8; # W/m^2/K^4, Stefanova-Boltzmannova konstanta L_Sun = 3.83e26; # W, zarivy vykon Slunce T(r,A) = ((1-A)*L_Sun/(16*pi*sigma * r**2))**(1./4.) ######################################################################## set xl "vzdalenost od Slunce r / AU" set yl "geometricke albedo A" set zl "teplota T / K" set xr [0.5:5.2] set yr [0:1] set ticslevel 0 set isosamples 9,11 sp T(x*AU,y) w l,\ "