Hi

As i asked at title; i want to determine natural frequency of a single or multi degree of freedom vibration system using MapleSim. but i can't and i don't know how i can do it. can every one help me plz?؟؟؟

Say I have a 3x3 matrix with real entries (like 3.4) and want to compute, say, A^(11). How can I do it and only have 5 digits in the output? I set Digits:=5; and used MatrixPower -- but got 10 digits, which is way too many for the classroom presentation.

Please see the attached worksheet for the error at the end, any idea what happens?

Besides, anyone know how to make plot array with each element has a plot and scatterplot in it?  I made my plot by making table first and then type the cmd individually.

Fitting_worksheet.mw

I also used a lot for for loop in this worksheet, can anyone suggest another way of doing it? And I think I initiated many Array(1..3) that is of no use.  But if I do not initiate and do the following

A[1]:=a : 

A[2]:=b:

A[3]:=c:

A;

output: A

Only if I initiate A:=Array(1..3) and then define the element I can see the complete A array at the end.

I am trying to write a procedure that have the following input:  expr, xdata,ydata. [dependent variables list]

I want a procedure to calculate Rsquare following the method described in 

http://www.mapleprimes.com/questions/120014-How-To-Calculate-Coefficient-Of-Determination

But I want the proc to take any abitary function and be able to identify the  dependent variable...how can I approach it? what function can I use to match [dependent variable list] in expr??

Thanks

Is it possible to set maple to fill in a header or footer for every new document I make?

Now i have to do it fot each new file by hand: 

View => Header footer ..=> 

 - select custom  - goto tab cusom headers put in the left textfield "&[File] pag. &[Page]/&[Pages]"

Hi,

I want to calculate a range.

eg. j=1...n-1 will use this equation

P(j,n):={((m/(2-1))(m/(2))^(j-1))/(2(m/(2))^(n)-1))} 

j=n will use other equation.

How to create the command? I keep receive errors.

Thank you.

Lina

what is the best maple book? could every one introduce books which have studied about maple !?? which is the best due to everyones opinion ? tnx for answering.

Hi,

I'm getting an error when I try to

plot(sin(x), x = -1..1, color = darkblue,tickmarks=[spacing(0.25),subticks=2]);

The error is
Error, (in plot) expecting option tickmarks to be of type list({"default", "piticks", "decimalticks", "pidefault", nonnegint, list(realcons), list(realcons = anything), spacing(realcons), spacing(realcons, realcons)}) but received [spacing(.25), subticks = 2]

What exactly am I doing wrong? After all, [...] is a list.

Even more confusingly, tickmarks=[[spacing(0.25),subticks=2],default] doesn't work either.

Any help appreciated!

Dear All,

I would like perform a symbolic integration to the following integrand with Heaviside function:

The integral I defined is as follows:

The output I got is

Could any one explain me how to avoid the undefined in the output. I tried to define the assumtions to constrain the solution. However, I failed to supress the undefined.

I expect a result of the integration should also be piecewise function.

Appreciate you constructive inputs.

Thanks!

restart:with(plots):

TSol:= f(x,y,Q); # The actual funtion is in the attached maple sheet because of its size.

animate( plot3d, [TSol,y=-1..1,x=0..1], Q=0..2,contours=65,style=patchnogrid, grid=[80,80], orientation=[-360,-180], lightmodel=light4, shading=zhue, transparency=0.3,style=contour);

plot3d(subs(Q=2,TSol),y=-1..1,x=0..1,contours=65,style=patchnogrid, grid=[80,80], orientation=[-360,-180], lightmodel=light4, shading=zhue, transparency=0.3,style=contour);

For Q = 2, why there is a big difference between these two ouputs?

 3d_anim_try.mw (If the maple sheet is not accessible then let me know, I will copy past the function (f(x,y,Q))?)

I have a question regarding maple reader.

It seems like in maple reader, one cannot interact with the output of the Explore function.... which mean I am forced to make my own component if I want to creat application for my co-worker to use..this is some what annoying.  Is there another way around?

please open the attached file in maple player and you will see the output created by Explore gives error message.  mapleplayer_test.mw

Download mapleplayer_test.mw

.  

Hello everyone, 

I have a problem solving with ODE's system solving. I have 2 equation and 4 initial  conditions. When i calculate like that,u can look this file.  beamsolvingfortim.mw. It is working it is giving me T1 and T2 equations depends on time. 

 In 2. system which i have a problem i want to calculate this equations depends on x(displacement). I have again 2 equation and 4 boundry conditions. it is solving the ODE'S system without boundry conditions. (It is giving with C1 C2 C3 and I)Problem is when i want to find its values with boundry conditions it is not giving a result. Is there a problem with complex number(I) or boundry conditions are not enough?   

> evalf(dsolve({sys}));

{X1s(x) = -3.060206320 + _C1 exp(-0.3487988669 x)

+ _C2 exp(0.3487988669 x) + _C3 exp(-0.3563227426 I x), X2s(x) =

-0.6321326989 _C1 exp(-0.3487988669 x)

+ 0.6321326989 _C2 exp(0.3487988669 x)

- 0.6053448484 I _C3 exp(-0.3563227426 I x)}

When i write like that, nothing is happening. I also upload the files. If u can help , i would be really appreciate. 

eq2 := dsolve({A, B, bc}, [X1s(x), X2s(x)])

beamsolving(thick.mw

Hello everybody,

I want to use the physics package in the context of special relativity.  I am mostly interrested with the Lorenzts transformations, not having to unprotect gamma all the time, the 4-vectors for space, time and momentum in (ct,x,y).  For exemple, to be able to calculate the invariant (ct)^2-(r^2).  A small document (or worksheet) would be very nice of you.

I understande that someone could say that I want to use a gun to kill a fly.  But this is a process that will lead me in using it in general relativity and using tensors in my calculations.  It would be very interresting to have a Student,Physics package.  Don't you think so?

Thank you in advance for your trouble and comprehension.

--------------------------------------
Mario Lemelin

Maple 18 Ubuntu 13.10 - 64 bits
Maple 18 Win 7 -  64 bits
messagerie : mario.lemelin@cgocable.ca
téléphone :  (819) 376-0987

It is well known that  operating in  geometry  package when setting objects symbolically some inconveniences arise. The error occurs if you do not set limits on the parameters. I wrote the procedure for finding center and radius of the circumscribed circle, devoid of drawbacks:

restart;

Circumcircle:=proc(a,b,c)

local x1, y1, x2, y2, x3, y3, n;

uses geometry;

x1,y1:=op(a); x2,y2:=op(b); x3,y3:=op(c);

n:=x2*y3-x3*y2+x3*y1-x1*y3+x1*y2-x2*y1;

point(A,a), point(B,b), point(C,c);

if type(n,realcons) then triangle(T,[A,B,C]); circumcircle(cc,T,'centername'=OO);

return [coordinates(OO),radius(cc)]  else assume(n<>0);

triangle(T,[A,B,C]); circumcircle(cc,T,'centername'=OO);

simplify([coordinates(OO),radius(cc)]); fi;

end proc:

At initial startup procedure works correctly:

Circumcircle([x1,y1], [x2,y2], [x3,y3]);

But if you change the arguments of the procedure, an error occurs:

Circumcircle([x1,1], [x2,y2], [x3,y3]);

But if you run the procedure again with the same arguments, the error disappears:

Circumcircle([x1,1], [x2,y2], [x3,y3]);

What is the reason?

Circumcircle.mws

Hi, the following Maple code produces an error for the value x=0: 0 instead of 1.

> for x from 0 to 5 do 

>print(evalf(sum((-1)^k*x^(2*k)/factorial(2*k), k = 0 .. 10)), evalf(cos(x))) 

> end do; 

0., 1. 
0.5403023059, 0.5403023059 
-0.4161468365, -0.4161468365 
-0.9899924966, -0.9899924966 
-0.6536436057, -0.6536436209 
0.2836642141, 0.2836621855 

Why such an error?

My solution to obtain the good value is:

> for x from 0 to 5 do 
> 
> print(evalf(subs(y=x,sum((-1)^k*y^(2*k)/factorial(2*k), k = 0 .. 10))), evalf(cos(x))) 
> 
> end do; 
                                1., 1. 


                      0.5403023059, 0.5403023059 


                     -0.4161468365, -0.4161468365 


                     -0.9899924966, -0.9899924966 


                     -0.6536436057, -0.6536436209 


                      0.2836642141, 0.2836621855 

Thanks.