Hi I was wondering if theres a way to smoothen a signal in maplesim.
I was trying to create a simple moving average but it doesnt seem to be able to create a feed back loop?
Much appreciate any help, or if there are other ways within maplesim to do it.

thanks
pjf

I am writing a package to define a few Records and collect some quantities and formulae to do some arithmetic with.

At some point I started running into a severe problem: Each time I want to use one of the exports form the module, the Maple kernel dies with a "Stack limit exceeded" message. The issue is: This seems to happen only on Maple 2016. I had some issues in Maple 2015 but that got sorted by rewriting some stuff. I never had any issue like that with Maple 15. I went as far as commenting out most of my .mapleinit code to make sure that did not cause issues. The code works as intended in Maple 15 and in Maple 2015.

I am posting this here and am asking others to see if they can reproduce this. Needless to say I will SCR this. But I actually need to get this to work.

The error can be triggered by simply entering

Booster;

after loading the package. Needless to say, it should just echo the name "Booster", to get more one would need to enter eval(Booster) (due to the last-name evaluation rules).

Thanks,

M.D.

BooMBA.mw

I solved two inequalities using

[solve]({f(theta) >0, f(theta) <=Pi}, theta)

for some complicated function f of theta, an solve gives me a list of intervals such as

[{theta <= Pi-arccos(-3/4+(1/4)*sqrt(13+16*sqrt(2))), arccos(3/4) < theta},

{arccos(-3/4+(1/4)*sqrt(13+16*sqrt(2)))+Pi <= theta, theta < 2*Pi-arccos(3/4)},

{theta <= 3*Pi-arccos(-3/4+(1/4)*sqrt(13+16*sqrt(2))), 2*Pi+arccos(3/4) < theta}]

I would like to convert each of the intervals to the range form a..b so I can plot the function over this range.

Is there a comand to convert something like {theta <=b, a < theta} to
a..b (assuming a and b are real numbers in some form)?  I can think of a rather complicated procedure to do it, but there should be something simpler.

I have computed the eigenfunction expansion for f(x)=x on 0<x<1 in terms of the eigenfunctions exp(-x/2)*sin(n*Pi*x).

I wish to calculate the weighted L2 error in this expansion (the weight function is w(x)=exp(x)).

Specifically, I want to determine how many terms in the eigenfunction expansion are necessary for the error to be less than say 0.3.

Here is the code:

f := x -> x
w := x -> exp(x)
assume('n', integer);
y :=  (n, x) -> exp(-x/2) sin(n Pi x)       
c := n-> (int(f(x)*y(n, x)*w(x), x = 0 .. 1))/(int(y(n, x)^2*w(x), x = 0 .. 1))
Fourierf := (n, x) -> sum(c(j)*y(j, x), j = 1 .. n)

fsolve(Lerror(n) = 0.3, n);

This seems to run forever without giving a value of n.  I know this is a large computation, but it seems that Maple should be able to handle it.  Does anyone have any suggestions?

Heather

Hello everybody,

I want to solve this pde.the desire solution is V(r,z). three boundar conditions are written that two of them are related to rhe radial and one is related to the longitudinal coordinate.

I attached the solution for you. but this solution is derived by Matlab. Now, I just want to resolve it by Maple, but I couldn't reach it. Please let me know the correct way asap.

Thanks a lot.PDE.mw
 

Download PDE.mwPDE.mw

For remainder of division of a (multivariable) polynomial to several polynomials at a same time one can use NormalForm in Maple. It is easy to write a procedure to also show the division but I wonder if there is any determined command such as NormalForm for this aim?

Hi,

do anybody know, how I can see how a function is programmed? I want to know more about dsolve function and what is hidden behind.

Thank you all

Hej all,

I have a problem. When i have a 9-1 vector "DD" how can i get maple to solve DD^2 my equation is A=(pi/4)*DD^2 should be simple.

Also when i have two vectors one is 9-1 and the other is 5-1 i want to divide the 5-1 / 9-1 but i want all possible soluitons? not just 5 results but 45 results.(the zip fuction only gives me 5 results)

Pls help 

Kind regards

Howdy!!

Just thought I'd put a quick one in. I'm trying to solve the Schrodinger equation given as:

under the constraint

I'm trying to plot this for a range of values of the parameters a,b and fixed p=0,1.

I'm aiming to get a solution of the form

and plotting it.

Thanks in advance

I wish to compare five letter words in a list of 5 letter words with a string like "*e*d*".  The asterisks represent unknown letters; the e and d, known letters.  I wish to compare this string with each word in the list and discard any words which do not match.  I would end up with a list of matching words eg ["fends", "sends", . . .]   Is there a string function to compare partial matches?  I've looked at the compare and match commands,  but they don't seem applicable.

   Any suggestions would be most appreciated.

Regards, David 

I copy a program from a book and try to run it, but when i try to run it return:

too many levels of recursion

Here is my program:

Hope every one can help me solve this problem. Thank!

I am trying to solve in Maple the system 
d{x}/dt =(1  1 \\1  -1).{x}
with initial conditions {x(0)} =(1 \\ -1),
where {x} is a vector.

This is what I did up to now

with(DEtools):
M1:=<1,1;1,-1>;
matrixDE(M1,t);
FM:=matrixDE(M1,t)[1];

but I am not sure how to correctly input the initial conditions to find a solution to the problem. I would appreciate any help provided.
 

In a trajectory. How can i return the range if i know the position-function. ... return x when y=0?

Recently I have used a FilterForge filter to generate hyperbolic tiles from images (Lua code with Möbius transform etc. available). If an input image (1 see below) is selected a tile (2) is generated if the hyperbolic parameters (5,3,2) were chosen.  How could this be made with Maple?

The origin of this problem is that I want to bound a norm below. After som assigning of variables etc it boils down to minimizing;

S(x,y,z,w,A,B,E,F,G,H):=A*x^2+A*y^2+B*x^2+B*y^2+(-E*G+F*H+1)*x*z/(1-2*E*G+2*F*H+(F*H)^2+(E*G)^2+(E*H)^2+(F*G)^2-4*E*F*G*H)+(-E*G+F*H+1)*y*w/(1-2*E*G+2*F*H+(F*H)^2+(E*G)^2+(E*H)^2+(F*G)^2-4*E*F*G*H)-(-E*H+F*G)*y*z/(1-2*E*G+2*F*H+(F*H)^2+(E*G)^2+(E*H)^2+(F*G)^2-4*E*F*G*H)-(-E*H+F*G)*x*z/(1-2*E*G+2*F*H+(F*H)^2+(E*G)^2+(E*H)^2+(F*G)^2-4*E*F*G*H)-6*E*x+6*F*y-6*G*z+6*H*w

I first treat the x,y,z,w as coefficents and optimize over them( these are real and imaginary parts of complex parameters of my norm which on paper is a sum).
I use standard way to get the minia w.r.t these 4 variables via "solve". Then I put the solution into the a new function N(A,B,E,F,G,H) and try to optimize that with NLP and get one solution, but this is only local and it has the following message attached to it;

"Warning, no iterations performed as initial point satisfies first-order conditions"

I want a global minima for my polynomial N in 6 variables and I have some constarints on them aswell. Furthermore I put all variables to be real.

What command or package can I use to get this?