From 3fe78dbd78a3fb66c86b55b6c752010b677d5406 Mon Sep 17 00:00:00 2001
From: MasaKato <NlGG@users.noreply.github.com>
Date: Sun, 7 Jun 2015 17:06:41 +0900
Subject: [PATCH] Create ex01
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

とりあえず。
---
 ex01 | 33 +++++++++++++++++++++++++++++++++
 1 file changed, 33 insertions(+)
 create mode 100644 ex01

diff --git a/ex01 b/ex01
new file mode 100644
index 0000000..ce7a015
--- /dev/null
+++ b/ex01
@@ -0,0 +1,33 @@
+from __future__ import division
+import math
+import functools  #for python3
+from random import uniform, normalvariate
+import matplotlib.pyplot as plt
+import numpy as np
+
+def envelop():
+    fig, ax = plt.subplots()
+    y = np.linspace(0, 100, 100)
+    colormap = ['b', 'g', 'r', 'c', 'm', 'y', 'k']
+    for k in range(1, 200, 1):
+        color = colormap[k % 6]
+        k = k / 10.0
+        c3 = 1/3*   1/50*((1/2*y-2*k**3) ** 3+(2*k**3)**3)
+        c2 = 1/2*   1/5*k*((1/2*y-(5*k**3-5*k**2+k))**2-(5*k**3-5*k**2+k) **2)
+        c1 = 1*1/4* 1/k*y + 1/100*4*k**6*y + (5*k**3-5*k**2+k)
+        c0 = 1/3*(k-2)**3+5*k+5
+        plt.xlim(0, 70)
+        plt.ylim(0, 70)
+
+        ax.plot(y, c3+c2+c1+c0, color='c', linewidth=1)
+
+    #c = 2 * math.sqrt(6) / 3 * y
+    #ax.plot(y, c, 'b-', linewidth=2, label='$k = ' + str(k) + '$')
+
+    #ax.legend()
+
+    plt.show()
+
+
+envelop()
+