Unanswered Questions

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

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.???

Dear forum users/ admins, I have some questions regarding the wind turbine vibrations model based on the MapleSim gallery. I am a new user, went through some tutorials and am currently trying to understand how the system works.

 

I will try to deliver my questions based on the pictures shown below:

 

 

 

 1. After running the model, 2 results were shown. One is "Latest Results" and the other is "zeta=0.01". Based on the images above, there is a comment that the result was dated Feb 14 2013. Is this the result on which the day the model was created? If it is, is it programmed to be in the result file everytime the simulation runs?

 

 

2. If I untick "zeta=0.01", a red graph appears, which says Latest Result. But the graphs are not labelled. How can I label this graph in MapleSim? 

3. I would also like to ask, what type of vibrations are occuring in the model, and what does the Deflection-Time graph means. Why does the deflection occur highest at 120s and then drops? Is it because the tower becomes stabilized after a while?

 

 

4. I understand that the model is run by a signal and fed into a speed component. May I know what is the f(critical) is, and how it affects the system?

 

 

5. From the Ramp signal, there is a "height" value. The pre-set value is 10. After changing this value, to lets say 50, the maximum deflections occur at a an earlier time as shown in the figure below:

 

 

6. I would like to understand how does the ramp affect the maximum deflection time that occurs.

 

Thank you very much!

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

let aa = map (+1.2) [1,2]
let bb = map (+1.4) [2,3]
foldr (aa . bb) [3.0,4.0]

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:

First 185 186 187 188 189 190 191 Last Page 187 of 365