>PS := `minus`(indets(sys1), {A1, A2, A3, A4, A5, A6, A7, A8});

{Mp1, Mp2, Mx, My, Mz, b1, b2, b3, b4, b5, b6, c1, c2, c3, c4,

c5, c6, c7, c8, c9, d1, d2, d3, d4, d5, d6, e1, e2, e3, e4, e5,

e6, e7, e8, e9, f1, f2, f3, f4, f5, f6, f7, f8, f9, g11, g12,

g13, g21, g22, g23, g31, g32, h1, h2, h3, h4, k1, k2, k3, k4,

k5, k6, m1, m2, m3, m4, m5, m6, n1, n2, n3, n4, n5, n6, p1, p2,

p3, p4, p5, p6, q1, q2, q3, q4, q5, r1, r2, r3, r4, r5, r6, r7,

r8, r9, s1, s2, s3, s4, s5, s6, s7, s8, s9, t1, t2, t3, t4, t5,

t6, t7, t8, t9, u1, u2, u3, x1, x10, x2, x3, x4, x5, x6, x9,

y1, y2, y3, y4, y5, y6, z1, z2, z3, z4, J[xx], J[yy], J[zz],

m[L1], m[L2], m[b]}

#random evaluation of parameters

>VP := [seq(c = op(convert(Statistics:-Sample(Uniform(-10.0, 10.0), 1), list)), c = PS)]:

#substitution of evaluated parameters in sys1

>SYS1 := eval(sys1, VP);

>solve(SYS1);

{A1 = 652.5529210, A2 = -771.1126525, A3 = -625.5533892,

A4 = -220.3553829, A5 = -750.3731914, A6 = -1.425032878,

A7 = -981.2020778, A8 = -582.7520760}

Download example.mw