我用matlab 中的 contour画clayton的密度函数等高线，期望得到类似于年轮的等高线图，目前只画出了分布函数等高线，并测试根据clayton的参数alpha的不断变化，分布函数等高线由椭圆形趋于矩形，棱角逐渐鲜明。



错误的密度函数等高线

我并不是懒惰，但是几经测试依然无法得到圆环状的等高线，也再没有找出自己的问题所在，请明白的高手指点一二。

%年轮等高线测试
clear all
clc
n=700;
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密度函数表达式
contour(uu,vv,z)

应该出现的等高线

求助！如果画出，也请讲解一下，个别代码编写作用。致谢，并追加论坛币。