# Octave, (mayby mathlab) syntax

# Projecting a circle onto a cylindar;
# Usedful when joining a small pipe on a bigger pipe

r=160/2
rb=370/2
y=-r:r

axis ("auto", "equal");


plot(y,asin(sqrt(r^2-y.^2)/rb)*rb,y,-(asin(sqrt(r^2-y.^2)/rb)*rb))
print('foo.ps','-deps')
