@Markiyan Hirnyk   My method gives the complete list of all solutions in explicit form. In some situations, it may be important!

@Markiyan Hirnyk   In Maple 12, we immediately obtain the result

RealDomain[solve](r^6+r^4+r^2-10 = 0);

In Maple 16 additional command is needed:

RealDomain[solve](r^6+r^4+r^2-10 = 0);

allvalues([%]);

@Voltageist   Look closely at the output. You can see the exact answer to your problem:

@Preben Alsholm  Thanks for improvement. Another equivalent variant:

H:=G@@2;

@Markiyan Hirnyk   you're right! Cases  a=0  or  h=0  or  k=0  require special consideration.

@John 2014  Your problem is reduced to solving a system of linear equations (see the above link).

Solution for the base L1 (the remaining items can be solved similarly):

Poly:=x*(t^2+1)+y*(t-2)+z*(t+3);  Poly2:=a1*t^2+ a2*t+ a3;

Sys:={seq(coeff(Poly,t,n)=coeff(Poly2,t,n), n=0..2)}; #  The system of linear equations

solve(Sys, {x,y,z});

1) Use Maple syntax for imputs of formulas, for example,  a1*t^2+a2*t+a3 . Your formulas are incomprehensible.

2) See  the question

• Showing that a list is an ordered basis for a subspace

@mahmood180   Replace  Int  by  int

@Markiyan Hirnyk 

execution.mws

Carl's code is working:

Sol:=fsolve({F,Fw,Fk,Ft}, map(`=`,{w,k,ki,T},0..infinity));
evalf(eval({F,Fw,Fk,Ft}, Sol));

@apm   Since  limit  is weaker than  MultiSeries:-limit  then before  limit  you have to use a simplifying transformation as shown by Axel Vogt. In addition, instead of  combine, you can use  expand  command:

restart;

Expr:=sqrt(-2*cos(alpha)*cos(alpha+d)+2-2*sin(alpha+d)*sin(alpha))/d:

limit(expand(Expr), d = 0, left);  limit(expand(Expr), d = 0, right);

                                                         -1

                                                          1

@Markiyan Hirnyk  Ill-conditioned and ill-posed systems of linear algebraic equations are the different things.

@Carl Love 

1) Thanks for comment on zero entries.

2)  Transpose of M does not help.  

In Maple 12:

 convert(M^%T, list)[];

@Carl Love  Done in Maple 12 Classic.

# Return the angle of a point in interval [0, 2*Pi).

Angle:= proc(P::[realcons,realcons])

local arg:= evalf(arctan(P[2],P[1]));

     `if`(arg < 0, evalf(arg+2*Pi), arg)

end proc:

SortByAngle:= (A::Matrix)->

     Matrix(sort(convert(A^%T, listlist), (P,Q)-> Angle(P) < Angle(Q)))^%T:

M:= < -1,2,3,1,1,-1 ; 0,2,3,1,0,0 >:

SortByAngle(M);

@wingjammer   My example is of a general nature, you can write in the brackets the order of columns that you need.

Specify your problem more clearly. I do not understand what you want in the general case.