# ejemplo resolucion grafica del problema 5 practico 2
# usando gnuplot

set logscale xy

#constants 
h = 6.62606876e-34   # Planck's constant
c = 2.99792458e8   # speed of light
k = 1.3806503e-23   # Boltzmann's constant

# defino formula de planck en frecuencia x, temp t
Planck(x,t) = (2*h*x**3/c**2)/(exp(h*x/(k*t))-1)

#factor reflexion del asteroide, poner valor correcto
f = 1e-2

set xlabel "log frecuencia"
set ylabel "log flujo recibido"
set grid

plot [1e10:1e16][1e-18:1e-8] \
 Planck(x,5770)*f lt 1 lw 2 title "reflejada",  Planck(x,270) lt 3 lw 2 title "emitida"

 pause -1
 
 # va grafica en png
 set terminal pngcairo size 1600, 1200
 set output "prob5prac2.png"
 replot