function [times,sigmatrue] = maketimes(jacobian, exactslowmodel, ncols,nrows,method,seed,scale,corlen)
  
rand('state',seed)


// Add some noise to the travel times
if method == 'N'
  noise = scale * rand(nrows,1,'normal'); 
else
  noise = scale * rand(nrows,1); 
end

if corlen ~= 0
  times = smooth(noise,corlen)' +    jacobian * exactslowmodel';
else
  times = noise + jacobian * exactslowmodel';
end

noisemean = mean(noise);
sigmatrue = sqrt((noise - noisemean) .* (noise - noisemean))

// sigmatrue = sqrt(var(noise));