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: 

i got two density, the part of it is similar the pattern of gamma, how to convert it into an expression of gamma times something?

test2 := x*Diff(P,x$2)+(a+1-x)*Diff(P,x)+mu*P=0;
Density := int(exp(-X*(I*u-a*i))*(-i*X)^(-a), X = -infinity .. infinity);

test2 := x*Diff(P,x$2)+(a+b-x)*Diff(P,x)+mu*P=0;
Density := int(-exp(-X*(I*u-a*i-b*i+i))*(-i*X)^(-a-b)*i*X, X = -infinity .. infinity);

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

Hi, I wrote a set of codes that worked perfectly fine 2 months ago, but when I try to run all my files again, this statement came out :

"Error, Got internal error in Typesetting:-Parse:-Postprocess:"intercal error:invalid object""

Could someone please explain what is the meaning of this please? Thanks

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...

hi everyone,

I need to plot a graph for system of nonlinear equations with multiple variables. 

here is some example of system nonlinear

f1=3.65a^2-2ac-7ad+3.65b^2+7bc-2bd-1.65a+6.22b+0.5=0

f2=13.22a^2-7ac+2ad+13.22b^2-2bc-7bd-6.22a-1.65b+0.24=0

f3=2c^3-2ac-7bc+2d^2+7ad-2bd-0.2=0

f4=c^2+d^2-1=0

with initial condition (0.7, 0.2, 0.6, 0.3)^t. I solved this using Newton's method and here the iteration

I have a procedure which accepts arguments of a (created by myself) type which are records, like so:

Subs:=proc(eqn::seq(equation),elemt::Element,$);
...

Here "Element" is  type I have defined; it is a record with certain fields present.

My question: will this work with option remember? I have the procedure sort-of working without the remember option, but I need to make it work recursively (the Element record can have fields which contain names of other records of type Element...

For the following procedure how do I dynamically grow the array? Is an Array the right datatype for this task?

How to find the maximum of sqrt((x-1)*(y-x))+sqrt((7-y)*(1-x))+sqrt((x-y)*(y-7)) for x=-2..3,
y=0..11 with Maple?
Here are  a few unsuccessful attempts.
> Optimization:-Maximize(sqrt((x-1)*(y-x))+sqrt((7-y)*(1-x))+sqrt((x-y)*(y-7)),
 x = -2 .. 3, y = 0 .. 11);
Error, (in Optimization:-NLPSolve) complex value encountered
> maximize(sqrt((x-1)*(y-x))+sqrt((7-y)*(1-x))+sqrt((x-y)*(y-7)), x = -2 .. 3, y = 0 .. 11);