switchtodesign;
closeall;
clear;

Ntrial = 100;
mu = matrix(1, Ntrial);
sigma = mu;

for (ii = 1:Ntrial) {
    switchtodesign;
    setnamed("NOISE_1", "automatic seed", 0);
    setnamed("NOISE_1", "seed", 0);
    
    run;
    
    signal = getresult("OSC_1", "signal");
    signal = signal.getattribute(signal.getattribute);
    avg_sig = mean(signal);
    var_sig = sum((signal - avg_sig)^2)/(length(signal) - 1);
    sigma_sig = sqrt(var_sig);
    
    mu(ii) = avg_sig;
    sigma(ii) = sigma_sig;
}

histc(mu, 10);