数学建模社区-数学中国
标题:
Maple的for循环求助
[打印本页]
作者:
No-one
时间:
2020-3-27 14:25
标题:
Maple的for循环求助
本帖最后由 No-one 于 2020-3-27 14:26 编辑
: B0 H* u( x9 }+ v2 `
2 R8 J+ T% F4 |
想通过Maple的for循环来把数据依次写入到一个零矩阵里，但结果是产生了一个table,矩阵里的0没有被替换，请问大佬们该如何修改。谢谢啦
: S5 y" S7 `: \7 P, U! d4 g
#######################################
% j( x3 e0 ]+ u! C
Qlist := [seq(i, i = 1 .. 2)];
; C. x9 R5 w8 s; n* O
klist := [seq(.1*i, i = 1 .. 2)];
, m) E, F8 ]; S( Q
WL_Matrix := Matrix(nops(Qlist), nops(klist), 0);
+ N; P" b4 z! y9 k7 Z
n := nops(Qlist);
* }) P9 \) {' q4 v/ S F4 Z
$ c. X5 I; ~& z0 _( E
for i to n do
' o1 z! z! ]. F/ x
$ n4 w, V8 H9 K* R
for j to n do
% `7 r* h' ?, p( q3 y" z# m
" T- z4 l1 y. I* x* ~& K) Z% _# W
eq5 := [solve(subs([Q = Qlist
, k = klist[j]
], W^2*W0*k+Q^2) = 0, W)];
{; `! n0 Q* k& J% Z D2 [
# \& m0 T" z. D
WL_Marix[i, j] := eq5[1]; print(WL_Marix);
" R* K7 z, Q0 F9 }4 i' l% I" I
" @5 Y' v$ J7 b7 m' R" g
eq5 := assign['eq5']
( z$ S" k; G( j
3 z% |7 U+ C+ X1 p8 h" Z/ q* m
end do
1 l8 T) t6 j1 _$ Q5 X* ^( K5 u
2 N( L& l6 R- N6 n; b4 \
end do;
& T" Y8 d4 p! Q$ m0 L# ^
#################################
6 m9 t, [$ E" r
' m" a" b' e. v; M% i# D, x# r5 n# G
3 u3 H5 Q! e, y
' v! N; ]1 N! O C
7 o6 f& f' c3 f6 R3 i( x
* y6 `" Y" G( `/ B- R; F
" C. @. G# N2 }" }4 E' T2 j, v5 ]
: M6 W; M% a7 l) V
S, X. M u* R r$ C p z
Capture1.PNG
(9.48 KB, 下载次数: 444)
2020-3-27 14:24 上传
点击文件名下载附件
Maple18运行结果
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
