我的分形函数编程哪里有问题吗，和论文上得出的结果不一样

结果对不上，求大神看看是程序哪里有问题吗

@clear all
clc
d=2.1;
BN=unifrnd(0,2pi,26);
AN=unifrnd(0,2pi,26);
c=normrnd(0,1,26);
x=linspace(1,2,26);
y=linspace(1,2,26);
[X,Y]=meshgrid(x,y);
g=0;
for n=1:100
g= g+0.01*1.5.^((d-3).*n). *sin((X.*cos(BN)+Y.*sin(BN))*1.5^n + AN);
end
subplot(2,1,1)
surf(X,Y,g)
legend(‘D=2.1’)


下面的是我算的结果

希望有做过的大神指导