9.3.10 “Isocaps” sample

C++ code

mglData a(61,51,40);
a.Modify("-2*((2*x-1)^2 + (2*y-1)^2 + (2*z-1)^4 - (2*z-1)^2 - 0.1)");
gr->Rotate(40,60);
gr->Light(true);
gr->CutMin = mglPoint(0,-1,-1); gr->CutMax = mglPoint(1,0,1.1);
gr->Box();
gr->Surf3(-1,a);
gr->ContF3(a,'x',-1); gr->ContF3(a,'y',-1);
gr->ContF3(a,'z', 0); gr->ContF3(a,'z',39);

MGL code

new a 61 51 40
modify a '-2*((2*x-1)^2 + (2*y-1)^2 + (2*z-1)^4 - (2*z-1)^2 - 0.1)'
rotate 40 60
light on
cut 0 -1 -1 1 0 1.1
box
surf3 a -1
contf3 a 'x' -1
contf3 a 'y' -1
contf3 a 'z' 0
contf3 a 'z' 39

Pure C code

HMDT a = mgl_create_data_size(61,51,40);
mgl_data_modify(a,"-2*((2*x-1)^2 + (2*y-1)^2 + (2*z-1)^4 - (2*z-1)^2 - 0.1)",0);
mgl_rotate(gr,40.,60.,0.);
mgl_set_light(gr,1);
mgl_set_cut_box(gr,0.,-1.,-1.,1.,0.,1.1);
mgl_box(gr,1);
mgl_surf3_val(gr,-1.,a,0);
mgl_contf3(gr,a,'x',-1, "", 7);
mgl_contf3(gr,a,'y',-1, "", 7);
mgl_contf3(gr,a,'z', 0, "", 7);
mgl_contf3(gr,a,'z',39, "", 7);
mgl_delete_data(a);

Fortran code

a = mgl_create_data_size(61,51,40);
call mgl_data_modify(a,"-2*((2*x-1)^2 + (2*y-1)^2 + (2*z-1)^4 - (2*z-1)^2 - 0.1)",0);
call mgl_rotate(gr,40.,60.,0.)
call mgl_set_light(gr,1)
call mgl_set_cut_box(gr,0.,-1.,-1.,1.,0.,1.1);
call mgl_box(gr,1)
call mgl_surf3_val(gr,-1.,a,'')
call mgl_contf3(gr,a,'x',-1, '', 7);
call mgl_contf3(gr,a,'y',-1, '', 7);
call mgl_contf3(gr,a,'z', 0, '', 7);
call mgl_contf3(gr,a,'z',39, '', 7);
call mgl_delete_data(a)

Python

a = mglData(61,51,40);
a.Modify("-2*((2*x-1)^2 + (2*y-1)^2 + (2*z-1)^4 - (2*z-1)^2 - 0.1)");
gr.Rotate(40,60);   gr.Light(True);
gr.SetCutBox(0,-1,-1,1,0,1.1);    gr.Box();
gr.Surf3(-1,a);
gr.ContF3(a,"x",-1);    gr.ContF3(a,"y",-1);
gr.ContF3(a,"z",0);     gr.ContF3(a,"z",39);