Hallo,

I am rather a beginner in Maple, so I will be grateful for understanding.

I try to solve the equation with one unknown. To create this equation I firstly had to solve the set of four equations so it is very complicated. Maple is able to create the plot of it but command fsolve is not sufficient to solve it. Here is my procedure:

restart:
z1:=H/mg*cosh(mg/H*(x+C11))-C21:
z2:=H/mg*cosh(mg/H*(x+C12))-C22:
tg_alpha1:=diff(z1,x):
tg_alpha2:=diff(z2,x):
r1:=H*subs(x=xP,tg_alpha2)-H*subs(x=xP,tg_alpha1)-P:
r2:=subs(x=0,z1):
r3:=subs(x=l,z2)-h:
r4:=subs(x=xP,z1)-subs(x=xP,z2):
mg:=2:l:=20:h:=5:xP:=l/5:P:=10:
C := solve({r1, r2, r3, r4}, {C11, C12, C21, C22}):
s1 := int(sqrt(1+(diff(subs(C, z1), x))^2), x):
s2 := int(sqrt(1+(diff(subs(C, z2), x))^2), x):
s := subs(x = xP, s1)-subs(x = 0, s1)+s2-subs(x = xP, s2):
eq1 := l0-subs(x = l, s):
l0 := 40:
plot(eq1, H=10 .. 11, colour = red);
H := fsolve(eq1 = 0, H = 10.7 .. 10.8);

I would be grateful for any help and indicating my mistakes.

Thanks

Iza

It seems if I copy some text into one of the cells of a maple spreadsheet I get an error.

Specifically the error only occurs if a space exists in the text.  Also entering text with spaces produces the error as well.

hello there - typed into a cell produces   Error, missing operator or ';'

Similarily    5 4 we receive Error unexpected number

There is only two options for the cell properties under type - Symbolic and Floating Point.  Another option to fix this issue would be to have a Label option?  

Hello all!

I did

with(plottools):

display(pieslice([0,0],5,Pi/4..7Pi/4,color=yellow))

but I don't know how to do "Animating Pacman". Can you help me? I read it on Mapleprimes but I don't really understand. Can you explain it to me? Thank you so much! Sorry, My English is not very good

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?