Items tagged with maple

Feed

Hi,

I'm not sure that I mean datatable component corectly.

I also consider that I was done somthing wrong

Thank you for advanced for any help.

restart

with(DocumentTools):

Oryginaly DataTable was inserted as a 3 x 3. I will traing to push maple to obtain 4 x 4 with specific row and column name.

``

SetProperty("DataTable0", visibleRows, 4);

DocumentTools:-SetProperty("DataTable0", visibleColumns, 4);

DocumentTools:-SetProperty("DataTable0", columnWidths, [20, 40, 80, 80]);

DocumentTools:-SetProperty("DataTable0", rowNames, [r1, r2, r3, r4]);

DocumentTools:-SetProperty("DataTable0", columnNames, [c1, c2, c3, c4]);

DocumentTools:-SetProperty("DataTable0", update)

``

``

``


wzel

Download datatable_problem.mw

The attached worksheet shows a small selection of new and improved results in integration for Maple 2016. Note that integration is a vast topic, so there will always be more improvements that can be made, but be sure that we are working on them.

Maple2016_Integration.mw

A selection of new and improved integration results for Maple 2016

New answers in Maple 2016

 

 

Indefinite integrals:

 

int(sqrt(1+sqrt(z-1)), z);

(4/5)*(1+(z-1)^(1/2))^(5/2)-(4/3)*(1+(z-1)^(1/2))^(3/2)

(1.1)

int(arctan((-1+sec(x))^(1/2))*sin(x), x);

-arctan((-(1/sec(x)-1)*sec(x))^(1/2))/sec(x)+(1/2)*(-1+sec(x))^(1/2)/sec(x)+(1/2)*arctan((-1+sec(x))^(1/2))

(1.2)

int(((1+exp(I*x))^2+(1+exp(-I*x))^2)/(1-2*c*cos(x)+c^2), x);

-x-2*x/c-x/c^2+I*exp(I*x)/c-I*exp(-I*x)/c-I*c*ln(exp(I*x)-1/c)/(c-1)-I*ln(exp(I*x)-1/c)/(c-1)-I*ln(exp(I*x)-1/c)/(c*(c-1))-I*ln(exp(I*x)-1/c)/(c^2*(c-1))+I*c*ln(-c+exp(I*x))/(c-1)+I*ln(-c+exp(I*x))/(c-1)+I*ln(-c+exp(I*x))/(c*(c-1))+I*ln(-c+exp(I*x))/(c^2*(c-1))

(1.3)

int(x^4/arccos(x)^(3/2),x);

(1/4)*(-x^2+1)^(1/2)/arccos(x)^(1/2)-(1/4)*2^(1/2)*Pi^(1/2)*FresnelC(2^(1/2)*arccos(x)^(1/2)/Pi^(1/2))+(3/8)*sin(3*arccos(x))/arccos(x)^(1/2)-(3/8)*2^(1/2)*Pi^(1/2)*3^(1/2)*FresnelC(2^(1/2)*3^(1/2)*arccos(x)^(1/2)/Pi^(1/2))+(1/8)*sin(5*arccos(x))/arccos(x)^(1/2)-(1/8)*2^(1/2)*Pi^(1/2)*5^(1/2)*FresnelC(2^(1/2)*5^(1/2)*arccos(x)^(1/2)/Pi^(1/2))

(1.4)

 

Definite integrals:

int(arcsin(sin(z)), z=0..1);

1/2

(1.5)

int(sqrt(1 - sqrt(1+z)), z=0..1);

((4/5)*I)*(2^(1/2)-1)^(3/2)*2^(1/2)+((8/15)*I)*(2^(1/2)-1)^(3/2)

(1.6)

int(z/(exp(2*z)+4*exp(z)+10),z = 0 .. infinity);

(1/20)*dilog((I*6^(1/2)-3)/(-2+I*6^(1/2)))-((1/60)*I)*6^(1/2)*dilog((I*6^(1/2)-3)/(-2+I*6^(1/2)))+(1/20)*dilog((I*6^(1/2)+3)/(2+I*6^(1/2)))+((1/60)*I)*6^(1/2)*dilog((I*6^(1/2)+3)/(2+I*6^(1/2)))+((1/120)*I)*6^(1/2)*ln(2+I*6^(1/2))^2-((1/120)*I)*6^(1/2)*ln(2-I*6^(1/2))^2+(1/40)*ln(2+I*6^(1/2))^2+(1/40)*ln(2-I*6^(1/2))^2+(1/60)*Pi^2

(1.7)

simplify(int(sinh(a*abs(x-y)), y=0..c, 'method'='FTOC'));

(1/2)*(piecewise(x < 0, 0, 0 <= x, 2*exp(-a*x))+piecewise(x < 0, 0, 0 <= x, -4)+2*piecewise(c <= x, -cosh(a*(-x+c))/a, x < c, (cosh(a*(-x+c))-2)/a)*a-exp(-a*x)+piecewise(x < 0, 0, 0 <= x, 2*exp(a*x))+4-exp(a*x))/a

(1.8)

int(ln(x+y)/(x^2+y), [x=0..infinity, y=0..infinity]);

infinity

(1.9)


Definite integrals with assumptions on the parameters:

int(x^(-ln(x)),x=0..b) assuming b > 0;

(1/2)*erf(ln(b)-1/2)*Pi^(1/2)*exp(1/4)+(1/2)*Pi^(1/2)*exp(1/4)

(1.10)

int(exp(-z)*exp(-I*n*z)*cos(n*z),z = -infinity .. infinity) assuming n::integer;

undefined

(1.11)


Integral of symbolic integer powers of sin(x) or cos(x):

int(sin(x)^n,x) assuming n::integer;

` piecewise`(0 < n, -(Sum((Product(1+1/(n-2*j), j = 1 .. i))*sin(x)^(n-2*i-1), i = 0 .. ceil((1/2)*n)-1))*cos(x)/n+(Product(1-1/(n-2*j), j = 0 .. ceil((1/2)*n)-1))*x, n < 0, (Sum((Product(1-1/(n+2*j+1), j = 0 .. i))*sin(x)^(n+2*i+1), i = 0 .. -ceil((1/2)*n)-1))*cos(x)/n+(Product(1+1/(n+2*j-1), j = 1 .. -ceil((1/2)*n)))*ln(csc(x)-cot(x)), x)

(1.12)

int(cos(x)^n,x) assuming n::negint;

-(Sum((Product(1-1/(n+2*j+1), j = 0 .. i))*cos(x)^(n+2*i+1), i = 0 .. -ceil((1/2)*n)-1))*sin(x)/n+(Product(1+1/(n+2*j-1), j = 1 .. -ceil((1/2)*n)))*ln(sec(x)+tan(x))

(1.13)

int(cos(x)^n,x) assuming n::posint;

(Sum((Product(1+1/(n-2*j), j = 1 .. i))*cos(x)^(n-2*i-1), i = 0 .. ceil((1/2)*n)-1))*sin(x)/n+(Product(1-1/(n-2*j), j = 0 .. ceil((1/2)*n)-1))*x

(1.14)

Improved answers in Maple 2016

 

int(sqrt(1+sqrt(x)), x);

(4/5)*(1+x^(1/2))^(5/2)-(4/3)*(1+x^(1/2))^(3/2)

(2.1)

int(sqrt(1+sqrt(1+z)), z= 0..1);

-(8/15)*2^(1/2)-(8/15)*(1+2^(1/2))^(3/2)+(4/5)*(1+2^(1/2))^(3/2)*2^(1/2)

(2.2)

int(signum(z^k)*exp(-z^2), z=-infinity..infinity) assuming k::real;

(1/2)*(-1)^k*Pi^(1/2)+(1/2)*Pi^(1/2)

(2.3)

int(2*abs(sin(x*p)*sin(x)), x = 0 .. Pi) assuming p> 1;

-2*(sin(Pi*p)*signum(sin(Pi*p))*cos(Pi/p)-p*sin(Pi/p)*cos(Pi*(floor(p)+1)/p)+sin(Pi*(floor(p)+1)/p)*cos(Pi/p)*p-sin(Pi*p)*signum(sin(Pi*p))-sin(Pi*(floor(p)+1)/p)*p+sin(Pi/p)*p)/((cos(Pi/p)-1)*(p^2-1))

(2.4)

int(1/(x^4-x+1), x = 0 .. infinity);

-(sum(ln(-_R)/(4*_R^3-1), _R = RootOf(_Z^4-_Z+1)))

(2.5)


In Maple 2016, this multiple integral is computed over 3 times faster than it was in Maple 2015.

int(exp(abs(x1-x2))*exp(abs(x1-x3))*exp(abs(x3-x4))*exp(abs(x4-x2)), [x1=0..R, x2=0..R, x3=0..R, x4=0..R], AllSolutions) assuming R>0;

(1/8)*exp(4*R)-29/8+(7/2)*exp(2*R)-5*R*exp(2*R)+2*exp(2*R)*R^2-(5/2)*R

(2.6)

Austin Roche
Mathematical Software, Maplesoft


with(PDEtools, casesplit, declare)
``

L := 1651.12; m := 3205.12; r1 := .1875; r2 := 2; z1 := 0; z2 := 12; ld := 4.5

NULL

declare(u(r, z), w(r, z))``

with(DEtools, gensys)

rr := (L+2*m)*(diff(u(r, z), r))+L*(diff(w(r, z), z))+L*u(r, z)/r

zz := L*(diff(u(r, z), r))+(L+2*m)*(diff(w(r, z), z))+L*u(r, z)/r

rz := m*(diff(u(r, z), z))+m*(diff(w(r, z), r))

BCS := {rr(r1, ld) = 0, rz(r1, z) = T, w(r, 0) = 0, zz(r, z2) = 0}

{3205.12*(diff(u(r, z), z))(.1875, z)+3205.12*(diff(w(r, z), r))(.1875, z) = T, 8061.36*(diff(u(r, z), r))(.1875, 4.5)+1651.12*(diff(w(r, z), z))(.1875, 4.5)+1651.12*(u(r, z))(.1875, 4.5)/r(.1875, 4.5) = 0, 1651.12*(diff(u(r, z), r))(r, 12)+8061.36*(diff(w(r, z), z))(r, 12)+1651.12*(u(r, z))(r, 12)/r(r, 12) = 0, w(r, 0) = 0}

(1)

``

NULL

sys3 := [(L+2*m)*(diff(u(r, z), r, r))+(L+m)*(diff(w(r, z), r, z))+(L+2*m)*(diff(u(r, z), r))/r-(L+2*m)*u(r, z)/r^2+m*(diff(u(r, z), z, z)) = 0, (L+m)*(diff(u(r, z), r, z))+m*(diff(w(r, z), r, r))+(L+2*m)*(diff(w(r, z), z, z))+(L+m)*(diff(u(r, z), z))/r+m*(diff(w(r, z), r))/r = 0]

pdsolve(sys3, BCS, numeric)

 

 

``

``


Download PDE_equation2.mw

Hi all,

I have the following PDE, is it solveable by Maple or not. Do I need a boundary condition and how many or I can get a general solution? I am new to Maple. Any help will be appreciated.

Thank you.

 

 

 

Can we get it in MapleSim, not in exactly this form, but in substance? (Not in Maple)
The line of intersection of surfaces:
(x1-0.5) ^ 4 + x2 ^ 4 + x 3 ^ 4-1. ^ 2 = 0.;
x1 ^ 2 + (x2-0.25) ^ 2 + x3 ^ 2-1. ^ 2 = 0.;
(Red) rotates about an axis oX3. During rotation, the line intersects with the fixed sphere ((0., 1.5, 0 .5); R = 1.725). One of the points of intersection is drawn in green. Green Dot and the center of the sphere connected to the blue segment.  In the sphere  of  fixed  trajectory of  the green point.
In other words, the geometric model  3d  cam mechanism and its kinematics.


Here is a simple program.

t:=(a/b)^2*x*y;

st:=cat(``,`t:=`,t,`;`);

fprintf(`outfile`,"%A\n",st);

Running this on Maple 2015 gives a correct Maple statement for t in outfile which can be read into another Maple session.

Running on another computer with Maple 17 gives

t:= || (a^2/b^2*x*y) || ;

in outfile which is not a valid statement.

I know that I can avoid this problem by using the save statement but I want to understand why the code gives the wrong result in Maple 17 and how to change it so it works on both versions.

Thanks.

Hello All.

Why Maple can’t do this Simple indefinite integral?

I'm have a integral :

 

Please compare to Mathematica:

Thanks in advance for your help.

 

I_Mariusz

test.mw

This is really weird.

restart;
eq:=phi(f(x,y(x),diff(y(x),x)),g(x,y(x),diff(y(x),x)));

Does any one see a "||" in the above? I do not. Then why Maple shows this:

r := cat("dsolve(",eq,",y(x))");

       r:="dsolve("||phi(f(x,y(x),y'(x)),g(x,y(x),y'(x)))||",y(x))"

This is Maple 2015 on windows. Running in worksheet mode

Here is screen shot

 

This problem does not show up all the time. For example, this works ok:

 

No "||" added. The problem with all the above, is that I can't make it all as a string in the first example above. I am trying to make a string of   "dsolve(" + ode + "),y(x)"

Here is the display options

 

thank you

Hello,

 

I'm trying to solve the integral u/(1-u) with Maple and noticed that it returned a result that doesn't accord to the solution I found by hand or the solution from WolframAlpha. This is a screenshot of the weird behaviour:

Does Maple do any weird conversions? Or did I do something wrong or is Maple wrong?

Thanks in advance,

As the year draws to a close, we start looking forward to a new year and a new release of Maple. With every new release comes many new features and updates to explore.

We are looking for several new beta testers with a good working knowledge of Maple; We need your input, your ideas, and your experience with our products to help us improve the software and get it ready for general release.

There are many benefits to becoming a beta tester:

  • You’ll get to use the new software before anyone else does.
  • You’ll help us make our software better in ways that work for you.
  • Your suggestions could determine the future direction of the software.
  • You’ll get feedback right from the development team.

If you are interested in becoming a beta tester for the next version of Maple, please email: beta (at) maplesoft.com for more information.

hi,

 

having some issue using the codegeneration function in maple. I have a script written in Maple that I would like to do in Python 2.7. I knew to maple and not sure if my syntax is incorrect or if the program doesnt support the conversion, as an example,

 

q[Fick] := proc (t) options operator, arrow; -hm0(t)*(ps(T[C5](t))-p[inf]) end proceq2 := m[C5](t)*Cp[C5l](T[C5](t))*(diff(T[C5](t), t)) = q[heat](t)

 

thanks

Hi

Hope somebody can help me! I installed Maple on my computer, but everything is very small , for example the bar in the top with all the small icons. It is the same thing when I make plots. The values on the axes are very small. How can I change this?

I have taken a photo of the problem, don't know if it helps.

The Maple 2015.2 update is now available for Windows and Linux users. It provides official support for Windows 10, connectivity to the latest MATLAB release, support for MapleSim 2015.2, updates to Physics, and a variety of small improvements throughout the product.  It is available through the Check for Updates system, and is also available from our website on the Maple 2015.2 download page.

The Mac update will be released as soon as we have finished fixing the problems on Mac OS X 10.11. We’re working on it.

eithne

We are looking for enthusiastic Maple users to become Maple Ambassadors, to inspire and educate others about the benefits that Maple brings to education.

 

As an Ambassador, you will have the opportunity to influence the development of Maple through regular meetings with Maplesoft developers, get advance news of upcoming features and products, get assistance with Maple events on your campus, and more. In return, we ask that you do what you are probably already doing – sharing your experiences with Maple, answer questions on forums (like this one!), sharing your Maple applications, providing us with feedback, etc.

 

You can find more information and an application form at Maple Ambassador Program. We’re looking forward to hearing from you!

 

Daniel

Maple Product Management

Hi,

     I have a list of 603 integrals that I want to evaluate. Unfortunately, I can't get Maple to do most of them. Mathematica can do some that Maple can't, and returns an answer in terms of BesselJ functions. So my question is 2-fold

1) Is there a way to make Maple do this integral?
2) If not, is there a way to efficiently convert 603 expessions to Mathematica and back?

 

EXAMPLE INTEGRAL
restart;
assume(k1::real, k2::real, R::real, R>0);
a :=cos(x)*exp(I*(k1*R*sin(x)+k2*R*sin(x)-4*x))*sin(x):
int(a, x=-Pi/2..Pi/2) assuming real;


Thanks! 

restart;

assume(k1::real, k2::real, R::real, R>0);

a :=cos(x)*exp(I*(k1*R*sin(x)+k2*R*sin(x)-4*x))*sin(x)

cos(x)*exp(I*(k1*R*sin(x)+k2*R*sin(x)-4*x))*sin(x)

(1)

int(a, x=-Pi/2..Pi/2) assuming real;

int(cos(x)*exp(I*(k1*R*sin(x)+k2*R*sin(x)-4*x))*sin(x), x = -(1/2)*Pi .. (1/2)*Pi)

(2)

Mathematica Answer

ans := -(1/((k1 + k2)^6*R^6))*2*I*Pi*
(
10*(k1 + k2)^4*Pi*R^4*BesselJ(2, sqrt((k1 + k2)^2*R^2))
+ 2*Pi ((k1 + k2)^2*R^2)^(3/2) (-30 + (k1 + k2)^2*R^2) *BesselJ(3, sqrt((k1 + k2)^2*R^2))
- (k1 + k2)^4*R^4*(-(k1 + k2)*R*cos((k1 + k2)*R) + sin((k1 + k2)*R))
+ 8*(k1 + k2)^2*R^2*(-(k1 + k2)*R*(-6 + (k1 + k2)^2*R^2)*cos((k1 + k2)*R) + 3*(-2 + (k1 + k2)^2*R^2)*sin((k1 + k2)*R))
- 8*(-(k1 + k2)*R*(
120 - 20*k2^2*R^2 + k1^4*R^4 + 4*k1^3*k2*R^4 +

 k2^4*R^4 + 4*k1*k2*R^2*(-10 + k2^2*R^2) +

 k1^2*(-20*R^2 + 6*k2^2*R^4))*cos((k1 + k2)*R) +

 5*(24 - 12*k2^2*R^2 + k1^4*R^4 + 4*k1^3*k2*R^4 + k2^4*R^4 +

 4*k1*k2*R^2*(-6 + k2^2*R^2) +

 6*k1^2*R^2*(-2 + k2^2*R^2))*sin((k1 + k2)*R)
)
);

-(2*I)*Pi*(10*(k1+k2)^4*Pi*R^4*BesselJ(2, (k1+k2)*R)+2*Pi((k1+k2)^2*R^2)^(3/2)*BesselJ(3, (k1+k2)*R)-(k1+k2)^4*R^4*(-(k1+k2)*R*cos((k1+k2)*R)+sin((k1+k2)*R))+8*(k1+k2)^2*R^2*(-(k1+k2)*R*(-6+(k1+k2)^2*R^2)*cos((k1+k2)*R)+3*(-2+(k1+k2)^2*R^2)*sin((k1+k2)*R))+8*(k1+k2)*R*(120-20*R^2*k2^2+k1^4*R^4+4*k1^3*k2*R^4+k2^4*R^4+4*k1*k2*R^2*(R^2*k2^2-10)+k1^2*(6*R^4*k2^2-20*R^2))*cos((k1+k2)*R)-40*(24-12*R^2*k2^2+k1^4*R^4+4*k1^3*k2*R^4+k2^4*R^4+4*k1*k2*R^2*(R^2*k2^2-6)+6*k1^2*R^2*(R^2*k2^2-2))*sin((k1+k2)*R))/((k1+k2)^6*R^6)

(3)

 

 


Download ToughIntegral.mw

Hello, forgive me if I used bad english, I am not a native speaker.

 

Anyhow: I have to decide the set of coordinates looking at the grafs intersections.  (I hope you understand that)

 

This is my function: f(x)=1/4*x3-x2-x+4

 

When you plot the function you see that the intersections is (-2,0) (2,0) and (4,0)

 

BUT I have to use a command to find these three intersections, plzz help me!

 

THANK YOU indeed.

 

 

[point of] intersection

1 2 3 4 5 6 7 Last Page 3 of 250