matpower opf problem

[anshu dahiya]

sir, in my 5 bus system the power is running but not the optimal power
flow. please check what is the problem in program , I have attached my
matpower file

function mpc = add_userfuc(casem5)
mpc.version = '2';
mpc.baseMVA = 100;
%No. Bus Pmin Pmax Marginal Cost

%      bus_i type Pd Qd Vm Va Vmax Vmin gs bs area basekva zone 
mpc.bus =[ 1 3 0 0 1 0 1.06 0.8 1 0 1 100 0.44;
         2 1 300 0 1.05 0 1 0.9 1 0 1 100 0.44;
         3 2 300 0 1.02 0 1 0.8 1 0 1 100 0.44;
         4 2 0 0 1 0 1.02 0.9 1 0 1 100 0.44;
         5 2 0 0 1 0 1.02 0.9 1 0 1 100 0.44];
 mpc.area=[1; 1; 1; 1; 1];
     %busno. pg qg qmax qmin vg mbase status pmax pmin pc1 pc2 qc1min qc1max 
qc2min qc2max rateA rateb ratec apf 
mpc.gen =[ 1 60 0 100 0 1 100 1 100 0 0 100 0 0 0 0 0 0 0 0 0;
           2 60 0 100 0 1 100 1 100 0 0 100 0 0 0 0 0 0 0 0 0;
           3 80 0 100 0 1 100 1 100 0 0 100 0 0 0 0 0 0 1 100 0;
           4 70 0 100 0 1 100 1 100 0 0 100 0 0 0 0 0 0 1 100 0;
           5 65 0 100 0 1 100 0.9 100 0 70 100 0 0 0 0 0 0 1 100 0];
%Frombus Tobus R  X B ratea rateb ratec ratio angle status angmin angmax
mpc.branch = [1 2 0.0281 0.281 0 300 350 355 1 0 1 -360 360;
2 3 0.0108 0.108 0 200 225 235 1 0 1 -360 360;
4 3 0.0297 0.297 0 200 225 235 1 0 1 -360 360;
4 5 0.0297 0.297 0 240 260 262 1 0 1 -360 360;
4 1 0.064 0.64 0 150 160 164 1 0 1 -360 360;
5 1 0.0304 0.304 0 200 230 240 1 0 1 -360 360];


%model stsrtup_cost shutdown_cost c2 c1 c0
mpc.gencost = [2 11.56 0 100 1.089 0.074;
    2 11.56 0 70 1.033 0.089;
    2 11.57 0 75 1.043 0.085;
    2 11.57 0 85 1.041 0.087;
    2 11.57 0 75 1.089 0.086];