Let all vertices of a cube be in an one-to-one correspondence with the elements of the set {seq(j^2,j=1..8)}. Let us consider the sum  of the products of the numbers at the vertices belonging to an edge over all edges of the cube. What is its minimum? How to find this with Maple?

In maple 16 I do the following.

rsolve({a(1) = 1+2/n, a(p) = 1+(p-1)*a(p-1)/n+2*b(p-1)/n, b(0) = 1, b(p) = 1+p*b(p-1)/n}, {a(p), b(p)});

subs(p = n-1, %);

I then copy the answer for a(n-1) and evaluate it at n=20, say.

evalf(subs(n = 20, (1/n)^(n-1)*n*GAMMA(n-1)*(Sum((n+2*(1/n)^p1*exp(n)*GAMMA(p1+1, n))/(n^2*(1/n)^(p1+1)*GAMMA(p1+1)), p1 = 1 .. n-2)+1+2/n)));

which gives

6.067606450

But if I plot a(n-1) with

How to programmatically insert the signs  +, −, ×, /, (, ) into 123456789  to make it equal to 100 ? Of course, with Maple. I understand more or less how to produce that by using formal grammars
(see http://en.wikipedia.org/wiki/Formal_grammar ).

This is a simple question,

How do I tell maple to give only one fraction instead of a sum of several with different denominators.

Example

if i sum 1/3 + x/5 maple simply gives me: 1/3 + x/5

I want maple to give me (5+3x)/15

Thanks

Hello,

I would like to do the following procedure from a basic concept, but I can't figure out what I'm doing wrong (sorry for that, I'm a beginner). I have tried to look up some more information about procedures on the Help menu with no success. Here you have what I'm planning to do:

Considering the following function (which approximates the cosine function):

aprcos:=(x,n)->add((-1)^i*(x^(2*i)/(2*i)!),i=0..n)

I would like...

I have simplified a fairly large expression derived by Maple and am chacking to see if the two expressions are the same.  Previously, when I did this, the result was pretty straightforward: when the output was 0 = 0, the expressions were the same; however, now when Maple simplifies the expression I can sometimes use a combination of the simplify, factor, expand, collect, and combine functions to get it to spit out 0 = 0, but sometimes I can't.  Here is an example which shows the problem.

how do i write a command of a length of an array?

for example, c:=[1 3 5 7 8 2 1]..the length is 7 right? so how do i write the command..

it's a pleasure if any can help..

restart;
a:= 2.75:  b:= 7.:  N:= 2^4:
X:= t-> add(sin(b^k*t)/a^k, k= 0..N):
Y:= t-> add(cos(b^k*t)/a^k, k= 0..N):
plot(
   [X, Y, 0..2*Pi], numpoints= 2^12
  ,scaling= constrained
  ,axes= none
  ,color= cyan
  ,caption= cat("Happy New Year ", StringTools:-FormatTime("%Y"))

Is it possible to insert a hyperlink to a pdf file ? I use links to other Maple worksheets and internet URLs, but cannot seem to find a way to create a link to a local pdf file. Thanks for any help.

Happy Holidays !

Laguerre
charc := -(-1/(-1+I*X))^alpha/(-1+I*X);
Density := int(-exp(-I*X*u)*(-1/(-1+I*X))^alpha/(-1+I*X), X = -infinity .. infinity);
can not evaluate

Meixner
charc := (1+mean/(-1+exp(I*X)*mean)-mean*exp(I*X)/(-1+exp(I*X)*mean))^gamma;
Density := int(exp(-I*X*u)*(1+mean/(-1+exp(I*X)*mean)-mean*exp(I*X)/(-1+exp(I*X)*mean))^gamma, X = -infinity .. infinity);

Bernoulli
charc := 1/(I*X/(-1+exp(I*X)))^n;
Density := int(exp(-I*X*u)/(I*X/(-1+exp(I*X)))^n, X = -infinity .. infinity);...

Custom distribution
restart;
with(Statistics):
Density := 2*exp(-t^2)*sqrt(Pi);
Dist := int(Density, t);

a.

RealDist := Distribution(Dist(t-1));
Error, (in Statistics:-Distribution) unexpected argument(s): Pi(t)*(erf(t))(t)

b.

RealDist := Distribution(PDF = (t->2*exp(-t^2)*sqrt(Pi)), CDF = (t->Dist(t)), Mean = 1);

RealDist := Distribution(PDF = (t->2*exp(-t^2)*sqrt(Pi)), CDF = (t->Dist(t)));
X:=RandomVariable(RealDist(t-1));

I'd like to be able to get a quicker low accuracy computation of the argument of the complex-valued JacobiTheta1 function.

I only need as much accuracy as to get a decent plot, which means... I don't know... about 1e-3 relative accuracy?

If I can get a formula or approach that is "evalhf'able" then I can construct a procedure which does this action and writes results inplace to a float[8] Matrix.

The naive approach in the code below takes about 6min...

Hi all, I am having two simultaneous equations , i dont know how to solve these two equations, any one help to solve these equations..

Im knowing the approximate/accurate answer for the equations,
These are the values for my equations

alpha=1.53222, beta= -0.555345, kappa= -0.399943,

need_help.mw

Thank you

Find all the positive integers ,such that