[問題] for迴圈中的矩陣問題已回收
最近寫了一個for迴圈 用來求兩個三元二次方程式的解
最後是用三維矩陣u表示解 程式如下
但是執行後 卻沒有反應(沒有當掉 也沒有出現錯誤)
想請大家幫我看一下 哪邊有錯誤
P.S if那邊的意思是,當條件成立,則有一三維矩陣u(i,j,k)
syms u;
%設邊界值
Lx=3;
Ly=2.5;
Lz=2;
h=0.1;
nx=Lx/h; ny=Ly/h; nz=Lz/h;
x=-Lx:h:Lx;
y=-Ly:h:Ly;
z=-Lz:h:Lz;
%開始算
for i=1:2*nx+1
for j=1:2*ny+1
for k=1:2*nz+1
if (((x(i))^2+(y(j))^2+(z(k))^2-1)*((x(i))^2+(y(j))^2+(z(k))^2-4))<=0
u(i,j,k)=1;
else
u(i,j,k)=0;
end
end
end
end
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 163.25.118.141
→
03/02 00:53, , 1F
03/02 00:53, 1F
→
03/02 00:55, , 2F
03/02 00:55, 2F
→
03/02 01:25, , 3F
03/02 01:25, 3F
→
03/02 02:59, , 4F
03/02 02:59, 4F
→
03/02 08:31, , 5F
03/02 08:31, 5F
→
03/02 08:32, , 6F
03/02 08:32, 6F
→
03/02 09:26, , 7F
03/02 09:26, 7F
推
03/02 11:09, , 8F
03/02 11:09, 8F
→
03/02 11:10, , 9F
03/02 11:10, 9F
→
03/02 11:21, , 10F
03/02 11:21, 10F
→
03/02 14:37, , 11F
03/02 14:37, 11F