## How to get procedure to handle f or f(x,y)...

I have a procedure to truncate algebraic equations. Found the basic answer hrere. I would like to make it handle equations and functions i.e f:=a x^2 +b y... and f(x,y):=c x y^2.... I can make it handle the first type but not the second.

 (1)

 (2)

 (3)

 (4)

 (5)

 (6)

 (7)

 (8)

## Arbitrary sign changes...

Could anyone tell me why Maple changes the sign of some expressions? For example...

How to prevent this changes?

Denis

## curious about dsolve for matrix equations...

I was pleased to observe that in Maple 2018 one can use dsolve to get solutions of matrix ODEs, e.g., X'(t)=A(t) X(t). The command is the obvious dsolve({diff(X(t),t)=A(t).X(t), X(0)=Xo}).

I was then surprised that Maple balked at the commands dsolve({diff(X(t),t)=A(t).X(t), X(0)=Xo},numeric) and dsolve({diff(X(t),t)=A(t).X(t), X(0)=Xo},numeric,{x(t),y(t)}).  Is there some error on my part?

## How can I linearize this dyamical system containin...

How to linearize this complex dynamical system around x1(t)=0, x2(t)=0, x3(t)=0. Basically i want to keep nonlinear part up to fifth order on the right hand side and later want to drop it. How can i do it. The maplesheet Dyn_system.mw is attached.

How can I generate a sequence of random numbers Xi and Yi in [-1,1] ?

## Why I couln't receive the value of Wn[0, 0, 1] in ...

Why I couln't receive the value of Wn[0, 0, 1] in the following loop?

> restart;

>for n from 1 to 3 do

Wn[0 .. 1, 0 .. 1, n] := n*ArrayTools:-Alias(Matrix(2, 2, [[3, 1], [1, 3]]), [0 .. 1, 0 .. 1])

end do:
>Wn[0, 0, 1];

Wn[0, 0, 1]

## Hard Integral numerically....

Hello

I want to compute integral with Maple,but returns unevaluated for me.

int((-5*ln(x)^4*Pi^4-20*ln(x)^2*Pi^4-8*Pi^4+120*MeijerG([[0, 0], [1, 1, 1]], [[0, 0, 0, 0, 0], []], x^Pi))/(120*Pi^4*(-1+x)^2), x = 4/10 .. 6/10, numeric);

# ???

Thanks.

## How do I get output diplayed with 4 significant di...

Hello all

I am teaching engineering and I use Maple to make handouts with solutions to assignments.

In Maple 2017 I used interface(displayprecision=4) to make maple it display all outputs with 4 significant digits.

Just to be clear:

If the actual value is 123.456789, then interface(displayprecision=4) makes it display 123.5 in Maple 2017. This is just how I want it.

If the actual value is 1.23456789, the interface(displayprecision=4) makes it diplay 1.235 in Maple 2017. Again, this is perfect, just how I (and most engineers) like it.

Now I have installed Maple 2018 and now interface(displayprecision=4) means 4 decimal points rather than significant digits.

So:

123.456789 => 123.4568 in Maple 2018

1.23456789 => 1.2346  in Maple 2018

This is NOT what I want.

Digits or, say, evalf[4] will not work as a fix, as I only want the DISPLAYED value to have 4 significant digits, I want the entire value used in calculations.

I am told (by Carl Love, thank you Carl) that displayprecision=4 means 4 decimal points in most versions of Maple, but I do need 4 significant digits.

Is there a fix in Maple 2018? Something that will make all output values display with 4 significant digits?

I am considering re-installing Maple 2017 and reverting to it, it turns out to be a huge work to format each output value.

I hope someone can help me, thank you in advance.

## how to find turning points, critical points and in...

Hi guys,

I have some information about function f(x)=x^2+y(1-x)^3/4 where x domain is [0,1] and y is a parameter but i can not obtained this information with maple.

I know function f(x) for y>1.4266 is monotically decreasing and there is no turning points,(actually i can compute infection points of above function: y=1.4266 and x=4/5 )

I know function for 1.4266>y>1.0144 has two turning points a minimum for x<4/5 and a maximum for x>4/5 , but i dont know how y=1.0144 is calculated and what is the nature of this point !!! and how this relates to x><4/5 !

I know fucntion for y<1.0144 has to two turning points that the global maximum of function happens at x<0.97702, but also i can not compute 1.0144 and how this relates to x<0.97702 !!

1.mw

Thanks

## How to convert to sum?...

I am trying to evaluate some functions when u=v so function goes to 0/0

Example given below of the expansion series I need to produce first. But I can't produce it.

 (1)

 (2)

 (3)

 (4)

 (5)

 (6)

 (7)

 (8)

 (9)

 (10)

 (11)

## Integration rule to Maple int(1/x,x) = ln(abs(x))...

Hello

I'd like Maple to return ln(abs(x)) for int(1/x,x) instead of ln(x).

I tried convert MMA code,but I failed.

realIntegrate[f_, x_Symbol] :=
Simplify[Integrate[f, x] /. Log[expr_] :> Log[Abs[expr]],
x \[Element] Reals]; Unprotect[Integrate];
Integrate[f_, x_Symbol] /; ! TrueQ[$flag] := Block[{$flag = True}, realIntegrate[f, x]]; Protect[Integrate];

My maple code works only for very simple cases.

REALINT := proc (f, x)

if typematch(int(f, x), ln(y::anything), 's') then

assuming([simplify(ln(abs(rhs(s[1]))))], [x in real])

else int(f, x)

end if

end proc;

REALINT(1/x, x);

#ln(abs(x))

REALINT(1/(x+1), x);

#ln(abs(x+1))

REALINT(1/x+1, x);

#x+ln(x) dosen't work.

Thanks.

## Simplifying with sqrt...

I have an expression involving a sqrt, and i'd like to simplify it, but simplify(expression,sqrt) doesn't seem to fully help.

Eq1 := [(1/2)*(R[b]*kh[a2]+R[m]*kh[a2]-Rh[m]*kh[a2]+sqrt(R[b]^2*kh[a2]^2+2*R[b]*R[m]*kh[a2]^2-2*R[b]*Rh[m]*kh[a2]^2+R[m]^2*kh[a2]^2-2*R[m]*Rh[m]*kh[a2]^2+Rh[m]^2*kh[a2]^2))/kh[a2] = 0, 0 = 0, -(1/2)*(R[b]*kh[a2]+R[m]*kh[a2]-Rh[m]*kh[a2]+sqrt(R[b]^2*kh[a2]^2+2*R[b]*R[m]*kh[a2]^2-2*R[b]*Rh[m]*kh[a2]^2+R[m]^2*kh[a2]^2-2*R[m]*Rh[m]*kh[a2]^2+Rh[m]^2*kh[a2]^2))/kh[a2] = 0]

simplify(Eq1, sqrt)

gives

[(1/2)*(R[b]*kh[a2]+R[m]*kh[a2]-Rh[m]*kh[a2]+sqrt(kh[a2]^2*(R[b]+R[m]-Rh[m])^2))/kh[a2] = 0, 0 = 0, -(1/2)*(R[b]*kh[a2]+R[m]*kh[a2]-Rh[m]*kh[a2]+sqrt(kh[a2]^2*(R[b]+R[m]-Rh[m])^2))/kh[a2] = 0]

which on paper simplifies to:

[R[b]+R[m]-Rh[m]=0,0=0,R[b]+R[m]-Rh[m]=0]

is there a way to get maple to show this?

[In part I am trying to better understand how to manipulate sqrt expressions in maple]

## Displaying two vector fields by its line elements ...

Please how do I display these two fields in the euclidian plane:

(X,Y)=(2e_1, e_2)=(2\partial/\partial x_1,\partial/partial x_2)

with their **elements representation** ? What is the appropriate command to display this?