## a multibody dynamic problem(solving an 2nd order O...

several beams are connected by gemel one after another, the 1st is fixed at one end on the roof, and the last is driven by a force(fx,fy vector summation) at the end point. i am using lagrange method for the virtual displacement of angle theta[i] , and momentum balance in 2 direction (x,and y)to build the model. in the last step, i got an ODE system of 2nd order, but i cant code in maple to solve it. i think it is solveable, (6 equations, 6 unknowns). Kindly let me know how to write the code to solve it! below is the code.

> restart:
interface(warnlevel=0):
with( plots ): with( plottools ):
with( DEtools ): with( PDEtools ):
g:=9.81:
x[0][2]:=0:
y[0][2]:=0:
#initialization function of each beam,in this function, we should give 3 variables, n=the number of this #beam, la=the length of this beam, ma=the  mass of this beam
beams:=proc(n,la,ma)
global l,m,x,y:
l[n]:=la:           #length of  beam n
m[n]:=ma:           #mass of beam n
#initialize the positions of each beam
x[n][2]:=0:        #the x position of the end point of each beam
y[n][2]:=0:        #the x position of the end point of each beam
x[n][1]:=0:        #the y position of the center point of each beam
y[n][1]:=0:        #the y position of the center point of each beam
end proc;
#define the forces that applied on the end of  beam n
Force:=proc(n,fm,fn)
global fx,fy:
fx[n]:=fm:
fy[n]:=fn:
end proc:

## how to solve implicit ODE system like this?...

i have this im plicit ODE system in modeling a multibody dynamic system. theta[1](t), theta[2](t), are the angle displacement of the beams, and they are time depentent, in these 6 equations, i have the term of cos(theta[1](t)), sin(theta[2](t)),cos(theta[2](t)), sin(theta[1](t)), multipling (diff(theta[1](t),t))^2, and diff(theta[1](t),t\$2), which make them implicit equations. there are 6 unknown variables fx1,fx0,fy1,fy0(which are the forces at the end of the beams),and theta[1](t), theta[2](t). i gave the initial condition as ics:={theta[1](0) = 1/6*Pi, theta[2](0) = 1/3*Pi, D(theta[1])(0) = 0, D(theta[2])(0) = 0} Kindly let me know how to write the code to solve it! Thanks. Still below is the code: >

## Symbolic Computations under Max-Plus Algebra?...

Is there a way to do a Symbolic Computations under Max-Plus Algebra? For instance, linear algebra computations under under Max-Plus Algebra. I would like to appreciate any comments on it.

## Primitive Polynomials...

I need to search for primitive polynomials of degree n over finite fields GF(q), where q is a PRIME POWER. I have major difficulties doing that, especially when I wish to fix one of the coefficients of the polynomial. On the other hand, polynomials over GF(p) pose no problem. Does anybody have any experience in this? PLEASE HELP!

## Maple 10 update...

I noticed on the maplesoft website that there was an update available for Maple10. Is this in response to a detected product flaw ? Also is getting this update an urgent issue for a novice like myself still learning ? What I'm trying to ask is would I be missing out on some necessary feature of Maple10 by not having this update ? v/r

## Auto- and cross-correlation with Maple...

How do I perform auto- and cross-correlations in Maple?

## Sorting through multiple solutions...

My students are sometimes asked to use Maple commands that give multiple answers (solve and dsolve, most commonly) and then are asked to select a certain solution out of those returned. E.G. dsolve(diff(y(x),x) = 1/2*(6*x^3+2*sin(8*x)^2)/(y(x)^2)); returns 1 real solution and 2 complex solutions. The students are asked to select the real solution and assign it to a variable. I have making sure that Maple will select the correct solution out of the set (or list or sequence) for me by using complicated if statements. E.G. The first part of the if statement I use for the dsolve() command above is if evalf[50](subs({x=1,_C1=1},Re(rhs(solsde[1])))) =evalf[50](subs({x=1,_C1=1},rhs(solsde[1]))) then real_solution:=rhs(solsde[1]); fi; I leave off the two elif parts after the first if ... then for the sake of brevity.

## and...

how do i incorporate an and clause for my if statement? like if a=1 and b<>3 then

## how do i output text?...

i want to output text. if n=b then; n may not be completely factored end if; how do i do this?

## Plotting the Dirac Function...

Hi there everyone. I have recently been trying to plot fourier transforms. After taking the fourier transform of a particular function of the time domain for example a sinusoid. I get a sum of terms each with a translated Dirac function. When I try and plot the series of Dirac functions it doesnt work in the sense that there is nothing displayed. The plot initializes but there are no Dirac functions. Is there a way to plot Dirac Functions? thanks

## interpolation or approximation with two variables...

Dear Maple Users, I have the following problem: How to find the best-fitted function P(x,y) = ax+by+cxy+d to the following data: i= 0, 1, 2, 3, 4 , 5 , 6 , 7 xi= 150, 152, 154, 156, 156, 157, 158, 158 yi= 134, 136, 106, 61, 90, 124, 41, 135 I know that in Maple functions that can solve interpolation and approximation problems with one variable are implemented. But how to solve the problem with two variables? Best regards, Wieslaw Kotarski kotarski@gate.math.us.edu.pl

## loops...

if i have a procedure that has a bunch of loops, but when i execute it seems to be taking a long time. is there a way to figure out which loop the procedure is 'stuck' in?