function [kW] = PowerCurve( Wind )
 
x = [3,4,5,6,7,8,9,10,11,12,13]';
y = [0,53,496,1096,1882,2857,3968,5102,6100,6803,7200]';
B = [x.^3, x.^2, x, x.^0];
A = inv(B'*B)*B'*y;
 
kW = 0*Wind;
for i=1:length(Wind)
    if(Wind(i) < 3)
        kW(i) = 0;
    elseif(Wind(i) > 13)
        kW(i) = 7200;
    else
        kW(i) = [Wind(i)^3, Wind(i)^2, Wind(i), 1]*A;
    end
plot(x,y,'r',Wind,kW,'r.')
end