Items tagged with mathematica mathematica Tagged Items Feed

I would like to announce a new unofficial record computation of the MRB constant that was finished on Sun 21 Sep 2014 18:35:06.

I really would like to see someone beat it with Maple!

It took 1 month 27 days 2 hours 45 minutes 15 seconds. I computed 3,014,991 digits of the MRB constant, (confirming my previous 2,00,000 or more digit computation was actually accurate to 2,009,993 digits), with Mathematica 10.0. I Used my version of Richard Crandall's code:

 

____________________________________________________________________________

(*Fastest (at MRB's end) as of 25 Jul 2014.*)

DateString[]

prec = 3000000;(*Number of required decimals.*)ClearSystemCache[];

T0 = SessionTime[];

expM[pre_] := 

  Module[{a, d, s, k, bb, c, n, end, iprec, xvals, x, pc, cores = 12, 

    tsize = 2^7, chunksize, start = 1, ll, ctab, 

    pr = Floor[1.005 pre]}, chunksize = cores*tsize;

   n = Floor[1.32 pr];

   end = Ceiling[n/chunksize];

   Print["Iterations required: ", n];

   Print["end ", end];

   Print[end*chunksize]; d = ChebyshevT[n, 3];

   {b, c, s} = {SetPrecision[-1, 1.1*n], -d, 0};

   iprec = Ceiling[pr/27];

   Do[xvals = Flatten[ParallelTable[Table[ll = start + j*tsize + l;

        x = N[E^(Log[ll]/(ll)), iprec];

        pc = iprec;

        While[pc < pr, pc = Min[3 pc, pr];

         x = SetPrecision[x, pc];

         y = x^ll - ll;

         x = x (1 - 2 y/((ll + 1) y + 2 ll ll));];(*N[Exp[Log[ll]/ll],

        pr]*)x, {l, 0, tsize - 1}], {j, 0, cores - 1}, 

       Method -> "EvaluationsPerKernel" -> 4]];

    ctab = ParallelTable[Table[c = b - c;

       ll = start + l - 2;

       b *= 2 (ll + n) (ll - n)/((ll + 1) (2 ll + 1));

       c, {l, chunksize}], Method -> "EvaluationsPerKernel" -> 2];

    s += ctab.(xvals - 1);

    start += chunksize;

    Print["done iter ", k*chunksize, " ", SessionTime[] - T0];, {k, 0,

      end - 1}];

   N[-s/d, pr]];

t2 = Timing[MRBtest2 = expM[prec];]; DateString[]

Print[MRBtest2]

MRBtest2 - MRBtest2M

_________________________________________________________________________.

 

I used a six core Intel(R) Core(TM) i7-3930K CPU @ 3.20 GHz 3.20 GHz with 64 GB of RAM of which only 16 GB was used.

t2 From the computation was {1.961004112059*10^6, Null}.

 

 

 

Hi all

I have a mathematical problem and I asked it in various sites but the answers till yet are not correct.

Assume that we have:

T[m]:=t->t^m:
b[n,m]:=unapply(piecewise(t>=(n-1)*tj/N and t<n*tj/N, T[m](N*t-(n-1)*tj), 0), t):

where n,N,tj are known constants. furthermore assume that we want to comute the following integral:

for following approximations:

I have written the following code but it seems to be incorrect:

V1:=Vector([seq(seq(b[n,m](t),m=0..1),n=1..3)]);
V:=evalf(V1.Transpose(V1));

the original program is :

taaylor.mws

I will be so grateful if any one can help me to solve it by maple

Mahmood   Dadkhah

Ph.D Candidate

Applied Mathematics Department

hi every one,i want to know what are the maple strentgh compraed to matlab and mathematica ? why someone should choose maple instead of these both ? i searched the net,but nothing useful has been found , tnx in advance .

Hi All,

I would like to request information the representation of the following result from Mathematica :

Mathematica result:    MeijerG[{{0, 1/2}, {}}, {{0, 1}, {-1, -1}}, a, 1/2]

Maple is able to take: MeijerG[{{0, 1/2}, {}}, {{0, 1}, {-1, -1}}, a] which is represented as
                                MeijerG([[0, 1/2], []], [[0, 1], [-1, -1]], a)

Could you please advise me on how to implement this function in maple. I would be grateful if you also include matlab in the discussion.

The matiematica output is a result of the following integration:

Integrate[r*(BesselI[1, al*r]* BesselK[1, al*r]), r]

 

Looking forward to your reply. Thanks in advance.

Regards

Raj

HelloI have Mathematica code(notebook) I want to convert it to Maple code (worksheet)
 Thanks 

Does anyone know if Maple has a counterpart to the Wolfram Lightweight Grid Manager that can be used with Maple 14 or Maple 16?

The latest Maple 17 bashing from Mathematica appears on their comparison webpage.

The first thing that paid my attention in their "comparison" was their comparison of the Mobius project to the Mathematica demonstrations project.  To my knowledge the Mobius project is an extension to application center which I thought was the same as Mathematica demonstration project anyways ???

Anyways anyone have thoughts on this recent comparison?

hi dear

i am a student of m.phil in maths.

i am doing my research in numerical solution of fractional derivatives but

i am facing a little problem in writing the code for fractional derivatives by hpm

in mathematica software if you are having any code kindly sent it to me

thanks

I came across this post http://www.stiglerdiet.com/2012/02/10/mathematica-a-love-story/

and in the commented section near the end someone mentions this can not be elegantly done in Maple.  I would like to disprove this comment.  He says "There’s NO WAY any numeric language like R or Maple can try to do this."

I don't believe it, and I welcome anyone to disprove his...

Hello,

I've been trying to force maple to find a solution of this differential equation:

ode := -((D@@2)(y))(r)-2*y(r)/r+y(r) = (2-V(r))*u[0](r)

dsolve(ode, y(r));

gives

y(r) = r exp(-r) _C2 + (exp(r) + 2 exp(-r) Ei(1, -2 r) r) _C1 - 4 

/ / /
|1 | |
exp(-3 r) |- | |
|2 | |
\ \/

/ 2 \ \

Hi,

I'm new in Maple and i've a project to do. I've to find at what time is 90% of the steady-state concentration of C achieved. And the eq is 

z'=k(1/2-z)(1/3-z), where k=0.3. 

So i had worked in Mathematica before and i know that it would look something like that in there: 

eqn1 = 0 == k (1/2 - z[t]) (1/3 - z[t])

DSolve[Evaluate[{eqn1}], z[t]].

But i have no idea how to find this 90%. I hope that someone will...

Hi,

I am trying to do somthing like this (I actually have a 5D Array, this is a simple example) (mathematica code)

a = { {{1,2,3},{4,5,6}},{{1,2,3},{4,5,6}},{{1,2,3},{4,5,6}}};

Transpose[a, {1,3,2}]

 (* this {1,3,2} dictates what gets transposed e.g. {1,2,3}, would leave the Array the same, and the example transposes the last two dimentions*)

in maple.

After some research it seems it makes sense in maple to use an...

Mathematica 9 was just released Nov 28 which caused me, of course, to have a look to see what wolfram has done in way of comparing itself to maple. 

Perhaps we can discuss, compare, debunk, clarify some of their claims or perhaps dismiss the results altogether

Here is one speed comparison chart I came across.  Is it relavent?  It shows Maple 16 as the worst performer.

 

Dear friends,

I was working on the integral of cosh(ax)/cosh(x) from zero to infinity with "a" being a rational number p/q where p<q and p-q being odd. I was quite shocked to discover that Maple 15 (X86 64 LINUX)cannot do this integral while Mathematica gets it right every time. I certainly hope this will be rectified soon.

For those who might be interested the calculation of the integral

there is a link

1 2 3 Page 1 of 3