mech_problem.mw

Hi

I solve a mechanical exercise but i had a problem.

I know M (mass) and S (stifness) matrices (6x6).

I want to solve the (λ²M+S)v=0  eigenvalue problem, where λ are the eigenvalues and

Hi, this  must have been asked numerous time, but I cannot find it in MaplePrimes:

How does one "collect" sqrt terms? collect(expression,sqrt) does not work because sqrt is actually represented as ^(1/2) in  Maple. collect on '^' does not work either... ??

Mac Dude

Hey guys..

I need some help on determining the correlation coefficient between my data points for an Influenza epidemic and the SIR model, with my estimated parameters Beta and Alpha. The initial conditions for the equations are defined at the beginning of my worksheet:

SIR_model.mw

Hope someone can help me ;)

How to successfully deploy Matlab function [B, IX] = sort (A, ...)
IX is a permutation vector of the corresponding column of A

M := Matrix(6, 3, [2, 3, 3, 5, 7, 8, 12, 5, 9, -3, 4.1, 7, 7, 7, -3, 9, 3, 8]):
M;
L := <(seq(op(convert(sort(M[() .. (), j]), list)), j = 1 .. op(M)[2]))>;
L := ArrayTools[Alias](L, [op(M)[1 .. 2]], Fortran_order);
P := Matrix(op(M)[1 .. 2],  (i, j)-> ListTools[Search](L[i, j], M[() .. (), j]));
# from Matlab

I had MAPLE 12 on WINDOWS VISTA machine and the harddrive crashed.  Now I have a newer WINDOWS 7 machine & would lke to put my individually licensed MAPLE 12 on this machine.  I do not have a DVD because the original installation occurred via MAPLESOFT website if I recall correctly.

Can anyone direct me as to what I need to do?

THANKS

Hi,

I'm having some problems with the generation of the typical stiffness matrices seen in FEA formulations. I have got an elementary matrix that is looped through the entire stiffness matrix to generate it. As there is juxtaposition between the elementary matrices, there is some recursive assignment involved that I can't handle properly. Here is the code for a 6x6 elementary matrix with a 3x3 juxtaposition. n is the nombre of recursions needed to fill the stiffness matrix. M is the elementary matrix

On my keyboard, Cntrl = action for evaluating an expression inline does not work. The reason is that = is the upper level action (with the shift key) on a key with 0 at the lower level. Doing Cntrl Shift = does not work either. How do I get this done on my keyboard?

What is the number of all the solutions of the equation frac(x*floor(x)) = 1/2 belonging to RealRange(1,100)? How to count it with Maple?

Hi,

is the series expansion for the following expression bugged in maple or am I missing some crucial thing of series?

series((1/epsilon-log(4)),epsilon=0);

results in epsilon^-1 -2 -2*ln(2)

for xmaple 16.01. The console version does not add the wrong term "-2".

Apparently it works here: 

How is it possible to let multiple cores compute a loop faster?

For exaple:

for i = 1 to 1000 do

*computation(i)*

end do

The goal is to let a quad core cpu split the 1000 independant computations in 4 packages of 250 each to speed up the computation by the factor 4 so that core1 computes i from 1 to 250, core2 from 251 to 500 and so on, at the same time. Is this possible?

A triangle is given.The points

Hi, i have the following question: i've been asked to investigate the following family of functions:

f(x,y)=(be^a(x)^2)+y^2

the question is: How does the shape of each curve depend of a and b.

Thank you

Hello every one,

I want to plot two parameters against each other from a system of ode.

The system is

restart:

with(plots):with(LinearAlgebra):

Eq1:=diff(f(eta),eta$3)+f(eta)*diff(f(eta),eta$2)-diff(f(eta),eta$1)^2-M^2*diff(f(eta),eta$1)+M^2*b+b^2=0;

Eq2:=diff(g(eta),eta$3)+f(eta)*diff(g(eta),eta$2)+(Lambda-2*diff(f(eta),eta$1)-M^2)*diff(g(eta),eta$1)+diff(f(eta),eta$2)*g(eta)=0;

N:=10: M:=0.1:

bc:=f(0)=0,D(f)(0...

Dear friends,

I ran into another problem while using Maple to do residue calculus. The following call

residue(1/x^(1/3)/(x^2+2*x*cos(phi)+1), x=-exp(I*phi));

returns zero -- it does not recognize the residue. On the other hand, if I do a subsitution like this

subs(x=-exp(I*phi), 1/x^(1/3)/(x^2+2*x*cos(phi)+1));

followed by

convert(%, trig); simplify(%);

then I get a divide by zero error, which shows that the...

Hello all,

I recently helped someone with a contour integration problem, you may find my calculation at the following link. It occurred to me to test this integral with Maple and Mathematica.

For example,

int(1/x^(1/3)/(x^2+2*x*cos(Pi/7)+1), x=0..infinity);

produces an unusable output from Mathematica. Maple is a little better,...