Continue to Site

Eng-Tips is the largest engineering community on the Internet

Intelligent Work Forums for Engineering Professionals

  • Congratulations pierreick on being selected by the Eng-Tips community for having the most helpful posts in the forums last week. Way to Go!

spdiags

Status
Not open for further replies.

calledol

Civil/Environmental
Mar 5, 2011
1
Hi,

I'm trying to create a 40000x40000 sparse matrix with 2 diagonals. This is my code:


n = 200;
s = n*n;

for j=1:2
for i=1:(s-1)
if mod(i,n) == 0
u(i,j) = 0;
else
u(i,j) = 1;
end
end
end

c = [-1;1];

Y = spdiags(u,c,s,s);

When I try to run it i get the following message:

??? Index exceeds matrix dimensions.

Error in ==> spdiags at 114
a((len(k)+1):len(k+1),:) = [i i+d(k) B(i+(m>=n)*d(k),k)];

Error in ==> uppg3 at 26
Y = spdiags(u,c,s,s);



I don't see where the problem is. Would really appreciate some help.

Best regards
Calle
 
Replies continue below

Recommended for you

Status
Not open for further replies.

Part and Inventory Search

Sponsor