二元正太积分   积分函数已经写好了  怎么样能把图象话出来啊  ？？？求大神指点library(pracma)
f1<-function(x,y)
   {
    (1/(2*pi*sqrt(1-c^2)))*exp(-(x^2-2*c*x*y+y^2)/(2*(1-c^2)))}
c=0.9
f2<-dblquad(f1,-Inf,1,-Inf,1#来自pracma包

请大神解答  小弟急求