Wonder if this can be accomplished in Maple.

so I have a list of 100 items labeled {1..100} of various value {$100, $160, $220, ......  , }

the task is to distribute these items among 3 people A,B,C so they get an approximately equal share.

Adding the values and dividing by 3 gives the dollar total to aim for. 

This post has C.Love procedure for evenly sized groups

 http://www.mapleprimes.com/questions/200480-Product-Grouping

but what i want is a method for different sized groups. ie 25 items for A, 35 for B and 40 for C (user defined).

additionally there is a fixed constraint: A has been bequeathed items 1,4,8; B items 2 and 20; C item 50.

restart:
S:= {3, 4, 5, 6, 8, 9, 28, 30, 35}:
SL:= [A,B,C,D,E,F,G,H,I]:
assign(Labels ~ (S) =~ SL); #Create remember table.
AllP:= [seq(P, P= Iterator:-SetPartitions(S, [[3,3]], compile= false))]:
lnp:= evalf(ln((`*`(S[]))^(1/3))):

Var:= proc(P::({list,set}(set)))
local r:= evalf(`+`(map(b-> abs(ln(`*`(b[]))-lnp), P)[]));
end proc:

Min:= proc(S::{list,set}, P::procedure)
local M:= infinity, X:= (), x, v;
     for x in S do
          v:= P(x);
          if v < M then  M:= v;  X:= x  end if
     end do;
     X
end proc:

ans:= Min(AllP, Var);
              [{3, 9, 35}, {4, 8, 28}, {5, 6, 30}]
subsindets(ans, posint, Labels);
               [{I, A, F}, {B, E, G}, {C, D, H}]

Hi

When I run a maple file it uses 25% of CPU.

When I run 2 files, the half of CPU is used.

How can I change the preferences for using most of CPU

Hello,

i need help to translate a system which is given below to a for loop.

Other wise i am writing it with myself. 

instead of doing it like that 

with(LinearAlgebra);

sys := [galerkin_funcs[1], galerkin_funcs[2], galerkin_funcs[3], galerkin_funcs[4], galerkin_funcs[5], galerkin_funcs[6], galerkin_funcs[7], galerkin_funcs[8], galerkin_funcs[9], galerkin_funcs[10]];

var := [w[1], w[2], w[3], w[4], w[5], w[6]];

Kmat, Fmat := GenerateMatrix(sys, var);

i want to do it like that.

N:=10

for i to N+4 do

sys(1,i) := galerkin_funcs[i] 

end do

for i to N do

var(1,i) := w[i] 

end do

After that i will generate matrix with this comman Kmat, Fmat := GenerateMatrix(sys, var);

But this for loop i wrote is not doing the i want to do.

Thanks for your help.

Error, (in fsolve/polynom) Digits cannot exceed 38654705646

I am using fsolve to find numerical approximations to the roots of many fairly large polynomials (degrees up to ~80).  I often get this error message and I'm not sure why.  Is there any workaround?  Any help is much appreciated.

Is there a bug for resizing the spreadsheet?  I am unable to drag the size to the full worksheet.  If it's a bug is there a workaround?

How can I generate some n-tuple random list of integers s.t. any component is between -30 and 50? For example if n=5 then four random of such 5-tuple are

[-1,2,8,7,9] , [0,-9,2,-3,-5] , [4,5,3,-8,-1] , [12, -5, 0, 6,8]

I am trying to obtain the splitting field of New_polyq. evala@AFactor did not complete. Applying splitting sequentially produced independent extensions from the first 2 (3?) factors. evala@Indep did not complete for the union of all 4 extensions.

What libraries would handle this better?

restart; _EnvExplicit:=false;interface(labelwidth=200);
Rho_polys:=rho[3,1]^3-2, rho[3,2]^2+rho[3,2]*rho[3,1]+rho[3,1]^2, 2*rho[6,1]^3+rho[6,1]^6-2, rho[12,1]^2+rho[6,1]^2-1, 2*rho[12,2]^2-rho[6,1]^2*rho[3,2]*rho[3,1]^2-2*rho[6,1]^2-2;
New_poly:=1/16*(-rho[6,1]^4*rho[3,2]*rho[3,1]-2-rho[3,1]^2*rho[6,1]^4-2*rho[6,1]*rho[3,2]*rho[3,1]-2*rho[3,1]^2*rho[6,1]+2*lambda^2)*(rho[6,1]^4*rho[3,2]*rho[3,1]+2*rho[6,1]*rho[3,2]*rho[3,1]-2+2*lambda^2)*(-2+2*rho[3,1]^2*rho[6,1]+rho[3,1]^2*rho[6,1]^4+2*lambda^2)*(-2+rho[6,1]^2*rho[3,2]*rho[3,1]^2+2*lambda^2);
sol:=solve({Rho_polys});
alias(op(sol));
New_polyq:=subs(sol,New_poly);

hi all

i have a plot and i want  make it symmetry into axes y.

restart;

I5:=Int((x^2)/((x^4-2*x0*x^2+x0^2+1)^(3/2)),x=0..infinity);
I6:=(1/2)*Int(1/(x^4-2*x0*x^2+x0^2+1)^(1/2),x=0..infinity);
f:=((-4/Pi)*((x0*I6-I5)/(x0*I5+I6)^(1/3)));
A:=1/(x0*I5+I6)^(2/3): #A:=Delta/eF

plot(A,x0=-3..5);

Hi,

I am trying to find the roots of Hankel function H1(2, z)?

j := 1;
for i from 0 to 10 do z0 := i*step+z_min; x[j] := fsolve(f = 0, z = z0, complex) end do;

for p to 10 do for j from p+1 to 9 do if `and`(Re(x[j])-Re(x[p]) < 0.1e-4, Im(x[j])-Im(x[p]) < 0.1e-4) then for i from j to 10 do x[i] := x[i+1] end do; p := p-1; break end if end do end do;

the first of these two processes works fine however the second does not. The second on is to get rid of same value solutions! I am not sure if I have missed anything and also is there a way to determine a max value in the complex domain and use it in a for loop?

any help would be great 

hi all

i have tow plots and i want to display them together, but i cant, plz help me

restart;

I5:=Int((x^2)/((x^4-2*x0*x^2+x0^2+1)^(3/2)),x=0..infinity);
I6:=(1/2)*Int(1/(x^4-2*x0*x^2+x0^2+1)^(1/2),x=0..infinity);
f:=((-4/Pi)*((x0*I6-I5)/(x0*I5+I6)^(1/3)));
B:=x0/(x0*I5+I6)^(2/3): #B:=mu/eF
C:=B/((16/Pi^2)*((x0*I6-I5)^2/(x0*I5+I6)^(2/3))): #C:=mu/epsilon0/2

p1:=plot(B,x0=0..4):
p2:=plot(C,x0=-3..0):
display({p1,p2},x0=-3..4);

After reading 

http://math.stackexchange.com/questions/674535/avoid-dividing-by-zero-with-just-variables-and-basic-operators 

I tried to check this 'trick' in maple, but see what happens in Maple 2015:

 what is the correct solution for 0^0?

Here the code for those who hate typing:

x := 0; 0^(x^2);
0
1
for x from -1 to 1 do print(x, evalf(0^(x^2), 20)) end do;
-1, 0.
0, 1.
1, 0.
for x from -1 by .1 to 1 do print(x, evalf(0^(x^2), 20)) end do;
-1, 0.
-0.9, 0.
-0.8, 0.
-0.7, 0.
-0.6, 0.
-0.5, 0.
-0.4, 0.
-0.3, 0.
-0.2, 0.
-0.1, 0.
0., Float(undefined)
0.1, 0.
0.2, 0.
0.3, 0.

Harry

Hello all! How to create this plot in Maple? 

hi all

how can i solve this program?

restart;
I5:=evalf(int((x^2)/((x^4-2*x0*x^2+x0^2+1)^(3/2)),x=0..infinity)):
I6:=evalf((1/2)*int(1/(x^4-2*x0*x^2+x0^2+1)^(1/2),x=0..infinity)):
A:=f-evalf((-4/Pi)*((x0*I6-I5)/(x0*I5+I6)^(1/3))):
B:=solve(A,x0);

I am attempting to write a Gaussian elimination routine to solve a system Ax = b using loops, but I have been having trouble.  Any help would be mcuh appreciated.  Thanks!

Hello,

I'm writing a code and I seem to have an issue when trying to implement a procedure. Here is the code:

Any help would be greatly appreciated.

Gambia Man