## Short-code way to obtain the 5 lowest orders in pa...

Asked by:

Hi all,

I have an expression of the form .

Is there a shorter way to obtain let's say the 5 lowest orders than  ?

Cheers

NOh

## NonlinearFit error estimation...

Asked by:

Hi everyone!

I wrote simple program which fits experimental result. I know that experimental points should be fitted with this function — Te=exp(a*X+b), where "a" and "b" are parameters. "X" — vector of independent variable, "Te" — vector of dependent variable.

My question is: can maple estimate error of my fitting? I want to obtain "a±Δa" and "b±Δb".

My program: Profile.mw

P.S. Excuse me, for my english. It isn't my native language.

## after updating to OS X Yosemite the License data a...

Asked by:

Der Group,

after updating to OS X Yosemite the License data of Maple 18 are invalid. What can I do?

Yours

Berthold

## find area between curves...

Asked by:

Hi Maple friends.

If I have y=-x^2+9 and y=x+3, how can I find the area that is bounded by the functions?

Easy enough to calculate by hand, but I need to use Maple to check my answers. Most preferably I'd like to use the Maple context menu to get the result, but any solution will be appreciated.

Thanks in advance.

## Insert matrices into expressions of quantum operat...

Asked by:

Hi again,

I am lost... again. I worked with the Physics package and have obtained my result in form of a lengthy equation of commutators of non-commuting objects. So basically my problem is solved, except that I now want to compute some examples to see, if everything is correct. For that purpose I want to replace now the non-commuting objects with explicit matrices, compute the result and see if it agrees with numerical studies. However, I have no clue which command to use, because if I try it with subs, then the multiplication is of the wrong type and Maple won't multiply my matrices. I tried to use eval but wasn't successful.

Here is an exampe to see what I mean.

Test2.mw

Thank you for your time

NOh

## trivial Commutator expression is not simplified in...

Asked by:

The following MWE shows what I mean:

with(Physics):Setup(mathematical=true):

Setup(noncommutativeprefix={MX,MY,MZ});

test:=proc()

local eq;

eq:=-Commutator(MX,MY)-Commutator(MZ,MY);

eq:=simplify(subs(MX=-MZ,eq));

return eq;

end proc:

test();  # yields -[-MZ,MY] - [MZ,MY]

%  # yields 0

Any ideas how I can solve this? I would like to return the simplified version.

## function call within sum, summation index is not e...

Asked by:

Hi guys

I have an if statement inside a function. This function then get's called in a sum. However, the result is not correct as the if statement checks wether the argument equals 0, but since the sum passes only the summation index name, it doesn't work as intended.

How can I tell Maple, that it has to evaluate the index before passing it on? MWE appended.

Cheers

NOh

 (1)

 (2)

 (3)

## Updating Maple Integration Database...

Asked by:

I have an integral that maple can not solve but I can solve it by hand. How can I add this to maple integration database?

f:=int(r^2*BesselJ(0,a*r)*BesselI(1,b*r),r)

Please see the file below.

Integral.mw

Many many thanks! :)

## Substitute a function in a matrix...

Asked by:

Hi

How can I substitute a function in a matrix? This is what I tried but it's not working at all:

lign1:=x-y+2*z=1:
lign2:=(-2)*x+y+z=0:
lign3:=(-4)*x+y+7*z=2:
lign4:=3*x-2*y+z=1:
T := GenerateMatrix([lign1, lign2, lign3, lign4],[x, y, z], augmented);

Then I'd like to substitute another function with lign1:

lign5:=x+y+z=6:
T2:=subs(lign1=lign5, T);

But T2 = T and that wasn't really the point.

Sorry if the answer is obvious no amount of searching has helped me so far. I hope this question isn't too ridiculous.

Thanks

## Simplifiying exponential expression (susb,algsubs,...

Asked by:

kappa := Vector(7, [1,w[1]*(1-phi+phi*(1-1/(1+exp(-mu[p]-tau[p3]))))+(1-w[1])*
(1-phi+phi*(1-1/(1+exp(-mu[p]-tau[p3]-eta[p2])))),w[1]*phi/(1+exp(-mu[p]-tau[
p3]))+(1-w[1])*phi/(1+exp(-mu[p]-tau[p3]-eta[p2])),w[1]*(1-phi+phi*(1-1/(1+exp
(-mu[p])))*(1-phi)+phi^2*(1-1/(1+exp(-mu[p])))*(1-1/(1+exp(-mu[p]-tau[p3]))))+
(1-w[1])*(1-phi+phi*(1-1/(1+exp(-mu[p]-eta[p2])))*(1-phi)+phi^2*(1-1/(1+exp(-
mu[p]-eta[p2])))*(1-1/(1+exp(-mu[p]-tau[p3]-eta[p2])))),w[1]*phi^2*(1-1/(1+exp
(-mu[p])))/(1+exp(-mu[p]-tau[p3]))+(1-w[1])*phi^2*(1-1/(1+exp(-mu[p]-eta[p2]))
)/(1+exp(-mu[p]-tau[p3]-eta[p2])),w[1]*(phi/(1+exp(-mu[p]))*(1-phi)+phi^2/(1+
exp(-mu[p]))*(1-1/(1+exp(-mu[p]-tau[p3]))))+(1-w[1])*(phi/(1+exp(-mu[p]-eta[p2
]))*(1-phi)+phi^2/(1+exp(-mu[p]-eta[p2]))*(1-1/(1+exp(-mu[p]-tau[p3]-eta[p2]))
)),w[1]*phi^2/(1+exp(-mu[p]))/(1+exp(-mu[p]-tau[p3]))+(1-w[1])*phi^2/(1+exp(-
mu[p]-eta[p2]))/(1+exp(-mu[p]-tau[p3]-eta[p2]))]);

Download kappa.txt

Here is the expression, I am trying to simplify, given a set of rules. NEW_Cole.mw

I have tried different substitutions, using simplify with side rules, applyrule, eval, subs, algsubs.

But none seem to be working as the way I want them to be.

Is there a better way?

Thanks!

## Solutions of Differential Equations with Embedded...

Maple 18

The Embedded Components are containers that currently use industries for modeling complex systems to find viable solutions in real time and thus avoid huge wait times and overload our computer; by this paper should show you how to implement a dynamic worksheet through Embedded Components in Maple; it goes from finding solutions to ordinary differential equations partial; which interact with the researcher using different parameters.
Using graphical programming will find immediate solutions to selected problems in science and engineering criteria of variability and boundary conditions evolving development with buttons on multiple actions.

cimac_2014.pdf

(in spanish)

Solutions_of_Differential_Equations_with_Embedded_Components.mw

Lenin Araujo Castillo

Physics Pure

Computer Science

## Quantum perturbation theory: energy denominators a...

Asked by:

Hi!

I seem to run into problems with (quantum) perturbation theory. In the following minimal working example, an energy denominator, as occuring in perturbation theory is not evaluated if I previously assume that the quantum number is a positive integer. It's supposed to return an error, as the energy denominator is 0.

 >
 >
 >
 (1)
 >
 >
 (2)
 >
 >
 (3)
 >
 (4)
 >

Download energy_denominator.mw

Best regards,

Sören

## spacecurve not working on Maple 17/18...

Asked by:

Hi,

I was helping for a computing session today.

Maple 17 and 18 were not plotting the following:

restart:

with(plots):
spacecurve([x,cos(x), sin(x)], x = 0..10);

On some of the machines, with campus site license, windows 7, x64.

The plot came up with a "frame" only. And nothing was plotted.

We really can't figure out why?

Something to do with the settings(options)?

Had anyone had a similar experience?

Thanks,

casper

## Simple code for MonteCarlo approximation of Pi...

Asked by:

Hello everybody,

I find a answer on that subject at this link: http://www.mapleprimes.com/maplesoftblog/95508-Monte-Carlo-Pi

But for starting easy, I try this:

with(RandomTools[MersenneTwister]);
directpi := proc (n)
local i, x, y, N;
N := 0;
Digits := 2;
for i to n do
x := GenerateFloat();
y := GenerateFloat();
if x^2+y^2 < 1 then
N := N+1:
end if:
end do:
N;
end proc;

evalf(directpi(4000)/(1000.0), 4);

3.176

I would like to know if this code is appropriate and if it is possible to have a random number between -1 to 1 so it will be easier to plot the points in the square with the circle inside.

Thank you.

```--------------------------------------
Mario Lemelin```
```Maple 18 Ubuntu 13.10 - 64 bitsMaple 18 Win 7 -  64 bits
messagerie : mario.lemelin@cgocable.ca
téléphone :  (819) 376-0987```

## Solving system of ODEs by Laplace transform...

Asked by:

I'm trying to solve this system of ODEs by Laplace transform.

> de1 := d^2*y(t)/dt^2 = y(t)+3*x(t)

> de2 := d^2*x(t)/dt^2 = 4*y(t)-4*exp(t)

with initial conditions

> ICs := y(0) = 2, (D(y))(0) = 3, x(0) = 1, (D(x))(0) = 2

Using

> deqns := de1, de2

and

> var := y(t), x(t)

I need to solve it for both y(t) and x(t), I have tried this by:

> dsolve({ICs, deqns}, var, method = laplace)

And

> dsolve({ICs, deqns}, y(t), method = laplace)

> dsolve({ICs, deqns}, x(t), method = laplace)

However I get this error message:

Error, (in dsolve/process_input) invalid initial condition

Any help is appreciated

 First 73 74 75 76 77 78 79 Last Page 75 of 87
﻿