﻿rm(list = ls())  # ワークスペースをクリア
# 同時密度関数
f <- function(x, y) {
  z <- 0.0
  if ((x >= 0.0) && (x < 1.0) && (y > 0.0) && (y <= 2.0 - 2.0 * x)) { 
    z <- (1.0 / (2.0 * (1.0 - x)))
  }
  return(z)
}
vf <- Vectorize(f)  # スカラー関数をベクトル化
# 3次元グラフを描画
x <- seq(0.0, 2.0, length = 61)  # 0から2までの等間隔の系列：長さ=61
y <- seq(0.0, 2.0, length = 61)  # 0から2までの等間隔の系列：長さ=61
z <- outer(x, y, vf)
persp(x, y, z, expand = 0.8, shade = 0.2, phi = 16, theta = 340)