## Approximate state-space realization of a given imp...

Hello,

I'm trying to calculate the convolution of the radiation forve in the time domain. I need to do a state space representation of a function K(t), where :

t=[0; 1; 2; 3];

k=[0.1; 0.2; 0.3; 0.4];

If you could help me with finding the appropriate command that give the matrix A, B C and D

just a remarque that in matlab used imp2ss (based on the Hankel SVD method)

thanks

## How do I solve Base Excitation Problem of any SDOF...

Is it possible to develop a code for base excitation problem where excitation data is of earthquake acceleration?

p.s. There are 12000 data available in EQ file.

## Code Editor Syntax and Symbolic Solutions, Weird L...

Hello Everyone,

I am currently trying to tidy up my worksheet and hide all of the script/code that I wrote in the Maple 2D-Math within the execution script of buttons etc.
Was going well so far until I started transferring code where I was using a lot of symbolic math (if I understood correctly)
There is one instance where I have a variable chi which equals x divided by a.
a being a number, e.g. 5000, while x is and remains a variable.
No problem in 2D-Math as the unknown x is dragged through all following equations.

Not so in the code editor, I think. Doing the calculation results not in x/5000, but in this:

It keeps getting worse because this term will be derived, assembled into a matrix, etc. where the expressions are getting ridiculous:

Without getting into the details of what my script should do, I believe this is a very general issue, not specific to my code.

Is the code editor / "Code when Used/Clicked" not capable of symbolic math? Oder is there some setting that could potentially help me?

Best Regards,
Lennart

## Trigonometric circle...

Hi,

I try to build a trigonometric circle with particular values Pi/6,Pi/4,Pi/3, but I find that the coordplot command with its options do not react .

Ideas ? Thanks

CercletrigonométriqueMaple.mw

## Maple UI is tiny on a high dpi screen in Linux...

I'm using Maple 2020 in Linux (Ubuntu 20.04). The screen resolution is 2560×1440, so I have system-wide setting of 200% scaling for all apps, but unfortunately Maple doesn't adopt the system setting and UI elements are tiny. To make it tolerable, in "Tools - Options - Interface", I've set Default Zoom to be 200% and checked the box "Large toolbar icons". However, menus (File, Edit, View etc.) and dialog boxes are still tiny.

Am I right that the Maple GUI is written in Java / Swing? I've heard that Java 9 has better support for high dpi scaling. Does the default Maple installation use Java 9?

## Can I solve one equation in one variable in parall...

Dear All,
I want to solve a highly nonlinear equation in one variable. The equation includes some undetermined parameters. The running time is long in serial programming. Can I help me to execute the following command in parallel? I know that the ‘solve’ command is not thread-safe.

solve(Eq1, beta);

Let me know solve the problem, if the number of equations is larger than one, as below:
solve({eq1,eq2,eq3},{beta1,beta2,beta3});
Best wishes

## how to solve equation iteratively?...

Hello experts,I need your help to obtain the solution of the mentioned equation in the attached picture by using the newton method by supposing the random value of involved constant like D  etc. I need an algorithm for the newton method for the mentioned equation. note h=1+x^2/2

## Bug in modp1(('Embed')(...))...

Hi all,

There is issue with modp1(('Embed')(...)) function.

bug_embed.mw

As you can see in lines 1, 3, 5 Embed function in combining with Constant cutting several digits from argument.

Thank you.

## how to guess 1 or 2 when rsolve is wrong at some n...

second has rsolve solution

compare original sequence

accsum(ListTools[Reverse]([seq(patterngroup[k], k=(mm)..(mm+lengthofsequence))])),

with rsolve solution , it start from third number wrong, some number are wrong but correct at next number,

but i just want to guess next number is 1 or 2

[seq(subs(k=kk, sol), kk=0..(lengthofsequence+1))],

with(SumTools):
with(gfun):
accsum := proc(xx)
local result,mm:
result := []:
for mm from 1 to nops(xx) do
temp := Summation(xx[k], k = 1 .. mm):
result := [op(result), temp]:
od:
return result:
end proc:

patterngroup := [2, 1, 1, 1, 2, 2, 1, 1, 2, 2, 2, 1]: # left is latest, 2,1,1..
lengthofsequence := 5:
for mm from 1 to 12 do
if guessgf(accsum(Reverse([seq(patterngroup[k], k=(mm)..(mm+lengthofsequence))])), x) = FAIL then
print(patterngroup[mm]):
else
sol := rsolve(listtorec(accsum(ListTools[Reverse]([seq(patterngroup[k], k=(mm)..(mm+lengthofsequence))])), u(n))[1], u(k)):
acxxc := accsum(ListTools[Reverse]([seq(patterngroup[k], k=(mm)..(mm+lengthofsequence))])):
olist := [seq(patterngroup[k], k=(mm)..(mm+5))]:
#start from 0
print(patterngroup[mm], olist, accsum(ListTools[Reverse]([seq(patterngroup[k], k=(mm)..(mm+lengthofsequence))])), [seq(subs(k=kk, sol), kk=0..(lengthofsequence+1))], sol, subs(k=(lengthofsequence+1), sol), subs(k=(lengthofsequence+1), sol)-acxxc[nops(acxxc)]):
end if:
od:

## Solving an optimal control problem...

Hello all

I wanna solve an optimal control problem and I have searched the Internet but I could not find any tutorial or video course on how to solve it with the Pontryagin maximum principle method. It is my first time that I want to use MAPLE for solving an optimal control problem and I would be thankful if someone can help me.

$$\max \int_{0}^{1} x_{2} [u(t)-u(t)^2] dt$$

$$\dot{x}_{0} = -(1-u(t)) x_{0}(t)+2 x_{1}(t)$$

$$\dot{x}_{1}(t) = (1-u_{t}) x_{0}(t) +2 x_{2}(t) -[3-u(t)]x_{1}(t)$$

$$\dot{x}_{2}(t) = (1-u(t))x_{1}(t) -2 x_{2}(t)$$

$$0 \le u_{t} \le \frac{1+t^2}{1+t}$$

Thanks

## Given a lagrangian in general relativity , how ca...

Given a lagrangian  in general relativity , how can I calculate equation of motion using euler lagrangian equation in maple software ?

the lagrangian is

L = 1/2 g[mu,nu] diff(x[~mu], lambda) diff(x[~nu],lambda)

here lambda is the affine parameter.  and the metric is

metric coordinate = (t,u,x,y)

## How can I solve two equations?...

How I can find the coefficient an, and bn according to the following solution?

the coefficients an and bn can be found by solving the
two linear equations that come from V = V[0] at eta=eta[0] and
V = V[1] at eta=eta[1], and comparing with following Eq in each
case.

Best

2.mw

## how to plot the Profiles of optimal control functi...

I an unable to prot the optimal  control function and the controls individualy for the diterministic mathematical model with an two or more optimal control. May  any one can give me a sample program on it. Such us the following figure

You may give me direction even.

## Bug in modp1(('Rem')(...))...

Hi there.

It seems like a bug in modp1(('Rem')(...)) with large polynomials with large coefficients.

bug_rem.mw

It needs the file polys.m there:

https://dropmefiles.com/9fATR

File polys.m too big for this forum so I used dropmefiles.

Polys.m contain two polynomials x and f_t with large degrees: degree(f_t) = m = 50021, degree(x) = 2*m - 2 = 100040 and large coefficients up to 2^N, where N = ceil(m / 2)+2 = 25013.

I just compute rem(x,f_t) mod 2^N.

As you can see in the first part of doc bug_rem.mw I decreased coefficients of polynomials by additional mod 2^N (with Embed function), where N = floor(N / 2) = 12506. WIth these decreased polynomials and decreased N modp1(('Rem')(...)) function works well and use maximum about 2.5 Gb of RAM.

But in the second part of doc bug_rem.mw with original polynomials and N = 25013 modp1(('Rem')(...)) use maximum about 3.5 Gb of RAM and crash with error:

Error, Maple was unable to allocate enough memory to complete this computation.  Please see ?alloc

This is strange and looks like a bug considering that the test server has 48 Gb of RAM.

Is it a bug or modp1(('Rem')(...)) just need more than 48 Gb of RAM?

How many RAM it needs for this computation?

Thank you