Advanced Search

Question: problem with fsolve

Question:problem with fsolve

gatsh 4 Maple
plot
April 05 2006
0
I'm trying to generate a list of solutions by changing two parameters as follows (T,x[1,a]): A:=0.2:Q:=273.15: > J:=10: > K:=283.15: > U:=0.4: > V:=0.1: > W:=0.8: > N:=(K-Q)/J+1: > M:=(W-U)/V+1: > data:=[]: > List:=[]: > forget(fsolve); > for T from Q by J to K do > eq[1]:=x[1,a]*alpha[1,a]=x[1,b]*alpha[1,b]: > eq[2]:=x[2,a]*alpha[2,a]=x[2,b]*alpha[2,b]: > eq[3]:=x[3,a]*alpha[3,a]=x[3,b]*alpha[3,b]: > eq[4]:=x[1,a]+x[2,a]+x[3,a]=1: > eq[5]:=x[1,b]+x[2,b]+x[3,b]=1: > forget(fsolve); > for x[1,a] from U by V to W do > solsf:=fsolve({eq[1],eq[2],eq[3],eq[4],eq[5]},{x[2,a],x[3,a],x[1,b],x[2,b],x[3,b]},{x[2,a]=0..1,x[3,a]=0..1,x[1,b]=0..0.1,x[2,b]=0..1,x[3,b]=0..1}): > L:=(eval([T,x[1,a],x[2,a],x[3,a],x[1,b],x[2,b],x[3,b]],solsf)): > print(L); > data:=[op(data),[L[4],L[3],L[1]],[L[7],L[6],L[1]]]: > List:=[op(List),L]: > od:od: > plot3d(data,X[3]=0..1,X[2]=0..0.3): [273.15, 0.4, 0.1987979621, 0.4012020379, 0.04790224960, 0.1129806320, 0.8391171184] [273.15, 0.5, 0.1707540969, 0.3292459031, 0.03723084836, 0.09304292907, 0.8697262226] [273.15, 0.6, 0.1360529023, 0.2639470977, 0.03044172656, 0.07806939740, 0.8914888760] [273.15, 0.7, 0.09675664427, 0.2032433557, 0.02529230202, 0.06519637414, 0.9095113238] [273.15, 0.8, 0.05339038964, 0.1466096104, 0.01937157403, 0.04795170845, 0.9326767175] [283.15, 0.4, 0.004544266383, 0.09545573362, 0.006838577972, 0.003132486469, 0.9900289356] [283.15, 0.5, 0.004544266383, 0.09545573362, 0.006838577972, 0.003132486469, 0.9900289356] [283.15, 0.6, 0.004544266383, 0.09545573362, 0.006838577972, 0.003132486469, 0.9900289356] [283.15, 0.7, 0.004544266383, 0.09545573362, 0.006838577972, 0.003132486469, 0.9900289356] [283.15, 0.8, 0.004544266383, 0.09545573362, 0.006838577972, 0.003132486469, 0.9900289356] for some reason which I don't understand the first solution with a set T is o.k but the second one even though it changes my T and even x[1,a] I'm getting the same result (5 times). what do is wrong with my code? please help me!!
Please Wait...