MaplePrimes Questions

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.

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

First 1096 1097 1098 1099 1100 1101 1102 Last Page 1098 of 2434