l=[100:10:40000];
i=planck(l);
l=l*1.e-9;
I=trapz(l,i);
x=6.96*1.e8/(149.6*1.e9);  %stosunek promienia slonca do odleglosci Ziemia-Slonce
I=I*x^2;