i try find some part of solution of this kind of pde but i can't get results my openion is maybe this pde is wronge when i defined 

pde.mw

Let eps be a real number greater than zero. Calculate the limit lim(x-->oo)[x^(1-eps)*Integral(from x to x+1) sin(t^2)dt] .

my solution is a little bit long which when i click on pdetest command i wait at least more than a hour but still is runing and i don't get any result and not give me error , which i don't know my result is true or not so How i can find that my pde by this solution it will be zero or not ?

Download test.mw

Hi,

I'm looking to animate the Riemann sum only on the interval [0,2], while keeping a global view of the curve (G). I tried the background option, but without success. Any ideas? Thanks!

S7_Riemann_Animation_Aire.mw

Hi. How to plot the 2d output like this? I don't know how to declare "a" and "b" random constant?

Download x=a_y=b_plot.mw

In a procedure what is the way to only allow specific values in a variable?

ie/

change:=proc(a::integer, b={10,20,30})#only values 10, 20 and 30 allowed for variable b
  a*b;
end proc:

Trying to find a formula such that when n is negative the value is +1 and when n is positive the value is 0, without using an if statement. 

The square grid is given in the Euclidean plane. All grid points therefore have integer coordinates. Investigate whether there are three points in this grid plane that form an equilateral triangle.

Consider the following:

MyTableElement1 := proc(L::list(nonnegint))
  ## L will have only two elements
  local M, x, y;
  M:=L;
  x:=convert(M[1],string);
  y:=convert(M[2],string);

  return cat("\\begin{tabular}{c} ",x," \\\\ ",y," \\end{tabular}")

end;

MyTableElement1([2,3]);

will output

"\begin{tabular}{c} 2 \\ 3 \end{tabular}"

as desired. However, if one inserts an \hline into the table,

MyTableElement2 := proc(L::list(nonnegint))
  ## L will have only two elements
  local M, x, y;
  M:=L;
  x:=convert(M[1],string);
  y:=convert(M[2],string);

  return cat("\\begin{tabular}{c} ",x," \\\\ \\hline ",y," \\end{tabular}")

end;

will output

"\begin{tabular}{c} 2 \\ \hline 3 \end{tabular}"

again as desired. However, if I copy and paste this int a document, I get

"\begin{tabular}{c} 2 \\ hline 3 \end{tabular}".

Note that \hline is now just hline. Putting "\\\\" in front of the hline outputs

"\begin{tabular}{c} 2 \\ \\hline 3 \end{tabular}",

but this does not compile properly. How can I get a proper \hline command to appear in the table? Thank uou for your consideration.

After a long chat with ChatGPT I finally received a fully working code for a proc for a generalized Woodbury Identity for the inversion of the sum of two or more positive definite matrices.
I was inspired by a family member who is a trained professional programmer, who told me that in his professional work he uses ChatGTP for an initial draft of his program.  
I did find out that ChatGTP makes errors: from simple ones like writing 'Simplify' instead of 'simplify' to serious conceptual errors, for example in recursive loops. However, ChatGTP seems to 'understand' the error after given specific feedback. Although, this does not mean that the next proposal does not contain the same logical error. But after a long chat I received a nice proc that seems to work. 
My second surprise was that Gemini suggested a formula for the generalized Woodbury lemma that was unknown to me, and I was unable to find on Scholar Google or https://math.stackexchange.com. Based on a special case of that formula, I was able to write the second proc myself. 
In conclusion, to start working it can be helpful to collaborate with AI friend, a little patience may help, AI may not be astute as someone on Mapleprimes wrote, but neither am I. I am now retired, and it is fun to play with Maple and AI. 
By the way, the search term Woodbury did not give a single hit on Mapleprimes.With_a_little_help_from_my_friends.mw
kind regards,Harry 

 

 Can I solve the Tolman-Oppenheimer-Volkoff equation with Maple ?  I'm having trouble with Einstein's equation with the energy tensor as the second member

Let's suppose I define a procedure P which uses a Maple built-in procedure BP.
BP is a procedure which uses optional parameters, let's say a and b.

I would like to call P with the same optional parameter names than BP, that is to write something like this

P := proc(....., {a::a_type:=a:-default, b::b-type:=b_default)
       ....
       BP(...., 'a'=a, 'b'=b)
     end proc.

The execution of P(..., a=1) leads to the error unexpected option: 1=1
 

A workaround is to define P this way

P := proc(....., {P_a::a_type:=a_default, P_b::b-type:=b_default)
       ....
       BP(...., a = P_a, b = P_b)
     end proc.

# example

P(..., P_a=1, P_b="yes")

This works but it would be easier for the users to use the same option names in the definition of P than in the definition of BP.

Question: Is that possible?

Here is an illustrative example where I would prefer to write

ecdf(S, color="red", thickness=3)

instead of

ecdf(S, ecdf_color="red", ecdf_thickness=3)

Download ecdf.mw

Thanks in advance

Hello,

I tried to apply Berlekamp's theory to the polynomial X^9+X^6-X+1 in F3.
The matrix of the linear application whose kernel dimension is sought is:

 The dimension should be equal to 2 according to the theory since this polynomial decomposes into 2 irreducible polynomials.
However, the nullspace instruction of maple does not give a result consistent with the theory. Please help me find what is happening.

I was trying to find the best way to determine what powers of 10 a number is. 

trunc(log10(a)) #a being any number

It works great for integers, however it fails when the number has is some integer with 7 nine's in the decimal.
trunc(log10(99.999999))
                                                  1

trunc(log10(99.9999998))
                                                 1

trunc(log10(99.9999999))
                                                 2

To get around that fail I found length to be another way. 
     length(trunc(a))-1  #a is any number

Any other ways or mayber there's a better way?

Download evcalc.mw