function energia

k=1.48;%Constantes de la funcion Weibull%
c=7.7;

v=1:25;%Vector de valores de velocidad en m/s%

for i=1:25%Vector P de valores de la funcion Weibull%
p(i)=(k/c)*(v(i)/c)^(k-1)*exp(-(v(i)/c)^k);
end

plot(v,p)%Representacion funcion Weibull%

for v=1:25%Creacion tabla con valores de velocidad y de frecuencia de la funcion Weibull%
wei='(k/c)*(v/c)^(k-1)*exp(-(v/c)^k)';
tabla(v,1)=v;
tabla(v,2)=eval(wei);
end

%Declaracion curvas de velocidad para cada modelo%
g87=[0,0,0,83.1,182.2,335,544.3,818.9,1154.5,1501,1769.8,1918.3,1977,1994.7,1999,1999.8,2000,2000,2000,2000,2000,2000,2000,2000,2000]*0.001;
s77=[0,0,0,44,129,241,396,594,846,1100,1318,1467,1500,1500,1500,1500,1500,1500,1500,1500]*0.001;
g58=[0,0,9.7,33.5,78.5,145.3,236.6,358.8,513.1,669.7,778.3,828.1,844.4,848.7,849.7,849.9,850,850,850,850,850]*.001;
e80=[0,0,8,59,140,264,439,670,965,1291,1630,1670,1670,1670,1670,1670,1670,1670,1670,1670]*0.001;

potencia_g87=0;%Inicializacion vectores potencia%
potencia_s77=0;
potencia_g58=0;
potencia_e80=0;

for i=4:24%Bucle calculo energia generada%
    integral_g87=g87(i)+(g87(i+1)-g87(i))*((0.5*g87(i+1)^2-(i+1))-(0.5*g87(i)^2-i));%integral curva potencia%
    pot_g87(i)=integral_g87*tabla(i,2);%potencia calculada%
    potencia_g87=potencia_g87+pot_g87(i); %sumatorio con la anterior%
end

for i=4:19%Bucle calculo energia generada%
    integral_s77=s77(i)+(s77(i+1)-s77(i))*((0.5*s77(i+1)^2-(i+1))-(0.5*s77(i)^2-i));%integral curva potencia%
    pot_s77(i)=integral_s77*tabla(i,2);%potencia calculada%
    potencia_s77=potencia_s77+pot_s77(i); %sumatorio con la anterior%
end

for i=3:20%Bucle calculo energia generada%
    integral_g58=g58(i)+(g58(i+1)-g58(i))*((0.5*g58(i+1)^2-(i+1))-(0.5*g58(i)^2-i));%integral curva potencia%
    pot_g58(i)=integral_g58*tabla(i,2);%potencia calculada%
    potencia_g58=potencia_g58+pot_g58(i); %sumatorio con la anterior%
    
end
for i=3:19%Bucle calculo energia generada%
    integral_e80=e80(i)+(e80(i+1)-e80(i))*((0.5*e80(i+1)^2-(i+1))-(0.5*e80(i)^2-i));%integral curva potencia%
    pot_e80(i)=integral_e80*tabla(i,2);%potencia calculada%
    potencia_e80=potencia_e80+pot_e80(i); %sumatorio con la anterior%
end

%energia anual generada, resultado de multiplicar potencia por tiempo%
energia_g87=potencia_g87*24*365
energia_s77=potencia_s77*24*365
energia_g58=potencia_g58*24*365
energia_e80=potencia_e80*24*365

