MaplePrimes Questions

i was wondering anyone had any advice for interpreting conformal plots vs standard plots in the context of symmetry exhibited in the conformal plot like the one attached.

Hi everybody

In the following attached file, I try to evaluate an integration, K_fL. Unfortunately, Maple does not evaluate it and just puts integration symbol and its bounds. I want to have final integration value. Is there any solution to this integration or maybe Maple can not solve this integration because of the complexity of integrand?

Thanks in advance

Q1.mw

Sorry when referencing something written in  explianatory notes in the function advisor / help pages, who do i quote as being  the author?

suppose one has a solution to a system of equations that returns results as RootOf(X) expressions. I would like to know if there is a procedure for extracting the "X" from the RootOf for the purpose of further algebraic analysis.

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

 

i have runed below system of PDES as:



c := 3*10^8;
                           
hbar := 0.105457148e-33;
                             
kB := 0.13806503e-22;
                                
epsilon0 := 0.885418782e-11;
                                                        
timeunit := 1;
                             
g := 1*timeunit;
                            
t0 := 0.10e-4/g;
                           
k := 5/(0.1e10*g);
                                                          
td := 0.10e-4/g;
B := 1;
                            
L := 0.1e-2;
                           
OD := 0.1e7;
              
eta := g*OD/(2*L);
                                 
                       
Omegap0 := .1*g;
      

Omegac(z):=(1e7*B*g)/(sqrt(1+((z-L)/(0.2*L))^(2)));


                 
PDE1 := diff(rho31(t, z), t) = -(1/2)*g*rho31(t, z)+(.5*I)*Omegac(z)*rho21(t, z)+(.5*I)*Omegap(t, z);
                        
PDE2 := diff(rho21(t, z), t) = (.5*I)*Omegac(z)*rho31(t, z);
                                                          
         
PDE3 := diff(Omegap(t, z), t) = -c*(diff(Omegap(t, z), z))+I*c*eta*rho31(t, z);
       

IBC := {Omegap(0, z) = 0, Omegap(t, 0) = Omegap0*exp(-((t-t0)/k)^2), rho21(0, z) = 0, rho31(0, z) = 0};
 
                              
pds := pdsolve({PDE1, PDE2, PDE3}, IBC, numeric, time = t, range = 0 .. 0.1e-2);
                      

 

 

every things is ok but I don't know how can i plot rho31, rho21 and Omegap in 2D or 3D plot.

I have tried

but i see the bellow error:

Error, (in pdsolve/numeric/plot) unable to compute solution for t>INFO["failtime"]:
unable to store -37500000000000000*I when datatype=float[8]

Please help me to solve the error?

Thanks

Hi,

I would like to thank everyone who takes their time to respond to posts on this page. I have another question.

I have about 11 trillion lists that are of the form A shown below. The lists ofcourse have more elements than A (about 50 elements). What I want is to compute this value I call f, and the moment this value is found to be 1, we stop, and if all the values of f are 0, then we print(A). Here is my code, which has a problem and doesn't give me the result I want. How do I modify this so that if any value of f=1, then we stop, elif no 1 is found, we print just one value for A and not for every 0 produced?

with(ListTools);
A := [[1, 2], [1, 7], [5, 6], [1, 6], [1, 9], [6, 5], [9, 1], [2, 1], [7, 1], [6, 1], [5, 6]]; n := 9;
for i to n do for j to n do for k to nops(A) do if [i, j] = A[k] then a := Search([i, j], A); b := Search([j, i], A); if a < b then f := `mod`(b-a, 2); if f = 1 then break else print(A) end if end if end if end do end do end do;

Presently what I am doing is to Append all values of f into an Array and check that 1 is not an element of that Array. This is however very inconvenient since my lists are huge, and they are many, and we do not have to cpompute all the values of f once we come across a 1. Other than this, is there a way I can make my program run faster? Any suggestions are welcomed.

Thanks,

VIC.

How do I define this sequence a(n) in Maple?

a(n) = 1/n for n=odd, a(n) = -1/n^2 for n = even.

Thanks!

hi

please help me for dsolve this equation...

error_length.mw

dsys3 := {(63/20000000)*(diff(u(x), x, x, x, x))-120000000000*(diff(u(x), x, x))-60000000000*(diff(w(x), x, x)), -2.100000000*10^(-30)*(diff(w(x), x, x, x, x, x, x))-4.860000000*10^(-14)*(diff(w(x), x, x, x, x))-(2400*(diff(u(x), x, x)+(diff(w(x), x))*(diff(w(x), x, x))))*(diff(w(x), x))-(1/50000000)*(diff(w(x), x, x))*(120000000000*(diff(u(x), x))+60000000000*(diff(w(x), x))^2-120000000000), u(0) = 0, u(1/2000000) = 0, w(0) = 0, w(1/2000000) = 0, (D(u))(0) = 0, (D(u))(1/2000000) = 0, ((D@@2)(w))(0) = 0, ((D@@2)(w))(1/2000000) = 0, ((D@@4)(w))(0) = 0, ((D@@4)(w))(1/2000000) = 0}

{(63/20000000)*(diff(diff(diff(diff(u(x), x), x), x), x))-120000000000*(diff(diff(u(x), x), x))-60000000000*(diff(diff(w(x), x), x)), -0.2100000000e-29*(diff(diff(diff(diff(diff(diff(w(x), x), x), x), x), x), x))-0.4860000000e-13*(diff(diff(diff(diff(w(x), x), x), x), x))-2400*(diff(diff(u(x), x), x)+(diff(w(x), x))*(diff(diff(w(x), x), x)))*(diff(w(x), x))-(1/50000000)*(diff(diff(w(x), x), x))*(120000000000*(diff(u(x), x))+60000000000*(diff(w(x), x))^2-120000000000), u(0) = 0, u(1/2000000) = 0, w(0) = 0, w(1/2000000) = 0, (D(u))(0) = 0, (D(u))(1/2000000) = 0, ((D@@2)(w))(0) = 0, ((D@@2)(w))(1/2000000) = 0, ((D@@4)(w))(0) = 0, ((D@@4)(w))(1/2000000) = 0}

(1)

dsolve(dsys3, numeric, initmesh = 3024, abserr = 0.1e-4)

`[Length of output exceeds limit of 1000000]`

(2)

``


Download error_length.mw

 

hello!, How can i to plot the next data in a 2D plot?. it is the result of a loop for.... do....while. I had tried with pointplot anda dataplot, but must be something wrong!

with(plots);
for i from 1.73205080756887853 by 0.5e-2 while i < 2.87500000000000000 do roll := i; f1max := NLPSolve(f1, {f2 = i}, x1 = 0 .. 2, x2 = 0 .. 3, method = sqp, maximize = false); c := op([1], f1max), print(c, roll) end do; plot(roll, c);
Warning, limiting number of major iterations has been reached
            2.86923108976435204, 1.73205080756887853
                2.40562374977021154, 1.737050808
                2.36135298525774395, 1.742050808
                2.34703023482192563, 1.747050808
                2.33762023747274306, 1.752050808
                2.33074752039182975, 1.757050808
                2.32549269761476607, 1.762050808
                2.32138461334482216, 1.767050808
                2.31814359109128576, 1.772050808
                2.31558754335936889, 1.777050808
                2.31359008686494683, 1.782050808
                2.31205927196994354, 1.787050808
                2.31092573645557353, 1.792050808
                2.31013563350152307, 1.797050808
                2.30964617400063910, 1.802050808
                2.30942269107110354, 1.807050808
                2.30943663600124838, 1.812050808
                2.30966416809946962, 1.817050808

I have a simple algebraic problem, but Maple can't eliminate the exp(3P) in each term. Please help.

I get the following error:

Error, (in solve) cannot solve for an unknown function with other operations in its arguments


restart

R13eqn := -2*H*Ybar3*Zbar-H*Z1-H1*Z+H1*Zbar+H4*Ybar3-H41 = 0

-2*H*Ybar3*Zbar-H*Z1-H1*Z+H1*Zbar+H4*Ybar3-H41 = 0

(1)

H := exp(3*P)*(Z+Zbar)

exp(3*P)*(Z+Zbar)

(2)

H1 := 3*P1*exp(3*P)*(Z+Zbar)+exp(3*P)*(Z1+Zbar1)

3*P1*exp(3*P)*(Z+Zbar)+exp(3*P)*(Z1+Zbar1)

(3)

H4 := H*(Z4+Zbar4)/(Z+Zbar)

exp(3*P)*(Z4+Zbar4)

(4)

H41 := ((H1*(Z4+Zbar4)+H*(Z41+Zbar41))(Z+Zbar)-H*(Z4+Zbar4)(Z1+Zbar1))/(Z+Zbar)^2

((3*P1(Z+Zbar)*(exp(3*P))(Z+Zbar)*(Z(Z+Zbar)+Zbar(Z+Zbar))+(exp(3*P))(Z+Zbar)*(Z1(Z+Zbar)+Zbar1(Z+Zbar)))*(Z4(Z+Zbar)+Zbar4(Z+Zbar))+(exp(3*P))(Z+Zbar)*(Z(Z+Zbar)+Zbar(Z+Zbar))*(Z41(Z+Zbar)+Zbar41(Z+Zbar))-exp(3*P)*(Z+Zbar)*(Z4(Z1+Zbar1)+Zbar4(Z1+Zbar1)))/(Z+Zbar)^2

(5)

simplify(R13eqn)

(((-3*Z(Z+Zbar)*P1(Z+Zbar)-3*P1(Z+Zbar)*Zbar(Z+Zbar)-Z1(Z+Zbar)-Zbar1(Z+Zbar))*Zbar4(Z+Zbar)+(-3*P1(Z+Zbar)*Z4(Z+Zbar)-Z41(Z+Zbar)-Zbar41(Z+Zbar))*Zbar(Z+Zbar)-3*P1(Z+Zbar)*Z(Z+Zbar)*Z4(Z+Zbar)+(-Z1(Z+Zbar)-Zbar1(Z+Zbar))*Z4(Z+Zbar)-Z(Z+Zbar)*(Z41(Z+Zbar)+Zbar41(Z+Zbar)))*(exp(3*P))(Z+Zbar)+exp(3*P)*(Z+Zbar)*(Zbar4(Z1+Zbar1)+Z4(Z1+Zbar1)+(3*P1-2*Ybar3)*Zbar^3+((3*P1-4*Ybar3)*Z+Zbar1)*Zbar^2+((-3*P1-2*Ybar3)*Z^2-2*Z1*Z+(Z4+Zbar4)*Ybar3)*Zbar-3*P1*Z^3+(-2*Z1-Zbar1)*Z^2+(Z4+Zbar4)*Ybar3*Z))/(Z+Zbar)^2 = 0

(6)

Zbar41 := -2*Zbar*Zbar1

-2*Zbar*Zbar1

(7)

Z41 := -2*Z1*Z

-2*Z1*Z``

(8)

Z4 := -Z^2

-Z^2

(9)

Zbar4 := -Zbar^2

-Zbar^2

(10)

simplify(R13eqn)

((3*P1(Z+Zbar)*Zbar(Z+Zbar)^3+(3*Z(Z+Zbar)*P1(Z+Zbar)+Z1(Z+Zbar)+3*Zbar1(Z+Zbar))*Zbar(Z+Zbar)^2+3*(Z(Z+Zbar)*P1(Z+Zbar)+(2/3)*Z1(Z+Zbar)+(2/3)*Zbar1(Z+Zbar))*Z(Z+Zbar)*Zbar(Z+Zbar)+3*(Z(Z+Zbar)*P1(Z+Zbar)+Z1(Z+Zbar)+(1/3)*Zbar1(Z+Zbar))*Z(Z+Zbar)^2)*(exp(3*P))(Z+Zbar)-exp(3*P)*(Z+Zbar)*(Zbar(Z1+Zbar1)^2+Z(Z1+Zbar1)^2+3*(Z+Zbar)*((P1+(1/3)*Ybar3)*Z^2+((2/3)*Zbar*Ybar3+(2/3)*Z1+(1/3)*Zbar1)*Z-Zbar*((P1-Ybar3)*Zbar+(1/3)*Zbar1))))/(Z+Zbar)^2 = 0

(11)

expand(((3*P1(Z+Zbar)*Zbar(Z+Zbar)^3+(3*Z(Z+Zbar)*P1(Z+Zbar)+Z1(Z+Zbar)+3*Zbar1(Z+Zbar))*Zbar(Z+Zbar)^2+3*(Z(Z+Zbar)*P1(Z+Zbar)+(2/3)*Z1(Z+Zbar)+(2/3)*Zbar1(Z+Zbar))*Z(Z+Zbar)*Zbar(Z+Zbar)+3*(Z(Z+Zbar)*P1(Z+Zbar)+Z1(Z+Zbar)+(1/3)*Zbar1(Z+Zbar))*Z(Z+Zbar)^2)*(exp(3*P))(Z+Zbar)-exp(3*P)*(Z+Zbar)*(Zbar(Z1+Zbar1)^2+Z(Z1+Zbar1)^2+3*(Z+Zbar)*((P1+(1/3)*Ybar3)*Z^2+((2/3)*Zbar*Ybar3+(2/3)*Z1+(1/3)*Zbar1)*Z-Zbar*((P1-Ybar3)*Zbar+(1/3)*Zbar1))))/(Z+Zbar)^2 = 0)

-3*(exp(P))^3*Ybar3*Zbar^4/(Z+Zbar)^2-2*(exp(P))^3*Z^3*Z1/(Z+Zbar)^2-3*(exp(P))^3*P1*Z^4/(Z+Zbar)^2-(exp(P))^3*Z^3*Zbar1/(Z+Zbar)^2+3*(exp(P))^3*P1*Zbar^4/(Z+Zbar)^2+(exp(P))^3*Zbar^3*Zbar1/(Z+Zbar)^2-(exp(P))^3*Z*Zbar(Z1+Zbar1)^2/(Z+Zbar)^2-(exp(P))^3*Z*Z(Z1+Zbar1)^2/(Z+Zbar)^2-(exp(P))^3*Z^4*Ybar3/(Z+Zbar)^2-(exp(P))^3*Zbar*Zbar(Z1+Zbar1)^2/(Z+Zbar)^2-(exp(P))^3*Zbar*Z(Z1+Zbar1)^2/(Z+Zbar)^2+3*P1(Z+Zbar)*(exp(3*P))(Z+Zbar)*Z(Z+Zbar)^3/(Z+Zbar)^2+3*P1(Z+Zbar)*(exp(3*P))(Z+Zbar)*Zbar(Z+Zbar)^3/(Z+Zbar)^2+3*(exp(3*P))(Z+Zbar)*Z(Z+Zbar)^2*Z1(Z+Zbar)/(Z+Zbar)^2+(exp(3*P))(Z+Zbar)*Z(Z+Zbar)^2*Zbar1(Z+Zbar)/(Z+Zbar)^2+(exp(3*P))(Z+Zbar)*Zbar(Z+Zbar)^2*Z1(Z+Zbar)/(Z+Zbar)^2+3*(exp(3*P))(Z+Zbar)*Zbar(Z+Zbar)^2*Zbar1(Z+Zbar)/(Z+Zbar)^2-4*(exp(P))^3*Ybar3*Z^3*Zbar/(Z+Zbar)^2-8*(exp(P))^3*Ybar3*Z^2*Zbar^2/(Z+Zbar)^2-8*(exp(P))^3*Ybar3*Z*Zbar^3/(Z+Zbar)^2-4*(exp(P))^3*Z^2*Z1*Zbar/(Z+Zbar)^2-2*(exp(P))^3*Z*Z1*Zbar^2/(Z+Zbar)^2-6*(exp(P))^3*P1*Z^3*Zbar/(Z+Zbar)^2+6*(exp(P))^3*P1*Z*Zbar^3/(Z+Zbar)^2-(exp(P))^3*Z^2*Zbar*Zbar1/(Z+Zbar)^2+(exp(P))^3*Z*Zbar^2*Zbar1/(Z+Zbar)^2+3*P1(Z+Zbar)*(exp(3*P))(Z+Zbar)*Z(Z+Zbar)^2*Zbar(Z+Zbar)/(Z+Zbar)^2+3*P1(Z+Zbar)*(exp(3*P))(Z+Zbar)*Z(Z+Zbar)*Zbar(Z+Zbar)^2/(Z+Zbar)^2+2*(exp(3*P))(Z+Zbar)*Z(Z+Zbar)*Zbar(Z+Zbar)*Z1(Z+Zbar)/(Z+Zbar)^2+2*(exp(3*P))(Z+Zbar)*Z(Z+Zbar)*Zbar(Z+Zbar)*Zbar1(Z+Zbar)/(Z+Zbar)^2 = 0

(12)

solve(-(exp(P))^3*Z^2*Zbar*Zbar1/(Z+Zbar)^2+(exp(P))^3*Z*Zbar^2*Zbar1/(Z+Zbar)^2-4*(exp(P))^3*Ybar3*Z^3*Zbar/(Z+Zbar)^2-8*(exp(P))^3*Ybar3*Z^2*Zbar^2/(Z+Zbar)^2-8*(exp(P))^3*Ybar3*Z*Zbar^3/(Z+Zbar)^2-4*(exp(P))^3*Z^2*Z1*Zbar/(Z+Zbar)^2-2*(exp(P))^3*Z*Z1*Zbar^2/(Z+Zbar)^2-6*(exp(P))^3*P1*Z^3*Zbar/(Z+Zbar)^2+6*(exp(P))^3*P1*Z*Zbar^3/(Z+Zbar)^2+3*P1(Z+Zbar)*(exp(3*P))(Z+Zbar)*Z(Z+Zbar)^2*Zbar(Z+Zbar)/(Z+Zbar)^2+3*P1(Z+Zbar)*(exp(3*P))(Z+Zbar)*Z(Z+Zbar)*Zbar(Z+Zbar)^2/(Z+Zbar)^2+2*(exp(3*P))(Z+Zbar)*Z(Z+Zbar)*Zbar(Z+Zbar)*Z1(Z+Zbar)/(Z+Zbar)^2+2*(exp(3*P))(Z+Zbar)*Z(Z+Zbar)*Zbar(Z+Zbar)*Zbar1(Z+Zbar)/(Z+Zbar)^2-(exp(P))^3*Z^3*Zbar1/(Z+Zbar)^2+(exp(P))^3*Zbar^3*Zbar1/(Z+Zbar)^2-(exp(P))^3*Z*Zbar(Z1+Zbar1)^2/(Z+Zbar)^2-(exp(P))^3*Z*Z(Z1+Zbar1)^2/(Z+Zbar)^2-(exp(P))^3*Z^4*Ybar3/(Z+Zbar)^2-(exp(P))^3*Zbar*Zbar(Z1+Zbar1)^2/(Z+Zbar)^2-(exp(P))^3*Zbar*Z(Z1+Zbar1)^2/(Z+Zbar)^2+(exp(3*P))(Z+Zbar)*Z(Z+Zbar)^2*Zbar1(Z+Zbar)/(Z+Zbar)^2+(exp(3*P))(Z+Zbar)*Zbar(Z+Zbar)^2*Z1(Z+Zbar)/(Z+Zbar)^2-3*(exp(P))^3*Ybar3*Zbar^4/(Z+Zbar)^2-2*(exp(P))^3*Z^3*Z1/(Z+Zbar)^2-3*(exp(P))^3*P1*Z^4/(Z+Zbar)^2+3*(exp(P))^3*P1*Zbar^4/(Z+Zbar)^2+3*P1(Z+Zbar)*(exp(3*P))(Z+Zbar)*Z(Z+Zbar)^3/(Z+Zbar)^2+3*P1(Z+Zbar)*(exp(3*P))(Z+Zbar)*Zbar(Z+Zbar)^3/(Z+Zbar)^2+3*(exp(3*P))(Z+Zbar)*Z(Z+Zbar)^2*Z1(Z+Zbar)/(Z+Zbar)^2+3*(exp(3*P))(Z+Zbar)*Zbar(Z+Zbar)^2*Zbar1(Z+Zbar)/(Z+Zbar)^2 = 0, P1)

Error, (in solve) cannot solve for an unknown function with other operations in its arguments

 

NULL

``


Download Help_Maple_divide_an_Exp_on_both_sides.mwHelp_Maple_divide_an_Exp_on_both_sides.mw

hello, i would like to know what am i doing wrong, if you can help me, ill be very grateful

ex

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!

First 970 971 972 973 974 975 976 Last Page 972 of 2308