The algorithm that I need to replicate is as follows:

real function f(x,y)

integer n; real a,b,c,x,y

f<-max(|x|,|y|)

a<-min(|x|,|y|)

for n=1 to 3 do

b<-(a/f)^2

c<-b/(4+b)

f<-f+2*c*f

a<-ca

end for

end function f

How can I define f,a as  functions that I am later using as variables(in f=f+2cf,b=(a/f)^2)? also, is n just a variable for iteration? 

Heun functions arise in the solutions of various differential equations, for instance for the Schroedinger equation for the hydrogen atom in physics, which is also of chemical interest.  Although they have been nominally included in Maple for several years, they are still in a primitive state; despite their obscurity and intractable nature, there seems not to exist much possibility, within Maple, to convert these functions into better known and characterised functions.  A similar condition holds for Lame and spheroidal functions that are invaluable in the solution of differential equations in physics but are not even mentioned in Maple. 

The compilation of mathematical functions by Abramowitz and Stegun was published half a century ago, but there are still important functions explained therein that are lacking from Maple, not to mention the successor in the NIST Digital Library of Mathematical Functions.

Integral equations are another weak component of Maple; the present content relies on a basis of work of ProfessorCorless and his student submitted to the 'Maple Share Library' -- decades ago.  Forty years ago, David Stoutemyer generated some procedures to solve non-linear integral equations in Reduce, but forty years later Maple has no benefit from that knowledge.

We can only hope that Maple 19 will remedy some of these gross deficiencies.  The teaching, learning and practice of physics will benefit from their implementation.

At the internet site of The Heun Project, a strong declaration is made that only Maple incorporates Heun functions, which arise in the solution of differential equations that are extremely important in physics, such as the solution of Schroedinger's equation for the hydrogen atom.  Indeed solutions appear in Heun functions, which are highly obscure and complicated to use because of their five or six arguments, but when one tries to convert them to another function, nothing seems to work.  For instance, if one inquires of FunctionAdvisor(display, HeunG), the resulting list contains

"The location of the "branch cuts" for HeunG are [sic, is] unknown ..." followed by several other "unknown" and an "unable". Such a solution of a differential equation is hollow.

Incidentally, Maple's treatment of integral equations is very weak -- only linear equations with simple solutions, although procedures by David Stoutemyer from 40 years ago are available to enhance this capability.

When can we expect these aspects of Maple to work properly, for applications in physics?

hi, is there a way to collect all commands in one place which have been used during a clickable math session of a document? in fact this would help in creating automated tasks.

How I can get a list of all label refernces of a worksheet. In fact I want to write a small code which would produce the LaTex output (written to a text file) for my indicated labels. thanx

Dear All

I have a question:

How can I computet the laplace inverse of a function that is below!

My version of Maple is 13:

It should be noted that r>0, \rho>0.5, y\in R.

> restart;
> with(MTM);
print(??); # input placeholder
> ilaplace(r^(rho-.5)*abs(y)^(rho-1)*exp(-r*y)/GAMMA(rho-.5), s);

I am new to Maple. I am working though the manual and in chapter two I tried to get the derivative of ln(x^2+1). I am getting something completely different that the actual derivative. I tried a simple derivatie (the derivative of 2x and get 0). I do not know what I am doing wrong. Any help will be appreciated.

Dear mapleprimes users,

I have a problem concerning this function:

Naief := proc(A::integer,B::integer, p::posint)
local x, a:= A mod p, b:= B mod p;

 for x from 1 to p-1 do
    if a^x mod p = b then return x end if
end do;
print(¨No solution.¨)

end proc

It works fine for what it should do, finding x for a^x = b mod p by inserting x from 1 to p-1 until it finds an
apprioperate x.

My problem is concerning its computation time, which I like to calculate with:

Codetools:- Usage(Naief(a,b,p), output = [cputime], quiet, iterations= 2^12)

The problem is that keeps repeating

alot even when a,b and p are immensely huge. I don`t know how to fix it, because I need real CPUtimes which increase
a,b and p increase. I get the idea that the values of the CPUtime are not realistic.

Thanks for the help!

we always have subscript variable in the math book, but how could this be natral done in maple

I want to get a seq a₁ a₂ a₃

seq(a[i],i=1..3)

but how could I get a  a_ij

seq(a[i_j],i=1..3);

and

seq(a[ij],i=1..3);  both was not right

I am currently running into an issue where the numerical solution to an equation (involving an integral, yes, but the value I am solving for is simly a constant in the integral) is taking significantly longer than I would hope it would to solve. I am solving a similar equation (with a simpler expression) and it is significantly easier to solve, and I am hoping for that kind of speed.

On the last two lines in the attachment, I have two expressions. The penultimate expression is the baseline speed that I would like to match. The last expression is the fsolve I would like to speed up.

Is there any way to numerically speed up the process? I found that when I did tracelast after halting the process, there were HUGE numbers being added and subtracted, multiplied and divided. Not only did this significantly slow down the proecss but it also adds much numerical instability, which I would also like to avoid.

All help would be appreciated.

pole-dragging-mapleprimes.mw

how to compute non-reduced hilbert series?

hilbertseries return a simplified version of hilbert series,

how to output non-simplified version of hilbert series?

Hi all,

I google and found a program using C# connect with Maple. The Maple file is mla file - a pakage library type of Maple. I want to review data structure and all interfaces funtions to understand the way to implement this features.

Please help me the way to read the original Maple code. I uploaded the .mla file into mediafire if you want to review it. Link http://www.mediafire.com/download/abd9kpk3q6oq8lb/CHEMISTRY.mla 

Regards,

Quan Nguyen

Why the disparity in the solution of P1 and P2? P1 uses finite difference while P2 uses midrich.

See them here P1.mw     and   P2.mw

I want to put the solution obtained from ShowSolution() command of Student[Calculus1] into my LaTeX file. How I can achieve that?

I succeeded in getting latex of ExpandSteps because it was a module, but ShowSolution() is not a module.

It's Alan's birthday and he's having a

party. Seven other people will attend.

Everyone will sit around the dining table.

The seating arrangement must meet the

following conditions:

* Amy and Alan sit together. Brad and

Beth sit together.

* Charles sits next to either Debbie or

Emily.

* Frances sits next to Debbie.

* Amy and Alan do not sit next to either

Brad or Beth.

* Brad does not sit next to Charles or

Frances.

* Debbie and Emily do not sit next to

each other.

* Alan does not sit next to either Debbie

or Emily.

* Amy does not sit next to Charles.

Who should sit where?