code snapshot我无法运行解决u [k]的一组方程的eqs表/矩阵,如附图和代码所示
Nn=2;h=pi/sqrt(Nn)
u[x_] := ([([Sum]), (i = (-Nn)), (Nn)]((u[i]) (cnqff[i]))
delta0[k_, j_] := If[k == j, 1, 0]
delta2[k_, j_] := If[k == j, -[Pi]^2/(3 h^2),(-2(-1)^(j-k))/(h^2(j- k)^2)]
eqs = Table[1/2[([Sum]), (k =(-Nn)), (Nn)](u[k] delta2[k, j])) + 1/2
qq[[j]] [([Sum]),(k = (-Nn)),(Nn)](u[k] delta0[k, j])) == [Lambda]([([Sum]),
(k = (-Nn)),(Nn)](u[k] delta0[k, j])) , {j, -Nn, Nn}]
请注意,sum是从-Nn到Nn的sigma符号,在此代码中Nn设置为2
有人可以帮忙解决这个问题吗?
这听起来很傻但是在我昨天关闭了一圈之前它跑了并且正在执行3个稀疏矩阵,第二个是我感兴趣的那个是11乘9矩阵,第2列和11个零能够提取一个9乘9系数矩阵,以便在我的计算中进行..
为什么它昨天运行而今天不运行..我重新启动了内核和整个系统,但不幸的是它没有运行。
非常感谢enter image description here和Br,