How to make the harmonic complex tones sound comfortable?
In Octave: plot(y1);

audiowrite will effectively clip floating point sample values to range -1 ... 1. This will sound unpleasant. You can divide your data by its maximum absolute value, before writing the audio file, to normalize it to that range.