# パッケージがインストール済なら次の行は不要
# import Pkg; Pkg.add("Plots")
using Plots  # グラフ描画のためのモジュール
function f(x, y)  # 同時密度関数
    cond = (
        (x .>= 0.0) .&
        (x .< 1.0) .&
        (y .> 0.0) .&
        (y .<= 2.0 .- 2.0 .* x)
    )
    return ifelse.(cond, 1.0 ./ (2.0 .* (1.0 .- x)), 0.0)
end
# 3次元グラフを描画
x = LinRange(0, 2, 61)
y = LinRange(0, 2, 61)
X = [x for y in y, x in x]
Y = [y for y in y, x in x]
Z = f(X, Y)
fig = surface(X, Y, Z, color=:steelblue1, alpha=0.7, 
              linecolor=:black, zlims=(0,8), camera=(40,25))
display(fig)