1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
pkg load mathgl
mathgl
dat = mathgl.mglData(); # data to for plotting
dat.Create(30,40);
for i=0:30
for j= 0:40
dat.SetVal(1/(1+(i-15)*(i-15)/9.+(j-20)*(j-20)/16.),i,j);
endfor
endfor
gr = mathgl.mglGraph(); # class for plot drawing
gr.SetRanges(0,2,0,2,0,1) # ranges of coordinates
gr.Rotate(50,60) # rotate axis
gr.Light(true) # enable lighting
gr.Surf(dat) # plot surface
gr.Cont(dat,"y") # plot yellow contour lines
gr.Axis() # draw axis
gr.Puts(mglPoint(1,1,1.2),"\\i f = \\dfrac{1}{1+(5x-5)^2+(5y-5)^2}")
gr.WriteFrame("sample.png") # save it
|
