geom3d

deff("[z]=surface(x,y)","z=sin(x)*cos(y)")

t=%pi*(-10:10)/10;

// 曲面の3次元プロット

fplot3d(t,t,surface,35,45,"X@Y@Z")

// ここで, (t,t,sin(t).*cos(t))は曲面上の曲線で,

// geom3d および xpolyにより描画される

[x,y]=geom3d(%pi/2,0,surface(%pi/2,0))

my_plot_desc          = "plot3d and use of geom3d";

  r = (%pi):-0.01:0;

x = r.*cos(10*r);

y = r.*sin(10*r);

deff("[z]=Surf(x,y)","z=sin(x)*cos(y)");

t=%pi*(-10:10)/10;

// 曲面をプロット

fplot3d(t,t,Surf,35,45,"X@Y@Z",[19,2,3]);

// 3次元グラフィック上に2次元グラフィックを追加

z=(sin(x).*cos(y));

[x1,y1] = geom3d(x,y,z);

xpoly(x1,y1,"lines");

// いくつかのプロット用パラメータを調整

BackgroundColorId           = color(70,174,255);

current_axe                 = gca();

plot_3d                     = current_axe.children(2);

plot_3d.hiddencolor         = 32;

polyline                    = current_axe.children(1)';;

polyline.foreground         = 8;

current_axe.rotation_angles = [70,47];

current_axe.background      = BackgroundColorId;

// 2番目の2次元グラフィック

[x1,y1] = geom3d([0,0],[0,0],[5,0]);

xsegs(x1,y1);

xstring(x1(1),y1(1),"The point (0,0,0)");

xtitle(my_plot_desc," "," "," ");