There is a enought complex function H(y) of one variable (includding definite integrals and 3 parameters). One needs to solve equations like that H'(y)=H''(y)=0 and build surface in parametrical space with respect to such solutions. Straightforward method leads to 'Warning, solutions may have been lost' for solve command. How to reformulate problem to become it solvable?

Hello,

I have got a problem using LPSolve with the option "method=interiorpoint". If I do not add this option it works fine but if I add this option I get the message "Error, (in Optimization:-LPSolve) unrecognized storage format, storage = []"

Could you please help me to deal with that? Thank you.

opoly := (1-x^2)*diff(f(x), x$2) - x*diff(f(x), x) + n^2*f(x) = 0;
sol := op(2,dsolve(opoly));
gen := sum(cos(n*arcsin(x))*y^n/n!, n=0..7);

with(gfun):
seriestoratpoly(series(gen, z, 8), ['egf']);

error not a series, why?

seriestoratpoly(gen, ['egf']);

refere to the following link, it is possible to use a few terms to calculate a generating function, but where is the function generating? how do it  work in gfun?

http://www.google.com.hk/url?sa=t&...

opoly := (1-x^2)*diff(f(x), x$2) - x*diff(f(x), x) + n^2*f(x) = 0;
sol := op(2,dsolve(opoly));
P := sum(cos(n*arcsin(x))*z^n/n!, n=0..7);
with(gfun):
algeqtoseries(P, x, z, 7, true);

Error, invalid input: gfun:-algeqtoseries expects its 1st argument, Pol, to be of type polynom(anything, [x, y]), but received 1+(1-x^2)^(1/2)*z+(1/2)*cos(2*arcsin(x))*z^2+(1/6)*cos(3*arcsin(x))*z^3+(1/24)*cos(4*arcsin(x))*z^4+(1/120)*cos(5*arcsin(x))*z^5+(1/720)*cos(6*arcsin(x))*z^6+(1/5040)*cos(7*arcsin(x))*z^7...

Hi,

How to solve this equation?

Any Idea?

Thx

1.778895759*Sigma-1831241.099/(76553.66445-.576e-5*Sigma^2)+6600.970252*Sigma/(76553.66445-.576e-5*Sigma^2)+.5739576533e-1*Sigma^2/(76553.66445-.576e-5*Sigma^2)+.4735119433e-4*exp(.7618258041e-2*Sigma)*Sigma^2/exp(.9051395693e-5*Sigma)/(76553.66445-.576e-5*Sigma^2)-39332.76308*exp(.7618258041e-2*Sigma)*(41.17+1/(1-exp(-160)))/exp(.9051395693e-5*Sigma)/(-69.17083220+Sigma)/(-69.17083220-Sigma)-629324.2088*exp(.7618258041e-2*Sigma...

Hello!

 I'm new to maple and I'm having some difficulties with coding.

I'll explain the problem

I have 2 vectors, a1 and a2, and I want to make a new set of vectors R(n,m) where R(n,m) = n*a1+m*a2 for n and m from 1 to 100 (as an example). This will create n*m vectors (in this case 10000) that will be marked as R(0,0), R(1,0), R(0,1), R(1,1), R(2,1), ...

It seems as though this would be simple but I haven't found a way to get the repetition statement ...

Whenever I try to plot a simple/small function in 2d, Maple literally takes forever to calculate it whenever the frame rate is too high.

For instance, the following command (from a tutorial, so it should work):

animate( plot, [exp(-x/5)*sin(x),x=0..t], t=0..20, frames=100 );

Will not compute, unless i change the number of frames to 50 or below.

I tried this on a friends computer (he has a mac, but has way weaker performance) and it worked just fine.

Hello,

To label plot axes as Greek letters, one should use the "labelfont=[SYMBOL]" parameter, but this option makes both axes in Greek letters. Is it possible to make only the horizontal axe as a Greek letter, while the vertical axe remains as an English word?

Thanks in advance.

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...