How do I get some sampling points from a 4 dimentional function at LHS(Latin Hypercube Sampling)?

function scores = rosenbrockfcn(x)
    scores = 0;
    n = size(x, 2);
    assert(n >= 1, 'Given input X cannot be empty');
    a = 1;
    b = 100;
    for i = 1 : (n-1)
        scores = scores + (b * ((x(:, i+1) - (x(:, i).^2)) .^ 2)) + ((a - x(:, i)) .^ 2);
    end

This is Rosenbrock function.

I want to get some samples from a 4 dimentional rosenbrock function.

like [x1 x2 x3 x4] and f(x)

Boundary conditions are [-10 10]

Please let me know. Thank you!