function res = RD_c_steps(num_trials,err)
%

num_err = length(err);
res = zeros(num_err,4);
del_theta = 6*pi/180;
theta = [0:del_theta:2*pi-del_theta];
locs = 9.554*[cos(theta)',sin(theta)'];

for e = 1:num_err
    locs_n = locs + err(e)*(0.5-rand(60,2));
    c_vals = zeros(num_trials,1);
    for t = 1:num_trials
        [e num_err t num_trials]
        [a,b,c] = RD_linear_dist(60,0.0001,1200,4,4,0.25,0.0625,0.05,...
            1,locs);
        c_vals(t) = c;
    end
    res(e,1) = mean(c_vals);
    res(e,2) = var(c_vals);
    res(e,3) = min(c_vals);
    res(e,4) = max(c_vals);
end
