3D Plots

* Declarations

SETS

planes Scenarios /Uno,Dos,Helix/

u /u1*u50/

v /v1*v20/

w /w1*w100/

axisitem /x,y,z/

;

PARAMETER

uval(u)

vval(v)

wval(w)

InterlockingTori(planes,u,v,axisitem)

HelixPlot(planes,w,v,axisitem)

;

* Tori Calculation

uval(u) = -3.14159 + ((ord(u)-1)/(card(u)-1))*2*3.14159;

vval(v) = -3.14159 + ((ord(v)-1)/(card(v)-1))*2*3.14159;

InterlockingTori("uno",u,v,"x")

= cos(uval(u))+.5*cos(uval(u))*cos(vval(v));

InterlockingTori("uno",u,v,"y")

= sin(uval(u))+.5*sin(uval(u))*cos(vval(v));

InterlockingTori("uno",u,v,"z")

= .5*sin(vval(v));

InterlockingTori("dos",u,v,"x")

= 1+cos(uval(u))+.5*cos(uval(u))*cos(vval(v));

InterlockingTori("dos",u,v,"y")

= .5*sin(uval(u));

InterlockingTori("dos",u,v,"z")

= sin(uval(u))+.5*sin(uval(u))*cos(vval(v));

* Tori Plotting

$setglobal gp_key "no"

$setglobal gp_style "lines"

$setglobal gp_grid "no"

$setglobal gp_xlabel "Earth"

$setglobal gp_ylabel "Water"

$setglobal gp_zlabel "Air"

$setglobal gp_title "Fire"

$libinclude gnuplotxyz InterlockingTori x y z

$setglobal gp_xlabel "no"

$setglobal gp_ylabel "no"

$setglobal gp_zlabel "no"

* Helix Calculation

wval(w) = ((ord(w)-1)/(card(w)-1))*31.4159;

vval(v) = ((ord(v)-1)/(card(v)-1))*6.28319;

HelixPlot("Helix",w,v,"x")

= (1-0.1*cos(vval(v)))*cos(wval(w));

HelixPlot("Helix",w,v,"y")

= (1-0.1*cos(vval(v)))*sin(wval(w));

HelixPlot("Helix",w,v,"z")

= 0.1*(sin(vval(v))+wval(w)/1.7-10);

* Helix Plotting

$setglobal gp_style "lines"

$setglobal gp_title "Plot of Helix"

$libinclude gnuplotxyz HelixPlot x y z