Open the model and define a vector of Mu values.

open_system('vdp');
Mu_Values = [0.5:0.25:5];
MuVal_length = length(Mu_Values);

Using Mu_Values, initialize an array of Simulink.SimulationInput objects. To preallocate the array, a loop index is made to start from the largest value.

for i = MuVal_length:-1:1
    in(i) = Simulink.SimulationInput('vdp');
    in(i) = in(i).setBlockParameter('vdp/Mu',...
        'Gain',num2str(Mu_Values(i)));
end

Simulate the model using parsim. Set to 'RunInBackground' to enable the use of command prompt while the simulations are running and to create an array of Simulink.Simulation.Future objects.

Future = parsim(in,'RunInBackground','on');

Use the wait method on Future to block the execution.

Ok = wait(Future)