range rx = [Au_NPs_Sc+Ta_commercial_Copy.txt]Au_NPs_Sc+Ta_commercial_Copy!col(A);  // Column A range
range ry = [Au_NPs_Sc+Ta_commercial_Copy.txt]Au_NPs_Sc+Ta_commercial_Copy!col(F);  // Column B range

sample_rate = 2000;  // Sampling frequency
seconds = 1;  // Number of seconds
num_samples = sample_rate * seconds;  // Number of samples
f = 10;  // Fundamental frequency
nonlin_deg = 0.25;  // Degree of non-linearity
nonlin_phi = -pi / 4;  // Skew in non-linearity

rx = data(0, seconds, seconds / (num_samples-1));  // Generate the time vector, and put it to column A

// The following is the equation 9 in [1]
factor = sqrt(1-nonlin_deg*nonlin_deg); 
num = nonlin_deg * sin(nonlin_phi) / 1 + factor;
numerator = num + sin(2*pi*f*rx);
denominator = 1 - nonlin_deg * cos(2*pi*f*rx + nonlin_phi);
ry = factor * (numerator / denominator);

ry = ry  - sin(2*pi*0.1*rx);  // Add part of cycle as a trend
white_noise ry;  // Add some noise