plot([sqrt(x^2-5*x+6)/log[10]((x+10)^2), [-11, t, t = -15 .. 25], [-9, t, t = -15 .. 25]], x = -30 .. 20, -15 .. 25, color = [red, blue, blue], linestyle = [1, 2, 2], tickmarks = [[-30, -20, -11, -9, 2, 3, 20], default], scaling = constrained, discont, numpoints = 10000);

The points of discontinuity of the function are  -11, -10, -9 . The point  -10  is the point of removable discontinuity. In fact, the range is  (-infinity, infinity) . It follows from the values of limits and the continuity of the function in the corresponding intervals:

Range(limit(sqrt(x^2-5*x+6)/log[10]((x+10)^2), x = -11, right),  Open(limit(sqrt(x^2-5*x+6)/log[10]((x+10)^2), x = -10))) ; 

Range( limit(sqrt(x^2-5*x+6)/log[10]((x+10)^2), x = 2),  limit(sqrt(x^2-5*x+6)/log[10]((x+10)^2), x = -9, right)); 

                                                                  Range(- infinity, Open(0))
                                                                   Range(0,  infinity)

f := x->add(sin(x/(k+1))/k, k = 1 .. 49)+sum(sin(x/(k+1))/k, k = 50 .. infinity):

plot(f, 0..50);

A:=f(x)+g(y)+x^2+y^3+65:

select(has, A, x);

select(has, A, y);

select(not has, A, {x,y});

a:=Matrix([[0,0,0,0,0,0,0,1,1,1,1], [0,0,0,1,1,1,1,0,0,0,0], [0,1,1,0,0,1,1,0,0,1,1], [1,0,1,0,1,0,1,0,1,0,1]]):

Seq:=seq(seq(a[i,j], j=1..11), i=1..4);

b:=cat(Seq);  # or

c:=[Seq];

If we want to use the  Basis  command, without conversion to  R^n  can not do. But this process can be automated with a simple procedure. Formal argument  S - a set or a list of any polynoms or any matrices (all matrices must be of the identical size):

Basises:=proc(S)

local N, n, T, T0, m;

uses LinearAlgebra;

N:=nops(S);

if type(S[1], polynom) then

n:=max(seq(degree(S[k]), k=1..N));

T:=[seq([seq(coeff(S[k], x, m), m=n..0, -1)], k=1..N)];

map(convert, T, Vector);

T0:=Basis(%);

print(convert([seq(add(T0[k][i]*x^(n-i+1), i=1..n+1), k=1..nops(T0))], whattype(S)));

fi;

if type(S[1], Matrix) then

m:=RowDimension(S[1]);  n:=ColumnDimension(S[1]);

T:=[seq([seq(seq(S[k][i,j], j=1..n), i=1..m)], k=1..N)];

map(convert, T, Vector);

T0:=Basis(%);

convert([seq(Matrix(m, n, convert(T0[k], list)), k=1..nops(T0))], whattype(S));

fi;

end proc:

Examples:

Basises({x^2+x+4,  x+3,  2*x^2-x-5,  5*x^2+x-7});

Basises({Matrix([[2,3],[5,6]]),  Matrix([[3,2],[0,1]]),  Matrix([[1,1],[0,5]]),  Matrix([[4,3],[0,6]])});

Your system has a lot of solutions, even in a narrow ranges near the origin. See the plot:

h := 1: lamda := 2:

sys:=[100*abs(((1-nc+I*kc)/(1+nc-I*kc)+(nc-I*kc-.518+I*0.23e-1)*exp(-3.14*(I*4)*h*(nc-I*kc)/lamda)/(nc-I*kc+.518-I*0.23e-1))/(1+(1-nc+I*kc)*(nc-I*kc-.518+I*0.23e-1)*exp(-(1/581)*(3.14*(I*4)*2500)*(nc-I*kc))/((1+nc-I*kc)*(nc-I*kc+.518-I*0.23e-1))))^2 = 43.12, 100*abs(((1-nc+I*kc)/(1+nc-I*kc)+(nc-I*kc-.176+I*0.23e-1)*exp(-3.14*(I*4)*h*(nc-I*kc)/lamda)/(nc-I*kc+.176-I*0.23e-1))/(1+(1-nc+I*kc)*(nc-I*kc-.176+I*0.23e-1)*exp(-(1/720)*(3.14*(I*4)*2500)*(nc-I*kc))/((1+nc-I*kc)*(nc-I*kc+.176-I*0.23e-1))))^2 = 44.37]:

plots[implicitplot](sys, kc = -0.1 .. 0.1, nc = -0.5 .. 0.5, color = [red, blue], numpoints = 10000);

So, first of all, you need to decide - what the root of you interested. If you interested the root closest to the origin, we can further narrow the range:

plots[implicitplot](sys, kc = 0 .. 0.01, nc = -0.07 .. -0.12, color = [red, blue], numpoints = 10000);

Using DirectSearch with the  initialpoint  option, you can specify the value of the root.

The last lines of your code should be

E[j] := simplify(fnormal(LinearAlgebra[Eigenvalues](C[j])));

if E[j][1] > E[j][2] then lambda[j] := E[j][2]/(E[j][1]+E[j][2]) else lambda[j] := E[j][1]/(E[j][1]+E[j][2])  end if;

lambda[j];

 end do;

y := x -> 3*x/(x-2):

plots[display](plottools[disk]([3, 9], 0.15, color = red), plot([y(x), 9+(D(y))(3)*(x-3)], x = 2 .. 9, -1 .. 15, color = [blue, green], thickness = 2, linestyle = [3, 1], scaling = constrained));

eval(implicitdiff(x^2*y-3*y^3*x=0, y, x), {x=3, y=1});

                                             1/6

A:=plot(sin(x), x=-Pi..Pi, thickness=2, scaling=constrained):

plottools[reflect](A, [[0,0], [1,1]]);

dd := ["0101101", "0000001", "0011001", "0000001"];
map(convert, parse~(dd), decimal, binary)

                         [45, 1, 25, 1]

See  http://www.mapleprimes.com/questions/153200-How-Do-I-Write-A-Function-That-Returns#answer153208

The result can be improved:

with(DirectSearch): with(plots):

 xy := Matrix([[0.2e-1, .158], [0.2e-1, .159], [0.3e-1, .161], [0.3e-1, .164], [0.3e-1, .166], [0.4e-1, .169], [0.6e-1, .173], [0.8e-1, .178], [.1, .185], [.11, .187], [.14, .193], [.19, .2], [.28, .21], [.38, .223], [.44, .233], [.58, .244], [.82, .256], [1.4, .278], [1.71, .281], [1.78, .282], [1.78, .282], [1.81, .282]]):

 X := xy[() .. (), 1]:

 Y := xy[() .. (), 2]:

 f := a+b*x^c+d*x^e:

 sol := DataFit(f, X, Y, x, evaluationlimit=20000):

 f1 := eval(f, sol[2]);

 p1, p2 := pointplot(X, Y), plot(f1, x = 0 .. 2):

 display(p1, p2);

f:= n-> a^n + b^n + c^n:

simplify(f(5), {f(1)=1, f(2)=2, f(3)=3});

                          6

L := [[2, 0], [6, 1], [8, 0]]:   eq := f(x) = a*x^2+b*x+c:

subs(solve({seq(eval(eq, {x = L[i, 1], f(x) = L[i, 2]}), i = 1 .. 3)}), eq);

Edit: another variant

L:=[[2,0],[6,1],[8,0]]:   f:= x->a*x^2+b*x+c:

'f'(x)=subs(solve({'f(L[i,1])=L[i,2]' $ i=1..3}), f(x));