Items tagged with constant constant Tagged Items Feed

For all real a, the partial sums sn= sum((-1)^k (k^(1/k) -a), k=1..n) are bounded so that their limit points form an interval [-1.+  the MRB constant +a, MRB constant] of length 1-a, where the MRB constant is limit(sum((-1)^k*(k^(1/k)), k = 1 ..2*N),N=infinity).

For all complex z, the upper limit point of  sn= sum((-1)^k (k^(1/k) -z), k=1..n) is the  the MRB constant.

We see that maple knows the basics of this because when we enter sum((-1)^k*(k^(1/k)-z), k = 1 .. n) 

maple gives

sum((-1)^k*(k^(1/k)-z), k = 1 .. n)

 

marvinrayburns.com

I have some data for a model in MapleSim that I would like to use a time look up table with.  I've found that the two options for interpolation are linear and 1st derivative, but the data was intended to be interpretted as piecewise constant.  Is there any way to acheive this option in MapleSim?

Hello guys ...

I used a numerically method to solve couple differential equation that it has some boundary conditions. My problem is that some range of answers has 50% error . Do you know things for improving our answers in maple ?

my problem is :

a*Φ''''(x)+b*Φ''(x)+c*Φ(x)+d*Ψ''(x)+e*Ψ(x):=0

d*Φ''(x)+e*Φ(x)+j*Ψ''(x)+h*Ψ(x):=0

suggestion method by preben Alsholm:

a,b,c,d,e,j,h are constants.suppose some numbers for these constants . I used this code:


VR22:=0.1178*diff(phi(x),x,x,x,x)-0.2167*diff(phi(x),x,x)+0.0156*diff(psi(x),x,x)+0.2852*phi(x)+0.0804*psi(x);
VS22:=0.3668*diff(psi(x),x,x)-0.0156*diff(phi(x),x,x)-0.8043*psi(x)-0.80400*phi(x);
bok:=evalf(dsolve({VR22=0,VS22=0}));

PHI,PSI:=op(subs(bok,[phi(x),psi(x)]));
Eqs:={eval(PHI,x=1.366)=1,eval(diff(PHI,x),x=1.366)=0,eval(PHI,x=-1.366)=1,eval(diff(PHI,x),x=-1.366)=0,
eval(PSI,x=1.366)=1,eval(PSI,x=1.366)=1};
C:=fsolve(Eqs,indets(%,name));
eval(bok,C);
SOL:=fnormal(evalc(%));


I used digits for my code at the first of writting.

please help me ... what should i do?

I need to get rid of the type "constant" for gamma. In Maple, gamma is defined as Euler's constant by default. While it is easy enough to unprotect(gamma) and then get rid of its value, Maple will refuse to solve an equation for gamma, as it remains of type constant even after deassigning it. So I need to regain gamma as a variable.

Some may feel this is an unwise thing to do. But it actually is not: I am writing a document involving physics, and gamma is the accepted symbol for the relativistic energy. I cannot avoid using that, lest mass confusion ensues (this involves students). I really don't want to write gammar instead. Euler's constant, otoh, does not figure at all in my document.

Note that I need a solution that works in Maple 15 and later as I am working in a heterogeneous environment as far as Maple versions are concerned.

Thanks in advance,

Mac Dude

 

Hello.

 

Could not evaluate numerical integral with constant in it. I use method = _cuhre.  Maple print solution like this:

 

Int(Int(Int(max(0., (0.9483573506e-3*(-1.*sin(a)*cos(w)-1.*cos(a)*sin(w)*sin(b)))*cos(a)*cos(b)^2*(-58.5*signum(cos(b)*sin(w)*sin(b))*kk+200.*cos(b)*sin(w)*sin(b))*Heaviside(-58.5+200.*cos(b)*sin(w)*sin(b)*signum(cos(b)*sin(w)*sin(b))*kk)*Heaviside(1.-.98*cos(b)^2)/sqrt(1.-.98*cos(b)^2)), a = 0. .. 6.283185308), b = 0. .. 1.570796327), w = 0. .. 6.283185308)+Int(Int(Int(min(0., (0.9483573506e-3*(-1.*sin(a)*cos(w)-1.*cos(a)*sin(w)*sin(b)))*cos(a)*cos(b)^2*(-58.5*signum(cos(b)*sin(w)*sin(b))*kk+200.*cos(b)*sin(w)*sin(b))*Heaviside(-58.5+200.*cos(b)*sin(w)*sin(b)*signum(cos(b)*sin(w)*sin(b))*kk)*Heaviside(1.-.98*cos(b)^2)/sqrt(1.-.98*cos(b)^2)), a = 0. .. 6.283185308), b = 0. .. 1.570796327), w = 0. .. 6.283185308)

How could it be taken.

Thank!!!

 

Hi!

Say, I got an expression that depends on two variables, x and y. How can I tell Maple, that y is actually just a (real) constant, so y does not depend on x?

Because when I apply a differentiation with the "D" - command, it would always also write out expressions, where y is differentiated w.r.t. x.

Thanks!

Hey, I have a system of ODE and I can't draw it's phase curve. I tried to use DEplot and phaseportrait, but it doesn't work. Here is my system:

 dx/dt=x

dy/dt=ky              (k is a constant)

 

Here is my piece of code:

DE := [diff(x(t), t) = x(t)];

DF := [diff(y(t), t) = k*y(t)];

with(DEtools);

phaseportrait([DE, DF], [y, x], t = -5 .. 5, y = -5 .. 5, x = -5 .. 5, k...

Hello, I have a problem with plotting graphs with a constant.. I have one system of equations and one another equation:

dx/dt=x(t)

dy/dt=ky(t)

and 

|y|=C|x|^k

 

I need to compare theit curves.

k and C are constants. I thought it is alright if I would remove constants, BUT there is one in the power, so I have no idea what to do.. 

I just tried to plot the system, but it doesn't work either:

Thanks a lot Markiyan, and I am sorry I may have overlooked this already answered question. But, now I have just one more thing to know. On using AllSolutions = true in solve I get the answer as x = _z1~ ∏ / k, I then try to subs my fav symbol N instead of _Z1~, but Maple doesn't take it! It keeps it as it is. What kind of constant is this _Z1~, I had thought it to be the same as _C1, which is used by Maple while solving any ODE.

Reason why I want to...

As many of you know now the MRB constant = sum((-1)^n*(n^(1/n)-1),n=1..infinity).

Here are some equations involving various forms of that summation.

The first one involves convergent series and is too obvious. The others involve divergent series.

The last two, however, are new!

 

Let c=MRB constant and a, c~, x, and y = any number.

 

sum((-1)^n*(c~*n^(1/n)-c~),n=1..infinity)= c*c~.

evalf(sum((-1)^n*(n^(1/n)-a),n=1..infinity)) gives c-1/2*(1-a).

evalf(sum((-1)^n*(x*n^(1/n)+y*n),n=1..infinity)) gives (c-1/2)*x-1/4*y.

And it appears that

evalf(sum((-1)^n*(x*n^(1/n)-a),n=1..infinity)) gives (c - 1/2)*x + 1/2*a.

> f := a+D(x); 

> convert(f(t), diff);

 

Is there a way to tell Maple that "a" is a constant, so (t) won't attach to a? thanks, kyle

The MRB constant can be computed in Maple by evalf(sum((-1)^n*(n^(1/n)-1),n=1..infinity)).

On my laptop restart; st := time(); evalf(sum((-1)^n*(n^(1/n)-1), n = 1 .. infinity), 500); time()-st gives a timming of 37.908 seconds.

Using the procedure posted at the bottom of this message st := time(); A037077(500); time()-st gives a much faster timing of 1.903 seconds.

My fastest timing for 500 digits of MRB comes from my...

A better approximation gives more digits of accuracy in the result per digit of precision used in the computation than a good approximation does.  I was wondering if anyone could come up with a better approximation to the MRB constant than 31/165, 

The MRB constant Z will probably have several parts.

The following example is from the Maple help pages
> with(GraphTheory);
> with(SpecialGraphs);
> H := HypercubeGraph(3);
DrawGraph(H)
 
 


What I would like to do in the MRB constant z,  MRB constant z part2, and etc. is to draw a series of graphs that show the some of the geometry of the MRB constant.

See http://math-blog.com/2010/11/21/the-geometry-of-the-mrb-constant/. I would like to draw a tesseract of 4 units^4, a penteract of 5 units^5, etc and take an edge from each and line the edges up as in Diagram 3:

`` 

 

 

 

As usual I'm asking for your help.

``

``

 

Download May262012.mw

 

This post can be downloaded here:  Download May202012.mw

Below we have approximations involving the MRB constant. The MRB constant plus a fraction is saved as P while a combination of another constant is saved as Q. We then subtract Q from P and always have a very small result!

1 2 3 Page 1 of 3