function sample_p = lhs_setup(N , lb, ub)
    n=length(ub);
    d=ub-lb;
    sample_p = lhsdesign(N,n);

    for i=1:n
        sample_p(:,i) = sample_p(:,i)*d(i)+lb(i);
    
    end
end
