Plot 创建三维弯管

我有一个问题，创造一个弯曲的空心管。 我看到了下面的帖子，但是程序代码有问题

在这方面有人能帮忙吗？ 我希望创建的曲线有一个很好的曲率作为一个环形结构。在上述情况下，弯曲圆柱体的后曲率在其域中不相同

% Parameters
r=5; l=5; nTheta=100, nL = 20;

theta = linspace(0,2*pi,nTheta+1);
x = r * cos(theta);
y = r * sin(theta);

z = linspace(0,l,nL)';
xshift = repmat( sin(z), 1, nTheta+1); %this is a function of z

X = repmat(x,nL,1) + xshift;
Y = repmat(y,nL,1);
Z = repmat(z, 1, nTheta+1);

% Plot
surf(X,Y,Z)