global r

fob=[];
rr=[];
for i=1:20
    r = 10.^((i-10)/6);
    [xo,Ot,nS]=powell('penal',[0 0 0],0,0,[-5 -5 -5],[5 5 5]);
    fob=[fob; Ot];
    rr=[rr; r];
end

plot(rr,fob);
