Hello all,

I have a square singular matrix E (12x12), and a matrix A (12x12). I need to find a matrix R (12x12) such that A = E. R. R^T .E^T 

How can I use Maple to find the matrix R?

Thank you so much,

STHence

Hi, I need to update a procedure prcLeastSquareLA so that it includes error estimate information as output in it. So far for the procedure I have this:

prcLeastSquareLA:=proc(data,degree)

local vars,y,A,V,k,e,i,j,v,c,lstVars;  

vars := seq(c_i, i = 0 .. degree);

y := unapply(`+`(seq((c_i)*t^i, i = 0 .. degree)), t);

 for k to nops(data) do 

e_k := y(data[k][1]) = data[k][2];

v_k := [coeffs(lhs(e_k))];

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.

We have just released a free collection of Maple and MapleSim classroom materials that helps bring modern technology to any introductory course in control design. This collection includes interactive classroom demonstrations that illustrate key concepts; lab projects and assignment questions; and example models ready to be explored.  It’s been designed primarily for instructors, but students should find much of it useful, too.

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?