Write a procedure using the variable args that will take an unspecified finite number of numbers, delete the smallest and the largest, and return the average of the rest as a floating point number. If there are fewer than 3 arguments have an error message say: "There are not enough arguments. There should be at least three." 

You should have no input parameters in the definition of the procedure. You may write it directly or you may use the Maple command sort as a part of your program. Do ?sort to see how sort works. Test your procedure with each of the following argument sequences:  

   50,40,40, 40, 40, 10

   1,2

   seq(100 - i, i = 1..100)

   seq(modp(n,111), n=1..1000).

If n people (numbered 1 to n) stand in a circle and someone starts going around the circle and eliminating every other person till only one person is left, the number J(n) of the person left at the end is given by 

    J(n) = 1                           if n = 1
    J(n) = 2*J(n/2) - 1          if n > 1 and n is even
    J(n) = 2*J((n-1)/2) + 1   if  n > 1 and n is odd

(i) Write a recursive procedure to compute J. [As a check the first 16 values (starting with 1) of J(n) are 1,1,3,1,3,5,7,1,3,5,7,9,11,13,15,1]. 
(ii)Compute the value of J(10000). 
(iii) Can you explain why this is so much faster than our recursive procedure to compute the n-th Fibonacci number?

Has anyone installed and run maple under ubuntu installed the Windows Subsystem for Linux?   We are having trouble doing this: specifically running programs with text input files, etc.

I have a say sum(F(k) ,k=1..n)  It fails unless n is an actual integer.  But sum(F(k), k=1..a) works. Then then eval(%,a=n) completes it. That is probably correct but not necessarily a valid assumption. Details in attached. Also it is hard to understand the error message.
 

Download Summation_problem.mw

Dear Users!

Hope you would be fine. I want to construct system of equations by comparing the likes powers of x^i*y^j*t^k1*exp(k2*eta) for an expression H1 present in attached file. Please see the fix my problem. I shall be very thankful for your kind help. 

Help.mw

From  time to tiime Malel emits a sort of "Boiing" noise when I make a keyboard error. I am not sure exaclty when it does this, but it is loud. Is it possible to turn it off? I have looked quite extensively and cannot find out anything.

How to compute inverse of a square matrix in maple?

Like below matrix:

Matrix(6, 6, [[1/EE__1, -nu__21/EE__2, -nu__31/EE__3, 0, 0, 0.], [-nu__12/EE__1, `#mn("1")`/EE__2, -nu__32/EE__3, 0, 0, 0.], [-nu__13/EE__1, -nu__23/EE__2, `#mn("1")`/EE__3, 0, 0, 0], [0, 0, 0, 1/GG__23, 0., 0], [0, 0, 0, 0., 1/GG__13, 0], [0., 0., 0, 0, 0, 1/GG__12]])

Maple help say there is a mode n number, but what actually is it?

In the following problem at two example are given. For Z=2 the sum is converging whereas at Z=4 it is not converging. Thank you

PROBLEM.mw

I have a problem writing a program for the numerical solution of nonlinear volterra integral equation using the method of reproducing kernel space. I have my algorithm as well as the program I tried to write, though they are full of error messages. Please could anyone give me a clue on how to go about my challenges. The algorithm is as follows:

Step 1. Fix 𝑎 ≤ 𝑥 and 𝑡 ≤ 𝑏.
If 𝑡 ≤ 𝑥, set 𝑅𝑥(𝑡) = 1 − 𝑎 + 𝑡.
Else set 𝑅𝑥(𝑡) = 1 − 𝑎 + 𝑥.
Step 2. For 𝑖 = 1, 2, . . . , 𝑚 set 𝑥_i = (𝑖 − 1)/(𝑚 − 1).

Set 𝜓_i(𝑥) = 𝐿_t𝑅𝑥(𝑡)|𝑡=𝑥_i .
Step 3. Set 𝑢₀(𝑥₁) = 𝑢(𝑥₁).
Step 4. For 𝑖 = 1, 2, . . . , 𝑚 set 𝛾_ij = [𝜓^-1]_ij.
Step 5. 𝑛 = 1.
Step 6. Set S_n = Σ𝑛
𝑘=1 𝛾_nk𝑢_k-1(𝑥_k).
Step 7. Set 𝑢_n(𝑥) = Σ𝑛
𝑖=1 S_i𝜓_i(𝑥).
Step 8. If 𝑛 < 𝑚then set 𝑛 = 𝑛 + 1 and go to step 6.
Else stop.

So I have this system of equations with which I am not sure if the result is the same or not using "series" and "limit" or what is going on here.

I hope it is clear what I mean.

Download Mapleprimes_-_Ionproduct.mw

As an example, the two-dimensional Array created by

A := Array(triangular[upper], 1..100, 1..100);

has both indexing function and storage "triangular[upper]", which is fine. However, the attempt

B := Array(triangular[upper], 1..100, 1..100, 1..100);

to make a three-dimensional analogue did not work: It returns "Error, triangular[upper] indexing is only valid with 2 dimensions". A similar error message is returned when I replace "triangular[upper]" by "symmetric".

(For definiteness, by a higher-dimensional symmetric Array I would like to understand an Array with entries that are invariant under every permutation of their indices. Similarly, I would call the Array upper-triangular if only its entries with non-decreasing indices can be non-zero.)

For a first solution attempt, I mimic a higher dimensional upper-triangular Array by instead creating a multiply nested one-dimensional Array, where the one-dimensional subarrays become shorter and shorter. I did some preliminary testing with CodeTools[Usage] and the memory and timing results seem to compare favorably to naively using standard rectangular Arrays.

It seems more natural to write my own indexing function. However, I am not sure how to write a suitable corresponding storage function, as the documentation on that latter subject mentions only Vectors and Matrices. Is it possible and advisable to write my own storage function, or is there yet another more natural and memory-efficient way to store higher-dimensional structured Arrays (with symbolic data) in Maple? 

Thank you very much for any insights, particularly documentation pointers.

Sebastiaan Janssens.

is there example data that can verify maxwell equations?

How to calculate determinant of a cube matrix ?

is there function to calculate 3x 3x 3 cube determinant?

How to proof classical and quantum both p summation equal to one in maple?

with(linalg):
A := matrix([[1,0],[0,1]]);
AA := exponential(t*A);

A := matrix([[-i/h,0],[0,-i/h]]);
AA := exponential(t*A);
Error, (in evalf/matrixexp) cannot determine if this expression is true or false: 0.9482523555e34*abs(t*i) < 1

arctan(sqrt(3)/2/(1/2));
cos(Pi/3) + i*sin(Pi/3);

exp(i*Pi/3);

cos(Pi/3) + i*sin(Pi/3) - exp(i*Pi/3);

with(ScientificConstants):
Constant(hbar);
h := evalf(%);
dsolve(diff(s(t),t) + i/h*s(t)=0,s(t));
dsolve(diff(s(t),t) + i*s(t)=0,s(t));
 

after solve differential equation, there is no Pi/3, 

can s(t) be a matrix ?

where do Pi/3 come from?

then i try Ket

with(Physics);
K1 := Ket(&psi;t);
dsolve(diff(K1,t) + i/h*K1=0,K1);

Error, invalid neutral operator;
Ket is matrix or complex number ?

and how to solve?