作用dworldD(下限,上限)
S=1000;
Q = 15;
黑色= 0.25;
裸= 0.5;
白色= 0.75;
死亡= 0.1;
lowgrowthK = 278;
highgrowthK = 313;
startarea = 0.01;
数= 100;
长度=(上下)*数目;
如果长度〜=圆(长)
长度= ROUND(长度);
J = ABS(长度);
别的J = ABS(长度);
结尾
如果长度> 0
跳=(上部* S - 低级* S)/(长度-1);
E(1,1)=低级* S;
为了a=2:长度
E(A,1)= E((A-1),1)+跳跃;
结尾
结尾
如果长度<0
跳=(低级* S-上部* S)/(J + 1);
E(1,1)=低级* S;
为了A = 2:J-
E(a,1)=E((a-1),1)-跳跃;
结尾
结尾
如果长度<0
长度= ABS(长度);
结尾
SB = 5.669 * 10 ^ -8;
Ablack = startarea;
一株白色= startarea;
为了一个= 1:长度
ABARE = 1-Ablack-一株白色;
Ap=Abare*裸露+白色*白色+黑色*黑色;
Temission =((E(A,1)/ SB)×(1-AP))^ 0.25;
Trock =((E(A,1)/ SB)*(1-裸))^ 0.25;
Tblack=q*(Ap黑色)+淡紫色;
Twhite = Q *(AP-白色)+ Temission;
betab =β(lowgrowthK,highgrowthK,Tblack);
betaw=贝塔(低生长厚度、高生长厚度、Twhite);
growthb = Ablack *(betab * ABARE死亡);
growthw = *一株白色(betaw * ABARE死亡);
Ablack = Ablack + growthb;
Awhite=Awhite+growthw;
如果Ablack
Ablack = startarea;
结尾
如果Ablack>1
Ablack=1;
结尾
如果阿维特
一株白色= startarea;
结尾
如果一株白色> 1
一株白色= 1;
结尾
Temissionresults(1,A)= Temission;
Ablackresults(1,A)= Ablack;
Awhiteresults(1,a)=Awhite;
Trock结果(1,a)=Trock;
Apresults(1,A)=鸭;
growthbresults(1,A)= growthb;
growthwresults(1,A)= growthw;
结尾