-2,1);i=ones(ny-1,1); Ax=spdiags(e*[1 -2 1],-1:1,nx-2,nx-2); Ay=spdiags(i*[1 -2 1],-1:1,ny-1,ny -1); Ay(1,1)=-3;Ay(end,end)=-3; A=kron(Ay/dy^2,speye(nx-2))+kron(speye(ny-1),Ax/dx^2); Dv=speye for the PPE e=ones(nx-2,1);i=ones(ny-2,1); Ax=spdiags(e*[1 -2 1],-1:1,nx-2,nx-2); Ay=spdiags(i -2,1);i=ones(ny-2,1); Tx=spdiags(e*[1 -2 1],-1:1,nx-2,nx-2); Ty=spdiags(i*[1 -2 1],-1:1,ny-2,ny -2); Tx(1,1)=-1; Tx(end,end)=-1; Tt=kron(Ty/dy^2,speye(nx-2))+kron(speye(ny-2),Tx/dx^2); Tt=speye((nx
calculate sparse matrix (AUTO)J_a = 2*(1/dx^2+1/dy^2);J_b = -1/dy^2;J_c = -1/dx^2; J=spalloc(Nx*Ny,Nx*Ny,(Nx
证明 f_{n}(x)=x^{n}+nx-2 ( n 为正整数)在 (0,+\infty) 上只有唯一正根 a_{n} ;(2)计算 \underset{n\rightarrow \infty}{\lim