I have a function that is a double summation and I cannot get maple to compute it. (I'm new to maple)

The function is as follows:

F(x) := sum from a=1 to x of [sum from b=0 to a-1 of [cos(2 pi x b / a)]]

for integer x, F should produce integers, but I am only getting 0's

naturally F is undefined for non integers. but maple is producing non 0 numbers for non integer entries. So, i have no idea what it is calculating.

Thanks in advance for any help. BC

I wrote the following code but am stuck on how to Update prcLeastSquareLA, so it does what prcLeastSquare does.

I believe you have to add error estimate information in prcLeastSquareLA but I have tried to no avail.

prcLeastSquareLA:= proc(data,degree)
local vars,y,A,V,k,e,i,j,v,c,1stVars;

vars:= seq(ci, i=0..degree);
y:=unapply(‘+’(seq(ci*t^j, i=0..degree)),t);
for k

How to find the minimum of (3*a/sqrt(1-u^2)+b/sqrt(1-t^2))/c under the constraints
{a^2+2*b*c*u >= b^2+c^2, a*t+b*u <= c, b^2*(t^2-u^2)/(t^2-1)+c^2 <= 2*b*c*u}, assuming all variables nonnegative, with Maple? Both  DirectSearch:-GlobalOptima and Optimization:-Minimize  fail in it.

I am writing to request some general opinion/support. I’m new to Maple and I’m interested in analysing a pile foundation (long steel pipe driven into the ground). The pile is normally modelled as a series of masses interconnected by straight elements and the soil surrounding the pile is modelled as a series of non-linear springs acting on each of the masses.

I would like to know if this is something that can be done in Maple or if MapleSim is needed.

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.