I'm trying to print a nice looking array (6 columns with various numbers of rows) from an array that storing the data with 5 indices.  Here's the code:

P[i,j,p,q,k]:=(W[1,i].W[1,j].W[1,p].W[1,q]=W[1,k]);

Hello Maple experts,

I am preparing a little package for Maple for which I would like to set the Digits variable to the maximum possible value with hardware floating point evaluation. I know that this is done with

Digits := trunc(evalhf(Digits));

However, as advertised in the help pages, a change of the Digits variable won't persist when done inside a procedure. When preparing my package I can specify some 'setup' commands that are executed when loading the package but changing the Digits setting seems to be impossible...

Hi,

I wondered if you would be able to help.  I've set up a for-loop that substitutes values of b into a given function (dependent only on b).

When substituting the values of b into the function f the values of f decrease to a point and then the values of f begin to increase.  I want to establish the interval in which this change occurs, i.e. the minimum of the function.

Below I have included a code similar to our situation:

Hello, I would like to perform some formal computations with non-commutative

variables in maple (Lie-bracket and further) in order to define Lie algebras structures;

any idea about how I could "detype" the commutativity assumed upon the variables?

Thank you for your help!

from an fsolve where I solve 2 variables (fout and u2)  I get get this anwer.

ddmax := {fout = 0.3873538368e-1, u2 = 9.708916648}

However they are not allways in the same order, somethimes u2 is listed first. So how can I determine which value is listed first.

Thank you for any solutions to this.

M.

A rather fundamental I expect but I have tried reading the help file and cannot understand it - apologies if that makes me rather dim.

I have some data points (not a function) which I would like to plot.

If I simply use plot then the points are joined by straight lines - I'd like it to be a curve. Having looked at this forum I found a post which pointed to the curvefitting help page but I've looked there, tried a few things but failed to plot anything.

I have these three functions:

Phi := proc (z) options operator, arrow; z^(-(1/2)*alpha)*(BesselJ(alpha, 2*sqrt(lambda*alpha*z))-BesselY(alpha, 2*sqrt(lambda*alpha*z))*BesselJ(alpha, 2*sqrt(lambda*alpha*a))/BesselY(alpha, 2*sqrt(lambda*alpha*a))) end proc

w := proc (z) options operator, arrow; alpha*z^alpha end proc

Hi

I founda something strange here:

diff(y(x),x$2)+4*y(x)=cos(2*x);

dsolve(%,y(x));

y(x) = sin(2*x)*_C2+cos(2*x)*_C1+1/8*cos(2*x)+1/4*sin(2*x)*x

that's the answer maple give me when i run it for the first time

take a look at the two cos(2*x) tern?

shouldnt that be just ONE constant times it?

ok then if you run it again

somtimes it can collect them together give you single cos(2*x) term with ONE constant

EI {dy(x,t)/dx } +CA{d y(x,t)/dx}= summation of (F_o (sin(5nvt/r)sin(nwt)))  from n=1 to n= infinity

the derivative here are partial not ordinary

where EI, CA ,F_o , v ,r and w are constants 

Hi guys,

What is the exterior product of

       1  0  3  0

 A=(                 )

      -2  3  1  4

       1  4  -2  2

 B=(                 )

      3  2  0  4

written as A^B?

Hi there,

does anyone know how to write standard multicolumn numerical data sets generated by a "for" loop into an array of data files with names "dataset%.txt", where % is iteratively given by the value of the loop index? Thanks a lot. 

hi,

how do i solve the following?

restart:

ode:=diff(y(t),t$2)+y(t)=piecewise(t in [0,Pi],t, t>Pi,Pi);

ics := y(0)=0, D(y)(0)=0;

dsolve({ode,ics});

Thanks

Hi everybody,

I'm not really sure of where sould I post my question.

Let me explain my problem :

I've got a convolution product of my function f and my filter called G (let say a gaussian for example, or a top hat filter). It is as follows :

filtre:=proc(g);
return x->(1/(2*Delta)*int(g(t),t=(x-Delta)..(x+Delta)));
end;

filtre:=proc(g) local a,G;
G:=x->sqrt(gamma/(Pi*Delta^2))*exp(-gamma*x^2/(Delta^2));
return x->int(g(t)*G(x-t),t=-infinity..infinity);
end;

Hi, i have a dgl i want to transform:

y''-7y'+10y = (sin(2x))/x

The laplace command works, but how can i insert initial conditions ? The help page says:

Hi,

Is there a command to extract the number of (significant) digits for each of the elements in a list?
Thanks for helping me.

Hans van Walen