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] .

 Why you delete my question  most of my equation are same but demand of the questions are different, i can make a 100 account and each time i asked by one of them, and right now most of my question taged like they dublicated but they are don't ، i am not jobles and sick to post a dublicate question, You started a riot.

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.

 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