dyneld.m

tic
format long;
n=length(data(:,1));
     Aeq=ones(1,n);
     a=data(:,1)';
          b=data(:,2');
               c=data(:,3)';
                    l=data(:,4)';
                         u=data(:,5)';
                         m=length(Pd1);
                         A=ones(1,n);
                          K=zeros(1,n);
                         for i=1:m
                             L((i-1)*10+1:i*10)=l;
                               U((i-1)*10+1:i*10)=u;
                                Aeq((i-1)*10+1:i*10)=K;
                                 A1((i-1)*10+1:i*10)=a;
                                  B1((i-1)*10+1:i*10)=b;
                              CC(i,(i-1)*10+1:i*10)=A;
                         end
Beq=0;
H=2*diag(A1);
P1=quadprog(H,B1,Aeq,Beq,CC,Pd1,L,U);
for i=1:m
    PP(i,:)=P1((i-1)*10+1:i*10);
end
for i=1:m*n
       C(i)=A1(i)*P1(i)^2+B1(i)*P1(i);
       end
    Tcost=sum(C)+m*sum(c);
   toc;