Unanswered Questions

This page lists MaplePrimes questions that have not yet received an answer

Hi Mapleprimes,

I wrote some quick maple code -

for a from 1 to 10 do
if isprime(a)=false then print(ifactor(a)) else print('1',a) end do;
end do

I want to use the factors command to reproduce oeis.org/A027750/

Can anyone help me with this?

Regards,

Matt

 

Hi Mapleprimes,

I wrote some quick maple code -

for a from 1 to 10 do
if isprime(a)=false then print(ifactor(a)) else print('1',a) end do;
end do

I want to use the factors command to reproduce oeis.org/A027750/

Can anyone help me with this?

Regards,

Matt

 

Hello,

In the maple.ini, i defined some folders where i have put my packages.
It seems that, when I made some altercations in the maple.ini, they are not taken into account immediately.
Do you have ideas about when the update of the maple.ini is taken into account ?
And how can I force the update ? Should I restart my computer so that the update of the file is taken into account ?
Thanks a lot for your feedback.

 

Find the set of solutions of each of the linear congruence:

a) x≡3x≡3 (mod 5).

b) 2x≡52x≡5 (mod 9).

Hello I have the following small piece of code.

XMLTools[Print](MathML[Export]('sin(theta)'=0.25));

which exports sin(theta)=.25, how do I force the 0 to display.

 

Also

why does the following fail:

sol:=solve([cot(x)=2,x>=0,x<2*Pi],x,AllSolutions, Explicit);

i tried to solve a nonlinear ode with numerical method but maple can't solve it and this error occur:

Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging

my maple codes are attached below:

numeriacal_sol.mw

can any help me?

Hello,
I need to formulate the follow relationship in proper math symbols:

Differential in A (last days value minus today's) has a tendency to reach the Differential in B (today's)

I though this could be expressed with

AΔ -> BΔ

BUt I guess there are more elegant and mathematically correct ways to do this in Maple?

thank you!
Dave

Hi

I have the transition matrix used in Markov chain

A := Matrix([[alpha, beta, gamma], [delta, epsilon, zeta], [eta, theta, mu]])

I would like to write a system of equations that can be solved  to get a Markov chain  irreducible and aperiodic

All the entries of the transition matrix are in the interval [0,1)

 

Many tanks for any help

for example

func1 := proc(system1)

for i from 1 to 100 do

solve([system1[1], system1[2]],[x,y]);

od:

end proc:

 

func1([diff(y,t) = data[i+t+1], diff(x,t) = data[i+t+1]])

i is depend on the for loop inside a function, but woud like to pass this system into a function with i

this will cause error

how to write better for passing a system as parameter using variable inside a function?

Hi guys,

I've had only a little experience with Maple, but I decided to use it for preliminary frequecy response calculations. The funny thing is, I have already the solution in some way, but I'm too stupid to get it working. I can't see the mistake, however, it should have something to do with the H_n(f) function and the other normalized functions.

NASA has published a nice paper which explains the calculations, however, they have used MathCAD. Anyway, I don't think this should be a problem. Here is the documentation: https://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/20070016024.pdf 

The error message is for example for the H_n(f) function:

Warning, expecting only range variable f in expression H__n(f) to be plotted but found name H__n

 

Please help me.

Maple file: calculations.mw

i use optimization package with constraint hello >= 0

Minimize(xx=0, {hello >= 0})

but solution only return the case when hello = 0

how about hello > 0?

i would like to find all possible set of solutions using this constraint

do i need to set upper bound, such as {hello <= 7, hello >=0}

can it return solution when hello = 1.1, 1.2, ...2, 2.1, 2.2, 2.3, ....7

I am looking for a numerical solver for a parabolic PDE (up to 2nd order derivatives but no mixed ones) on the spatio-temporal domain [X x Y x T], either as an external package or as MAPLE code.  

I have coded the method of lines on the domain [X x T] and indeed also used pdsolve as a check for that case. However, pdsolve (numerical) cannot solve the PDEs on the domain [X x Y x T].  The run times and memory requirements for the latter case would of course be significantly greater.  

I am about to code up the method of lines (in MAPLE) on the domain [X x Y x T], but am wondering whether there exist external FORTRAN or C code packages that would be faster if called up in MAPLE and whose results would then be post-pocessed in MAPLE.

Does anyone have any suggestions?

MRB

 

hi--how i can solve following equation?

thanks

 

Eq.mw


Maple Worksheet - Error

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

Download Eq.mw

 

I've implemented the optimal taxation model proposed in this paper using Maple.

But it never stops running and get stuck in the last line for integral computation. Any idea of what's wrong with that?

This is the last line:

Here is the full code.


Please help me on this :

restart; with(PDETools), with(plots)

n := .3:

Eq1 := (1-n)*(diff(f(x, y), `$`(y, 3)))+(1+x*cot(x))*f(x, y)*(diff(f(x, y), `$`(y, 2)))-(diff(f(x, y), y))/Da+(diff(f(x, y), y))^2+n*We*(diff(f(x, y), `$`(y, 2)))*(diff(f(x, y), `$`(y, 3)))+sin(x)*(theta(x, y)+phi(x, y))/x = x*((diff(f(x, y), y))*(diff(f(x, y), y, x))+(diff(f(x, y), `$`(y, 2)))*(diff(f(x, y), x))):

Eq2 := (diff(theta(x, y), `$`(y, 2)))/Pr+Nt*(diff(theta(x, y), y))^2/Pr+Nb*(diff(phi(x, y), y))*(diff(theta(x, y), y))/Pr+(1+x*cot(x))*f(x, y)*(diff(theta(x, y), y)) = x*((diff(f(x, y), y))*(diff(theta(x, y), x))+(diff(theta(x, y), y))*(diff(f(x, y), x))):

Eq3 := Nb*(diff(phi(x, y), `$`(y, 2)))/(tau*Pr)+Nt*(diff(theta(x, y), `$`(y, 2)))/(tau*Pr)+(1+x*cot(x))*f(x, y)*(diff(phi(x, y), y)) = x*((diff(f(x, y), y))*(diff(phi(x, y), x))+(diff(phi(x, y), y))*(diff(f(x, y), x))):

ValWe := [0, 5, 10]:

bcs := {Nb*(D[2](phi))(x, 0)+Nt*(D[2](theta))(x, 0) = 0, f(0, y) = ((1/12)*y)^2*(6-8*((1/12)*y)+3*((1/12)*y)^2), f(x, 0) = 0, phi(0, y) = -.5*y, phi(x, 12) = 0, theta(0, y) = (1-(1/12)*y)^2, theta(x, 0) = 1, theta(x, 12) = 0, (D[2](f))(x, 0) = Da^(1/2)*(D[2, 2](f))(x, 0)+Da*(D[2, 2, 2](f))(x, 0), (D[2](f))(x, 12) = 0}:

pdsys := {Eq1, Eq2, Eq3}:

p1 := ans[1]:-plot(theta(x, y), x = 1, color = blue):

plots[display]({p1, p2, p3})

 

``


 

Download untitle_2_(1).mw

1 2 3 4 5 6 7 Last Page 1 of 200