clear vars;

 I = 9;
 n = 2;
 theta = 1;

 alpha = 0.5;

 locPoint = ((I - 1)/theta) +1;
 index = 1;

 for i = 1:locPoint
     xTemp2 = [];
     for n1 = 1+alpha:n+1+alpha
        n12 = n1 - alpha;
        d1 = gamma(n12);
        ind = n12 -1;
        yData = i^ind;
        x1 = (gamma(n12)/d1)*yData;
        xTemp2 = [xTemp2, x1];
     end
     bmat(index,  1:n+1) = xTemp2;
     index = index + 1;
 end

 for k1 = 1: n+1
    for k2 = 1: I
        xValTemp = 1+((k2 - 1)*theta);
        xTemp(k2 , 1) = xValTemp^(k1 - 1);
    end
    finalMat(1:length(xTemp), k1) = xTemp;
 end

 [locr, locc] = find(finalMat == I);
 xMat = finalMat(1:locr, :);


 weightTemp = xMat;
 finalWeight = (weightTemp'*weightTemp);
 finalWeight = finalWeightweightTemp';
 finalMovingCoeff = bmat*finalWeight;