Hi,

I am trying to realize the following calculation in Maple.

$
  \left[\sum_{i=0}^n y_i(x) \partial_x^i , \sum_{j=0}^m z_j(x) \partial_x^j \right]  \\
=   \sum_{i=0}^n \sum_{j=0}^m \sum_{l=0}^i  \binom il y_i(x) \left( \partial_x^{i-l} z_j(x)\right) \partial_x^{l+j} \\
- \sum_{j=0}^m \sum_{i=0}^n \sum_{l=0}^j  \binom jl z_j(x) \left( \partial_x^{j-l} z_i(x)\right) \partial_x^{l+i} \ .

$

Is there a way to make maple understand d/dx as a differential opperator and calculate with it? When i for example try to calculate diff(d/dx, x) it should give me d^2/dx^2 as a result. Unfortunately i don't know how to realize this.

Basic problem is i don't know how to realize operator expressions in maple like for example:

f(x) d/dx      ( f(x) is a smooth function of x here )

where when applied to a function h(x) it should result in f(x) d/dx h(x) .

Is that possible?

Thank you very much in advance.

What I tried is at

http://apfp.sourceforge.net/maple_interval_arith.pdf

I would like to be able to do comparisons of intervals.

At least take min or max - I thought I saw that on help pages but didn't get it to work.

Would also like to use with complex numbers - tried evalrC - but did not do what I expected.

Hello,

please explain how to write a code to calculate and output the actual area using integration for y=X^3 over range (0,2) using left-hand rule and 200 subdivisions?

Thank you 

I currently have a function quadsum(n) that determines the [x,y] solutions of the above equation for an integer n. :

quadsum:= proc(n::nonnegint)
local
k:= 0, mylist:= table(),
x:= isqrt(iquo(n,2)), y:= x, x2:= x^2, y2:= y^2;
if 2*x2 <> n then x:= x+1; x2:= x2+2*x-1; y:= x; y2:= x2; end if;
while x2 <= n do
y:= isqrt(n-x2); y2:= y^2;
if x2+y2 = n then k:= k+1; mylist[k]:= [x,y] end if;
x:= x+1; x2:= x2+2*x-1;
end do;
convert(mylist, list)
end proc:

How would I alter this so that I get [x,y] for n= (5^a).(13^b).(17^c)(29^d) for non-negative integers a,b,c,d?

please is there any one can help me to find a solution of a sytem of 3 non linear equations each with 3 variable and with more than 30 unknown coefficients

this is the system

solve({EEE_x(x, y, z) = 0, EEE_y(x, y, z) = 0, EEE_z(x, y, z) = 0}, {x, y, z})

where x,y,r are the unknowns

and the three equations are simply the partial derivative with respect to x,y and z repectively

EEE_x(x,y,z):=(&DifferentialD;)/(&DifferentialD; x) EE(x,y,z)

EEE_y(x,y,z):=(&DifferentialD;)/(&DifferentialD; y) EE(x,y,z)

EEE_z(x,y,z):=(&DifferentialD;)/(&DifferentialD;z)EE(x,y,z)

the main equation is EE where (it has 3 variables and more than 30 qunknowns coefficients

(x, y, z) ->

1
----------------------------------------------------------------
2
/ 2 2 2\
hh \ii + jj x + ll z + mm y + 100. y + nn y z + oo x + pp z /

/ 2 2 2 2 3 2
\p z y + q z y + l z x + g z x + o z y + n z x + m y x

2 2 2 3 2 2 2
+ j y x + k y x + i z y + d z y + f z x + h z y

2 2 4 3 2 3
+ e y x + u z y x + v z y x + a + b x + c x + r x + s z

2 2 4 3 4 \
+ t z + bb z + cc y + dd y + ee y + ff y + gg z + aa x/

1.  a procedure quadsumstats whose input is an integer n. This procedure should return a list of length 

n whose kth  entry is the number of solutions to
x^2 + y^2 = k 
for
1 <= k and k <= n

I am sort of confused as to how to construct that list of length n and how to obtain integer solutions to the equation in maple.

2.

a procedure firstCount(k) that finds the first integer
n
with
k
representations as
"x^2+y^2= n." What does it mean for an integer to have k representations?

This is a question I have also submitted to the technical support, I am worried that it is a bit too technical for them, however:)

I am debugging a C program which links against the OpenMaple API library (under Linux and with Maple 17 and 18). I am using valgrind memcheck, because I am experiencing strange behavior which could be due to writes beyond allocated blocks of memory.  

The first thing which jumps to my eye, are many errors of the types

Use of uninitialised value of size (4/8/16)

Invalid read of size (4/8/16)

Conditional jump or move depends on uninitialised value(s) 

The same errors are also printed when I use the examples that ship with Maple. For instance, I compile "simple.c" with

gcc  -Wl,--no-as-needed -lmaple -lmaplec -lrt -L /usr/lib -L $MAPLEDIR/bin.X86_64_LINUX -I $MAPLEDIR/extern/include -o simple simple.c

and run valgrind as

valgrind --tool=memcheck --error-limit=no --log-file=memcheck.log ./simple

memcheck.txt 

Some, but not all, of the errors occur in __intel_sse2_strcpy or __intel_sse2_strlen. Furthermore, according to valgrind there are definite memory leaks. which appear in the library. 

Practically this makes it hard for me to identify my potential own errors. I am a bit surprised to see so many warnings because I tend to fix my own programs until memcheck does not print these anymore (before I give it away at least). The question is: Can I consider these errors as safe to ignore? How would I distinguish real errors which may appear in my application?

Hello all,

Suppose I have a vector valued function f (dimensions of vectors is 2). I can use fieldplot to show me the function's behaviour over some window.

I'd like to do the same thing, but I want to compact the space down to the unit circle. Basically, given a vector <x,y>, I'd like to "fieldplot" f(<x,y>) but with each arrow f(<x,y>) appearing at <x,y>/(1 + Norm(<x,y>,2)) instead of <x,y>.

Is this possible?

If I have the following system of first order diff eq's:

x'(t)=2x(t)+3y(t)

y(t)=-3x(t)-2y(t)

then can I consider the coefficient matrix A=<<2,-3>,<3,-2>> and compute the eigenvalues of A and infer as follows:

if the eigenvalues are of the same sign- eq point is a node

if they are of opposite signs- eq point is a saddle

if they are pure imaginary- eq point is a center

if they are complex conjugates- eq. point is a spiral

I've been given these conditions but my text says for a linear system of the form x'=Ax, the eigenvalues of A can be used to identify the nature of the eq. point. I am confused as to whether this applies to the given system as well; I have obtained 5 different trajectories and drawn the phase diagram for the system

I'm trying to solve a system of four pdes and I know that the Newton method won't converge.

Are there other numerical methods that I can use?

Any help would be greatly appreciated!

Thanks,

Eve

Hi,

I'm an industrial engineering student who's running into problems with solving a simple system of non-linear equations.

Why do I get this error when issuing the solve command? I'm pretty sure I'm only passing two lists into the function?

Thanks a lot in advance,

Joshua

Hi!

Is there any way to remove the empty space that comes under images when printing the project, while using document mode?

Regards

Nicolai

I have a matrix A for which the basis of the left null space using NullSpace() is the empty set {}  while the column space is {e1,e2,e3}. By definition, we need every vector in col space . every vector in basis of left null space =0 but how would I show that in this case? Can I determine another basis for the left null space?

Can anyone coax Maple to solve this reccurence relation? It seems harmless enough but Maple is struggng a bit with "hypergeomsols."

f := c -> (2*n-c)*f(c-1) - (c-1)*n*f(c-2);

f(0) := 1;

f(1) := 2*n-1;

I've got the following matrix :

A:=[<a,a-1,-b>|<a-1,a,-b>,<b,b,2a-1>] where <> are the column elements of A, a is  a real number defined on [0,1] and b^2=2a(1-a) 

a) to show A is an orthogonal matrix, I understand that I need A.Transpose(A)=Identity(3*3) but is there a way in which I can let a take a random real numbered value between 0 and 1? The rand() only returns an integer within a range. Directly multiplying A and Transpose(A) will return an expression in a, so what's the right approach?

b) from a) we can infer that A is a matrix that describes a rotation in e1,e2,e3 where these are the standard bases vectors in R3. How can I determine the rotation axis? The hint I've been given says I need to consider the Eigenvalues and eigen vectors but I don't quite understand how.