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))];

I need to solve a quartic for on variable in terms of others. I am getting this error: 

Error, (in Engine:-Dispatch) should not happen: Rename expects the input to contain unknown functions

I am not sure how to fix this problem

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 ...