Question: trouble creating proper sequence for solving equations

hello ,

I have figured out how to solve for corresponding values but when i try to solve the set of equations by using for example,

a :=[10 values]

b:=[10 values]

w:= [10 values]

##here i want to solve the set of equations for  [a1 b1 w1] ,[a2 b2 w1] , [a3 b3 w1]....like i want to solve the equations for all correponding values of a,b @ each w ##

set of equations (function of a,b,c)

i cant make it work ,i end up taking one value at a time for a , b and solve ...is there a way to make the solver solve for all the values



k1vals:=[169.7662,1481.894,9260.411]:

k2vals:=[1.0889,0.7873,0.60432]:

h1vals:=[22.8947,25.9618,29.0708]:

h2vals:=[-86.3961,-83.02,-79.5941]:

h3vals:=[-356.109,-350.556,-344.905]:

h4vals:=[-212.375,-350.556,-344.905]:

h5vals:=[-31.1687,-23.6689,-15.7705]:

h6vals:=[29.2084,32.5816,36.0019]:

wvals:=[0.5324,0.46 ,0.5945]:

    S := {seq({k1=k1vals[i], k2=k2vals[i], h1=h1vals[i],h2=h2vals[i],h3=h3vals[i],h4=h4vals[i],h5=h5vals[i],h6=h6vals[i],w=wvals[i]}, i=1..3)};

(1)

 

    sys:={a2+a3+a5-1, 2*a1+2*a4+4*a5-1.512-(2*w),a2+2*a3+a4-0.66-w-(2*m), (k2*a2*a4)-(a1*a3), (k1*((a1+a2+a3+a4+a5+(3.76*m))^2)*(a5*a4))-(a1^3*a2),-127.7395+(-221.862*w)-(h1*a1)-(h2*a2)-(h3*a3)-(h4*a4)-(h5*a5)-(m*3.76*h6)};

(2)

    sols := {seq(seq([op(i),op(j)], j=RootFinding:-Isolate(subs(i,sys),[a1,a2,a3,a4,a5,m])),i=S)};

(3)

 



Download paper_model_expvalue.mw

 

 

Please Wait...