Dear all

How can I finish the attached code and I plot   the set of points (x[i],y[j], U[i,j]) 

U[i,j] represent the value of U at the position  (x[i],y[j])

Many thanks for your help

Plotsetofpoint.mw

Before Maple 2015, units were enlosed in [[ double brackets ]] when displaying 2D input and output. In later versions these brackets are displayed only when editing input, and never in output, which in my opinion detracts from readability. Does anybody know if it is possible to revert to the old behaviour?

Here is my code:

with(DEtools): with(plots): with(linearAlgebra):

DE1 := diff(y(x), x) = d^3*y(x)/dx^3+3*d^2*y(x)/dx^2+4*d*y(x)/dx+12*y(x);

DEplot(DE1, [y(x), x = 0 .. 5][[y(0) = -3, (D(y))(0) = 0, ((D^(2)(y))(0) = 0, ((D^(3)(y))(0)]]);

I get the Error message: Error, (in DEtools/DEplot) called with too few arguments

Any help would be greatly appreciated, thanks.

Hi I have a function i d live to plot and integrate but maple tells me there is a probleme with the range when i want to plot and will not give me a numerical value of the integral :

What am i doing wrong ??

I am trying to repeat the 'first example' in Stephani & MacCullum, Differential equations see chapter 16 and (16.5).  The differential equation is,

pde := b1(x, y)*(diff(u(x, y), x))+b2(x, y)*(diff(u(x, y), y)) = 0

and I would like to compute DeterminingPDE.  The textbook answer is linear in b1 and b2 (as my 'by hand calculation' is) but Maple's,

DeterminingPDE(pde)

contains quadratic terms such as b1^2, b1*b2 etc.  I don't understand the appearance of the quadratic terms.  Is it possible for Maple to return an answer which is linear in b1 and b2?

In the link below to my worksheet I have attempted to construct the Fourier series for the function 1-x.  I should mention that the I often interchnge the variables x & t (probably a bad habit).  It appears that I am getting something close, but my term for n=0 appears to be undefined.  Given the plot I generated the term for n=0 should be approximately 1/4 to get the base of the triangle to intersect the t or x axis.

So I thought if I can get another set of eyes on my work maybe someone can point out what I have incorrect.

Appreciate any assistance

asymmetric_triangle_fourier_series.mw

I'm learning Maple with the provided user manual, and it clearly says to use ctrl+= to display inline results. It simply doesn't work on  my system. Running Maple 2017.3. Is this a bug?

deq := diff(f(z), z$2)-(2*z^2+z-a-1)/z*diff(f(z), z)-((2*a+4)*z+a+1)/(2*z)*f(z):

DETools:-formal_sol(subs(a = .1, deq), f(z), order = 3)[1];
                        1+.5000000000*z+.6845238095*z^2+O(z^3)

evalf[20](DETools:-formal_sol(subs(a = .1, deq), f(z), order = 3)[1]);
                        1.+.50000000000500000000*z+.68452380953750000000*z^2+O(z^3)

The second output is padded to 20 digits, but only the first 10 digits are correct, apparently taking the cached values from the first computation. I think that's really wrong.

Numerical evaluation of HeunBPrime fails if abs(z)>1:

HeunBPrime(.2, .3, .4, .5, 1.5);
Warning, breaking the computation of HeunBPrime after 20000 terms, the series is not converging
                      9.604689581*10^15896

while this gives the correct value:

subs(z = 1.5, convert(series(HeunBPrime(.2, .3, .4, .5, z), z = 0, 40), polynom));
                          15.37195056

How can solve system of equations by cramer's rule or LUDecomposition method in maple

I need to solve a system of differential-integral equations (see them on the picture). The T(x) function contains an integral with variable upper limit of the wanted function r(x), while the r(x) contains T(x). I don't know how to load the integral into a system, what syntaxix should I use. The boundary conditions and constants are written in the file, as is the error that appears when I try to enter the T(x) function. In the end I need a graphic of r(x) and T(x).

1_с_Y(theta).mw

1. This seems wrong:

modp1(ConvertIn((x^2+1)^5, x), 2);
Error, modp1: invalid arguments

ConvertIn works for other unexpanded inputs, and the documentation doesn't say that the polynomial should be expanded.

2. It's not clear what the correct way to test for equality is:

p := modp1(Monomial(1, x), 2):
q := modp1(Multiply(p, One(x)), 2):

is(p=q);
Error, (in property/ConvertRelation) invalid terms in sum: modp1(ConvertIn(x, x), 2)

evalb(p = q);
                              true

There is no Equal function, and IsZero(Subtract(p, q)) is rather clumsy.

3. In the documentation, with "Display Examples with 2D math input" checked, the second input looks like this:

a := (x^4-x^2+2) mod p

which doesn't make any sense.

4. In a 2D Input cell, using a label to refer to a previous output which is zppoly gives an error:

Error, '_Inert_ZPPOLY' is not a valid inert form

The same label in a 1D Input cell works.

5. It's not quite consistent that sometimes one-argument modp1(zppoly) works -- because zppoly stores the modulus -- but sometimes it doesn't:

modp1(IsZero(Subtract(p, q)));
Error, modp1: invalid arguments to function Subtract

6. There is no documentation for modp1/Embed.

 I want to learn about maple software. How to additte any text in maple work sheet? Kindly guide me in such a cases.

Regarded Nadeem Abbas (PhD scholar Mathematics PK)

My question is within the worksheet.

assignation.mw

I must be missing something in my Fourier integral.  My understanding is that the sinc function is the transform of a square wave.  In the link below I am getting something slightly different.  I have the parameter tau to define relative to the period, T, to vary the width aspect ratio of the wave,  If tau=T I do get the sinc function.

What am I missing or is what I have correct?

Sq_wave_Fourier_transform.mw