These are Posts and Questions associated with the product, Maple 15

Hello,

I'm rediscovering maple, and i want to find the zeroes of an expression by numerical calculation. But it seems that the function fsolve can not handle the monster i created...

(The monster in question here is a determinant of a 12x12 matrix with one variable an trigonometric functions that i don't want to copy pas (NDA))

The problem might comes from the extrem values it gives for a numerical calculation (slope around 1e6)...

what is the difference between using DEtools forexample bernoullisol and dsolve directly?

I have been solved a optimization problem by

When defining a polynomial as follows

p:=x->x^(r+1)-(r+1)*Sum((-1)^(r+k)*(r+k)!*(r+1)!*r!/((k!)^2*(r-k+1)!*(2*r+1)!)*x^k,k=0..r);

the result of

simplify((eval(p(1),r=1)));

is given as 1 whereas hand evaluation and also plotting the polynomial clearly shows that the result is 0. How can this behaviour be explained?

with(PIR):RPP:='PIR/homalg';

'homalg/default':=RPP;

Error, illegal use of an object as a name

var:=[I];PIR['Pvar'](var);

test1 := (3*f-f^2)*(Diff(g, x))+(-f+3*f^2-2)*(Diff(g, y))+(2-f^2)*(Diff(g, z));

test1r := solve(test1=0,f);

how to solve this in terms of Diff(g,x), Diff(g,y) etc.

after solve, it is like a bug, because result is Diff(1,x) etc where is g?

How can I install package "gravitational" for system Win7 (maple 15 and maple16)?

Hi,

I am trying to define a system of units so that I get results in the units I want.

Specifically; I want to use energy in MeV (or eV, if I have to) and masses*c^2 in MeV (or eV) as well. Related, I need ScientificConstants like Planck's h in MeVs etc.

So I say

AddSystem('Accelerator','MeV','meter','second'); UseSystem('Accelerator');

which is accepted. I even get energies converted to MeV if I enter one in a different unit and use

I am dealing with second order linear differeantial equations for RLC and forced RLC circiuts.

I am having trouble with plotting a solution curve in the phase plane.

> restart;

with(plots); with(DEtools);

L := 10; R := 20; C := 1/100;

de1 := diff(x(t), t) = y(t);

de2 := L*(diff(y(t), t))+R*y(t)+x(t)/C = 0;

DEplot({de1, de2}, [x(t), y(t)], t = 0 .. 5, [[x(0) = 10, y(0) = 0]], stepsize = 0.5,

I have a problem solving two equations with trigonometric functions. They are as follows:

eq1:= 9180pi - 16295.070arccos(-1.0 + 0.001356852103xu) + 22.1100sin(-1.0pi + 1.0arccos(-1.0 + 0.001356852103xu))(737.0-xu)+7115.0700arccos(-1.513347022 + 0.002053388090xu) - 14.6100sin(-1.0pi + 1.0arccos(-1.513347022 + 0.002053388090xu))(737.0-xu) + (0.1050As(xu-125)/(1349-xu)) - 0.1050As = 1688

eq2:= 9180pi- 16295.070arccos(-1.0 + 0.001356852103xu) + 22.1100sin...

I ran into an interesting situation with numeric integration.

For a simulation task I need to evaluate an overlap integral of a pair of 3-d ellipsoids which have a certain aspect ratio and the axes of which are tilted against each other. The distribution in each plane is Gaussian, but with largely differing sizes. While I have an analytic expression for the result; the integration (along a certain direction) cannot be done analytically and so I integrate numerically....

One of those vexing things that should be easy:

I am reading in a data file using ImportMatrix. Since the file has one column with texual data I cannot use datatype=float since that bombs. I >can< ImportMatrix with datatype=string or anything.

I then extract columns using the appropriate command and want to do something numeric with them (only with those columns that have correct floating point numbers). That fails since my data are now not floats (even when using anything...

Hi

I have problem with solve polynomial...maple can't manage with it.

My polynomial :

1-r+4*r^12*x^12+16*r^6*x^4-8*r^7*x^6+2*r^8*x^8+24*r^6*x^2-36*r^7*x^4+24*r^8*x^6-6*r^9*x^8-24*r^7*x^2+60*r^8*x^4-80*r^9*x^6+60*r^10*x^8-24*r^11*x^10-x+8*r^4*x^2-4*r^5*x^4-16*r^5*x^2+6*r^4+2*r^2-6*r^5-r^7+4*r^6+8*r^8*x^2-28*r^9*x^4+56*r^10*x^6-70*r^11*x^8+56*r^12*x^10-28*r^13*x^12+8*r^14*x^14-r^15*x^16-5*r^3

I'm getting for *good roots* and something "1/4(RootOf(.......

Hi all,

I am looking for a package to calculate the power mean (or Hölder mean) defined like

Mean_p(x_n)=(1/N Sum(x_n^p,n=1..N))^(1/p)

in a hopefully obvious notation. I can obviously implement this definition directly in Maple, but there are subtleties (like the limit of p=0) which maybe a mature implementation already covers (?).

I am still mostly on Maple 15 which does not have such a Mean. If 16 or 17 have it I would not becessarily know.

k:=<1,2,3>:

mu:=<4,5,6>:

sigma:=<7,8,9>:

solutions:=<1,2,3,4>

product(add(1/((I*.5)*exp(solutions[i]*sqrt(2)*sigma[j]+mu[j]-ln(k[j]))+1)^k[j], i = 1 .. 4), j = 2 .. 3);

returns "Error, bad index into Vector"

You must be logged into your Facebook account in order to share via Facebook.

Click the button below to share this on Google+. A new window will open.

You must be logged in to your Twitter account in order to share. Click the button below to login (a new window will open.)

Please log-in to your MaplePrimes account.

Wrong Email/Password. Please try again.