結合密度関数
x=seq(-15,15,length=50) y=seq(-15,15,length=50) pdf=function(x,y){ return((1/(2*pi*70))*exp(-(1/2)*(((x-(-1))^2)/49))*exp((-1/2)*(((y-3)^2)/10))) } z=matrix(0,50,50) #z for(i in 1:50){ for(j in 1:50){ z[i,j]=pdf(x[i],y[j]) } } persp(x,y,z,theta = 50, phi = 30, expand = 0.5, col = "lightblue",zlim=c(0,0.002))