consider:

assume(kf1>0,kf2>kf1,q>kf1,q<kf2);

1/2*q*Pi^2*(ln(kf1)*kf1^2-ln(kf2)*kf2^2-ln(q)*kf2^2-arctanh(1/(kf1^2+q^2)^2*(kf1^4+q^4-6*q^2*kf1^2))*kf1^2+ln(q)*kf1^2-ln(2)*kf1^2+ln(-kf1+q)*kf1^2+ln(kf1^2-q^2)*q^2-ln(kf1-q)*q^2+2*kf1^2-2*q^2-1/2*ln((q-kf2)^2)*kf2^2+3*ln(kf2+q)*kf2^2+4*q*kf2-3*q^2*ln(kf2+q)+ln(2)*kf2^2+arctanh((kf2^4+q^4-6*kf2^2*q^2)/(kf2^2+q^2)^2)*kf2^2+ln(-q+kf2)*q^2-ln(kf2^2-q^2)*q^2-4*q*kf1-3*ln(q+kf1)*kf1^2+2*arctanh((-kf1^2+q^2...

Hi all! 
Is there any possibility in Maple to perform symbolic operations with arbitrary matrices (possibly given only the size of the matrix or it's structure (block matrices, diagonal matrix and so on)). For example, I would like to obtain the following results (input -> output). 

Possible code: 

assume(A, Matrix(n)): assume(B, Matrix): assume(E,IdentityMatrix(n)): 

I tried to use multiple light sources with transparency for 3D plots and could not get it to work.  There are indications in the plot,structure help that this is supposed to be possible.  I have attached a worksheet that demonstrates some of the problems with the "light" and "LIGHT" options for 3D plotting.  Much of the output has been removed since the graphics take so much space.  The same plots are generated in both Windows and Linux.

Does anyone...

In the LinearAlgebra packet, I use the command MatrixFunction to calculate a function with respect to a matrix J as follow:

> restart; with(LinearAlgebra);

> J := Matrix(3, 3, {(1, 1) = 2, (1, 2) = 1, (1, 3) = 0, (2, 1) = 0,
 (2, 2) = 2, (2, 3) = 1, (3, 1) = 0, (3, 2) = 0, (3, 3) = 2});

> M := MatrixFunction(J, ln(1+x), x);

Then I got the answer:

Matrix(3, 3, {(1, 1) = ln(3), (1, 2) = 1/3, (1, 3) = -1/18, (2, 1) = 0, (2, 2) = ln(3...

I would like to pay attention to http://www.scientificcomputing.com/atlas-032408.aspx . It seems to be a powerful tool to research, to teach, and to learn.

I am jsut a beginner!

I was trying to get the numerical value of e^10,

but Maple did not show me the approximate commond. How to get the numerical value of this ?

Hello,

(sorry if this appeared elsewhere, but it seems my first try to ask this failed)

I tried to use Maple last night to do some heavy computational work, but it turns out (I let him print time()) that he was only active for about an hour.  He did not do anything when my Macbook Pro was on stand-by.

This morning, when he woke up, he was unwilling to stop!  I pushed the "Interrupt the current operation" button...

On Monday, August 6 at 1:31 a.m. EDT, NASA will attempt the landing of a new planetary rover, named Curiosity, on the surface of Mars.  The Mars Science Laboratory project is managed by the NASA Jet Propulsion Laboratory (JPL) in Pasadena, California, a world-renowned center for robotic space exploration and advanced science and engineering.  JPL recently began a widespread adoption of Maplesoft technology, and Maplesoft’s products are expected to help JPL save...

With the addition of ten new Clickable-Calculus examples to the Teaching Concepts with Maple section of the Maplesoft website, we've now posted 63 of the 154 solved problems in my data-base of syntax-free calculations. Once again, these examples and associated videos illustrate point-and-click computations, but more important, they embody the

if my eqsn is x^4+x^3+x^2 . i can take a common factor x^2 ,then the eqsn will be x^2(x^2+x+1) ...also i can take a common factor of x ...then the eqsn will be x(x^3+x^2+x) ...how can i direct maple about my common factor ...u think this is useless but i have just made my problem simple to ask a qsn so that everyone can understand what i really want to tell

Hi !

I am in desperate need of help to solve a system of non-linear equations. I guess it cant be solved analytically, so maybe numerically? The system is

 EQ1 := u*alpha+u*s*A*Ӫ^k = alpha

EQ2: c*s^(c-1) + A*Ӫ^k*[(w-b+s^c)/(beta+alpha+s*A*Ӫ^k)]=0

EQ3 := Ӫ = [(1-u)/u]*[T/(1-T)]

I would like to solve for Ӫ, u and s. All other parameters are constant.

Help pleaseeee!

Apart from the online description of this new Maple 16 feature here, there is also the help-page for subexpressionmenu.

I don't know of a complete listing of its current functionality, but the key thing is that it acts in context. By that I mean that the choice of displayed actions depends on the kind of subexpression that one has selected with the mouse cursor.

Apart from arithmetic operations, rearrangements and some normalizations of equations, and plot previews, one of the more interesting pieces of functionality is the various trigonometric substitutions. Some of the formulaic trig substitutions provide functionality that has otherwise been previously (I think) needed in Maple.

In Maple 16 it is now much easier to do some trigonometric identity solving, step by step.

Here is an example executed in a worksheet. (This was produced by merely selecting subexpressions of the output at each step, and waiting briefly for the new Smart Popup menus to appear automatically. I did not right-click and use the traditional context-sensitive menus. I did not have to type in any of the red input lines below: the GUI inserts them as a convenience, for reproduction. This is not a screen-grab movie, however, and doesn't visbily show my mouse cursor selections. See the 2D Math version further below for an alternate look and feel.)

The very first step above could also be done as a pair of simpler sin(x+y) reductions involving sin(2*a+a) and sin(a+a), depending on what one allows onself to use. There's room for improvement to this whole approach, but it looks like progress.

Download trigident1.mw

In a Document, rather than using 1D Maple notation in a Worksheet as above, the actions get documented in the more usual way, similar to context-menus, with annotated arrows between lines.

Download trigident2.mw

I am not quite sure what is the best way to try and get some of the trig handling in a more programmatic way, ie. by using the "names" of the various transformational formulas. But some experts here may discover such by examination of the code. Ie,

eval(SubexpressionMenu);

showstat(SubexpressionMenu::TrigHandler);

The above can leads to noticing the following (undocumented) difference, for example,

> trigsubs(sin(2*a));
              
                                 1       2 tan(a)
[-sin(-2 a), 2 sin(a) cos(a), --------, -----------,
                              csc(2 a)            2
                                        1 + tan(a)

    -1/2 I (exp(2 I a) - exp(-2 I a)), 2 sin(a) cos(a), 2 sin(a) cos(a)]

> trigsubs(sin(2*a),annotate=true);

["odd function" = -sin(-2 a), "double angle" = 2 sin(a) cos(a),

                               1                       2 tan(a)
    "reciprocal function" = --------, "Weierstrass" = -----------,
                            csc(2 a)                            2
                                                      1 + tan(a)

    "Euler" = -1/2 I (exp(2 I a) - exp(-2 I a)),

    "angle reduction" = 2 sin(a) cos(a),

    "full angle reduction" = 2 sin(a) cos(a)]

And that could lead one to try constructions such as,

> map(rhs,indets(trigsubs(sin(a),annotate=true),
>                identical("double angle")=anything));

                             {2 sin(a/2) cos(a/2)}

Since the `annotate=true` option for `trigsubs` is not documented in Maple 16 there is more potential here for useful functionality.

Hi all,

Im trying to use the solve command to solve the eqn:= 40 = x arccosh(56/x) for x. The command solve(eqn, x) only provides the left solution(for low values of x) while I know there must be a solution as well for higher values of x. How can I get Maple to look at the other solution?

Thanks in advance!

Hi guys, I'm trying to solve the following function:

>with(inttrans):

>invlaplace (15/(s³ + 6s² + 15s + 15), s, t)

but I'm not getting a viable answer - can maple answer this?

I've also simplified the above function a bit to: H(s) = 15/[(s + 2.3222)(s² + 3.6778s + 6.4594)] and then tried to find it's inverse laplace but still don't get an answer.

please a need help

i try to organize the matrix by indice , i.e 

B := Vector[row]([1, 4, 13]);
 ind:=indices(B);

A := Matrix(3, 3, [1, 3, 4, 2, 4, 1, 2, 3, 2]);
#i try to do that

A:=B(ind,..); # and didn´t work  

#this is what i wont,

A := [1, 1], [1, 2], [1, 3], [3, 1], [3, 2], [3, 3], [2, 1], [2, 2], [2, 3]

help