Unanswered Questions

This page lists MaplePrimes questions that have not yet received an answer

Hi,

 

Say I run a procedure, it could take about 3 minutes to run, or maybe hours, or longer.

Is it possible to trigger a sound, like "ding", after it's finished?

I may like to read something while the program is running, but I want to see the results immediately.

It would be nice to have such function.

 

Also, for procedures that runs for a "long" time, could I get 10 (or 15) minutes sound notice?

 

Thanks,

 

casper

I'd like to try the new Maple IDE

http://www.maplesoft.com/products/toolboxes/IDE/index.aspx

Maple IDE - Integrated Development Environment

Is there a way to obtain a trial version even for a week just to try it out? I do not want to spend all this money then find out it does not work well.

any one here uses this and can recommend it or comment on how well it works? For example, can one run Maple code directly from it and have the code execute in Maple, or must one load the file manually from the editor into Maple to run it each time they make any changes to the code?

Hi all

I have a mathematical problem and I asked it in various sites but the answers till yet are not correct.

Assume that we have:

T[m]:=t->t^m:
b[n,m]:=unapply(piecewise(t>=(n-1)*tj/N and t<n*tj/N, T[m](N*t-(n-1)*tj), 0), t):

where n,N,tj are known constants. furthermore assume that we want to comute the following integral:

for following approximations:

I have written the following code but it seems to be incorrect:

V1:=Vector([seq(seq(b[n,m](t),m=0..1),n=1..3)]);
V:=evalf(V1.Transpose(V1));

the original program is :

taaylor.mws

I will be so grateful if any one can help me to solve it by maple

Mahmood   Dadkhah

Ph.D Candidate

Applied Mathematics Department

Hi all

 

I am trying to maximize a function f(x,y,z,w) in terms of x. (Only x is treated as a variable, and the others are treated as parameters).

However, all I know is that y,z,w they are parameters and they are non-negative. I have already tried with the "optmization help page" from maplesoft's website, and it looks like it will search the range of x,y,z,w, and it will return numerical values at which this function is maximized. 

 

 

However, what I want is instead a close-form solution of x=g(y,z,w) that will maximize the function.   In other words, I would like to keep the parameter in symbolic forms. 

 

Can Maple do that?

How can I get the Automorphism Group of a group G as a Permutationgroup?

 

Best regards

 

Kurt Ewald

Hello,

this is the second time I'm writing.

I posted this question in June http://www.mapleprimes.com/questions/201781-System-Of-Parametric-Equations.

This time I have  a similar problem because I'm trying to find a solution for a parametric system of equations but the number of equations and parameters is much bigger and using the tips you gave me last time I couldn't reach any result.

Here is the system:

1) alpha[1]=v*a*u*b ;
2) alpha[2]=v*a*u*(1-b);
3) alpha[3]= v*z*c*(1-a) ;
4) alpha[4]=v*z*(1-a)*(1-c) ;
5) alpha[11]=1/2*v*a* u* b* (-p*u*b+p*u*b*a+b*g-g);
6) alpha[22]=1/2*v*a*u*(1-b)* (p u b-p u b a-b g-p u+p u a);
7) alpha[33] =1/2*v*c*z*(1-a)* (c* (-z*p*a+q)-q);
8) alpha[44]=1/2*v*z*((1-a)*(1-c)* (c*z*p*a-z*p*a-q*c);
9) alpha[12]=v*a*u*b*(1- b)*(-p*u+p*u*a+g) ;
10) alpha[13]=v*a*u*b*z*c*p*(1-a) ;
11) alpha[14]=a*u*b*z*(1-a)*(1-c) ;
12) alpha[23]=a*u*z*c*(1-a)*(1-b);
13) alpha[24]=v*a*u*z*p*(1-a)*(1-b)*(1-c);
14) alpha[34]= v*c*z*(1-a)*(1-c)*(-z*p*a+q);

 

I have 14 equations/unknowns and 8 parameters (a, b, c, u, v, z, p, q).

I would like to write this system only in terms of alphas. In order to do so, I usually try to find the value for the parameters and the substitute them into the equations (and I have already found b,c,g,q using this technique) but I couldn't manage to find all of them. 

Howveer, as you suggested me, with Maple there is the command "eliminate" that implement exactly what I'm looking for but I can't make it work.

This is my code:

> sys := {alpha[1] = v*a*u*(1-b), alpha[2] = v*a*u*b, alpha[3] = v*z*c*(1-a), alpha[4] = v*z*(1-a)*(1-c), alpha[11] = (1/2)*v*a*u*(1-b)*(p*u*b-p*u*b*a-b*g-p*u+p*u*a), alpha[12] = v*a*u*b*(1-b)*(-p*u+p*u*a+g), alpha[13] =      z*c*a*u*(1-a)*(1-b), alpha[14] = v*z*a*u*p*(1-a)*(1-b)*(1-c), alpha[22] = (1/2)*v*a*u*b*(-p*u*b+p*u*b*a+b*g-g), alpha[23] = v*z*c*a*u*b*p*(1-a), alpha[24] = z*a*u*b*(1-a)*(1-c), alpha[33] = (1/2)*v*c*z*(1-a)*(c*(-z*p*a+q)-q), alpha[34] = v*c*z*(1-a)*(1-c)*(-z*p*a+q), alpha[44] = (1/2)*v*z*(1-a)*(1-c)*(c*z*p*a-z*p*a-q*c)};

> eliminate(sys, {a,b,c, p, q, u, v, z});

> simplify(%, size);

 

I also tries to substitute in the system the four parameters I already found but still I can't find a solution.

What am I doing wrong? Or the problem is that it is too complicated?

 

Thank you for your attention,

Elena

I have a Maple program that fits a nonlinear function with 9 parameters to a dataset that contains multiple records; each record gets fitted separately. The records are ordered by a certain experimental variable (an angle of a target). The initial parameters for each fit are the solutions of the previous fit; in this case I get relatively good initial estimates and not too much trouble fitting. I use Statistics:-NonLinearFit to do the actual fitting.

Since this is a bit time consuming and I need to do this a fair amount of times, I'd like to see if I can't speed his up by employing parallel programming.

One way to speed up a fit is to compute the sum of residuals in parallel; is there a way to do this in Maple?

Another way would be to run several fits at the same time. The problem I have with that in this case is that I use the results of the previous fit as starting values for the next fit. I might be able to increase the stride of this scheme; i.e. for 4 processes, use the result of fit 1 as starting point for fit 5; fit 2 for fit 6 etc. But before I embark on this non-trivial change; can Statistics:-NonLinearFit even be used in this way (i.e. is it thread safe)?? My fitting function is an actual function (as opposed to a procedure with local variables) so I don't think the function itself would cause trouble.

I could also contemplate the Grid package; except that I am concerned about the overhead moving the data over to each node (although it would all sit on one machine with 4 processors).

I realize that I may have to do some experimentation by myself, but if there is any experience out there I'd like to know about it...

Mac Dude

my int is in below.

my_int.mw

my int has hyperbolic and trigonal parts thus when they are multiplied , Maple is not able to integrate!

 

Now, How can I solve this integral?

thank you for help

my int is:

 

``

restart

eq1 := m*(diff(diff(w(x, t), t), t))+diff(diff(EIy*(diff(diff(w(x, t), x), x))+(EIz-EIy)*(diff(diff(w(x, t), x), x))*(theta(x)+phi(x, t))^2, x), x)-Pz = 0

l := 16; m := .75; EIy := 0.2e5; EIz := 0.400e7; GJ := 0.1e5; mj := .1; Pz := 5000; Mz := 0; theta := proc (x) options operator, arrow; 0 end proc

w := proc (x, t) options operator, arrow; q[1](t)*(cosh(1.8751*x/l)-cos(1.8751*x/l)+(-1)*.734096*(sinh(1.8751*x/l)-sin(1.8751*x/l)))+q[2](t)*(cosh(4.6941*x/l)-cos(4.6941*x/l)+(-1)*1.01847*(sinh(4.6941*x/l)-sin(4.6941*x/l)))+q[3](t)*(cosh(7.8548*x/l)-cos(7.8548*x/l)+(-1)*.999224*(sinh(7.8548*x/l)-sin(7.8548*x/l))) end proc

phi := proc (x, t) options operator, arrow; q[4](t)*sqrt(2)*sin(1.5708*x/l)+q[5](t)*sqrt(2)*sin(4.7124*x/l)+q[6](t)*sqrt(2)*sin(7.8540*x/l) end proc

NULL

f[1] := int(lhs(eq1)*(cosh(1.8751*x/l)-cos(1.8751*x/l)-.734096*(sinh(1.8751*x/l)-sin(1.8751*x/l))), x = 0 .. l)

``

restart;
eq:=diff(y(x),x) - y(x)^2 - x*y(x) - x + 1;
r:=dsolve(eq,y(x));
latex(r);

Gives an error "(in fprintf) string expected for string format". The problem seems when the math contains "erf"? Both examples below fail with result that contains this.

Maple 18, on windows 7.

Another example:

restart;
eq:=diff(y(x),x) + a*y(x)*(y(x)-x) - 1;
r:=dsolve(eq,y(x));
latex(r);

Another example, which seems to fail for different reason

restart;
eq:=diff(y(x),x) + a*x*y(x)^3+b*y(x)^2;
r:=dsolve(eq,y(x));
latex(r);

hello

we have an exam next week and I want to know

how I can write (fordo) in maple for numerical integration

in different methods such as trapezoid , newton cotes and so on.

thanks

This integral is defined as

Where

 

Now, How can I solve this integral?

tnx for your help.

 

 

Hi,

With Maplets, it is possible to create popup windows to tell or ask the user something, but only when it is appropriate. That is, it is not there all the time. A particular usage would be to give the user a facilitating question when they have incorrect entries.

Is that possible in some way with Möbius apps also? Or does one have to always keep a blank TextArea available, just a blank space if invisible perhaps but using up space in the presentation nonetheless, to give the student the information?

AOA. I want to plot the graph of the following function 

 

A new generalized complex representation of Euler gamma function in terms of Dirac delta function, which is

GAMMA(s) = 2*Pi*(sum((-1)^n*Dirac(s+n)/factorial(n), n = 0 .. infinity))

where*s = sigma+i*tau

for differenet values of parameters

 

Hi everyone,

I was trying to write a Maple sheet to calculate some perturbation theory expension as close to the "book notation" as possible. Lets for example consider a linearly perturbed harmonic oscillator with H = hω(n+1/2) and V = λ(a+a), this is one of the classical examples as it can be solved analytically by completing the square. In the Kato formulation (as used in the appended Maple worksheet), as well as in Rayleigh-Schrödinger perturbation theory one uses the projector onto the complement of the unperturbed state. I would like to do this by defining a projector in Maple. I've read the examples concerning projectors in Maple but it seems I have to use a rather ugly workaround by treating the states "below" the unperturbed state and the ones "above" separately. While this is a little annoying in one dimension it becomes a major nuisance in higher dimensions.

My question would be: How do I define the projector onto the complement of some state?

Cheers, Sören

restart; with(Physics); Setup(mathematicalnotation = true)

a := Annihilation(N, 1):

assume(`and`(`in`(m, nonnegint), m > 0)):

Physics:-Ket(N, m)

(1)

H := Physics:-`*`(Physics:-`*`(h, omega), n+1/2):

`&Delta;E__2` := simplify(value(Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(Dagger(psi), V), 1/(E(m)-H)), V), psi)))

-lambda^2/(h*omega)

(2)

`&Delta;E__4,111` := simplify(value(Dagger(psi).V.(1/(E(m)-H)).V.(1/(E(m)-H)).V.(1/(E(m)-H)).V.psi))

Error, (in Physics:-Dagger) numeric exception: division by zero

 

`&Delta;E__4,201` := -simplify(value(Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(Physics:-`*`(Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(Typesetting:-delayDotProduct(Dagger(psi), V), 1/(E(m)-H)^2), V), psi), Dagger(psi)), V), 1/(E(m)-H)), V), psi)))

lambda^4*(2*m+1)/(h^3*omega^3)

(3)

`&Delta;E__4` := `&Delta;E__4,111`+`&Delta;E__4,201`

`&Delta;E__4,111`+lambda^4*(2*m+1)/(h^3*omega^3)

(4)

``

Download kato_perturbation_theory.mw

First 256 257 258 259 260 261 262 Last Page 258 of 361