File:Moebius Surface 0.1.png

DescriptionMoebius Surface 0.1.png	A moebius strip parametrized by the following equations: ${\displaystyle x=\cos u+v\cos {\frac {nu}{2}}\cos u}$ ${\displaystyle y=\sin u+v\cos {\frac {nu}{2}}\sin u}$ ${\displaystyle z=v\sin {\frac {nu}{2}}}$, where n=0.1. This plot is part of a series depicting n for 0 to 1. See Möbius Strip for the rest of the series
Date	24 June 2007
Source	Self-made, with Mathematica 5.1   This diagram was created with Mathematica by n.
Author	Inductiveload
Permission (Reusing this file)	Public domainPublic domainfalsefalse I, the copyright holder of this work, release this work into the public domain. This applies worldwide. In some countries this may not be legally possible; if so: I grant anyone the right to use this work for any purpose, without any conditions, unless such conditions are required by law.
Other versions	Image for use by itself

Mathematical Function Plot
Description	Moebius Strip, 1 half-turn (n=0.1)
Equation	${\displaystyle x=\cos u+v\cos {\frac {nu}{2}}\cos u}$ ${\displaystyle y=\sin u+v\cos {\frac {nu}{2}}\sin u}$ ${\displaystyle z=v\sin {\frac {nu}{2}}}$
Co-ordinate System	Cartesian (Parametric Plot)
u Range	0 .. 4π
v Range	0 .. 0.3

Please be aware that at the time of uploading (15:15, 19 June 2007 (UTC)), this code may take a significant amount of time to execute on a consumer-level computer.

This uses Chris Hill's antialiasing code to average pixels and produce a less jagged image. The original code can be found here.

This code requires the following packages:

<<Graphics`Graphics`

MoebiusStrip[r_:1] =
    Function[
      {u, v, n},
      r {Cos[u] + v Cos[n u/2]Cos[u],
          Sin[u] + v Cos[n u/2]Sin[u],
          v Sin[n u/2],
          {EdgeForm[AbsoluteThickness[4]]}}];

aa[gr_] := Module[{siz, kersiz, ker, dat, as, ave, is, ar},
    is = ImageSize /. Options[gr, ImageSize];
    ar = AspectRatio /. Options[gr, AspectRatio];
    If[! NumberQ[is], is = 288];
    kersiz = 4;
    img = ImportString[ExportString[gr, "PNG", ImageSize -> (
      is kersiz)], "PNG"];
    siz = Reverse@Dimensions[img[[1, 1]]][[{1, 2}]];
    ker = Table[N[1/kersiz^2], {kersiz}, {kersiz}];
    dat = N[img[[1, 1]]];
    as = Dimensions[dat];
    ave = Partition[Transpose[Flatten[ListConvolve[ker, dat[[All, All, #]]]] \
& /@ Range[as[[3]]]], as[[2]] - kersiz + 1];
    ave = Take[ave, Sequence @@ ({1, Dimensions[ave][[#]], 
    kersiz} & /@ Range[Length[Dimensions[ave]] - 1])];
    Show[Graphics[Raster[ave, {{0, 0},
   siz/kersiz}, {0, 255}, ColorFunction -> RGBColor]], PlotRange -> {{0, siz[[
    1]]/kersiz}, {0, siz[[2]]/kersiz}}, ImageSize -> is, AspectRatio -> ar]
    ]

deg = 0.1;
gr = ParametricPlot3D[Evaluate[MoebiusStrip[][u, v, deg]],
      {u, 0, 4π},
      {v, 0, .3},
      PlotPoints -> {249, 5},
      PlotRange -> {{-1.3, 1.3}, {-1.3, 1.3}, {-0.7, 0.7}},
      Boxed -> False,
      Axes -> False,
      ImageSize -> 800,
      PlotRegion -> {{-0.22, 1.15}, {-0.5, 1.4}},
      DisplayFunction -> Identity
      ];
finalgraphic = aa[gr];

Export["Moebius Surface " <> ToString[deg] <> ".png", finalgraphic]