![]() |
XLPack for Matplotlib Experimental Version 0.1
例題集
|
三角格子上に等高線を描く.
Sub Test_Plt_Tricontour()
Const N = 256, Nlev = 7
Dim X(N - 1) As Double, Y(N - 1) As Double
Dim Z(N - 1) As Double, Zmin As Double, Zmax As Double
Dim Levels(Nlev - 1) As Double
Dim I As Integer, J As Integer
'-- Make data
For I = 0 To N - 1
X(I) = -3 + 6 * Rnd()
Y(I) = -3 + 6 * Rnd()
Z(I) = (1 - X(I) / 2 + X(I) ^ 5 + Y(I) ^ 3) * Exp(-X(I) ^ 2 - Y(I) ^ 2)
If Z(I) < Zmin Then Zmin = Z(I)
If Z(I) > Zmax Then Zmax = Z(I)
Next
For I = 0 To Nlev - 1
Levels(I) = Zmin + I * (Zmax - Zmin) / 6
Next
'-- Plot
Call Plt.Plot(N, X(), Y(), "o", "markersize=2, color='lightgrey'")
Call Plt.Tricontour(N, X(), Y(), Z(), Levels())
Call Plt.Show
End Sub