matlab2018b 中的双重积分最后出的结果很差 但不知道是哪里有问题

clc;

clear;

for n=1:100

  for u=1:100

    sum1=0;

    for p=1:n

      for q=1:u

        T=10;

        k=1.38*10^-23;

        f1(q,p)=exp((-k*T)*(p.*(p-1)/2)-p*q);

      end

    end

  end

  f2(n,u)=sum1+sum(f1(:));

end

f3=pinv(f2);

for n=1:100

  for u=1:100

    sum2=0;

    for c=1:n

      for d=1:u

        f4(c,d)=c.*exp((-k*T)*(c.*(c-1)/2-c*d));

      end

    end

  end

  f5(n,u)=sum2+sum(f4(:));

end

N=f3.*f5;

n=1:100;

u=1:100;

[n,u]=meshgrid(n,u)

mesh(n,u,N)

我想做一个N 随着n和u变化的三维图像，但是N的结果只是一个数值。在N（1，100）处