Ϊʲôһֱ³ö´íѽ¡£¡£¡£¡£ÄÇλ¸ßÊÖ°ï°ïæ¡£¡£
function =aaaa(n)ny=n;
nx=3*(ny-1)+1;
A=zeros(nx-1,ny-1);
for i=1:(nx-1)*(ny-1)
p=mod(i,nx-1);
q=fix(i./(nx-1))+1;
if 0<p<ny && p<q<ny
k=q+(p-1)*(nx-1);
A(i,i)=1;
A(i,k)=-1;
end
end
¸ßÊÖÔÚÄÄÀïѽ£¿£¿£¿ :sleepy::sleepy: ÄÇλ¸ßÊÖ°ï°ïæඡ£¡£¡£¡£:) ÌáʾÄÄÀïÓÐ´í£¿ ??? Attempted to access A(60,-58); index must be a positive integer or logical.
Error in ==> aaaa at 10
A(i,q+(p-1)*(nx-1))=-1;
ÔËÐоͳöÎÊÌâ¡£¡£¡£¡£:( ¸ßÊÖÔÚÄÄÀïѽ£¿¸ßÊÖÔÚÄÄÀ ¸ßÊÖÔÚÄÄÀïѽ£¿¸ßÊÖÔÚÄÄÀ ¸ßÊÖÔÚÄÄÀïѽ£¿¸ßÊÖÔÚÄÄÀ