我用matlab 中的 mesh画clayton的密度函数图，期望得到类似于帽子的密度函数图。
错误的密度函数图：

代码如下：
n=70;
x=linspace(-2,2,n);
y=linspace(-2,2,n);
u=normpdf(x');
v=normpdf(y');
[uu,vv]=meshgrid(u,v);
onen=ones(n);
r=1;
z=(1+r).*((uu.*vv).^(-r-1)).*((uu.^(-r))+(vv.^(-r))-onen).^((-2-r)/r); % Clayton copula密度函数表达式
mesh(uu,vv,z)


因为之前在论坛看到一位前辈画过Clayton Copula的 密度函数等高线图如下：

这是代码：
n=70;T = 100;
tic;
xx = (-2:4/(T-1):2)';
uu = normcdf(xx);
v = (0.02:0.03:0.2);
kappa = 1;
zz = normpdf(xx)*ones(1,T);% this the the part of the pdf from the X variable
tic;
for ii=1:T;
zz(:,ii) = zz(:,ii).*normpdf(xx(ii)).*clayton_pdf(uu,uu(ii),kappa);
x=linspace(-2,2,n);
y=linspace(-2,2,n);
end
u=normpdf(x');
v=normpdf(y');
[uu,vv]=meshgrid(u,v);
onen=ones(n);
r=1;
z=(1+r).*((uu.*vv).^(-r-1)).*((uu.^(-r))+(vv.^(-r))-onen).^((-2-r)/r); % Clayton copula密度函数表达式
contour(xx,xx,zz,v,'Color',[0 0.7 0]);
title('Clayton copula, \kappa = 1');


可以知道Clayton Copula是可以画出帽子那样形状的密度函数图的，大概形状应该是这样（用的二元正态分布作比较。。）：


求助！！！怎样画出Clayton Copula 帽子形状的密度函数图。致谢，并追加论坛币。