It took me a while to figure out what was causing this unusual sort result;

I loaded the MTM package because I want to convert numbers to integers for indexing. This sort result seem like it could be a bug.

Download MTMsort.mw

f=2*x*(int(sqrt(1+y^2*(t*x-1)^2/(1-(t*x-1)^2)), t = 0 .. 1))    ,0<x<2,   0<y<1,

I want to Determine the value of f for the various values of x,y For a domain So that x , y Changing with steps 0.1. for example:

Hello Everybody,

I was trying to apply the Newton-Raphon method("Newton" in maple) for the following problem to obtain the constants (c1...cM) without success.

eq[1] := -0.0139687 c[2] - 0.0132951 c[1] = 0
eq[2] := 24806.4 c[2] - 0.0139687 c[1] = 0

If anybody has experience with it or knows how to use it I would really appreciate your help.

Just to introduce the previous steps of calculations leading to the problem:To calculate the critical buckling force N and the shape of a rectangular uniformly loaded plate the governing diff. equation is the following

D11*w''''(x)+(2*D12+4*D66)*w''(x)''(y)+D22*w''''(y)+N*w''(x)=0

For the follwing solution of a boundary value problem(Boundary conditions:clamped-clamped: w(x=0,x=a)=0 & w'(x=0,x=a)=0) I applied the Ritz method:

w:=sum(c[i]*(cos(2*i*Pi*x/a)-1)),i=1..M)

Thus the Potential Energy P is:

P:=16537.6*c[2]^2-0.0139687*c[2]*c[1]+1033,60*c[1]^2-9.86960*N*c[2]^2-2.46740*N*c[1]^2

 deriving P to each constant c and setting them =0 leads to:

eq[1] := -0.0139687 c[2] + 2067.20 c[1] - 4.93480 N c[1] = 0
eq[2] := 33075.2 c[2] - 0.0139687 c[1] - 19.7392 N c[2] = 0

After calculating N=418,902  and feeding eq1 and eq 2,the follwing equations if two terms are considered:

eq[1] := -0.0139687 c[2] - 0.0132951 c[1] = 0
eq[2] := 24806.4 c[2] - 0.0139687 c[1] = 0

Everything I tried resulted in any c1 = c2 = 0 which is not  realistic. Maybe I made a mistake earlier.

Thanks a lot in advance.

Sam

 Plate_Buckling.mw

Hello,

Does anyone know how to find the repeated part of a fraction. I know the command cfrac, but I have

heard there should be another command which give the repeated part of a fraction as a result.

Who can help me?

 3/11+3*sqrt(3) = [5, 2, 7, 1, 1, 6, 1, 627, 1, 6, 1, 1, 7, 2, 9, 1, 12, 16, 1, 10, 1, 4, 3, 1, 1, 2, 1, 14, 2, 3, 10, 9, 1, 1, 139, 5, 5, 3, 2, 3, 6, 1, 1, 1, 9, 1, 2, 1, 6, 1, 2, `...`]

Hi guys,

I can not convert a series with label.

If anyone can confirm this?

Maybe-Bug.mw

Another_way_to_calculate_the_definite_integral.mw

p_tin the file I upload is a complicated expression containing θ and r. Will the accuracy of the result be affected if I use the second method in the file to calculate the definite integral?

Hi guys,

I can not execute any command in worksheet. I even can't set a variable to a number. When I write some lines and hit enter, it does not respond in anyway. Basicly I can't enter anything it doesn't do anything. How cold I solve this?

Dear Maple users

I know how to fit a function with some parameters to some data, but how can it be done if the data is 2-dimensional? I mean: I have some time array T, some X array and some Y array. How do I fit a function with certain parameters to the data: (x(t),y(t)) to fit (X(T), y(T)), ...

Erik

Hi

I've this equation ,what is displacement,

U := sum(C[i]*lambda[i]*h^2*M*(C[13]+C[55])*exp(Z*lambda[i]*M)*cos(M*x)/C[11], i = 1 .. 4);
 But when I used the data ,I do not get the right result

the datas are

h := .25;
M := Pi/L; L := 1;
E[3] := 1; E[1] := 40; E[2] := 1;
1
V[13] := .25; V[31] := .25; V[12] := .25; V[23] := .25; V[32] := .25; V[21] := .25;

G[44] := .2; G[55] := .2; G[66] := .2;

G[13] := .5;

C[55] := G[13];

Omega := -2*V[13]*V[21]*V[32]-V[12]*V[21]-V[13]*V[31]-V[23]*V[32]+1;

C[11] := (-V[23]*V[32]+1)*E[1]/Omega;

C[12] := (V[13]*V[32]+V[12])*E[2]/Omega;


C[33] := (1-V[12]*V[12])*E[3]/Omega;
 

C[1] := -.2596702449; C[2] := -2.411545729; C[3] := -0.6051254981e-1; C[4] := .4085627399;
 

lambda[1] := .6891428810; lambda[2] := -.6891428810; lambda[3] := 9.177422410; lambda[4] := -9.177422410;
 

could you help me plase

I'm trying to plot a piecewise function whose return value has units, but I get an error saying that it found a "name" (i.e., the unit symbol).

Example:

> f := t -> 2*t;

> fp := piecewise(f(tt) <= 4, 1*Unit(W), f(tt) > 4 and f(tt) <= 8, 2*Unit(W), f(tt) > 8, 3*Unit(W));

> plot(fp, tt = 0..10);

I've tried using the option useunits, but no success. Any help will be appreciated.

Limitation of Maple spreadsheet is only 100 x 52

I don't know why there's a limitation on the size but there's no way of increasing it either.  Is there?

Just to be sure that some outside operation wouldn't nudge the spreadsheet size I tried

with(Spread):
a:=CreateSpreadsheet():
b:=Matrix(100,100):

SetMatrix(a,b)
                           Error, (in Spread:-SetMatrix) spreadsheets cannot exceed 52 columns

So unfortunately until the feature gets built into Maple I don't believe there is a workaround.

Download pseudo.mw

I read research paper about elasticity, and I don't not how can get the results for displacment ant shaer stress,, could you help me please.

please help me

Kind Regards

I have a matrix A and a matrix B defined as B:=A. When an entry of B is changed in my Maple code, the matrix A is also changed. 

restart:
A := Matrix([1, 1]):
B := A:
B[1, 1] := 2:
print(A):

The result is

[2,1] but as you see, I did not change the entries of A during the code.  

This never happen for numbers. 

Please help me.

http://tinypic.com/r/5a048l/9

I have a multivariate polynomial in x and y. How can I firstly collect the terms with respect to the powers of x and y and then simplifying their coefficients.

Download A.mw