Hi Mr Preben Alsholm .

i before asked you about a problem.now i have another question...

for determine unknown parameter oemega(or eigenvalue) according to line ''newsys2 := subs(omega^2 = omega2*10^5,newsys)'' and also line ''approxsoln = [omega2 = 1, f(x) = x^2*(1-x)^2]''

power (5) in term'' 10^5''  in  (omega2*10^5,newsys) and initial guess (1)  (omega2 = 1) are very effective on final result eigenvalues .sometimes gain different eigenvalues ,which it is not impossible recognize that

which one is correct and convege.however i need Positive

 eigenvalue that is minimum between them,If and only if , it converge after
some iteration in two section of maple file which is attached as
Eigenvalue.mw

  According to my code converge occuerd until  eigenvalue gained in first section (before line `11` := ((1+6*alpha2)*(1/12))*(int(fy11^2, x = 0 .. 1))...........) are equal with those obtain in second section .(please see below file for example for_example.mw)

it is necessary mention that between section 1 and 2 is relations.for example amont of first obtained eigenvalue which obtain from first section must be repleaced in line ''approxsoln = [omega2 = 0.661514014001420, h(theta) = theta^2*(1-theta)^2]...........''in second section,

and this procedure should be continued between section 1 and 2 until convergence with desired accuracy occurred.Another relate is that,the first ODE system  can be solved using the first set of boundary data to obtain first estimate for  second section (fy11). At the next step by repeating in the same manner, 

this time by obtained function (g) at the end of maple code (Solution of the second ODE system  with second set of boundary conditions leads to the first estimate of function g),that at this stage, the first iteration is completed.

Now by replacing omega2 which is determined in this section in to the first section ,fy11 is updated and gained Further.

Next, the updated function , by continuing the iterative procedure.

in matlab bvp4c rule is used for this purpose.is this impossible in maple software ?if not please help me for solve and gain correct omega..

thanks alot

When I right click on a plot and add a label name for the vertical axis, it is shown in horizontal way next to it, what should I do if I want to have it in a vertical display along the axis written upward?

Error, invalid piecewise

but

is OK.

can someone tell me how can I derive and solve d(g1(x))=0 a complicated function as

I tried with the instructions simplify and combine but without results.

Thanks in advance.

My goal is to plot the integral J with respect to t and as you can see J is a piecewise function.

This is my code.   hw2_numerical_2.mw

Actually it's a problem about adiabatic invariants.

If you want to know the backgroud please see this link.

www.mapleprimes.com/questions/206645-How-To--Numerically--Solve-It-

The title pretty much says it, but I wrote a bunch of code that operates on a matrix of given size. Now I want to increment that size with a loop, how can I put a loop around all the code? It won't let me just copy paste it into my loop body for some reason. Thanks!

I am using solve() to give me a two variable solutions (e.g. a = 3, b = 5). Very often, I have 2 solutions, so the output will be something like {a = .8025072621, b = -.8025072621}, {a = 1.326426441, b = 1.326426441}. I store this output in a matrix called Answers, so that looks like:

Answers:= solve(equation that needs solving). 

Later when I want to access these answers, I'll do Answers[1][1] to get the first a, and Answers[1][2] to get the first b. Overwhelmingly this all works fine, but occassionaly Answers only has one set of solutions, so Answers[1][1] becomes meaningless, for example, if solve() returned {a = 3, b = 4}, then Answers[1] just returns a, and Answers[2] returns b, but Answers[1][1] gives an error. As I cannot predict when solve() will give me one solutions as opposed to 2, how can I write code to reliably access the first solution for use later?

example of non linear fractional integraL Equation 

x(t)=(t^3)+1/5 * x(t) * I^(1/2)x^2

how input the integral fraction int(^1/2) for the any function in maple ? 

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

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.

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);