clear clc clf A = 50; w = 2*pi; t = 0:.01:5*pi; R = 100 + t; x = A .* sin(w*t); th = asin(x./R); y = - R .* cos(th); plot(x, y)