I am trying to calculate the volume, lateral surface, ..., of a cone when we have two of three: radius of base, height, generator of cone.
My code only find with a element of list. How can I print all the results of the list? My code

restart;
 L := [[3, 4, 5], [5, 12, 13], [7, 24, 25], [8, 15, 17], [9, 40, 41], [11, 60, 61], [12, 35, 37], [13, 84, 85], [15, 112, 113], [16, 63, 65], [17, 144, 145], [19, 180, 181], [20, 21, 29], [20, 99, 101], [21, 220, 221], [23, 264, 265], [24, 143, 145], [28, 45, 53], [28, 195, 197], [32, 255, 257], [33, 56, 65], [36, 77, 85], [39, 80, 89], [44, 117, 125], [48, 55, 73], [51, 140, 149], [52, 165, 173], [57, 176, 185], [60, 91, 109], [60, 221, 229], [65, 72, 97], [68, 285, 293], [69, 260, 269], [84, 187, 205], [85, 132, 157], [88, 105, 137], [95, 168, 193], [96, 247, 265], [104, 153, 185], [105, 208, 233], [115, 252, 277], [119, 120, 169], [120, 209, 241], [133, 156, 205], [140, 171, 221], [160, 231, 281], [161, 240, 289]];
R := op(1, L[1]);
h := op(2, L[1]);
 l := op(3, L[1]);
peremeter_base := 2*Pi*R;
S_base := Pi*R^2;
V := (1/3)*Pi*R^2*h;
S_lateral := Pi*R*l;
S_total:= Pi*R*l + Pi*R^2;

I am trying to calculate log[1/3](x)-log[sqrt(3)](x^2)+log[x](9) with log[3](x) = a.
I tried:

restart;
sol := solve(log[3](x) = a, x);
f :=x->log[1/3](x)-log[sqrt(3)](x^2)+log[x](9) ;
simplify(f(sol))

I don't get the answer (2-5a^2)/a.

How can I get that answer?

I want to make some equation and convert to LaTeX to have some Exercises "Solve the equation" 
My code:
restart; A := [1, 2, 3/2, 4, 5/9]; B := [2, 5, -7, 1, 8]; C := [-1, -2, -3, -4, -5/3]; for i to nops(A) do 'solve*the*equation'*A[i]*x^2+B[i]*x+C[i] = 0 end do
How can I get code in LaTeX, it's mean, I get "Solve the equation $x^2+2x-1=0$" and get roots (if have)

I posted my question at here https://math.stackexchange.com/questions/2314488/how-can-i-find-maximum-and-minimum-modulus-of-a-complex-number.
With Mathematica, I got min is 22/5. This result is different from my solution by hand.

Repeat my problem. Let be the number z so that $|z+1| + 4|z-1| = 25$. Find the greastest and the least of the modulus of $z$. How can I find greastest and the least of modulus of z with Maple.

I have a list
restart;
L := [10*x+y-6*z+66 = 0, 26*x-10*y-3*z+96 = 0, 2*x-y+6 = 0, 2*x-y+6 = 0, 4*x-5*y+3*z-6 = 0, 30*x-22*y+7*z+48 = 0, 2*x-2*y+z = 0, 26*x-19*y+6*z+42 = 0, 22*x-17*y+6*z+30 = 0, 14*x-30*y+23*z-96 = 0, 6*x-10*y+7*z-24 = 0, 2*x-7*y+6*z-30 = 0, 4*x-y-z+10 = 0, 42*x-2*y-19*z+20 = 0, 38*x-2*y-17*z+20 = 0, 30*x-y-14*z+10 = 0, 18*x+2*y-11*z-20 = 0, 6*x+13*y-16*z-130 = 0, 18*x+7*y-16*z-70 = 0, 6*x+2*y-5*z-20 = 0, 14*x+15*y-22*z-150 = 0, 18*x+13*y-22*z-130 = 0, 6*x+4*y-7*z-40 = 0, 26*x-3*y-10*z+110 = 0, 10*x-6*y+z+52 = 0, 2*x-z+8 = 0, 26*x+6*y-19*z+92 = 0, 22*x+6*y-17*z+76 = 0, 2*x+y-2*z+6 = 0, 14*x-18*y+11*z+92 = 0, 2*x-6*y+5*z+20 = 0, 2*x+6*y-7*z-4 = 0, 10*x-3*y-2*z+22 = 0, 26*x-7*y-6*z+66 = 0, 2*x-y = 0, y-z+11 = 0, 4*x-9*y+7*z-77 = 0, 6*x-5*y+2*z-22 = 0, 24*x-13*y+z-11 = 0, 2*x-y = 0, 2*x-y = 0, y-z+11 = 0, 4*x-5*y+3*z-33 = 0, 6*x-4*y+z-11 = 0, 4*x-3*y+z-11 = 0, 18*x-19*y+10*z-110 = 0, 14*x-13*y+6*z-66 = 0, 2*x-2*y+z-11 = 0, 26*x-6*y-7*z+68 = 0, 10*x-2*y-3*z+24 = 0, 2*x-z+2 = 0, 6*x+y-4*z-1 = 0, 4*x+y-3*z-3 = 0, y-z-7 = 0, 4*x+3*y-5*z-17 = 0, 18*x+10*y-19*z-52 = 0, 8*x+17*y-21*z-111 = 0, 4*x+5*y-7*z-31 = 0, 4*x-3*y+z-2 = 0, 10*x-6*y+z+4 = 0, 14*x-10*y+3*z-4 = 0, 26*x-15*y+2*z+14 = 0, 22*x-13*y+2*z+10 = 0, 18*x-26*y+17*z-84 = 0, 2*x-6*y+5*z-28 = 0, 38*x-6*y-13*z-16 = 0, 6*x-y-2*z-2 = 0, 18*x-2*y-7*z-16 = 0, 6*x+y-4*z-22 = 0, 6*x+y-4*z-22 = 0, 14*x+11*y-18*z-138 = 0, 2*x+y-2*z-14 = 0, 2*x+y-2*z-14 = 0, 22*x+2*y-13*z+40 = 0, 14*x+3*y-10*z+22 = 0, 2*x-2*y+z+8 = 0, 2*x-5*y+4*z+14 = 0, 2*x+5*y-6*z-6 = 0, 2*x-z+2 = 0, 6*x+y-4*z+8 = 0, 22*x-12*y+z-2 = 0, 8*x-9*y+5*z-10 = 0, y-z+2 = 0, y-z+11 = 0, 2*x-y-9 = 0, 8*x-3*y-z-25 = 0, 34*x-18*y+z-164 = 0, 18*x-10*y+z-92 = 0, 14*x-10*y+3*z-96 = 0, 28*x-17*y+3*z-159 = 0, 8*x-5*y+z-47 = 0, 2*x+3*y-4*z+29 = 0, 4*x-3*y+z-23 = 0, 32*x-15*y-z-103 = 0, 20*x-9*y-z-61 = 0, 10*x-6*y+z-44 = 0, 14*x-10*y+3*z-76 = 0, 2*x-2*y+z-16 = 0, 4*x-3*y+z-23 = 0, 12*x-19*y+13*z-89 = 0, 16*x-19*y+11*z-87 = 0, 20*x-3*y-7*z-73 = 0, 14*x-3*y-4*z-43 = 0, 4*x+3*y-5*z-47 = 0, 16*x-3*y-5*z-53 = 0, 22*x-2*y-9*z-92 = 0, 12*x+5*y-11*z-105 = 0, 16*x+3*y-11*z-107 = 0, 12*x+y-7*z-13 = 0, 6*x+y-4*z-7 = 0, 10*x-14*y+9*z+4 = 0, y-z-7 = 0, 2*x-z-7 = 0, 2*x+9*y-10*z-70 = 0, 18*x+y-10*z-70 = 0, 12*x+11*y-17*z-119 = 0, 2*x+y-2*z-14 = 0, 2*x-4*y+3*z+15 = 0, 4*x+y-3*z-15 = 0, 2*x-11*y+10*z+50 = 0, 14*x+3*y-10*z-50 = 0, 8*x+13*y-17*z-85 = 0, 4*x+13*y-15*z-75 = 0, 2*x+5*y-6*z-30 = 0, 12*x-7*y+z-101 = 0, 2*x-y-15 = 0, 6*x-4*y+z-56 = 0, 6*x-7*y+4*z-47 = 0, 8*x-7*y+3*z-51 = 0, 2*x-4*y+3*z-24 = 0, 6*x+y-4*z-46 = 0, 2*x+3*y-4*z-10 = 0, 2*x-z-16 = 0, y-z+2 = 0, 6*x-2*y-z-52 = 0, 10*x-2*y-3*z-84 = 0, 10*x+y-6*z-78 = 0, 2*x+y-2*z-14 = 0]

Now I want to write the equation of the planes passing two points A and B and perpendicular to all the equations of the above list. I can't to write a proc. With each plane of List, I tried
with(geom3d);
point(A, 1, 5, 3): 
point(B, -2, -4, 3):
with(geom3d): 
Equation(plane(P, L[5], [x, y, z])): 
line(d, [A, B], t):
 n1 := ParallelVector(d):
 n2 := NormalVector(P):
 with(linalg):
 n3 := crossprod(n1, n2):
 k := sort(primpart(lhs(Equation(plane(R, [A, n3], [x, y, z])))));
k*signum(lcoeff(k)) = 0

I got true answer.