In this procedure I roll x dies and request a sum of the dies, n.
I then specify the third argument, o, to be exactly, minimum or maximum, so that it generates random rolls until the sum matches the conditions of my third argument.
It then returns the roll in a list and the number of rolls needed to match the conditions.

optional_argument.mw
Can I somehow make the third argument optional, with exactly being the default, so that I can enter only sim(3,12) to get a roll that is exactly 12?
I've written through this maplesoft document, but don't quite seem to understand how to make an optional argument in  this regard: https://fr.maplesoft.com/support/help/maple/view.aspx?path=parameter_classes

Any help in understanding this is appreciatet :)

Dear Sir

There is error appears when we solve system of differential equations contains 13 equations ,

the form of error

Error, (in DEtools/convertsys) invalid terms in sum: -9*Y[8]+2 .. Y[3]+.3*Y[6]+.4*Y[7]+3.5*Y[11]+3.5*Y[13]

Hint we use dsolve, numerically

How to overcome on this error 

Thanks

Hello,

1 .

If A and B are non-commuting operators and I want to calculate

(A+B)^n for some integer n, how can I tell maple to not commute them?

2.

And if I have this operator construct A^3. How can I replace it with some defined function f, so A^3 becomes (f@@3)(...)

Why I maple dont send out 1????

I expcect Q be 1, but receive  .9261583073444211949, Why?

Download suuual.mw

Hi,

WARNING: This is a pretty silly question. I know it, but I've been on this for hours already...
I have this function rather simple but which depends on a natural number n. When I try a Fourier transform on it, it cannot evaluate the result for any n.

I have to define another function, and evaluate n=2, in order to get an explicit result.

h := piecewise(abs(t) < 2*Pi*n, cos(t), 0)
                    piecewise(abs(t) < 2*Pi*n, cos(t), 0)
   F.T.               
---------->         fourier(piecewise(abs(t) < 2*Pi*n, cos(t), 0), t, w)         # doesn't work for any n



h2 := eval(h, n = 2)
                      piecewise(abs(t) < 4*Pi, cos(t), 0)

   F.T.               2 w sin(4 Pi w)
---------->          ------------------                                          # works with fixed n=2
                      (w - 1) (w + 1)

How can I have an explicit result for any natural n ?

Thanks for your time

after i tryrring to instell another vergion of maple and it failed when i reinstell the older vergion of maple i get this error 
what is missing and i did wronge 

l := (n+m+sum(a[i], i = 1 .. n)+sum(b[j], j = 1 .. m))*ln(alpha)+n*ln(lambda[1])+m*ln(lambda[2])+lambda[1]*(sum(x[i], i = 1 .. n))+lambda[2]*(sum(y[j], j = 1 .. m))-(sum((2+a[i])*ln(exp(lambda[1]*x[i])-1+alpha), i = 1 .. n))-(sum((2+b[j])*ln(exp(lambda[2]*y[j])-1+alpha), j = 1 .. m));

Error, (in property/ConvertProperty) invalid input: PropRange uses a 2nd argument, b, which is missing

Given the vector A = 4.00i^+7.00j^ , find the magnitude of the vector.

ithprime(n)  gives the unique prime whose index is n, thus ithprime(10) returns 29. Question: is there an inverse function in Maple, which given prime p would return n such that ithprime(n) = p?

Or if no such function exists can anyone suggest a way such a function could be coded?

Thanks in advance

David.

These are the timings for various algorithms, using different starting points deriving surfaces of dimension 5, 4, 3, 2, 1

times3:=[[], [.140], [1.344, .891], [1.578, 1.312, 1.375, 1.437, 1.922, 2.625, 6.406], [2.188, 2.312, 1.687, 2.110, 2.047, 1.578, 8.953, 1.891, 1.875, 9.344, 2.203, 55.969, 2.266, 2.531, 81.078, 2.172, 50.641, 2.500, 3.141, 61.656, 3.406, 3.375]]

times1:=[[.718], [.766, 4.703], [.750, .797, 7.594, 3.938], [6.594, 7.718, 11.969, 8.485, 11.391, 130.583, 548.284, 974.435], [7.281, 8.515, 65.569, 7.016, 8.312, 9.500, 8.562, 9.766, 10.641, 12.609, 13.281, 17.453, 18.640, 1763.860, 2659.990, 7812.89, 8189.139]]

So far i can get a boxplot of either:
Statistics:-BoxPlot(`~`[`~`[log10]](times3));
Statistics:-BoxPlot(`~`[`~`[log10]](times1));

but what I'd like is a boxplot like this but i can't work out how to do this.
 

A := (a^6)^(1/3)*(-b^3)^(1/3)/a^3

I Can't get this to simplify to -b/a. 

suppose i have afunction F(x) and i want to draw a sample for X from F(x)

Even and odd complement each other

how to find other sequence which complement each other?

such as 3 sequences divided integers or 5 sequences divided whole integers

is there monomials creation method such that solve result about coefficient and power are integers when right side columns are sequences?

i find even multiplication numbers are always solved into integers coefficient and power.

is there any more other sequences?

if I choose six multiplication table sequence, 

what is this complement of six multiplication table sequence?

AJUSTEMENTIMAGE.mwAJUSTEMENTIMAGE.mwHi,

How do to insert image in Plot0 ?

I am working on a simple program (4 lines excluding header/footer and debug command) that calculates a Lie-derivatives of a function G, supplied by the user, with the rule:
 

derivative of x[i]=F[i]  (an input supplied by the user),
derivative of y[i]=y[i+1].

this calculation revolves around a dot product, and I can't seem to work out what has gone wrong. My impression is that is because nops of a Vector gives you  1 more than the length of the Vector, so the Vector handed to DotProduct is wrong- but i can't see why it doesn't work in the case of G2 and F2 but does for G and F.

lieDer2_problems_mwe.mw

I am developing an algorythm which returns some differential equation, which I want to simplify. Here is an example:

eqq:= k[t]*(`ℓ`^2)*(diff(q[3](tau), tau, tau)+(5*alpha-sigma+2*theta+1)*q[3](tau)+(-4*alpha+sigma-theta)*q[2](tau)+q[1](tau)*alpha) = -(sqrt(m*(1/k[t]))*`ℓ`*k[t]*`Δθ`*(q[3](tau)-q[2](tau))*sin(sqrt(Lambda*k[t]*(1/m))*sqrt(m*(1/k[t]))*tau)+2*xi*sqrt(lambda*k[t]*m)*(diff(q[3](tau), tau)))*`ℓ`*(1/sqrt(m*(1/k[t])))

I want the parameters to be associated to the the variables, , ,  and their derivattives. So, I have used "collect" command, as below:

vars:= {q[1](tau),q[2](tau),q[3](tau),diff(q[1](tau),tau),diff(q[2](tau),tau),diff(q[3](tau),tau),diff(q[1](tau),tau$2),diff(q[2](tau),tau$2),diff(q[3](tau),tau$2)}:
collect(eqq,vars);

The problem is that the equations remain with non-simplified terms, such as the terms inside the "sine" functions and the term "k[t]*ell^2". The command "simplify" does not have any effect. Ideally, I would like to have something like this:

(diff(q[3](tau), tau, tau))+alpha*q[1](tau)+(-4*alpha+sigma-theta)*q[2](tau)+(5*alpha-sigma+2*theta+1)*q[3](tau)+2*xi*sqrt(lambda)/`ℓ`*(diff(q[3](tau), tau))-`Δθ`*sin(sqrt(Lambda)*tau)*q[2](tau)+`Δθ`*q[3](tau)*sin(sqrt(Lambda)*tau) = 0;

Does anyone know how to solve that?