I have a procedure, lfGn, with two or three real arguments, and I wish to find the maximum value of the procedure.  I was trying the Optimization package, and found this to be extremely slow.  Optimization was called with the following code: the first few lines just set some variables, for use in the last line, which does the call.

max2lfGn := proc (X::Vector) local a, d, e, f, H, m, s;
e := 10.^(1-floor((1/2)*Digits));
a:= Statistics:-Mean(X);
d := Statistics:-StandardDeviation(X);...

Mazimize only gives me the first maximum it finds. I want to find all maximums in an interval, or just the largest NUMERICALLY. is there a workaround?

I should convert factors of a generic equation into standard form (scientific notation) :

 example_eq := 345.1234 x1(t) + 34.5691 u1(t) - 5423.8932 y1(t) + ... many other terms...

Is there some maple command or procedure to convert automatically every factor of a generic (very long)   equation ?

Manipulating a system of long equations maple gave in output the following error message:

I have a function in maple with one variable (say f=2*x). I would like to generate a nx2 matrix with the variable (x) in the first column from 0 to 100 and the function value (f) in the second column. I would then like to export this numerical data. How can I do this in Maple?

sc_vs_t.mw

sc_vs_t.pdf

I want the solution of an ODE for different values of d say d=3,4,5,6... etc (see the mw file). I want to plot a(t) vs t for these different values of d.

Now what is wrong with my code ?  How to fix it ? 

I also  added a Mathematica work of the same thing (see the pdf) ,...

Dear guys, Can anyone tell me why my program does not work for some initial values? All my equations are:

> alpha := (6*h^2)^(1-n)*(1-k)^(1-n)*((1-m+k)/(2*n-1-k));

> eq := z-> 1=(m*(1+z)^3-k*(1+z)^2)*h^2/(H^2)+(((2*n-1)-(k*(1+z)^2*h^2/H^2))*((1-m+k)/(2*n-1+k))*((((H^2/h^2)-k*(1+z)^2)/(1-k))^(n-1)));

> Y := z->if not type(z,numeric) then 'procname(z)' else (fsolve(eq(z), H=h)) end if;

> l := dsolve({D(L)(z) = L(z)/(1+z)+(1+z...

lets say I have the folowng

a:=1

b:=3

c:=5

is it posible to have it like this

h:=a+b*c=1+3*5=16

hi I'm new to maple and I want it t perform some filtering

I want it to select a nummber in array within given condiitions

> t := [1, 4, 9];

>

> for i from 1 to 3 do if t[i]>2 and t[i]<t[i+1] then  s:=t[i] else s:=t[i]+7end ;
print(`output redirected...`); # input placeholder
Error, unterminated loop
Typesetting:-mambiguous(for i from 1 to 3 do  if t(i)gt2  then s

I have a function with one variable which I would like to plot the Kurtosis (statistical term) of (for a specific range). The purpose is to produce a plot which will show the 'peakiness' of the curve.

Say y = x^3 + 2*x^2 + 3*x;

For the range, say x = 0..10, I would like to plot the kurtosis (Maple function is also called Kurtosis). I have been able to do this numerically with a matrix of data but am unsure how to do this in Maple.

Your assistance is much appreciated.

Hi, 

I am relatively new to Maple and I am seeking an efficient way to create a probability histogram for a discrete probability distribution. For example,

x:      1     3   5    7

p(x): 1/8 1/4 1/2 1/8

I would like rectangles with height corresponding to the probabilities with bases centered at each x value, allowing for the possibility of gaps in the distribtuion.

http://www.solitaryroad.com/c681.html

refer to this link, dsolve the differential equation can not get generating function.

would you mind demonstrating the steps with maple?

I'm having the problem reported here back in 2008 and more recently in 2010:

http://www.mapleprimes.com/questions/39313-Postscript-Plot-Sizing

http://www.mapleprimes.com/questions/96487-Graphics-Output#comment96501

I'm not sure when the axiswidth, axisheight options of plotsetup got broken for me. I'm not certain, but it looks like I was able to control the width and height of my plots not so long ago (version 14 or 13) or am I mistaken and this bug,...

hi,

I'm having trouble with the symbolic integration of a function with respect to time. The equation to be integrated does not only contain time but also some variables which are time dependent. The equation is:

M = z²₁ * sin( z₅ * t⁾²

this has to be integrated with respect to time however z₁ and z₅ are also time dependent as they are states of a state space model:

z₁

I am a Maple 12 user (under window 7) and now trying to solve a system of equations with PHCpack. Somehow it does not work. I do follow steps provided in http://www.math.uic.edu/~jan/. The error message shows up as follows;

> restart;

> setPHCloc("C:\\phcPACK");

> f := [x*2 - x*y + 2, x^2 + y^2 - 1]:

Please see attachment.

Suppose particle A moves in counter-clockwise circular motion about the origin with radius 3,

 completing one revolution every time interval.  Particle B is in counter-clockwise circular 

orbit  about particle A with radius 1/2, completing 12 revolutions every time interval.  Find

 a  possible  parameterization of the path of particle B for one complete time interval.  

Animate your solution.