Download q1stringtool.mw

I recieved the following error:

Error, (in ifactor/QuadraticSieve:-SieveCube) sieving failure

But when I review the procedure ifactor, it doesnt tell me anything about A Quadratic Sieve algorithm, and it's really long and looks dodgey and suspicious, like line 24 for example, why is it computing the greatest integer divisor of a local variable and a random enormous square free number? and then another with an additional factor a few lines later? 

Hello

I need to solve or reduce (similar to the command Reduce in Mathematica) sets of nonlinear equations.  One such example is shown below:

eqns := {-1+theta[3, 6] = 0, 1-theta[3, 6] = 0, alpha+rho-theta[2, 2]+theta[3, 3] = 0, -theta[3, 6]^2+1 = 0, theta[2, 2]*theta[3, 6]-alpha = 0, theta[2, 2]*theta[3, 6]^2-alpha = 0, -2*theta[3, 3]*theta[3, 6]-2*rho = 0, theta[1, 2]*theta[2, 1]*theta[3, 6]^2+1 = 0, -alpha^2+rho^2+theta[2, 2]^2-theta[3, 3]^2 = 0, -theta[2, 2]^2*theta[3, 6]+2*theta[2, 2]*theta[3, 3]*theta[3, 6]+alpha^2+2*alpha*rho = 0, -theta[1, 3]*theta[2, 2]^2*theta[3, 0]+theta[1, 3]*theta[2, 2]*theta[3, 0]*theta[3, 3]-alpha^2*beta-alpha*beta*rho = 0, -theta[1, 2]*theta[2, 1]*theta[2, 2]*theta[3, 6]+2*theta[1, 2]*theta[2, 1]*theta[3, 3]*theta[3, 6]-alpha-2*rho = 0, -theta[1, 2]*theta[2, 1]*theta[2, 2]*theta[3, 3]+theta[1, 2]*theta[2, 1]*theta[3, 3]^2+theta[1, 3]*theta[2, 2]*theta[3, 0]*theta[3, 6]+alpha*beta+alpha*rho+rho^2 = 0, -alpha^2*rho-alpha*rho^2+theta[1, 2]*theta[2, 1]*theta[2, 2]-theta[1, 2]*theta[2, 1]*theta[3, 3]+theta[1, 3]*theta[3, 0]*theta[3, 6]-theta[2, 2]^2*theta[3, 3]+theta[2, 2]*theta[3, 3]^2+alpha+beta+rho = 0}

 and the indeterminates are:

fc := {theta[1, 2], theta[1, 3], theta[2, 1], theta[2, 2], theta[3, 0], theta[3, 3], theta[3, 6]}

Since I do know the solution, I issued the following command to check for typos.

seq(subs(theta[1,2]=-1,theta[1,3]=-1,theta[2,1]=1,theta[2,2]=alpha,theta[3,0]=beta,theta[3,3]=-rho,theta[3,6]=1,eqns[i]),i=1..nops(eqns))

and the outcome is zero for all equations.

When I try the command solve as follows:

solve(eqns,fc);

the result is

{theta[1, 2] = theta[1, 2], theta[1, 3] = theta[1, 3], theta[2, 1] = -1/theta[1, 2], theta[2, 2] = alpha, theta[3, 0] = -beta/theta[1, 3], theta[3, 3] = -rho, theta[3, 6] = 1}

that should be right but it is not what I am expecting.  

How can maple return the solution needed?

Some sets of solutions do not have a solution as the one above.  Some indeterminates cannot be found, is there a way maple returns the solution of the ones that can be solved and reduced the set of equations into two parts, solved ones e non solved ones?  I can provide an example if needed.

Many thanks.

Ed

Non-Linear overdetermined equations - which is best method? with less number of iterations

When I use Jacobi it takes 25 iterations.

Any other method which takes less iterations?

I search saveplot and sleep function here

i had kept close filename function and move plot setup default to first line

but thread no sleep in maple 12

then I remove sleep and can save image file

but after I call it in a function and run a for loop for this function

the file is not updated after sleep 3 seconds 

how to run a for loop call it and it can refresh 

I tried to control z-axes by writing

plot3d(sech((1/2)*sqrt(3/5)*(x-2*t)), x = -10 .. 10, t = -10 .. 10, scaling = constrained, style = patchnogrid)

but I could not change the height of z axes. I want to change the dimensions like in the following figure.

example-graph.pdfexample-graph.pdf

Although the range of z-axix is small but it appears very clear. Can anyone help me, please.

The following message is displayed at the beginning of my worksheet when i reopen it, no data is lost because the nature of my code retrives it's optimal boundaries for the loops but i dont know what the hell this is telling me this for:

"Warning, .hdb help databases are deprecated, 'C:\Users\the_r\maple\toolbox\Syrup\lib\Syrup.hdb' will not be used, see ?HelpTools,Migrate help page for more information"

It appears in blue standard size courier font if that is help thank you

Is it asking me to migrate? My passport was stolen by clandestine pharmacists 
 

 Dear fellows
I have the problme with first seq command.
for i while i <= M1-2 do for j while j <= M1-2 do for k while k <= M1 do Eq[i, jk] := simplify(eval(R, [x = i/(M1-2.), y = j/(M1-2.), t = k/(M1-1.)])) = 0 end do end do end do:

Sol := fsolve({seq(BC1[m8, m9]$m8 = 1 .. M1, m9 = 1 .. M1), seq(BC2[m10, m11]$m10 = 1 .. M1, m11 = 1 .. M1), seq(BC3[m12, m13]$m12 = 1 .. M1, m13 = 1 .. M1), seq(BC4[m14, m15]$m14 = 1 .. M1, m15 = 1 .. M1), seq(IC1[m4, m5]$m4 = 0 .. M1, m5 = 1 .. M1), seq(IC2[m6, m7]$m6 = 0 .. M1, m7 = 1 .. M1), seq(Eq[m1, m2, m3]$m1 = 1 .. M1-2, m2 = 1 .. M1-2, m3 = 1 .. M1-1)});

Got the following error.

Error, invalid input: seq expects its 3rd argument, step, to be of type numeric, but received m3 = 1 .. 1

Please help me in this regard.

Error, (in Optimization:-NLPSolve) the objective gradients at the initial point are too small

So this error indicates that the initial point is a local minimum or very very close to it right?

I'm very new to Maple (and coding in general). I'd like to know a very large coefficient of a certain 𝑞-series involving the Dedekind 𝜂-function. For example

s:=series(subs(q=q^6, eta^4), q, N);

where eta := series(q^(1/24)*product(1 - q^n, n = 1 .. 100), q, N). I know the command coeftayl(s, q=0, N) but for certain 𝑞-series it doesn't recognize the numerical value of this coefficient when 𝑁 is too large (even after subtracting the principal part of 𝑠, if there is any) - instead it realizes the coefficient as a HUGE limit of something. Cheers.

I need the final line of the uploaded worksheet to be the "Select Action", preprogrammed prior to the execution of the insert xml command, how do i do this?

Maple Worksheet - Error

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

Download ASK.mw

I am trying to do matrix multiplication with Dirac gamma matrices, but am getting a recursion error. 

the final input line, which is causing the error, doesn't seem to be inserted, but it is:

Library:-PerformMatrixOperations((mt . (Dgamma[2])) . m)

can anybody tell me what i'm doing wrong? i don't understand why 

Dgamma[2] . m 

works but 

mw . Dgamma[2] . m 

gets the recursion error. 

then i saw something that made me think i needed to use 'matrix', instead of 'Matrix'. so i tried that(gamult2) - now the 2nd multiply works, and the first fails.

any help greatly appreciated!

thanks,

larry

Download gammult.mw

gammult2.mw

hello

why this pdsolve not be answer

pde := diff(T(x, y, t), x, x)+diff(T(x, y, t), y, y) = (1/0.44e-2)*(diff(T(x, y, t), t));
conds := T(0, y, t) = 1, T(1, y, t) = 1, T(x, 0, t) = 1, T(x, 1, t) = 1, T(x, y, 0) = 0;
`assuming`([pdsolve([pde, conds])], [0 <= x and x <= 1, 0 <= y and y <= 1, t >= 0]);
 

Dear all

I have a fourth-order PDEs, with some boundary condition, I would like to compute an explicit solution or plot the solution in the interval [0,2]

PDEsolve.mw

thinks

Hello everyone,

Say I have the differential equation:

x₁'(t) = 3C₁*cos(t) + 2C₂*sin(t)

x₂'(t) = 5C₁*cos(t) + 1C₂*sin(t)

How can one draw phase portraits in maple?