# -*- coding: utf-8 -*-

import matplotlib.pyplot as plt
import numpy as np
import scipy.integrate as sinteg

plt.figure(1)
plt.clf()

datos=np.genfromtxt('sun_AM0.dat',float)

longitud=datos[:,0]*10 #Angstroms
flujo=datos[:,1]*10**10

plt.plot(np.log10(longitud), np.log10(flujo))   

plt.xlabel('x log10[Longitud de Onda (Angstrom)]')
plt.ylabel('y log10[Flujo (CGS)]')
plt.title('Espectro del Sol (Flujo vs Longitud de Onda)')
plt.legend()

"""
Dado que x1-x0 = 2 = constante podemos considerarlo como un trazador lineal
por lo que podemos aplicar el metodo del trapecio
"""

contador=0
integral=0
trapecio=0
while contador<1696:
    trapecio=flujo[contador]+flujo[contador+1] #no multiplico porque delta x=2
    integral=integral+trapecio
    contador=contador+1
print integral

integralcalc=sinteg.trapz(flujo,longitud)

@@ -226,17 +226,17 @@
 1%5D*10**1
-0
+1
 %0A%0Aplt.pl
@@ -459,16 +459,17 @@
 1-x0 = 2
+0
  = const
@@ -640,16 +640,17 @@
 rapecio=
+(
 flujo%5Bco
@@ -678,39 +678,44 @@
 r+1%5D
- #no multiplico porque delta x=
+)*10 #dividiendo la altura (20) por 
 2%0A  
@@ -825,8 +825,156 @@
 ngitud)%0A
+%0A#superficie esf%C3%A9rica a la distancia en que se encuentra la Tierra 2.81e23 m2%0Actesolar=integral*10**-11 # en W/m2%0A%0Aluminosidadtot= ctesolar*2.81e23%0A

pregunta2 definitiva