Unanswered Questions

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

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:

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

I have created a model for a robot in Solidworks and have imported it into Maplesim using the CAD toolbox. The problem I have is that the robot has 3 arms that are supposed to come together on a central piece pictured below in figure 1, but attempting to simulate the model with all arms connected with a revolute joint as in figure 2 yields an error that says "The system is underdetermined" the location of the error is main.

For the purposes of the image below I only connected one of the arms, this allows Maplesim to run the file successfully.

figure 1 showing the central piece that the 3 arms are supposed to connect to.

 

Figure 2 showing the problem revolute joints circled in black, the error at the bottom and the setting of the revolute joint on the right.

 

Essentially my question is how do I get the model to work? I apologise if this problem is not terribly well demonstrated, this is my first post onto this forum so I am not sure of all the standards.

First 183 184 185 186 187 188 189 Last Page 185 of 363