How can I add bootstrapping to a script?

I have 387 electric vehicles and 36 households. I have scripts that calculate how much the maximum power demand increase with the deployment of EVs.

There will be A LOT of possible combinations, like 9.4*10^135. I have 36 different scripts. If a household is denoted X, each script calculates:

How does the total maximum power demand of X households change with the introduction of EVs.

The output from the scripts are a cell arrary (Ivalues), containing one matrix for each possible combination of households. The matrices in turn contain an I.value for all the possible combinations of EVs. Every vehicle has a unique ID-nr, so which EVs that are combined is listed in IDnr.mat (has the same structure as Ivalues.mat).

The example in Ivalues.mat considers 4 households and 4 EVs, 3 households has been combined (with 4 households there is 4 possible combinations).

I want to conduct a bootstrap estimation to make my scripts computationally feasible, I want to integrate it into H3.mat (I have H4, H5... H36) to make the scripts computationally feasible for larger numbers of vehicles and households.

My end-goal is to get samples that is representative for the data.