sample the waveform with nsample points for fft

mc0 = -0.5; jlr0 = -0.8; M = 0.99; alpharange = 6.5; %// degrees alpharange2 = 180; %// degrees

mcalpha = (mc0 - (1 / M) - 1) * cosd(alpharange) + jlr0 * sind(alpharange) + (1/M) + 1; jlralpha = ((-mc0) + (1 / M) + 1) * sind(alpharange) + jlr0 * cosd(alpharange);

alpha1 = linspace(0, alpharange); jlralph = ((-mc0) + (1 / M) + 1) * sind(alpha1) + jlr0 * cosd(alpha1);

alpha2 = linspace(alpharange, alpharange2); jlrgamma = ((-mcalpha) + (1 / M) - 1) * sind(alpha2 - alpharange) + jlralpha * cosd( alpha2 - alpharange); plot([alpha1, alpha2], [jlralph, jlrgamma])