Does anyone know if there is a simple way to write a falling factorial (pochhammer symbol) in Maple. It seems pochhammer is the rising factorial with no option to change to falling.

Thanks in advance.

I assigned

before an algebraic calculation so I would like to get  or have the program print the 70 digits of the answer and not just 10 digits. Because when I press ENTER, I get only 10 digits.

fieldplot is a wonderful tool for plotting vector fields. The option 'fieldstrength' is very useful to scale the arrows so that one can better visualize the field. I often use fieldstrength=log.

However, if one includes the RealDomain library, the fieldstrength=log option fails. I don't see why it should. Can someone enlighten me before I report this as a bug.

Hi all, 
I was only wondering if there is a way how to trunk a solution using a maple command.

I want to use it to give a truncation error which is equal to :(14/45)(D^5)(y)(0) h^5

The result I am getting is as followed: (it is correct, I just want to cut off the bit with the power of 6)
........
......
>error:=expand(Yx[i+1]-Yx[i-3]-(4*h/3)*(2*f[i-2]-f[i-1]+2*f[i]));
                          
          error :=  (14/45)(D^5)(y)(0) h^5  + 7/10 (D^6)(y)(0) h^6
                  
thank you in advance

Hi, 

I'm being unable to plot the solution to the equation

schro := {-(diff(psi(x), x, x))+(2*a*b*x^4+a^2*x^6+(b^2-a*(2*p+3))*x^2-(2*p+1)*b)*psi(x) = 0};

for the special cases of a=b=1 and p=0 and p=1

I've used dsolve and am getting Heun functions :-(

The claim is that the solutions come out to be exponentials of the form:

psi(x)=(x^p)*exp(-(a*(x^4))/4 - (b*(x^2))/2)

thanks in advance

I am trying to implement Subresultant p.r.s. algorithm for calculating greatest common divisor. The algorithm decribed in the book:

My code return the correct GCD, however the sub-resultant terms are different from the result of the built-in function. The last term a[i-1] is huge and involves fractions. I think my implementation is same as the algorithm described in the textbook.

I have attached the file. Could anybody spot anything wrong in the code? Why do fractions still appear? In my code, "lsr" is last subresultant term returned from the built-in function, the second one is my result.

Download subresultant.mw

I am having difficulty to solve the following ODE

ode:=diff(R(x),x$6)=Dirac(x-y);

with initial and boundary conditions

bcs := ((D@@3)(R))(1) = 0, R(0)-((D@@5)(R))(0) = 0, ((D@@2)(R))(0)-((D@@3)(R))(0) = 0, ((D@@5)(R))(1) = 0, (D(R))(0)+((D@@4)(R))(0) = 0, ((D@@4)(R))(1) = 0;

I want to solve the above ODE for two cases 1) x>y and 2) x<=y

Thanks

Let us consider 

MmaTranslator:-FromMma(" Table[0,{n=10},{m=2}]");

[seq([seq(0,i=1..(m := 2))],j=1..(n := 10))]

The above result is syntactically incorrect in Maple language. The translation should be 

Matrix(10,2)
#or
[seq([seq(0,i=1.. 2)],j=1..10)]

up to the Mmma's result

Table[0, {n = 10}, {m = 2}]

{{0, 0}, {0, 0}, {0, 0}, {0, 0}, {0, 0}, {0, 0}, {0, 0}, {0, 0}, {0,0}, {0, 0}}

Another bug is as follows.

MmaTranslator:-FromMma("Sinc[x]");

Sinc(x)

whereas the expected result is piecewise(x=0,1,sin(x)/x) up to http://reference.wolfram.com/language/ref/Sinc.html

In general, the MmaTranslator package is outdated. It often returns working Mma's commands as incorrect (Concrete examples are long and need a context. These may be exposed on demand.). The question arises about the quality of other Maple translations. 

Hello,

I want plot this equations with values independents x,A,B

But this is error.

Any idea?

Thanks

I resolved the coefficients to a 2nd order diff eq of the form:ay''+by'+cy=f(t)

I have included the .mw file for convenience at the link at the bottom of the page.  I resolved the coefficients in 2 different ways & they do not concur.  The 1st approach used the LaPlace transform & partial fraction decomposition.  The coefficient results are given by equations # 14 & 15.  The 2nd approach used undetermined coefficients where I assumed the particular solution and then applied the initial conditions to resolve the coefficients pertaining to the homogeneous solution which are given in the results listed in equation #23.  Noted in the 1st case the coeff's are A₃ & A₄ and for the 2nd approach the coeff's are A₁ & A₂.  I have worked this numerous times & do not understand why they do not concur.  So I thought I should get some fresh eyes on the problem to find where I may have gone wrong.

Any new perspective will be greatly apprecieated.

I had trouble uploading the .mw file so I have included an alternative link to retrieve the file if the code contents is illegible or you cannot dowlad the file drectly from the weblink  Download coeffs_of_homogen_soln_discrepancy.mw.  You should be able to download from the alternative link below once you paste the link into your browser.  If you cannot & wish for me to provide the file in some other fashion respond with some specific instructions & I will attempt to get the file to you.

https://unl.box.com/s/dywe90wwpy0t4ilkuxshkivz2z26mud8

Thanks 4 any help you can provide.

Download coeffs_of_homogen_soln_discrepancy.mw

I need to implement gaussian elimination with cross multiply. I followed the method decribed in the textbook, however the result returned a wrong determinant. I noticed that cross multiplication changes the determinant in every step. Does anybody know how to use the cross multiplication without altering the determinant?

Find the first 6 non-zero terms of MacLaurin series of the function erf(x)

Consider the expression ∑_n=1 (-1)ⁿ e^*n/n²

1. Find the symbolic value of this sum.

2. Find an approximation for the value when  = -2.

3. Build a function to calculate an approximation for the value of the given expression for

any value for 

Hi all, first poster just getting to grips with Maple here. I am having some problems with the 'verify' function in Maple - I can't understand why it is returning 'FAIL' An example of my printout is given below. 

> verify( 2*q/(q^4+q^2+1)-1/((q-1)^2*(q+1)), 0, 'greater_equal') assuming q>7;

                                             true

> verify( (2*q+1)/(q^4+q^2+1)-1/((q-1)^2*(q+1)), 0, 'greater_equal') assuming q>7;

                                             FAIL

Now it's th second line i'm interested in and it was only after trial and error that i found the first option worked. The fact the 1st does but 2nd doesn't makes no sense to me. Perhaps it's something stupid i'm missing - any help would be great. Thanks.

Dear all,

I'm investigating the vibration performance of timber beams. I have sample data from my test which shows the vibration of the beam. I want to determine the eigenfrequency from this data. The problem I face is that I'm not finding the probber eigenfrequency. I have two data rows; time and amplitude. I'm able to plot the amplitude with SignalPlot but not the time, therefore I have to adjust the samplerate. I have the same problem with the fourier analysis. Is it possible to include the time period as well?

Regards,

Maurits