>> [a,q]=solve('a*(1+q+q^2)-168=0','a^2*q^5-42=0','a,q')
a =
119.20689035758144068333445851826
156.36366254051205597017003904258-47.469747489656238958112044329905*i
204.03289228069722368816273169830-20.317000840188657242033100304375*i
204.03289228069722368816273169830+20.317000840188657242033100304375*i
156.36366254051205597017003904258+47.469747489656238958112044329905*i
q =
.31198183952016233717742279931626
.53626889262154021237010436362974e-1+.26972390191281752291610673377378*i
-.20887376140318757077810278840206+.13943509098092458998395590708750*i
-.20887376140318757077810278840206-.13943509098092458998395590708750*i
.53626889262154021237010436362974e-1-.26972390191281752291610673377378*i
>> q1=q(1);a1=a(1);
>> a6=a1*q1^5
a6 =
.35232862692763666587757313118912
>>
□□□□□□□□□
□□□□□□□■□
□□□□□□■□□
□■□□□■□□□
□□■□■□□□□
□□□■□□□□□
□□□□□□□□□