[Theone714]

首先向大神敬礼！

其中参数L，Vr，nmax和矩阵vn，NnL都已经给出。以下求V矩阵，运行结果为一系列代数式，而不是具体的数，求教如何解决。谢谢了！马上就要答辩了！
以下为程序：

%求矩阵V
syms r;
for n=1:nmax;
for m=1:nmax;
S=r^(2*L)*exp(-(vn(n)+vn(m))*r^2)*r^2*Vr;
SVr=int(S,r,0,inf);
V(n,m)=NnL(n)*NnL(m)*SVr;
end
end
V