天津大学张琪昌等编的<分岔与混沌理论及应用>中Mathematica程序1

wy19830711

一维:

1. aa = Table[a[i, j], {i, remainder}, {j, 2, korder}];
   
2. hx = Table[0, {remainder}];
   
3. Do[hx[] = Sum[a[i, j]
   
4. x[1]^j, {j, 2, korder}] + O[x[1]]^(korder + 1), {i, 1, remainder}];
   
5. my = Table[0, {remainder}];
   
6. Do[y = hx[], {i, 1, remainder}];
   
7. Do[nx = D[hx[], x[1]] dx[1] - dy == 0; my[] = LogicalExpand[nx], {i,
   
8.       1, remainder}];
   
9. myok = Flatten[my];
   
10. aij1 = Flatten[aa];
    
11. my2 = Solve[myok, aij1];
    
12. my2 = Flatten[my2];
    
13. dx[1] = Simplify[dx[1] /. my2];
    
14. dx[1] = Normal[dx[1]];
    
15. 
    
16. centermaniford = {}
    
17. Do[dx[jj] = dx[jj];
    
18.   dx[jj] = Expand[dx[jj]];
    
19.   len = Length[dx[jj]];
    
20.   eff = Table[0, {len}];
    
21.   Do[item = dx[jj][];
    
22.     ee = 0; Do[
    
23.       eii = Exponent[item, x[ii]];
    
24.       ee = ee + eii;
    
25.       , {ii, 1, centerdimension}];
    
26.     If[ee > korder, eff[] = eff[], eff[] = eff[] + item];
    
27.     , {i, 1, len}];
    
28.   dx[jj] = 0;
    
29.   len1 = Length[eff];
    
30.   Do[
    
31.     dx[jj] = dx[jj] + eff[];
    
32.     , {i, 1, len1}];
    
33.   centermaniford = Join[{dx[jj]}, centermaniford];
    
34.   , {jj, 1, centerdimension}]
    
35. 
    
36. centermaniford = Reverse[centermaniford]
    
37. centermaniford = Simplify[centermaniford]

复制代码

wy19830711

1. remainder = total - centerdimension;
   
2. 
   
3. ahij2 = {};
   
4. Do[
   
5.     Do[
   
6.       Do[
   
7.         i = l - j; ahij2 = Join[{a[h, i, j]}, ahij2]
   
8.         , {j, 0, l}]
   
9.       , {l, 2, korder}]
   
10.     , {h, 1, remainder}];
    
11. Do[
    
12.     hyx = 0;
    
13.     Sum[Do[i = l - j; hyx = hyx + a[h, i, j] x[1]^i x[2]^
    
14.         j, {j, 0, l}], {l, 2, korder}];
    
15.     y[h] = hyx
    
16.     , {h, 1, remainder}];
    
17. nx = Table[0, {remainder}];
    
18. coe = {};
    
19. Do[
    
20.     nx[[h]] = Sum[D[y[h], x[ii]] dx[ii], {ii, 1, centerdimension}] - dy[h]
    
21.     , {h, 1, remainder}];
    
22. Do[
    
23.     Do[
    
24.       Do[
    
25.         i = ip - j;
    
26.         mycoe = Coefficient[nx[[h]], x[1], i];
    
27.         mycoe = Coefficient[mycoe, x[2], j] /. {x[1] -> 0, x[2] -> 0};
    
28.         coe = Join[{mycoe == 0}, coe]
    
29.         , {j, 0, ip}]
    
30.       , {ip, 2, korder}]
    
31.     , {h, 1, remainder}];
    
32. me = Solve[coe, ahij2];
    
33. me = Flatten[me];
    
34. Do[
    
35.     y = y /. me
    
36.     , {i, 1, remainder}];
    
37. 
    
38. centermaniford = {}
    
39. Do[
    
40.   dx[jj] = dx[jj];
    
41.   dx[jj] = Expand[dx[jj]];
    
42.   len = Length[dx[jj]];
    
43.   eff = Table[0, {len}];
    
44.   Do[
    
45.     item = dx[jj][];
    
46.     ee = 0;
    
47.     Do[
    
48.       eii = Exponent[item, x[ii]];
    
49.       ee = ee + eii;
    
50.       , {ii, 1, centerdimension}];
    
51.     If[ee > korder, eff[] = eff[], eff[] = eff[] + item];
    
52.     , {i, 1, len}];
    
53.   dx[jj] = 0;
    
54.   len1 = Length[eff];
    
55.   Do[
    
56.     dx[jj] = dx[jj] + eff[];
    
57.     , {i, 1, len1}];
    
58.   centermaniford = Join[{dx[jj]}, centermaniford];
    
59.   , {jj, 1, centerdimension}]
    
60. 
    
61. centermaniford = Reverse[centermaniford]
    
62. centermaniford = Simplify[centermaniford]

复制代码

wy19830711

1. remainder = total - centerdimension;
   
2. 
   
3. ahij2 = {};
   
4. Do[
   
5.     Do[
   
6.       Do[
   
7.         i = l - j; ahij2 = Join[{a[h, i, j]}, ahij2]
   
8.         , {j, 0, l}]
   
9.       , {l, 2, korder}]
   
10.     , {h, 1, remainder}];
    
11. Do[
    
12.     hyx = 0;
    
13.     Sum[Do[i = l - j; hyx = hyx + a[h, i, j] x[1]^i x[2]^
    
14.         j, {j, 0, l}], {l, 2, korder}];
    
15.     y[h] = hyx
    
16.     , {h, 1, remainder}];
    
17. nx = Table[0, {remainder}];
    
18. coe = {};
    
19. Do[
    
20.     nx[[h]] = Sum[D[y[h], x[ii]] dx[ii], {ii, 1, centerdimension}] - dy[h]
    
21.     , {h, 1, remainder}];
    
22. Do[
    
23.     Do[
    
24.       Do[
    
25.         i = ip - j;
    
26.         mycoe = Coefficient[nx[[h]], x[1], i];
    
27.         mycoe = Coefficient[mycoe, x[2], j] /. {x[1] -> 0, x[2] -> 0};
    
28.         coe = Join[{mycoe == 0}, coe]
    
29.         , {j, 0, ip}]
    
30.       , {ip, 2, korder}]
    
31.     , {h, 1, remainder}];
    
32. me = Solve[coe, ahij2];
    
33. me = Flatten[me];
    
34. Do[
    
35.     y = y /. me
    
36.     , {i, 1, remainder}];
    
37. 
    
38. centermaniford = {}
    
39. Do[
    
40.   dx[jj] = dx[jj];
    
41.   dx[jj] = Expand[dx[jj]];
    
42.   len = Length[dx[jj]];
    
43.   eff = Table[0, {len}];
    
44.   Do[
    
45.     item = dx[jj][];
    
46.     ee = 0;
    
47.     Do[
    
48.       eii = Exponent[item, x[ii]];
    
49.       ee = ee + eii;
    
50.       , {ii, 1, centerdimension}];
    
51.     If[ee > korder, eff[] = eff[], eff[] = eff[] + item];
    
52.     , {i, 1, len}];
    
53.   dx[jj] = 0;
    
54.   len1 = Length[eff];
    
55.   Do[
    
56.     dx[jj] = dx[jj] + eff[];
    
57.     , {i, 1, len1}];
    
58.   centermaniford = Join[{dx[jj]}, centermaniford];
    
59.   , {jj, 1, centerdimension}]
    
60. 
    
61. centermaniford = Reverse[centermaniford]
    
62. centermaniford = Simplify[centermaniford]

复制代码

wy19830711

1. remainder = total - centerdimension;
   
2. 
   
3. ahij2 = {};
   
4. Do[
   
5.     Do[
   
6.       Do[
   
7.         i = l - j; ahij2 = Join[{a[h, i, j]}, ahij2]
   
8.         , {j, 0, l}]
   
9.       , {l, 2, korder}]
   
10.     , {h, 1, remainder}];
    
11. Do[
    
12.     hyx = 0;
    
13.     Sum[Do[i = l - j; hyx = hyx + a[h, i, j] x[1]^i x[2]^
    
14.         j, {j, 0, l}], {l, 2, korder}];
    
15.     y[h] = hyx
    
16.     , {h, 1, remainder}];
    
17. nx = Table[0, {remainder}];
    
18. coe = {};
    
19. Do[
    
20.     nx[[h]] = Sum[D[y[h], x[ii]] dx[ii], {ii, 1, centerdimension}] - dy[h]
    
21.     , {h, 1, remainder}];
    
22. Do[
    
23.     Do[
    
24.       Do[
    
25.         i = ip - j;
    
26.         mycoe = Coefficient[nx[[h]], x[1], i];
    
27.         mycoe = Coefficient[mycoe, x[2], j] /. {x[1] -> 0, x[2] -> 0};
    
28.         coe = Join[{mycoe == 0}, coe]
    
29.         , {j, 0, ip}]
    
30.       , {ip, 2, korder}]
    
31.     , {h, 1, remainder}];
    
32. me = Solve[coe, ahij2];
    
33. me = Flatten[me];
    
34. Do[
    
35.     y = y /. me
    
36.     , {i, 1, remainder}];
    
37. 
    
38. centermaniford = {}
    
39. Do[
    
40.   dx[jj] = dx[jj];
    
41.   dx[jj] = Expand[dx[jj]];
    
42.   len = Length[dx[jj]];
    
43.   eff = Table[0, {len}];
    
44.   Do[
    
45.     item = dx[jj][];
    
46.     ee = 0;
    
47.     Do[
    
48.       eii = Exponent[item, x[ii]];
    
49.       ee = ee + eii;
    
50.       , {ii, 1, centerdimension}];
    
51.     If[ee > korder, eff[] = eff[], eff[] = eff[] + item];
    
52.     , {i, 1, len}];
    
53.   dx[jj] = 0;
    
54.   len1 = Length[eff];
    
55.   Do[
    
56.     dx[jj] = dx[jj] + eff[];
    
57.     , {i, 1, len1}];
    
58.   centermaniford = Join[{dx[jj]}, centermaniford];
    
59.   , {jj, 1, centerdimension}]
    
60. 
    
61. centermaniford = Reverse[centermaniford]
    
62. centermaniford = Simplify[centermaniford]

复制代码

ieluxinhua

能不能介绍一下这三个程序？

gghhjj

原帖由 ieluxinhua 于 2006-11-14 15:54 发表
能不能介绍一下这三个程序？

自己找《分岔与混沌理论及应用》对照看一下就行了

yzsldj

将这个程序用来计算实际问题时，无法得到正确结果，我试过几个"centerdimension=2"的例子，不知是何原因。

看来这个程序只能计算张琪昌<分岔与混沌理论及应用>书中的例4.4.1，输出正确的结果，其它问题的计算均报错。

wy19830711

应该不会吧,我都试过啊

gghhjj

原帖由 yzsldj 于 2006-11-18 11:05 发表
将这个程序用来计算实际问题时，无法得到正确结果，我试过几个"centerdimension=2"的例子，不知是何原因。

看来这个程序只能计算张琪昌<分岔与混沌理论及应用>书中的例4.4.1，输出正确的结果， ...

可以把你的例子给一下，大家看看

yzsldj

我的问题是：

centerdimension=3;
dx[1]=y[1]+x[1]*z[1];
dy[1]=-x[1]+y[1]*z[1]-x[1]*w[1];
dz[1]=0;
dw[1]=-w[1]+alpha*x[1]^2;
korder=5;
total1=4;

要么出错，要么无输出结果。
我把代码也附上，请大家看看问题出在哪里。

1. This program is to compute the centermaniford
   
2. 
   
3. ReadList["e:\center1.txt"]
   
4. remainder=total1-centerdimension
   
5. "******************************"
   
6. "To compute the zero space of the stable maniford and y[i]"
   
7. "******************************"
   
8. If [centerdimension==1,
   
9. aa=Table[a[i,j],{i,remainder},{j,2,korder}];
   
10. hx=Table[0,{remainder}];
    
11. Do[
    
12. hx[[i]]=Sum[a[i,j] x[1]^j,{j,2,korder}]+O[x[1]]^(korder+1),{i,1, remainder }];
    
13. my=Table[0,{remainder}];
    
14. Do[
    
15. y[i]=hx[[i]],{i,1,remainder}];
    
16. Do[
    
17.   nx=D[hx[[i]],x[1]] dx[1]-dy[i]==0;
    
18.   my[[i]]=LogicalExpand[nx]
    
19. ,{i,1,remainder}];
    
20. myok=Flatten[my];
    
21. aij1=Flatten[aa];
    
22. my2=Solve[myok,aij1];
    
23. my2= Flatten [my2];
    
24. dx[1]=Simplify[dx[1]/.my2];
    
25. dx[1]=Normal[dx[1]];
    
26. ]
    
27. 
    
28. If[centerdimension==2,
    
29.   ahij2={};
    
30.   Do[
    
31.      Do[
    
32.         Do[
    
33.             i=l-j; ahij2=Join[{a[h,i,j]},ahij2]
    
34.         ,{j,0,l}]
    
35.      ,{l,2,korder}]
    
36.   ,{h,1,remainder}];
    
37.   Do[
    
38.     hyx=0;
    
39.     Sum[Do[i=l-j;hyx=hyx+a[h,i,j] x[1]^i x[2]^j,{j,0,l}],{l,2,korder}];
    
40.     y[h]=hyx
    
41.   ,{h,1,remainder}] ;
    
42.   nx=Table[0,{remainder}] ;
    
43. coe={};
    
44. Do[
    
45.     nx[[h]]=Sum[D[y[h],x[ii]] dx[ii],{ii,1,centerdimension}]-dy[h]
    
46.    ,{h,1,remainder}];
    
47.    Do[
    
48.      Do[
    
49.         Do[
    
50.             i=ip-j;
    
51.             mycoe=Coefficient[nx[[h]],x[1],i];
    
52.             mycoe=Coefficient[mycoe,x[2],j]/.{x[1]->0,x[2]->0} ;
    
53.             coe=Join[{mycoe==0},coe]
    
54.          ,{j,0,ip}]
    
55.        ,{ip,2,korder}]
    
56.        ,{h,1,remainder}];
    
57.       me=Solve[coe,ahij2];
    
58.       me=Flatten[me];
    
59.        Do[
    
60.           y[i]=y[i]/.me
    
61.        ,{i,1,remainder}];
    
62. ]
    
63. 
    
64. If[centerdimension==3,
    
65.   ahijp3={};
    
66.   Do[
    
67.      Do[
    
68.         Do[
    
69.            Do[
    
70.               p=ip-i-j;
    
71.               ahijp3=Join[{a[h,i,j,p]},ahijp3]
    
72.            ,{j,0,ip-i}]
    
73.          ,{i,0,ip}]
    
74.       ,{ip,2,korder}]
    
75.     ,{h,1,remainder}];
    
76. 
    
77.   Do[
    
78.     hyx=0;
    
79.        Do[
    
80.           Do[
    
81.              Do[
    
82.                  p=ip-i-j;
    
83.                  hyx=hyx+a[h,i,j,p] x[1]^i x[2]^j x[3]^p
    
84.               ,{j,0,ip-i}]
    
85.             ,{i,0,ip}]
    
86.           ,{ip,2,korder}];
    
87.         y[h]=hyx;
    
88.       ,{h,1,remainder}];
    
89.   nx=Table[0,{remainder}];
    
90. coe={};
    
91. Do[
    
92.    nx[[h]]=Sum[D[y[h],x[ii]] dx[ii],{ii,1,centerdimension}]-dy[h]
    
93. ,{h,1,remainder}];
    
94.   Do[
    
95.      Do[
    
96.         Do[
    
97.            Do[
    
98.                p=ip-i-j ;
    
99.                mycoe=Coefficient[nx[[h]],x[1],i];
    
100.                mycoe=Coefficient[mycoe,x[2],j] ;
     
101. mycoe=Coefficient[mycoe,x[3],p]/.{x[1]->0,x[2]->0,x[3]->0};
     
102.             coe=Join[{mycoe==0},coe]
     
103.          ,{j,0,ip-i}]
     
104.         ,{i,0,ip}]
     
105.        ,{ip,2,korder}]
     
106.       ,{h,1,remainder}];
     
107.   me=Solve[coe,ahijp3];
     
108.   me=Flatten[me];
     
109.   Do[y[i]=y[i]/.me,{i,1,remainder}];
     
110. ]
     
111. 
     
112. If[centerdimension==4,
     
113.   ahijpl4={};
     
114.   Do[
     
115.      Do[
     
116.         Do[
     
117.            Do[
     
118.                Do[         
     
119.      l=ip-i-j-p;
     
120.               ahijp14=Join[{a[h,i,j,p,l]},ahijpl4]
     
121.            ,{p,0,ip-i-j}]
     
122.          ,{j,0,ip-i}]
     
123.         {i,0,ip}]
     
124.       ,{ip,2,korder}]
     
125.     ,{h,1,remainder}];
     
126.   Do[
     
127.     hyx=0;
     
128.        Do[
     
129.           Do[
     
130.              Do[
     
131.                 Do[
     
132.                  l=ip-i-j-p;
     
133.                  hyx=hyx+a[h,i,j,p,l] x[1]^i x[2]^j x[3]^p x[4]^l;
     
134.               ,{p,0,ip-i-j}]
     
135.             ,{j,0,ip-i}]
     
136.            ,{i,0,ip}]
     
137.           ,{ip,2,korder}];
     
138.         y[h]=hyx;
     
139.      ,{h,1,remainder}];
     
140.     nx=Table[0,{remainder}] ;
     
141. coe={} ;
     
142. Do[
     
143.    nx[[h]]=Sum[D[y[h],x[ii]] dx[ii],{ii,1,centerdimension}]-dy[h]
     
144. ,{h,1,remainder}] ;
     
145.   Do[
     
146.      Do[
     
147.         Do[
     
148.            Do[
     
149.               Do[
     
150.                l=ip-i-j-p;
     
151.                mycoe=Coefficient[nx[[h]],x[1],i];
     
152.                mycoe=Coefficient[mycoe,x[2],j] ;
     
153.                mycoe=Coefficient[mycoe,x[3],p] ;
     
154. mycoe=Coefficient[mycoe,x[4],l]/.{x[1]->0,x[2]->0,x[3]->0,x[4]->0};
     
155.             coe=Join[{mycoe==0},coe]
     
156.          ,{p,0,ip-i-j}]
     
157.         ,{j,0,ip-i}]
     
158. ,{i,0,ip}]
     
159.     ,{ip,2,korder}]
     
160.    ,{h,1,remainder}] ;
     
161.   me=Solve[coe,ahijpl4] ;
     
162.   me=Flatten[me];
     
163.   Do[
     
164. y[i]=y[i]/.me
     
165. ,{i,1,remainder}];
     
166. ]
     
167. 
     
168. If[centerdimension==5,
     
169.   ahijplm5={};
     
170.   Do[
     
171.      Do[
     
172.         Do[
     
173.            Do[
     
174.                Do[
     
175.      Do[         
     
176.      m=ip-i-j-p-l;
     
177.               ahijplm5=Join[{a[h,i,j,p,l,m]},ahijplm5]
     
178.            ,{l,0,ip-i-j-p}]
     
179.          ,{p,0,ip-i-j}]
     
180.         ,{j,0,ip-i}]
     
181.       ,{i,0,ip}]
     
182.     ,{ip,2,korder}]
     
183.     ,{h,1,remainder}];
     
184.   Do[
     
185.     hyx=0;
     
186.        Do[
     
187.           Do[
     
188.              Do[  
     
189.             Do[
     
190.                 Do[
     
191.                  m=ip-i-j-p-l;
     
192.                  hyx=hyx+a[h,i,j,p,l,m] x[1]^i x[2]^j x[3]^p x[4]^l x[5]^m;
     
193.               ,{l,0,ip-i-j-p}]
     
194.             ,{p,0,ip-i-j}]
     
195.            ,{j,0,ip-i}]
     
196. ,{i,0,ip}]
     
197.           ,{ip,2,korder}];
     
198.         y[h]=hyx;
     
199.      ,{h,1,remainder}];
     
200.     nx=Table[0,{remainder}] ;
     
201. coe={} ;
     
202. Do[
     
203.    nx[[h]]=Sum[D[y[h],x[ii]] dx[ii],{ii,1,centerdimension}]-dy[h]
     
204. ,{h,1,remainder}] ;
     
205.   Do[
     
206.      Do[
     
207.         Do[
     
208.            Do[
     
209.                Do[
     
210.     Do[
     
211.                m=ip-i-j-p-l;
     
212.                mycoe=Coefficient[nx[[h]],x[1],i];
     
213.                mycoe=Coefficient[mycoe,x[2],j] ;
     
214.                mycoe=Coefficient[mycoe,x[3],p] ;
     
215. mycoe=Coefficient[mycoe,x[4],l] ;
     
216. mycoe=Coefficient[mycoe,x[5],m]/.{x[1]->0,x[2]->0,x[3]->0,x[4]->0,x[5]->0} ;
     
217.             coe=Join[{mycoe==0},coe]
     
218.          ,{l,0,ip-i-j-p}]
     
219.         ,{p,0,ip-i-j}]
     
220. ,{j,0,ip-i}]
     
221. ,{i,0,ip}]
     
222.     ,{ip,2,korder}]
     
223.    ,{h,1,remainder}] ;
     
224. Print[ "To solve the equation"];
     
225.   me=Solve[coe,ahijplm5] ;
     
226.   me=Flatten[me];
     
227. Print[me]
     
228.   Do[
     
229. y[i]=y[i]/.me,
     
230. {i,1,remainder}];
     
231. ]
     
232. 
     
233. If[centerdimension==6,
     
234.   ahijplmn6={};
     
235.   Do[
     
236.      Do[
     
237.         Do[
     
238.           Do[
     
239.            Do[
     
240.                Do[
     
241.      Do[         
     
242.      n=ip-i-j-p-l-m;
     
243.               ahijplmn6=Join[{a[h,i,j,p,l,m,n]},ahijplmn6];
     
244.            ,{m,0,ip-i-j-p-l}]
     
245.          ,{l,0,ip-i-j-p}]
     
246.         ,{p,0,ip-i-j}]
     
247.       ,{j,0,ip-i}]
     
248.         ,{i,0,ip}]
     
249.     ,{ip,2,korder}]
     
250.     ,{h,1,remainder}];
     
251.   Do[
     
252.     hyx=0;
     
253.        Do[
     
254.           Do[
     
255.              Do[  
     
256.             Do[
     
257.                Do[
     
258.                 Do[
     
259.                  n=ip-i-j-p-l-m;
     
260.                  hyx=hyx+a[h,i,j,p,l,m,n] x[1]^i x[2]^j x[3]^p x[4]^l x[5]^m x[6]^n;
     
261.               ,{m,0,ip-i-j-p-l}]
     
262.             ,{l,0,ip-i-j-p}]
     
263.            ,{p,0,ip-i-j}]
     
264. ,{j,0,ip-i}]
     
265. ,{i,0,ip}]
     
266.           ,{ip,2,korder}];
     
267.         y[h]=hyx;
     
268.      ,{h,1,remainder}];
     
269.     nx=Table[0,{remainder}] ;
     
270. coe={} ;
     
271. Do[
     
272.    nx[[h]]=Sum[D[y[h],x[ii]] dx[ii],{ii,1,centerdimension}]-dy[h]
     
273. ,{h,1,remainder}] ;
     
274.   Do[
     
275.      Do[
     
276.         Do[
     
277.            Do[
     
278.                Do[
     
279.     Do[
     
280.         Do[
     
281.                n=ip-i-j-p-l-m;
     
282.                mycoe=Coefficient[nx[[h]],x[1],i];
     
283.                mycoe=Coefficient[mycoe,x[2],j] ;
     
284.                mycoe=Coefficient[mycoe,x[3],p] ;
     
285. mycoe=Coefficient[mycoe,x[4],l] ;
     
286. mycoe=Coefficient[mycoe,x[5],m] ;
     
287. mycoe=Coefficient[mycoe,x[6],n]/.{x[1]->0,x[2]->0,x[3]->0,x[4]->0,x[5]->0,x[6]->0} ;
     
288.             coe=Join[{mycoe==0},coe]
     
289.          ,{m,0,ip-i-j-p-l}]
     
290.         ,{l,0,ip-i-j-p}]
     
291. ,{p,0,ip-i-j}]
     
292. ,{j,0,ip-i}]
     
293. ,{i,0,ip}]
     
294.     ,{ip,2,korder}]
     
295.    ,{h,1,remainder}] ;
     
296. Print[ "To solve the equation "];
     
297.   me=Solve[coe,ahijplmn6] ;
     
298.   me=Flatten[me];
     
299. Print[me]
     
300.   Do[
     
301. y[i]=y[i]/.me,{i,1,remainder}];
     
302. ]
     
303. "******************************"
     
304. "To compute the centermaniford and crop the higher order items"
     
305. "******************************"
     
306. centermaniford={}
     
307. Do[
     
308.    dx[jj]=dx[jj];
     
309.    dx[jj]=Expand[dx[jj]];
     
310. len=Length[dx[jj]];
     
311. eff=Table[0,{len}];
     
312.   Do[
     
313.      item=dx[jj][[i]];
     
314.      ee=0;
     
315.      Do[
     
316.        eii=Exponent[item,x[ii]];
     
317.        ee=ee+eii;
     
318.        ,{ii,1,centerdimension}];
     
319.        If[ee>korder,eff[[i]]=eff[[i]],eff[[i]]=eff[[i]]+item];
     
320.      ,{i,1,len}];
     
321.      dx[jj]=0;
     
322.      len1=Length[eff];
     
323.      Do[
     
324.         dx[jj]=dx[jj]+eff[[i]];
     
325.        ,{i,1,len1}];
     
326.     centermaniford=Join[{dx[jj]},centermaniford];
     
327.   ,{jj,1,centerdimension}]
     
328. 
     
329. centermaniford=Reverse[centermaniford]>>e:\center.out
     
330. centermaniford=Simplify[centermaniford]
     
331. Print["This is the end of the program"]

复制代码

无水1324

一个不错的东西，顶上来！希望需要的人参考。程序中确实存在问题。有兴趣的可以一起研究一下！

octopussheng

这本书我看了一两遍了，后面的程序还没有算过，我还打算用这些程序出点结果呢，不会真的有问题吧，我也拿例子试试看！

无水1324

试一下，才知道真正的问题在那里

octopussheng

我算了一下书上的例4.4.1，为什么出来的结果是这样子的啊？

gghhjj

原帖由 octopussheng 于 2007-6-28 21:05 发表

                               
登录/注册后可看大图

我算了一下书上的例4.4.1，为什么出来的结果是这样子的啊？

没有这本书，能够把详细内容贴一下？