gankin = 100000 riritsu = 0.05 yokin = gankin for nen in range(11): rishi = yokin * riritsu print(nen,"年目：預金額",yokin,"円：利子",rishi,"円") yokin = yokin + rishi

import matplotlib.pyplot as plt gankin = 100000 riritsu = 0.05 yokin = gankin x = [0] y = [yokin] for nen in range(11): rishi = yokin * riritsu print(nen,"年目：預金額",yokin,"円：利子",rishi,"円") yokin = yokin + rishi x.append(nen+1) y.append(yokin) plt.plot(x,y) plt.show()

chosuimax = 15000 chosuistart = 12000 hosuiwariai = 0.15 hosuiteishiwariai = 0.20 daymax = 30 hosuiteishi = chosuimax * hosuiteishiwariai chosui = chosuistart for day in range(0,daymax+1): hosui = chosui * hosuiwariai chosuikari = chosui - hosui if chosuikari < hosuiteishi: hosui = chosui - hosuiteishi print(day,"日：貯水量",chosui,"：放水量",hosui) chosui = chosui - hosui

import matplotlib.pyplot as plt chosuimax = 15000 chosuistart = 12000 hosuiwariai = 0.15 hosuiteishiwariai = 0.20 daymax = 30 x_day = [] y_chosui = [] y_hosui = [] hosuiteishi = chosuimax * hosuiteishiwariai chosui = chosuistart for day in range(0,daymax+1): hosui = chosui * hosuiwariai chosuikari = chosui - hosui if chosuikari < hosuiteishi: hosui = chosui - hosuiteishi print(day,"日：貯水量",chosui,"：放水量",hosui) x_day.append(day) y_chosui.append(chosui) y_hosui.append(hosui) chosui = chosui - hosui plt.plot(x_day,y_chosui) plt.plot(x_day,y_hosui) plt.show()

import random sosu = 100 kakuritsu = 0.7 coin100 = 0 coin500 = 0 coin100min = 300 for kyaku in range(1,sosu+1): if random.random() < kakuritsu: coin100 = coin100 + 2 else: coin500 = coin500 + 1 coin100 = coin100 - 3 if coin100 < coin100min: coin100min = coin100 print("客",kyaku,"：100円硬貨",coin100,"：500円硬貨",coin500) print("100円硬貨 最小枚数：",coin100min)

import random import matplotlib.pyplot as plt sosu = 100 kakuritsu = 0.7 coin100 = 0 coin500 = 0 Kyaku = [] Coin100 = [] Coin500 = [] coin100min = 300 for kyaku in range(1,sosu+1): if random.random() < kakuritsu: coin100 = coin100 + 2 else: coin500 = coin500 + 1 coin100 = coin100 - 3 if coin100 < coin100min: coin100min = coin100 print("客",kyaku,"：100円硬貨",coin100,"：500円硬貨",coin500) Kyaku.append(kyaku) Coin100.append(coin100) Coin500.append(coin500) print("100円硬貨 最小枚数：",coin100min) plt.plot(Kyaku,Coin100) plt.plot(Kyaku,Coin500) plt.show()

import random inside = 0 for point in range(1,1001): x = random.random() y = random.random() distance = x * x + y * y if distance <= 1: #点が四分円内にあるとき実行 inside = inside + 1 if point == point//100*100: #point の 100 とびに実行 pai = 4 * inside / point print("点の総数",point,"：円周率の近似値",pai)

import random import matplotlib.pyplot as plt x_out = [] ; y_out = [] x_in = [] ; y_in = [] x_pai = [] ; y_pai = [] inside = 0 for point in range(1,1001): x = random.random() y = random.random() distance = x * x + y * y if distance <= 1: #点が四分円内にあるとき実行 inside = inside + 1 x_in.append(x) y_in.append(y) else: x_out.append(x) y_out.append(y) if point == point//100*100: #point の 100 とびに実行 pai = 4 * inside / point print("点の総数",point,"：円周率の近似値",pai) x_pai.append(point) y_pai.append(pai) plt.scatter(x_in, y_in, color="red") plt.scatter(x_out, y_out, color="blue") plt.show() plt.plot(x_pai,y_pai) plt.show()

# 関数 jikan =========================== def jikan(dosu):     import random     kaikyusu = len(dosu)     kakuritsu = [0] * kaikyusu     ruiseki = [0] * kaikyusu     #度数の合計 -----------------------     sum = 0     for n in dosu:         sum = sum + n     #確率 -----------------------------     for i in range(kaikyusu):         kakuritsu[i] = dosu[i] / sum     #累積確率 -------------------------     ruiseki[0] = kakuritsu[0]     for i in range(1,kaikyusu):         ruiseki[i] = ruiseki[i-1] + kakuritsu[i]     #出現確率に応じた階級の取得 -------     r = random.random()     kaikyu = 0     while True:         if r < ruiseki[kaikyu]:             return kaikyu + 1         kaikyu = kaikyu + 1 # Main Routine =========================== #最初の客の設定 ----------------------- chakuji[0] = 0 kaishiji[0] = 0 shiyoji[0] = jikan(shiyobnp) shuryoji[0] = shiyoji[0] machiji[0] = 0 #２番目以降の客の設定 ----------------- for kyaku in range(1,kyakusu):     chakukan[kyaku] = jikan(chakubnp)     chakuji[kyaku] = chakuji[kyaku-1] + chakukan[kyaku]     if shuryoji[kyaku-1] > chakuji[kyaku]:         kaishiji[kyaku] = shuryoji[kyaku-1]     else:         kaishiji[kyaku] = chakuji[kyaku]     shiyoji[kyaku] = jikan(shiyobnp)     shuryoji[kyaku] = kaishiji[kyaku] + shiyoji[kyaku]     machiji[kyaku] = kaishiji[kyaku] - chakuji[kyaku] #状況の表示 --------------------------- for kyaku in range(kyakusu):     print("客",kyaku,              "：到着間隔",chakukan[kyaku],              "：到着時刻",chakuji[kyaku],              "：開始時刻",kaishiji[kyaku],              "：使用時間",shiyoji[kyaku],              "：終了時刻",shuryoji[kyaku],              "：待ち時間",machiji[kyaku])