2つの独立な正規分布

Rohit さんにはフリーランスの仕事が2つあります。各仕事の報酬は、互いに独立な2つの正規分布に従います。

最初の仕事の income1 は平均 $500、標準偏差 $50
2つ目の仕事の income2 は平均 $1,000、標準偏差 $200

Rohit さんは支出を適切に計画できるよう、収入のシミュレーションを手伝ってほしいと考えています。サンプリングを使って、2つの仕事を合わせた総収入の95%信頼区間を求めましょう。

ここでは、モンテカルロシミュレーションで最も重要な確率分布の一つである正規分布を使ってシミュレーションを行います。

次のモジュールはすでにインポート済みです：NumPy は np、SciPy の stats モジュールは st。

st.norm.rvs() を使って、それぞれ適切な平均と標準偏差を設定し、正規分布から1,000回サンプリングして income1 と income2 に代入します。
income1 と income2 を加算して total_income を近似します。