Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Recently, entered a engineering graduate program and will be taking a graduate engineering class in my retirement.  In our Into classes we reviewed engineering math through ODE. In the second class, we performed some numerical methods using MatLab with a small introduction.  My formal Math wass in the 1980s, but now have to fire it up for a Engineering Analysis using Michael Greenberg's Advanced Engineering Mathematics, 2nd which has Maple imbedded code and directions.

 

What would be my assignment for getting up to speed with using Maple to solve algebraic-symbol-programming problems? For instance, which tutorials in what order, or webinars? i think I once had Fortran 77 and remember zip.

 

Appreciate any good advice.

MathMan2

i tend to study a particular function each day, and noticed that the CAS literally ignores me when i as for an evaluation of the mellin transform for the lerch phi function. i just assume there to return an error if the function has properties that do not satisfy the criteria for a particular transform, in this case, null, it just pretends i didnt ask.

Assuming the maple kernel hasnt evolved conciousness and has decided to hate me im guessing someone else has run into similar circumstance.

Hi,

I am trying to download some MapleSim Robotic models from the Maplesoft website but I have 'Invalid File Format' error. Could you help please?

Best

Onder

how i can simplify

(f(x[n])/Df(x[n]));
in code

restart;
taylor(f(x), x = gamma, 8);
f(x[n]) := subs([x-gamma = e[n], f(gamma) = 0, seq(((D@@k)(f))(gamma) = factorial(k)*c[k]*(D(f))(gamma), k = 1 .. 1000)], %);

1 2
f(gamma) + D(f)(gamma) (x - gamma) + - @@(D, 2)(f)(gamma) (x - gamma)
2

1 3 1 4
+ - @@(D, 3)(f)(gamma) (x - gamma) + -- @@(D, 4)(f)(gamma) (x - gamma)
6 24

1 5 1 6
+ --- @@(D, 5)(f)(gamma) (x - gamma) + --- @@(D, 6)(f)(gamma) (x - gamma)
120 720

1 7 / 8\
+ ---- @@(D, 7)(f)(gamma) (x - gamma) + O\(x - gamma) /
5040
2 3
c[1] D(f)(gamma) e[n] + c[2] D(f)(gamma) e[n] + c[3] D(f)(gamma) e[n]

4 5 6
+ c[4] D(f)(gamma) e[n] + c[5] D(f)(gamma) e[n] + c[6] D(f)(gamma) e[n]

7 / 8\
+ c[7] D(f)(gamma) e[n] + O\e[n] /

taylor(D(f)(x), x = gamma, 8);
Df(x[n]) := subs([x-gamma = e[n], f(gamma) = 0, seq(((D@@k)(f))(gamma) = factorial(k)*c[k]*(D(f))(gamma), k = 2 .. 1000)], %);

D(f)(gamma) + @@(D, 2)(f)(gamma) (x - gamma)

1 2 1 3
+ - @@(D, 3)(f)(gamma) (x - gamma) + - @@(D, 4)(f)(gamma) (x - gamma)
2 6

1 4 1 5
+ -- @@(D, 5)(f)(gamma) (x - gamma) + --- @@(D, 6)(f)(gamma) (x - gamma)
24 120

1 6
+ --- @@(D, 7)(f)(gamma) (x - gamma)
720

1 7 / 8\
+ ---- @@(D, 8)(f)(gamma) (x - gamma) + O\(x - gamma) /
5040
2
D(f)(gamma) + 2 c[2] D(f)(gamma) e[n] + 3 c[3] D(f)(gamma) e[n]

3 4
+ 4 c[4] D(f)(gamma) e[n] + 5 c[5] D(f)(gamma) e[n]

5 6
+ 6 c[6] D(f)(gamma) e[n] + 7 c[7] D(f)(gamma) e[n]

7 / 8\
+ 8 c[8] D(f)(gamma) e[n] + O\e[n] /

(f(x[n])/Df(x[n]));
this last term did not use f(x[n]) value from above to solve it. plxx help if any one can solve it...

Let a non-planar non-self-intersecting closed polygon in three dimensions P  be given, say
with(plots): with(plottools):
P := polygon([[0, 1, 1], [1, -1, 2], [3, 0, 5], [1, 1, 1]]):
How to find the minimal surface  with the boundary P?
There is no chance to find the solution as a closed-form expression.
Thus, the numerical solution (or/and a triangulation which approximates the minimal surface up to the given accuracy) is required.

how i can find order of convergence of newton method by expanding taylor series?? plz send me code???

 

how we can solve an ode with boundaries by using RK4 and any other technique??

how we can solve an ode with exat soltuion n calculate the unknw constant which aries after solving an ode.???

I wish to study the trend of medical consultations each day during six years. Thus I expect near 2200 datas to analyse.

But some parameters are to consider :

- I don't have yet the datas per day, but the mean is about 2 consultations per day

- as it is difficult to do more than 3 or perhaps 4 consultations during one day (9h30 am - 13h pm), the others if they exist will probably be seen the next day (aso if the next day 3 news consultations occured)

- then, I don't know actually (as I expect the datas or each day but don't have now these datas) if the better distribution will be simply follow a Poisson' law, or exponential, or negative binomial, ..

- do someone have a clue for the better law given what i said ?

 

Further, I don't have a stastic program especially used for time trend, excepting Systran 13, but I don't believe that this program can be used with a theoric model of distribution, I recall that it does usual tasks, autocorrelations, saisonnal adjustments, .. but with continuous distributions I believe, and a linear model (removing the basic frequencies)

As such program (study of temporal series) is usually sold about 3000$ in France, that I don't expect to be a trader, with only one calculus to do, could anyone tell me how to adjust the better model to the 2200 datas that could be expected ?

Thx for your help, friendly yours;

Milos

I wonder why the thole procedure becomes ... when converting from 1-D math to 2-D math.

hi...please help me for solve this nonlinear equations with pdsolve

thanksoffcenter2.mw

La := .25; Lb := 0.1e-1

h := 0.4e-2

rho := 7900

E := 0.200e12

nu := .3

ve := 5

g := 9.8

M := .5

Z0 := 0.1e-2

K := 5/6

C := sqrt(E/rho)

NULL

 

PDE[1] := diff(u(x, t), x, x)+(diff(w(x, t), x))*(diff(w(x, t), x, x)) = (diff(u(x, t), t, t))/C^2

diff(diff(u(x, t), x), x)+(diff(w(x, t), x))*(diff(diff(w(x, t), x), x)) = 0.3949999999e-7*(diff(diff(u(x, t), t), t))

(1)

PDE[2] := K*(diff(phi(x, t), x)+diff(w(x, t), x, x))/(2*(1+nu))+(diff(w(x, t), x))*(diff(u(x, t), x, x))+(diff(u(x, t), x))*(diff(w(x, t), x, x))+(3/2)*(diff(w(x, t), x, x))*(diff(w(x, t), x))^2 = (diff(w(x, t), t, t))/C^2

.3205128205*(diff(phi(x, t), x))+.3205128205*(diff(diff(w(x, t), x), x))+(diff(w(x, t), x))*(diff(diff(u(x, t), x), x))+(diff(u(x, t), x))*(diff(diff(w(x, t), x), x))+(3/2)*(diff(diff(w(x, t), x), x))*(diff(w(x, t), x))^2 = 0.3949999999e-7*(diff(diff(w(x, t), t), t))

(2)

 

PDE[3] := diff(phi(x, t), x, x)-6*K*(diff(w(x, t), x)+phi(x, t))/(h^2*(1+nu)) = (diff(phi(x, t), t, t))/C^2

diff(diff(phi(x, t), x), x)-240384.6154*(diff(w(x, t), x))-240384.6154*phi(x, t) = 0.3949999999e-7*(diff(diff(phi(x, t), t), t))

(3)

 

 

#####################################

(4)

at x= La

PDE[a1] := diff(u(x, t), x)+(1/2)*(diff(w(x, t), x))^2-M*(g-(diff(u(x, t), t, t))-Z0*(diff(phi(x, t), t, t)))/(E*Lb*h) = 0

diff(u(x, t), x)+(1/2)*(diff(w(x, t), x))^2-0.6125000000e-6+0.6250000000e-7*(diff(diff(u(x, t), t), t))+0.6250000000e-10*(diff(diff(phi(x, t), t), t)) = 0

(5)

PDE[a2] := diff(phi(x, t), x)-12*M*Z0*(g-(diff(u(x, t), t, t))-Z0*(diff(phi(x, t), t, t)))/(E*Lb*h^3) = 0

diff(phi(x, t), x)-0.4593750000e-3+0.4687500000e-4*(diff(diff(u(x, t), t), t))+0.4687500000e-7*(diff(diff(phi(x, t), t), t)) = 0

(6)

PDE[a3] := w(x, t) = 0

w(x, t) = 0

(7)

NULL

############################################

``

at x=0 NULL

(8)

PDE[b1] := u(x, t) = 0 

PDE[b2] := w(x, t) = 0

PDE[b3] := diff(phi(x, t), x) = 0

diff(phi(x, t), x) = 0

(9)

################################################

at t=0 for x= [0,La]

u(x, t) = 0

u(x, t) = 0

(10)

w(x, t) = 0

w(x, t) = 0

(11)

phi(x, t) = 0

phi(x, t) = 0

(12)

diff(phi(x, t), t) = 0

diff(phi(x, t), t) = 0

(13)

diff(w(x, t), t) = 0

diff(w(x, t), t) = 0

(14)

diff(phi(x, t), t, t) = 0

diff(diff(phi(x, t), t), t) = 0

(15)

diff(w(x, t), t, t) = 0

diff(diff(w(x, t), t), t) = 0

(16)

######################################################

at t=0 for x= [0,La)

diff(u(x, t), t) = 0

diff(u(x, t), t) = 0

(17)

diff(u(x, t), t, t) = 0

diff(diff(u(x, t), t), t) = 0

(18)

###################################################

at t=0 for x=La

NULL

diff(u(x, t), t) = -ve

diff(u(x, t), t) = -5

(19)

diff(u(x, t), t, t) = g

diff(diff(u(x, t), t), t) = 9.8

(20)

NULL

NULL

 

Download offcenter2.mw

I tried to use a for loop to assign the bounds for creating an RTable in a Maple C wrapper (actually it does not matter what I do here). I checked that the for loop syntax is correct C syntax, for example even as simple as for (int i = 0; i < 2; i++){...} but the loop becomes an infinite loop and takes forever to run. Has anyone encountered the same problem as I did

 Hi everybody,

Until recently I was using Maple 2015.0 on this PC :

Windows XP
2 dual core proc Intel(R) Xeon(R) CPU E3-1225 V2 ; 3.19 GHz
64 bytes,
16 Go RAM

For two days now I am using Maple 2015.0 on this one

Windows 7
4 dual core proc Intel(R) Xeon(R) CPU 5-2637 V3 ; 3.50 GHz
64 bytes,
64 Go RAM


Comparisons of the running times for exactly the same code, consisting in 10000 independent Monte Carlo  simulations distributed over all the nodes (resp 4 and 8) give :

Windows XP :  504 sec  (with a variation of the order of  +/- 3s for different replicates)

Windows 7   :  343 sec  (with a variation of the order of  +/- 3s for different replicates)


The expected running time should be 504 * (3.2/3.5) * (4/8) = 230 s
(or 504 * (4/8) = 252 s if you neglect the acceleration due to the clock rate)
The realized running time (343 s) is thus at least  40% larger than the expected one.


Does anybody has already observed this kind of performance loss during XP -> 7 migration ?
Are there some comparisons of cpu times between Windows XP and 7 ?


Thank you all for sharing your own experience.



postscript : a NOTIONALexample of how the computations are distributed is given below

NbOfRuns := 10000:

# Data is a Matix(NbOfRuns, NbOfCols, …) constructed elsewhere

NbOfNodes := Grid:-NumNodes():
NbOfRunsPerNode := NbOfRuns / NbOfNodes:
for k from 0 to NbOfNodes do
   FirstData := 1 + k * NbOfRunsPerNode:
   LastData := (k+1) * NbOfRunsPerNode:
   Grid:-Run(k+1, MyCode, [Data[FirstData..LastData, ..], …])
end do:

I'm trying to view the source of the webpage fuelsonline.ca

I used HTTP[Get] like I did before which doesn't retrieve much of anything useful.  I thought of the sockets package which I can't seem to pull any information from.  Then I had an idea and viewed the page source which finally has information I can use if I can get Maple to read it in.  If I can get Maple to pull in the contents of the page source I should at least be able to carry it from there. 

Any ideas from anyone?

Hello, I run Maple to solve Binary Integer Programming problem which contain about 1340 constraint and its goal to maximize the objective function.

At first, it's running for 2 hours and said that the iteration limit was reached. So I try to add 'iterationlimit' at LPSolve opts and set it to 10000, but after 3 or 4 hours it said that the iteration limit was reached. So I set 'iterationlimit' to 100000000 and now Maple keep evaluating more than 12 hours.

I run Maple at my notebook with these spesification:

Processor: Intel Corei3-5005U 2.0 GHz

Memory: 4GB RAM

Windows 10

 

It is normal? Or I must run Maple in higher notebook spesification?

Thank you in advance.

 

Below is my Maple file, hope you can help me.

ISL_2017_FASE3.mw

First 137 138 139 140 141 142 143 Last Page 139 of 2097