Items tagged with maple

I've implemented the optimal taxation model proposed in this paper using Maple.

But it never stops running and get stuck in the last line for integral computation. Any idea of what's wrong with that?

This is the last line:

Here is the full code.

restart:
tf:=2:
n:=2:
X[1]:=t->t^2:
A[1]:=diff(X[1](t),t,t):
plot(A[1](t),t=0..tf,labels=["Time t","n times derivative of X[1](t)"]);

I have a code like above. I run it. Then,

in the labels; it was written "n times derivative of X[1](t)"] why not "2 times derivative of X[1](t)"]) ? and why was not the subscript appeared ? I mean that I want to see X1 instead of  X[1]?

Help me please.

Thank you.

 

When Maple 2017?. Are there anybody about the new features of this version?. I,m waiting for a special project.

In maple we have type/symmfunc check for a symmetric function. I want to write cycfunc function.

Example: cycfunc(a^2*b + b^2*c + c^2a, [a,b,c]) return true.

Thanks you very much.

This is simplified from a larger example.  I never saw Maple exit from a numercal calculation.

> 123456789 ^ 987654321:
Execution stopped: Stack limit reached.

Process Maple R2016 exited abnormally with code 158

> kernelopts(maxdigits);
                                  38654705646

> lprint(123456789. ^ 987654321.);
.4339566080e7991619731

kernelopts(memusage) shows about 2~3 meg bytes used.

The system has 16G Ram.  16G swap.

Seems like something maple should have caught.

restart;
eq1 := T - M*g = 0;
eq2 := T1*sin(theta1) + T2*sin(theta2) - T = 0;
eq3 := -T1*cos(theta1) + T2*cos(theta2) = 0;
soln := solve({eq1,eq2,eq3},{T,T1,T2});

limit(subs(theta2=theta1,soln[3]),theta1=0);

This returns a limit -> ... expression.  Since there is a sin(theta1) in the denominator, this should be undefined or infinity.

Hi all

I have 4 - 2a and (2x^2 - 2x + sqrt(2))^2. How to factor it become 2(2-a) and 2(sqrt(2)x^2 - sqrt(2)x + 1)^2 ?

Thanks you very much.

 

I have a PDE with boundary conditions, from a NASA paper.  I always seem to have problems expressing the bc.

PDE := diff(u(x, t), t) = (1/4)*exp(2)*exp(-u(x, t))*(diff(diff(u(x, t), x), x))/(x^2+2);

The initial/boundary conditions are

@t=0, u(x, t) = 2-2*ln(-x^2+2)

@x=0, diff(u(x,t),x)=0  ## this is the bc I have problem expressing

@x=1, u(x,t) = 2+ln(1+t)

The exact solution given in the paper:

2 + ln(1+t) - 2*ln(2-x^2)

I tried

ics := u(x, 0) = 2*(1-ln(2-x^2));

bcs := D[1](u(0,t))=0, u(1,t)=2+ln(1+t);

PDEtools[Solve]([PDE,ics,bcs]); ## no solution

How do I do this?

Tom Dean

We have just released an update to Maple.  It includes updates to the Maple Workbook, the video component, the Physics package, and many other small improvements throughout the product. It is available through Tools>Check for Updates in Maple, and is also available from our website on the Maple 2016.2 download page.

eithne

Hi all, i wan to write getCoeff() function get coefficient general.

Example

f := a^2*b^2*c^2 + 2*a^2*b^2 + 2*a^2*c^2 + 2*b^2*c^2 + a^2 - 6*a*b - 6*a*c + b^2 - 6*b*c + c^2 + 8

getCoeff(f,a^2*b^2*c^2) return 1.

getCoeff(f,a^2*b^2) return 2.

getCoeff(f,a*b) return -6.

getCoeff(f,a^2) return 1, ...

and how to get coefficients freedom ?

Thank you very much.

I'd like to pay attention of Maple community to the recent work by Alex Degtyarev in algebraic geometry done with Maple.

Bertini.zip

The Joint Mathematics Meetings are taking place this week (January 4 – 7) in Atlanta, Georgia, U.S.A. This will be the 100th annual winter meeting of the Mathematical Association of America (MAA) and the 123nd annual meeting of the American Mathematical Society (AMS).

Maplesoft will be exhibiting at booth #118 as well as in the networking area. Please stop by our booth or the networking area to chat with me and other members of the Maplesoft team, as well as to pick up some free Maplesoft swag or win some prizes.

There are also several interesting Maple-related talks and events happening this week:

 

Teaching Cryptology to Increase Interest in Mathematics for Students Majoring in Non-Technical Disciplines and High School Students

Wednesday, January 4, 0820, L401 & L402, Lobby Level, Marriott Marquis

Neil Sigmon, Radford University

 

Enigma: A Combinatorial Analysis and Maple Simulator

Wednesday, January 4, 0900, L401 & L402, Lobby Level, Marriott Marquis

Rick Klima, Appalachian State University

 

MYMathApps Calculus - Building on Maplets for Calculus

Thursday, January 5, 0800, Courtland, Conference Level, Hyatt Regency

Philip B. Yasskin, Texas A&M University 
Douglas B. Meade, University of South Carolina 
Andrew Crenwelge, Texas A&M University

 

Maple Software Technology as a Stimulant Tool for Dynamic Interactive Calculus Teaching and Learning

Thursday, January 5, 1000, Courtland, Conference Level, Hyatt Regency

Lina Wu, Borough of Manhattan Community College-The City University of New York 

 

Collaborative Research: Maplets for Calculus

Thursday, January 5, 1400, Marquis Ballroom, Marquis Level, Marriott Marquis

Philip Yasskin, Texas A&M University 
Douglas Meade, U of South Carolina

 

Digital Graphic Calculus Art Design in Maple Software

Thursday, January 5, 1420, International 7, International Level, Marriott Marquis

Lina Wu, Borough of Manhattan Community College-The City University of New York 

 

Maplesoft will also be hosting a catered reception and brief presentation on Teaching STEM Online: Challenges and Solutions, Thursday January 5th, from 6:00pm – 7:30pm, at the Hyatt Regency, Hanover AB, on the exhibitor level. Please RSVP at www.maplesoft.com/jmm or at Maplesoft booth #118.

 

If you are attending the Joint Math meetings this week and plan on presenting anything on Maple, please feel free to let me know and I'll update this list accordingly.


See you in Atlanta!

Daniel

Maple Product Manager

I'm trying to generate 100 equation to solve them numerically , buy maple won't let me , here's the code :

for i by 0 to 99 do
a := a[i]+a[i+1] = i^2
end do ;

Let L = {a1, a2, ..., an} (or L = [a1, a2, ..., an]) with is homogeneous polynomial. How to create function createDegree() return elements of degree <= n from L.

Expamle. L = {a-2b, b^2, (a+c)^2}.

createDegree(L,1) return {a - 2b}.

createDegree(L,2) return {b^2, (a - 2b)^2, (a+c)^2, a - 2b}.

createDegree(L,3) return {b^2, (a - 2b)^2, (a - 2b)^3, (a+c)^2, (a - 2b)b^2, (a - 2b)(a+c)^2, a - 2b}.

Thanks you very much.

 

 

Is this a bug?

hypergeom([1, -1, 1/2], [-12,-3], 1);
Error, (in hypergeom/check_parameters) function doesn't exist: missing appropriate negative integers in the first list of parameters to compensate the negatives integer(s): [-3], found in the second list.
 

Yet this hypergeometric series terminates and Maple should be able to handle it, at least according to the Maple help page (the second rule below applies, yet the numerator has a smaller absolute value, so the first rule below applies).

If some   n[i] is a non-positive integer, the series is finite (that is,   F(n, d, z)  is a polynomial in    z).
If some  d[j]  is a non-positive integer, the function is undefined for all non-zero  z, unless there is also a negative upper parameter of smaller absolute value, in which case the previous rule applies.
 

 

Interestingly, the Wolfram Mathematica app can evaluate this to 311/312.

 

 

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