Hi,
When I try to make to matrices equal (A:=B), whenever I make a change on B it changes my matrix A. This doesn't happen with normal variables. Is there a way to make to matrices equal (without changing every A[i,j] one by one) such that they don't interact.

Thanks ! 

I'm having an issue with the implicitplot3d command in Maple 16.  The sequence of commands I've listed below works just fine in Maple 15, but in Maple 16 I'm getting a problem.  Sometimes no graph will display;  other times, Maple may initially graph the surface correctly, but as soon as I rotate the graph, the surface disappears.  Furthermore, the following strange things happen:

1) Re-executing the command after rotating does not redraw the image ...

Hello!  Hoping that Maple 16 is now better with function compilation, I have compiled the rather simple function that returns the complex value.  As I understand from the code below, Maple automatically assumes the arguments to be complex.  That's fine.

The compiler gcc gives some errors:

/var/tmp/sergey-2181/_m71f6dec6b3c75a742676f2dcb3fa1864.c: In function ‘_c71f6dec6b3c75a742676f2dcb3fa1864’:

/var/tmp/sergey-2181/_m71f6dec6b3c75a742676f2dcb3fa1864...

I would like to plot a parametric equation

y = 7at -t^2

x= t^5-(3+a)t^3-3

and vary the value of (a) between -3 and 3

i tried to use:

plot([seq((7at -t^2, t^5-(3+a)t^3-3),a=-3..3)],t);

but I just can't get it to work

any ideas :S

Mark

My result is complicated, ant it was very large, when I simplified the question, it just as this program.

The main problem is that result in YY is not in order, so the converted result is also disorder, please help me to store the result in order. Thank you!

So, my code is this:

u := u1+(u2-u1)*JacobiSN((1/2)*phi*sqrt(s*(u3-u1)), k)^2;

r := 1/u 

where u1,u2, u3, s, and k, are just some constants

and here is where I think it's having problems 

Top := [r*cos(v), r*sin(v), 2*sqrt(s*(r-s))];
Bottom := [r*cos(v), r*sin(v)-2*sqrt(s*(r-s))]
figure(22) := plot3d([Top, Bottom], r = s .. 20, v = 0 .. 2*Pi)

It gives me this error: 

Error, (in plot3d...

I want to write a code to eliminate the noise in the data, but am unfamilar with how to do it.  This is the data:

bdata:=[[8,15],[0,18],[0,19],[-1,20],[-9,23],[0,19],[-2,20],[4,17],[-3,19,],[-1,20],[6,16],[6,17],[-2,20],[-8,23],[-2,19],[0,18],[-6,21],[-6,21],[7,16],[-9,22],[5,16],[-9,23],[-7,22],[-2,20],[4,17],[5,17],[-2,19],[-1,20],[-6,21],[5,16],[6,17],[2,17],[0,18],[5,16],[-2,20],[7,16],[-2,20],[1,18],[0,18],[-8,29],[-9,27],[-10,32],[-6,28],[7,11],[5,9],[8,7]]:

Hello all,

I have a singular square matrix E (12x12 with its rank of 10). I need to find 2 invertible matrices S and T such that S.E.T is a square diagonal matrix (in Smith Normal form). 

Using Maple with the following commands:
>Temp := SmithForm(E);
>Rank(Temp) 
Temp is a identity Matrix and Rank of (Temp)  is 12!!! or Rank of (S.E.T) = 12.

However, the determinant of E is less than 10^-29 and rank of E is 10.

Hello,

 I have a function that on two variables y1,y2. I know that the different variables, y2,y3 are orthogonal, i.e. y2*y3=0, hence I want to substitute this information into the function, but it does only work if only numbers apear in the expression.

#A minimal example with numbers as coefficients and parameters as coefficents

v(y2,y3):=y2^(2)+ y2*2*y3+y3^(2);  # numbers
u(y2,y3):=y2^(2)+ y2*h*y3+y3^(2);  # parameters

Hi there,

I am trying to solve a system that has 5 unkown variables. I have 4 equation and a if statment in order to solve the system as it can be seen in the figure.

The problem is that nch4(t) and nar(t) are the input variables. If i put the IF EXPRESSION before the equation i get this...

How to extract affine terms from a differential equation after linearization?

Vanderpol Eqn:

To be precise, I linearized...

My program is to compute the sequence of integrals of Jn = int(x^(2n)*sin(pi*x), x=0..1) for n=1...N to 10s.f where N=10.

This is my program but it keeps telling me it is missing an operator and I'm unsure why.

N:= 10:
In:= evalf(2/Pi):

for n from 0 by 1 to (N - 1) do

  In:=evalf((1/(Pi^2))*(Pi-((2n+1)*(2n+2)*In))):

   maple_In := evalf(Int((x^(2n))*sin(Pi*x), x = 0..1), 20):

   printf( "%d %10.10f %10.10f %10.10f \n" , n + 1 , In , maple_In , abs( In - maple_In ...

i use Infeld's The Factorization method and find this Hermit is type D. since not exactly match r i use R(x,m) = r(x,m) + f(m) and compare them, then find k = x but eigenfunction Y := C*exp(-int(x,x)); has constant C, range x should between negative infinity and infinity. C is 1/Pi^(1/4) after solved. i guess this is wrong, but i do not know where is wrong.

another problem is that i assume C = 1, after sum for generating function the solution is not the same as mentioned in book which is exp(2*x*z-z^2)...

Hi. Can anyone help me with getting all the prime numbers from 1 to 2000 in maple with "from 1 by 1 to 2000 do" and maybe use "if"? How do i make this?

I am trying to compare two matrices in a while loop. Here is part of my maple code with the comparison of Matrix B and Matrix E