## Monte Carlo integration...

I am trying to use Monte Carlo integration example given at: https://www.maplesoft.com/products/maple/new_features/maple15/examples/montecarlo.aspx.

After coding the procedure, the statement approxint(x^2, x = 1 .. 3)

gives error: Error, (in approxint) invalid input: `if` expects 3 arguments, but received 1

But I have used in exactly the same way as given in the page. What is the problem

## Solve several inequalities with 1 variable...

Dears,

Realted to this post I am trying to solve the following inequalities in the variable "t"

Cosenos := proc (m, x) local C, x1, x2, y, x0, i, j, j0, k, K0, Kaux, K1, L, R, S, t1, t2, t, r1, r2, d; C := NULL; K0 := NULL; Kaux := NULL; K1 := NULL; L := NULL; R := NULL; S := NULL; d := nops(x); for k to d-1 do if evalf(x[k]) = 1 then j0 := m else j0 := 1+floor(x[k]*m) end if; K0 := K0, [(j0-1)/m, j0/m] end do; K0 := [K0]; S := K0[1][1] <= t, t <= K0[1][2]; for k from 2 to d-1 do S := S, K0[k][1] <= 1/2-(1/2)*cos(Pi*m^(k-1)*t), 1/2-(1/2)*cos(Pi*m^(k-1)*t) <= K0[k][2] end do; K1 := solve({S}, t) end proc

However, this run very very slow....For instance, for Cosenos(50, [.3225, .25877, .325, 1])

we interrupt the opertion after more than 10 minutes...I am doing something wrong?

## How do I get Maple to show only one instance of ea...

I've asked Maple to solve an equation, and the result is difficult to read because it is separating single variables into multiple multiplicative instances of that variable raised to different powers. Is there a way to tell Maple to combine all instances of a variable into one instance raised to the appropriate power?

For example:

pceq := (2/3)*Pi*G*rho^2*R^2:
pnonreleq := (1/5)*(3*Pi^2)^(2/3)*hbar^2*n^(5/3)/m:
rhocomposition := rho = m/((4/3)*Pi*R^3), n = Z*rho/(A*m_p):
obj := subs(rhocomposition, rhocomposition, pceq = pnonreleq):
req := R = simplify(solve(obj, R)[1]);

The above code produces

In the first place, I can't see why Maple would show what should be (Z/A)^(5/3) as z^3 / (A * [A*Z^2]^(2/3)), but I'd like to somehow tell it to not break up the variables like this.

Is this possible? Thanks for the help.

## Unable to create a custom probability distribution...

I want to use Burr distribution and hence I am creating custom probability distribution. The pdf of Burr is not piecewise. It gives error.

burrpdf := alpha*k*((x-gamma)/beta)^(alpha-1)/beta(1+((x-gamma)/beta)^alpha)^(k+1)

burr_distribution := Statistics[Distribution](PDF = (alpha > 0, beta > 0, gamma > 0, k > 0), burrpdf)

Error, (in Statistics:-Distribution) invalid input: too many and/or wrong type of arguments passed to NewDistribution;

## Unable to compute integral...

I am unable to compute integral, it displays the function itself.

burrpdf := alpha*k*((x-gamma)/beta)^(alpha-1)/beta(1+((x-gamma)/beta)^alpha)^(k+1)

mean_burr_temp := simplify(x*burrpdf)

int(mean_burr_temp, x = gamma .. t)

## isolve for inequalites...

Hi!

Looking the Maple's help, I see that the command "isolve"  tries to solve an equations   over the integers. Then, given m>1 and t in the interval [0,1], How can used this command to find an integer j>=1 such that (j-1)/m<=t<j/m. That is, fin j such that t belongs to the interval [(j-1)/m,j/m].

## Maple incorrectly ignores fixed points of permutat...

Let G be a permutation group like S3. I have a loop in which x runs over the elements of this group.

When I ask the cycle type of x or the number of orbits of x, Maple gives an incorrect result in some cases. This is because when x=(12), for example, Maple considers x to be an element of the permutation group S2 instead of S3, i.e. it fails to understand that (12) should really be (12)(3) in this case.

Is there a way I can instruct Maple to consider the permutation (12) as a member of S3 instead of S2?

I am posting the code as requested. Something like this:

with(GroupTheory);

for x in Elements(S3) do print(x, numelems(Orbits(PermutationGroup(x)))) end do;

(I don't know how to format this, sorry)

This code is simply asking how many orbits exist in the action of each element of S3. The result will be 1 for x=(12), however, because Maple will assume that (12) is acting on the set {1,2}, even though I am considering (12) as an element of S3. I want it to consider (12)=(12)(3) and let it act on the set {1,2,3}, thereby producing 2 orbits, not 1.

## How to plot the following function containing an i...

How to draw the following function containing an imaginary number

 >
 >
 >
 >
 >
 >
 >
 >

## Determination of the angles of the manipulator...

by: Maple 17

On the example of a manipulator with three degrees of freedom.
A mathematical model is created that takes into account degrees of freedom of the manipulator and the trajectory of the movement from the initial point to the final one (in the figure, the ends of the red curve). In the text of the program, these are the equations fi, i = 1..5.
Obviously, the straight line could be the simplest trajectory, but we will consider a slightly different variant. The solution of the system of equations is the coordinates of the points of the manipulator (x1, x2, x3) and (x4, x5, x6) in all trajectory. After that, knowing the lengths of the links and the coordinates of the points at each moment of time, any angles of the manipulator are calculated. The same selected trajectory is reproduced from these angles. The possible angles are displayed by black color.
All the work on creating a mathematical model and calculating the angles can be done without the manipulator itself, is sufficient to have only the instruction with technical characteristics.
To display some angles, the procedure created by vv is used.
MAN_2.mw

## Kernel connection lost during summation...

Hello everybody.

Maybe I'm so lucky but I found again the problem connected with "Kernel connection lost"; error (see one of my previous post regarding this ). But now I used Intel-based processor.

Here the simple code of my program:

restart;
R0 := exp((2*Pi*I)*n^2*z);
R1 := sum(R0, n = -infinity .. infinity);
R1 := abs(R1)^2;
R2 := exp((2*Pi*I)*(n+1/2)^2*z);
R2 := evalf(sum(R2, n = -infinity .. infinity));
R2 := abs(R2)^2;
R := evalf(sqrt(Im(d))*(R1+R2));
plot(R, k = 1 .. 10);

After this I got message from Windows 7 (64 bit) that mserver.exe has stopped and then above mentioned error with Kernel connection lost. Anybody knows what is the problem?

P.S. Maple 17 (64 bit).

## Find an integer...

Dears,

Consider the problem min{|a+2*i| : i integer}, where a is a number of the form 2k*t for a fixed integer k>1 and t in [0,1]. A simple checking shows that i must belongs to the interval [-(1+a)/2,(1-a)/2] and as the legth of this interval is 1, there is two integers in this interval.

How can I compute this integer i with Maple? The commad "ceil" not seem very suitable. Some idea?

## using @ for multiple applies...

Say I want to use multiples commands on an expression

simplify(expand(eval(expr,x=y+1)))

Now similar I could use

(simplify@expand)(expr)

but how do I do it with sth like eval which needs extra arguments

(simplify@expand@eval(x->y, x=y+1))(f)

??

## plot functions in the same figure...

i want to plot these four functions in the same graphe f(x)=1/1+x^2 on the intervalle [0,1] also the second function the line curve that connects the two points(0,f(0)) and (1/3,f(1/3)), the third function the line curve that connects the two points (1/3,f(1/3)) and (2/3,f(2/3)), the fourth function the line curve that connects the two points (2/3,f(2/3)) and (1,f(1))

can anyone help me ?

## combine while excluding a variable...

I have this expression:

-sin(alpha)*(sin(theta1)*cos(theta)-cos(theta1)*sin(theta))

which can be combined using a trigonmetric formula...

However when using combine() it also tried to sneak alpha into the result...

Is it possible to freeze the variable alpha when combining?