Let us consider the definite integral

J:=int(abs(x-(-x^5+1)^(1/5)), x = 0 .. 1);

Maple fails with it, Mathematica 10.1 finds it in terms of  special functions. Let us look at the integrand:
plot(x-(-x^5+1)^(1/5), x = 0 .. 1);

We see the expression under the modulus changes its sign at the unique point of RealRange(0,1). Therefore

solve(x-(-x^5+1)^(1/5));

Then

J:= int(-x+(-x^5+1)^(1/5), x = 0 .. (1/2)*2^(4/5))+int(x-(-x^5+1)^(1/5), x = (1/2)*2^(4/5) .. 1);

which outputs a complicated expression

(1/8)*2^(4/5)*(4*hypergeom([-1/5, 1/5], [6/5], 1/2)-2^(4/5))+(1/2)*2^(4/5)*((1/2)*2^(1/5)-(1/4)*2^(4/5))-(1/25)*Pi*csc((1/5)*Pi)*(-(25/2)*sin((1/5)*Pi)*GAMMA(4/5)*2^(4/5)*hypergeom([-1/5, 1/5], [6/5], 1/2)/Pi+(5/4)*sec((3/10)*Pi)*cos((1/10)*Pi)*2^(3/5)*Pi^(1/2)*csc((3/10)*Pi)/GAMMA(7/10))/GAMMA(4/5).

At the same time we have

int(abs(x-(-x^5+1)^(1/5)), x = 0 .. 1, numeric);

                          0.5000000000

How to obtain 1/2 symbolically?

I am trying to integrate product of exp(t+s) and a piecewise polynomial but the result can not be read and not usefull. also I used numerical integration function "Quadrature" but the result did not change.

error.mwerror.mw

Hi, I'm having a problem regarding this equation:
                                                           sum(v(i), i = 1 .. 28)

where, v is a Data Table created under the Components Tab with 28 Rows, and 1 Column.

I keep on getting an error of: 

                                                     Error, unsupported type of index, i

I didn't use or declare 'i' in any part of the program.

Thank you for your help! 

On a related post here, I have posted some screen shots for fun. Today is a special day for me in particular, because I passed my PhD viva with minor corrections! If you look colsely on the date, by coincidence, it was TWO years ago exactly. So I feel like I maybe 'update' the old post with a bit more, to make it V2.

So here comes the new screent shots: (Bear with me for now, I will post some actual content about how I feel when using Maple later.)

Alright, enough with the screent shots. If you click to see the big picture, you would notice that I have managed to get a complete sets of Maple versions from Maple V to Maple 2015!

===============================================

I am originally from Shanghai and I have always done well in Maths and science. I first heard about Maple in 2004 from a Maths teacher. He introduced me to this software. I played just a couple of times with Maple 6 on his computer, to get a first impression. At that time, all I could say was, hmmm insteresting.

In 2006, I went to UK for a foundation course. That was the first time I was actually taught how to use Maple 9.5. So I had access to it on the university computers. I discovered a lot about Maple and used it for ALL my Maths homeworks. Yes, I am lazy. I hope that Maple can do everything!

Then I got fascinated with Maple. Unfortunately, Maple was not taught in the undergraduate course. But there are materials for self-learning. By that time, I had become quite good with most of the contents in those materials. So I went look for other things to try. I went to ask my foundation course classmates, to see if they have anything on Maple and if they needed the help. One of them was at Imperial College and the questions were a bit chanlleging (finally!). That was the time I first met MaplePrimes! Hello, how have you been?

After that I have never stopped using Maple or Maple Primes. I may have been quite on the forum, but I was there. My PhD research rely hugely on the ability to compute the symbolic rank of certain matrices as well as the ability to simplifiying complicated expressions using siderules.

I do not want to talk too much about the technical details when using Maple. But I do

THANK Maple and it creators, developers and all relavent staff around the globe.

THANK Maple Primes and all its users. Some users have been particulary helpful!

Lastly, I dont have any other words to say other than this:

It's been really fun and enjoyable. Thanks!

The model of fixed-bed adsorption column

Fluid phase:

PDE:= diff(U(x, tau),tau)+ psi*Theta*diff(U(x, tau),x)-(1/Pe)*psi*Theta*diff(U(x, tau),$(x, 2))=-3*psi*xi*(U(x, tau)-Q/K);

IBC:={U(x, 0) = 0,U(0, tau) = 1+(1/Pe)*(D[1](U))(0, tau),(D[1](U))(1, tau)=0};

Particle:

PDE:= diff(Q(r, tau), tau) = diff(Q(r, tau), $(r, 2))+(2/r)*diff(Q(r, tau),r);

IBC:={Q(r, 0) = 0,(D[1](Q))(0, tau) = 0,(1/K)*(D[1](Q))(1, tau)=xi*(U-Q(1, tau)/K)};

Pe:=0.01:

psi:=6780:

Theta:=3.0:

xi:=10000:

I will really appreciate your help. Thanks in anticipation.

Hi guys,

I'm trying to draw a phase portrait based on a system of differential equations, but executing the DEplot command gives me the response in the title. 

The command I entered is this one :

DEplot(sysdif,[u(t),nu(t)],t=0..50,{[u(0)=0.831,nu(0)=0.7]},linecolor=red,numpoints=1000, thickness=1, u=0.5...1, nu=0.6...1, color=black);

the system of equation is this one :

eqd1 := diff(u(t), t) = u(t)*[-gamma + (rho)/(mu - nu(t)) - delta*u(t) + pii - (tau)/(alpha + epsilon*exp(-beta*((omega-u(t))/(kappa-u(t)))))] ;
eqd2 := diff(nu(t), t) = nu(t)*[(omega-u(t))/(kappa-u(t)) - (tau)/(alpha+epsilon*exp(-beta*((omega-u(t))/(kappa-u(t)))))] ;

the model is calibrated. I understand that maple cannot store a kind of number but even changing the parameters won't help. I've been looking for people with same error message but using solutions provided by forum members don't work. Before getting that error message i did have the one with "vars must be declared as a list ..." so I did :

u:='u' ;

nu:='nu';

but now I have the error  "unable to store  '[HFloat(0.005711776872341132)]' when datatype=float[8]".

Does anybody have an idea of the solution to my problem ?

Thanks for your time,

best regards,

Louis

I'm just using it this way:

with(QDifferenceEquations):
QPochhammer(-1,5,10)

Error, (in QDifferenceEquations:-QPochhammer) wrong type of arguments

Am I doing something terribly wrong???

Hi,

     I'm trying to numerically solve a PDE in Maple for different boundary conditions, however I'm having trouble even getting Maple to numerically solve it for simple boundary conditions.

I have cylindrical coordinates, r, z, theta, and I treat r = r(z, theta) for convenience to plot my solution surface. The initial coundary condition is that at z = epsilon (z = 0 is singular) , r = constant and of course r is periodic in theta. This is just a circle, and the analytical solution is know to be a half-sphere  r = sqrt(R^2 - z^2). I entered my initial boundary conditions into Maple, but it doesn't like the periodic one

IBC := { r(epsilon, theta) = R - epsilon__r,
              r(z, 0) = r(z, 2*Pi) };

pdsolve(
  PDE,
  IBC,
  numeric,
  indepvars = [z, theta],
  time = z,
  range = 0..2*Pi);
Error, (in pdsolve/numeric/par_hyp) Incorrect number of boundary conditions, expected 2, got 1

I'm not sure how to make this work, and then generalize it to more arbitrary intial slices r(epsilon, theta) = f(theta).

Here's the attached worksheet, ForMaplePrimesSUbmission.mw

Any help is appreciated,

Thanks

I have written some statements in maple and when integrating a vector I encounter an error . whats the problem? the error is:

Error, (in evalf/Ei/taylor) unsupported type of index, 1.

F1 := ScalarMultiply(Phi^%T, f);
F := int(F1, x = 0 .. 1);

which f is exp(x-t) and Phi is a vectror.

r := sqrt(x^2+y^2+z^2);
divE := diff(e*x/(4*Pi*r^3), x)+diff(e*x/(4*Pi*r^3), y)+diff(e*x/(4*Pi*r^3),z);
simplify(divE);

of a square matrix having its entries 1,2,...16? See http://oeis.org/A085000 for info.

I want to plot this contour 4x^2+9y^2=160  onto the same f(x,y)=4x^2+9y^2.  The error I keep getting is

Error, (in plot3d) first argument must be either in standard or parametric form

Thanks

I can;t seem to get DSOLVE working.  Will some kindly soul please send me a complete example, starting with an f(x) & clearly documented?   Many thanks.

AB

Suppose I have a file, test.txt, which the content is

2
2 0.0 0.0 0.0
1 0.0 0.0 1.6

I would like to save 

Z:=2;

C[1,1]:=2;C[2,1]:=0.0;C[3,1]:=0.0;C[4,1]:=0.0;

C[1,2]:=1;C[2,2]:=0.0;C[3,2]:=0.0;C[4,2]:=1.6;

However, if I use readdata("...",4), seems it do not work.

Is there any flexible command in maple can readfile with different column? Thank you very much!

I am having difficulty with the contourplot3d.  When I hit enter it comes back with a blank plot.  Either Im doing somthing wrong or my machine can't handle it.