What modifications to be done to find the minimum

function Y = myFunc(P,P0)

N = 2;

Pa = 9;

k = floor((P0*N/Pa));

Y = F(P)+k*F(N*P0-P*k);

function Fp = F(P)

R = 3;

beta = 8;

Fp = 1 - exp(-((2^R-1)./P).^(beta/2));

%Step 2: Minimize it within the bounds:

P0 = 9;

[Pi, FVal] = fminbnd(@(x) myFunc(x,P0),0,1.1);

What modifications to be done such that P0 varies from 1 to 10 and P belongs to P1 = [0, Pb]∪ ((N*P0)/(k+1)).