In a recent post, the following inequality was proved with Maple:

(a,b,c,d >= 0).

Here is another direct proof attempt.

f:=(a+b+c+d)^2*(a*b+a*c+a*d+b*c+b*d+c*d)^2-144*(a^2+b^2+c^2+d^2)*a*b*c*d:
g:=expand(eval(f,d=1)):
s:=minimize(g, [a=0..infinity,b=0..infinity,c=0..infinity]):
length(s);   # huge
        304856
map(evalf@evalf[500],s);

map(evalf@evalf[1000],s);

So, Maple returns the expression

min(0,r1,r2,r3)

where r1,r2,r3 are huge expressions containing RootOfs. In order to evaluate them, several hundreds of digits are needed.
The solution seems to be correct, but the question is: may we (mathematically) accept it? What do you think?

Does anyone know how to do crossover between block of binary number?

Example :

If I have 4 blocks, [0011001],[0011100], [0101010], and [1001110].

the first block will crossover with the second block, the second block will crossover with the third block, the third block will crossover with the last block and the last block must crossover with the first block. 

I had tried the code as below. 

for i from 1 to 4 do

for j from 1 to 7 do

block[i][j]:=block[i+1][j];

end do;

new_block[i]:=block[i];

end do;

The problem is the crossover can be done from 1 with 2, 2 with 3 and 3 with 4. But this code unwork for crossover for last block with the first block. Can someone help me? thank you in advance.

Hi,

I have a linear mathematical model that is in deterministic form, but actually one of the input parameters has a stochastic nature and follows a log normal distribution function. and I need to include this in the formulation. Please help me with this questions:

1) Do I need to use stochastic programming methods to solve it? if yes, how?

2) How can I include the distribution function in my model?

3) Can I solve the these kinds of problems using Maple?

How to random generate 2 integer number with condition x must be less than y ?

I had tried the coding below. But the output shown not what i want.

x:=RandomTools[Generate](integer(range=1..7);

y:=RandomTools[Generate](integer(range=1..7);

if x>= y then x:=RandomTools[Generate](integer(range=1..7); y:=RandomTools[Generate](integer(range=1..7); end if;

Example:

I get x:= 5 , y:=1, then i need to regenerate the value of x and y using the last statemet above. but the problem is sometime i might still get the value of x>= y for example on second generate x:3, y=2. I need a code that x and y must be regerate until it get x<y. Can someone help ? Thank you in advance.

How difficult is it to simulate gravitational influences and perturbing effects on celestial orbits with Maple? Could this syntax http://www.maplesoft.com/applications/view.aspx?SID=4484&view=html be altered without excessive changes to consider these aspects?

Are there somewhere worksheets to take a look at as an introduction and to see how such goals would be approached and implemented?

Some time ago, I have used the Virtual 3D Solar System code to plot a few interesting figures with asteroids. It can be found here: http://www.maplesoft.com/applications/view.aspx?SID=4484&view=html . I plotted the asteroid belt as a closed, warped surface according to how the asteroids are inclined, with the arguments of perihelion and longitude of ascending node where they cross the equatorial plane to the other celestial hemisphere. The plot was only done with axes in units of distances, for convenience's sake in astronomical units of 1 AU = 150 million km.

But I am strongly wondering if it were somehow possible to alter or adjust the code such that there is a zylinder placed around such plot where one can see celestial coordinates (Rectascension, Declination) on the inside, so that if one wants to know what coordinates an orbit of a specific asteroid has that one can look in the chart and go out on the balcony and set the telescope to these coordinates. Of course the asteroid most likely is not there on his entire orbital path, but you have at least the location where the asteroid CAN be. Would this be possible to plot?

This probably doesn't work if a closed cylinder is placed around the Virtual Solar System. But how about scrolling the cylinder and the orbits flat on a 2D plot?

hello dear,

I was intersted in finding Weyl scalars for a given metric, would be really helpful if anyone can give me an example of how to find a Weyl Scalar for a given metric or set of null tetrads using Debever formalism in Tensor package(more specifically, how to define h in debever formalism for a given set of null tetrad or for a given metric).

Regards,
Suresh

hello dear,

I was intersted in finding Weyl scalars for a given metric, would be really helpful if anyone can give me an example of how to find a Weyl Scalar for a given metric or set of null tetrads using Debever formalism in Tensor package(more specifically, how to define h in debever formalism for a given set of null tetrad or for a given metric).

Regards,
Suresh

hello,
I try to draw together point (0,0,12) and the line  x=4t, y=-2t, z=2t on one graph 3D. Somebody has some pointers? Thank you 
Agnieszka

Dear all

I have such an expression in the file (test.mw). Obviously it is not a simple form. 10⁹ should be cancelled both in denominator and nominator. I have applied the function simplify(). However, it doesn't work.  I hope someone can do me a favor. That will help me a lot. Thanks

hi.please help me for solve algebric equations below?

solve or fsolve dos not any answer?

i think use Newton  Raphson, because solve or fsolve not work.

thankscomparision.mw

Download comparision.mw

if three particles set into 3 energy states ( each values of energy states 0 , e, 2e)

how to calculate P in each configuration by a particle that follows boson and f-b

thank you

I'm having trouble plotting a couple things. I have

eq := diff(y(x), x$3)+3*diff(y(x),x$2)+12*y(x);
soln := dsolve(eq, y(x));
soln := evalf(soln);
PartSoln1 := dsolve({eq, y(0) = a,y'(0) = 0,y''(0) = 0}, y(x));
curves := {seq(PartSoln1, a = -3 .. 3)};
Then when I try plot(curves, x = -1..5, y = -5..5); I get Warning, expecting only range variable x in expression PartSoln1 to be plotted but found name PartSoln1.

Also,
charEq := r^3+r+1 = 0;
soln := solve(charEq, r);
soln := [evalf(soln, 5)];
soln := map(Re, soln);

I tried a few things, but can't figure out how to plot charEq, including the real roots.

Thanks for any help,

Heather

Here is a code

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/EQ.mw .

Download EQ.mw

A fragment of code

for b in extra_bcs do try print(b = 10^(-2)); res[b] := dsolve(dsys4 union {b = 10^(-2)}, numeric, initmesh = 2024, output = listprocedure, approxsoln = [omega2 = 0.1e-2, s(x) = cosh(upsilon*x)-cos(upsilon*x)-(cosh(upsilon)+cos(upsilon))*(sinh(upsilon*x)-sin(upsilon*x))/(sinh(upsilon)+sin(upsilon)), g(x) = sin(((2*n+1)*(1/2))*Pi)], abserr = 0.1e-1) catch: print(lasterror) end try end do; indx := indices(res, nolist); nops([indx]); res[indx[i]]; seq(subs(res[indx[i]](1), omega2(1)), i = 1 .. nops([indx]))

Hello everybody,

I have an interesting math problem, here it is:

Add the following operators: +, -, * in the middle of the following string 123456789 to get the final result: 2016

The operators can be added randomly at any place, it isn't necessarily placed in between single numbers.

To convince your result, one should provide the code or explain explicitly your method.

Have fun!!!!