
Re: tubes program not working in version 9
Posted:
Dec 8, 2012 1:28 AM


On 12/7/2012 12:37 AM, Roger Bagula wrote: > The problem seems to be sometimes it works but mostly it doesn't. > The tubes program was written for an earlier version > ( works in version 5 I think) by Mark McClure > and has worked fine for literally years. > I haven't got a clue what has gone wrong. > ( technically these are called channel surfaces or the like). > > TubePlotFrenet[curve_List, {var_, min_, max_}, radius_, opts___] := > Module[{tangent, unitTangent, normal, unitNormal, biNormal}, > tangent = D[curve, t]; > unitTangent = tangent/Sqrt[tangent.tangent]; > normal = D[unitTangent, t]; > unitNormal = normal/Sqrt[normal.normal]; > biNormal = Cross[unitTangent, unitNormal]; > ParametricPlot3D[ > curve + radius Cos[s] unitNormal + radius Sin[s] biNormal // > Evaluate, {var, min, max}, {s, 0, 2 \[Pi]}, opts]] > > r = Cos[q0*t] + 2; > x = r* Cos[p0*t]; > y = r *Sin[p0*t]; > z = Sin[q0*t]; > w0 = 8*{x, y, z} > h = TubePlotFrenet[w0, {t, 0, 2 \[Pi]}, 1, Axes > None, > Boxed > False, PlotPoints > {64, 16}, > ColorFunction > "CandyColors", MeshFunctions > {#3 &}] > > My other tubes program which came from an answer in this group > still works. >
Noticed few problems:
tangent = D[curve, t];
when it should be
tangent = D[curve, var];
after all, var is the independent variable symbol being passed and is the one build in.
Also, there are no numerical values for the q0 and p0. When I wrote
p0=1; q0=1; and fixed the above issue with var vs. t, then it worked.
Nasser

