function [p,x] = xvsp(delta)
global P_VALUE

p = 1:delta:2;
n= max(size(p));

for i=1:n
  P_VALUE = p(i);
  x(i) = fzero('rho_prime',0);
end
