我是octave的新手,我使用的是版本4.2.1。,32位。当我尝试使用大于cca的矩阵进行操作时,我收到此错误。 5000x5000。使用较小的矩阵,它工作得很好。有没有办法增加这种允许的索引类型?
我的代码:
function [dt,F,u0,n1u]=solve_diff_sec_ord()
%n odd!!!!
n = input('Input number: '); #slice the interval[-1,1] to n+1 pieces
A = [];
for i=1:1:n
a = zeros(1,n+2);
a(1,i)=1;
a(1,i+1)=-2;
a(1,i+2)=1;
A = [A;a];
end
B = A(:,[2:n+1]);
dt=det(B);
F=norm(B,'fro');
b=zeros(1,n);
h=2/(n+1);
for i=1:1:n
if i <= (n+1)/2
b(1,i)= (-1+i*h^2)+1;
else
b(1,i)= -(-1+i*h^2)+1;
end
end
b=transpose(b);
x=linsolve(B,b,struct()); #basically solving D.E.
u0=x((n+1)/2,1);
n1u=norm(x,1);
- &GT;错误:对于Octave的索引类型
,内存或维度太大