#コブダグラス　Y=(LK)**(1/2) (w=5 r=20 固定費用100 ) Y=20x+100 規模に関して収穫一定
import matplotlib.pyplot as plt
import numpy as np
fig, ax = plt.subplots()
x = np.linspace(0, 60, 200)
plt.ylim(0,1500)
plt.title("constant returns to scale")
plt.xlabel("yield")
plt.ylabel("cost")
y=20*x+100
ax.plot(x,y,'r-',linewidth=5)
for i in range(1,100):
    y = 100+(5*(x**2/i))+i*20
    ax.plot(x, y, 'b-', linewidth=2)
plt.show()

#コブダグラス Y=(LK)**(1/4)( w=5 r=20 固定費用100 )Y=20x**2+100　規模に関して収穫逓減
import matplotlib.pyplot as plt
import numpy as np
fig, ax = plt.subplots()
x = np.linspace(0, 10, 200)
plt.ylim(0,1500)
plt.title("increasing returns to scale")
plt.xlabel("yield")
plt.ylabel("cost")
y=100+20*x**2
ax.plot(x, y, 'r-', linewidth=5)
for i in range(1,100):
    y = 100+(5*(x**4/i))+i*20
    ax.plot(x, y, 'b-', linewidth=2)
plt.show()

#コブダグラス Y=LK (w=5 r=20 固定100) Y=20x**(1/2)+100　規模に関して収穫逓増
import matplotlib.pyplot as plt
import numpy as np
fig, ax = plt.subplots()
x = np.linspace(0, 1000, 200)
plt.ylim(0,1500)
plt.title("decreasing returns to scale")
plt.xlabel("yield")
plt.ylabel("cost")
y=100+20*(x**(0.5))
ax.plot(x, y, 'r-', linewidth=5)
for i in range(1,100):
    y = 100+(5*(x/i))+i*20
    ax.plot(x, y, 'b-', linewidth=2)
plt.show()