[PYTHON] Zeichnen Sie eine "Boobs-Kurve" in einem 3D-Diagramm (2).

Es ist mir gelungen, das Beispiel vorerst mit "Zeichnen eines Diagramms mit der Programmiersprache Julia" zu verschieben, aber mit den von meinem eigenen Programm erstellten Daten (in der Sprache Julia) Ich beschloss herauszufinden, wie man ein Diagramm mit zeichnet.

Dieses Mal führte ich ein Experiment durch, um "Brüste" unter Bezugnahme auf "Längenberechnung durch Gleichung für gekrümmte Oberflächen von Brüsten" zu vergrößern. Ich versuchte zu. Das von Julia erhaltene Berechnungsergebnis wird von Matplotlib als Grafik angezeigt. (Bitte erstellen Sie eine Umgebung und verschieben Sie sie)

• Beim Speichern des Quellcodes sollte die Zeichenkodierung utf-8 sein.
• Um das Beispiel zu starten, doppelklicken Sie im Explorer auf Oppai2.jl (Julia-Quellcode).
• Durch Erhöhen von ALPHA im Quellcode werden die "Brüste" erhöht. Durch Erhöhen der BETA wird der "Nippel" verlängert.
• Eine Datei mit dem Namen "Oppai2.png " wird im aktuellen Verzeichnis erstellt.

Hinweis: Informationen zum Erstellen von Umgebungen finden Sie unter "Zeichnen von Diagrammen in der Programmiersprache Julia". Die Schritte 4, 6 und 11 wurden zum 04.07.2014 hinzugefügt. ** Wenn Sie die Umgebung anhand des obigen Artikels vor dem 04.07.2014 erstellt haben, führen Sie bitte die zusätzlichen Schritte 4, 6 und 11 aus. **Vielen Dank.

↓ Beispiel für $ ALPHA = 1 $. "Titten" sind klein. Ziehen Sie die linke Taste, um das 3D-Diagramm zu drehen. Sie können vergrößern / verkleinern, indem Sie die rechte Taste ziehen.

↓ Beispiel für $ ALPHA = 3 $. "Titten" sind groß. Ziehen Sie die linke Taste, um das 3D-Diagramm zu drehen. Sie können vergrößern / verkleinern, indem Sie die rechte Taste ziehen.

Oppai2.jl

# -*- coding: utf-8 -*-

using PyPlot

using PyCall
@pyimport matplotlib.cm as cm
@pyimport matplotlib.font_manager as fm
@pyimport matplotlib.pyplot as plt
@pyimport numpy as np

ALPHA = 1.0
BETA  = 1.0
GAMMA = 1.0 / ALPHA

function bust(x, y)
ALPHA * (6 * exp(-((2 / 3 * abs(x) - 1) ^ 2 + (2 / 3 * y) ^ 2) - 1 / 3 * (2 / 3 * y + 0.5) ^ 3) 
+ BETA * 2 / 3 * exp(-2.818 ^ 11 * ((abs(2 / 3 * x) - 1) ^ 2 + (2 / 3 * y) ^ 2) ^ 2) 
- GAMMA * (2 / 3 * x) ^ 4) / 8
end

fig = plt.figure()
ax = fig[:gca](projection="3d")

nx = 201
ny = 201
x = linspace(-3, 3, nx)
y = linspace(-3, 3, ny)
z = zeros(Float64, length(y), length(x))

for iy = 1:length(y)
    for ix = 1:length(x)
        z[iy, ix] = bust(x[ix], y[iy])
    end
end

x, y = np.meshgrid(x, y)
sf = ax[:plot_surface](x,y,z,rstride=4,cstride=4,cmap=cm.coolwarm,linewidth=0.1,antialiased=true)
fig[:colorbar](sf,shrink=0.5,aspect=5)

ax[:set_zlim](-4,3)
ax[:view_init](0,-90)

@windows? (
begin
    #fp = fm.FontProperties(fname="C:\\WINDOWS\\Fonts\\msgothic.ttc")
    fp = fm.FontProperties(fname="C:\\WINDOWS\\Fonts\\msmincho.ttc")
    plt.title("Julia+Matplotlib-Probe", fontproperties=fp, fontsize=25)
end
: begin
    plt.title("Julia+Matplotlib-Probe", fontsize=25)
end
)

plt.savefig("Oppai2")
plt.show(block=true)