Maple 2018 Questions and Posts

These are Posts and Questions associated with the product, Maple 2018

possible to solve following equation with unknown parameter omega.

parameter constant.

I see before for one dimension ode this type equation was solved.

Now for 2d equation is possible?

can consider or I can send again.

Best

2d-2

 

In worksheet mode when me is writing my code, running code, evaluating is and can not work well with maple?

How can inactive  action? Because have to wait many times and is very boring for me.

 

How to study this ellipse with LinearAlgebra without "geometry" eq := -185173378616457/6178315520000*x+86813215770519/24713262080000*(y^2)+126906272070543/24713262080000*(x^2)+256107247454961/6178315520000+(2514994832007/950510080000*x)*y-9123740375967/6178315520000*y = 0 Axis ? foci ? ...Thank you

Based on the equation at here https://www.mapleprimes.com/questions/209660-Problem-With-RealDomainsolve

I tried solve the equation (x-1)*sqrt(x^2 - 4)=0 in Real domain. My code
restart;
RealDomain:-solve((x-1)*sqrt(x^2-4) = 0, x);


I got there solutions are 1, 2, -2.  I think, If we solve the given in RealDomain, we only get two solutions -2 and 2.

My question is: How many solutions are there in the equation (x-1)*sqrt(x^2 - 4)=0 by RealDomain:-solve?

Two pictures by using Mathematica.

How I can calculate integral?

Thanks


 

"restart;    f[1,1](r,theta,phi):=r^4 sin(6 theta) sin(3 phi):    L(r, theta,phi):=(2.784615385 10^10 ((∂)^2)/(∂r^2) `f__11`(r,theta,phi)+(2.784615385 10^10 (2+2 r cos(theta)) ((∂)/(∂r) `f__11`(r,theta,phi)))/(r (2+r cos(theta)))-(0.1175000000 (((∂)^4)/(∂theta^4) `f__11`(r,theta,phi)))/(r^4)-(0.1175000000 (((∂)^4)/(∂phi^4) `f__11`(r,theta,phi)))/((2+r cos(theta))^4)-(0.1175000000 (((∂)^4)/(∂phi^2∂r^2) `f__11`(r,theta,phi)))/((2+r cos(theta))^2)-(0.1175000000 (((∂)^4)/(∂r^2∂theta^2) `f__11`(r,theta,phi)))/(r^2)-(0.2350000000 (((∂)^4)/(∂phi^2∂theta^2) `f__11`(r,theta,phi)))/(r^2 (2+r cos(theta))^2)+(0.1175000000 ((cos(theta))^2 r^2+4 (cos(theta))^2 r^4+16 cos(theta) r^3-4+17 r^2) (((∂)^2)/(∂theta^2) `f__11`(r,theta,phi)))/((2+r cos(theta))^2 r^4)+(0.1175000000 (2 (cos(theta))^2 r^2+4 (cos(theta))^2 r^4+16 cos(theta) r^3+4+12 r^2) (((∂)^2)/(∂phi^2) `f__11`(r,theta,phi)))/(r^2 (2+r cos(theta))^4)-(2.784615385 10^10 (2 (cos(theta))^2 r^2+4 r cos(theta)+4) `f__11`(r,theta,phi))/(r^2 (2+r cos(theta))^2)+(0.2350000000 (((∂)^3)/(∂r∂theta^2) `f__11`(r,theta,phi)))/(r^3 (2+r cos(theta)))-(0.2350000000 (((∂)^3)/(∂phi^2∂r) `f__11`(r,theta,phi)))/(r (2+r cos(theta))^3)+(0.2350000000 sin(theta) (((∂)^3)/(∂theta^3) `f__11`(r,theta,phi)))/(r^3 (2+r cos(theta)))-(0.2350000000 sin(theta) (((∂)^3)/(∂phi^2∂theta) `f__11`(r,theta,phi)))/(r (2+r cos(theta))^3)+(0.1175000000 sin(theta) (((∂)^3)/(∂r^2∂theta) `f__11`(r,theta,phi)))/(r (2+r cos(theta)))-(0.1175000000 (2 r cos(theta)+3) sin(theta) (2 r cos(theta)+5) ((∂)/(∂theta) `f__11`(r,theta,phi)))/(r (2+r cos(theta))^3)) r^4 sin(6 theta) sin(3 phi):"

with(Student[Calculus1]); K[rr, s] := evalf(ApproximateInt(L(r, theta, phi), r = .2 .. 1, method = simpson)); KK[rr, s] := evalf(ApproximateInt(K[rr, s], theta = 0 .. 2*Pi, method = simpson)); k2 := evalf(ApproximateInt(KK[rr, s], phi = 0 .. 2*Pi, method = simpson))

``


 

Download

X := [0, -3]; Y := [3, 1]; Z := [5, -2]; PP := proc (M::list, N::list, K::list) local A, B, P, distPH, H, lambda, mdist; A := M; B := N; P := K; H := simplify(expand((1-lambda)*A+lambda*B)); distPH := sqrt((H[1]-P[1])^2+(H[2]-P[2])^2); mdist := diff(distPH, lambda); lambda := solve(mdist, lambda); simplify(distPH); H end proc; debug(PP); PP(X, Y, Z);#the calculation isn't finished

Hello,

How I can separate three functions (u__ru__theta and u__phi) in the equation attached below.

Please also see attached figure. I want to calculate L11L22L33.

Thanks.

SEPARATE
 

restart

B := simplify((r*(R+r*cos(theta))^2*(mu+lambda)*(diff(`u__θ`(r, theta, phi), r, theta))+2*r^2*(R+r*cos(theta))^2*(mu+(1/2)*lambda)*(diff(u__r(r, theta, phi), r, r))+r^2*(mu+lambda)*(R+r*cos(theta))*(diff(`u__φ`(r, theta, phi), phi, r))+mu*(R+r*cos(theta))^2*(diff(u__r(r, theta, phi), theta, theta))+(diff(u__r(r, theta, phi), phi, phi))*mu*r^2-3*(R+r*cos(theta))^2*(mu+(1/3)*lambda)*(diff(`u__θ`(r, theta, phi), theta))+(2*(R+2*r*cos(theta)))*r*(R+r*cos(theta))*(mu+(1/2)*lambda)*(diff(u__r(r, theta, phi), r))-r^2*sin(theta)*(mu+lambda)*(R+r*cos(theta))*(diff(`u__θ`(r, theta, phi), r))-3*r^2*(mu+(1/3)*lambda)*cos(theta)*(diff(`u__φ`(r, theta, phi), phi))-r*mu*sin(theta)*(R+r*cos(theta))*(diff(u__r(r, theta, phi), theta))-(2*(2*cos(theta)^2*r^2+2*cos(theta)*R*r+R^2))*(mu+(1/2)*lambda)*u__r(r, theta, phi)+r*`u__θ`(r, theta, phi)*sin(theta)*(3*r*(mu+(1/3)*lambda)*cos(theta)+R*mu))/(r^2*(R+r*cos(theta))^2))

``

``

 


 

Download SEPARATE


 

Hello, I bring here a problem with maple that in result the  "".

I can remove '' in final results?

Thanks


 


  ;
  ):


  X= (R+r*costheta)*sinphi):
  = (R+r*costheta)*cosphi):
  Z= (r*sintheta)):

  
  r, theta, phi,
    X, Y, Z
  ):

 


   
  
    (1/2*Curl
      
        < u__rrtheta, phi),
            u__thetartheta, phi),
            u__phirtheta, phi)
          >,
          rtheta, phi
        )
      )
    )
  );

Vector(3, {(1) = (1/2)*csgn(r)*csgn(R+r*cos(theta))*(csgn(1, R+r*cos(theta))*(R+r*cos(theta))*`#msub(mi("u"),mi("&phi;",fontstyle = "normal"))`(r, theta, phi)+((R+r*cos(theta))*(diff(`#msub(mi("u"),mi("&phi;",fontstyle = "normal"))`(r, theta, phi), theta))-r*sin(theta)*`#msub(mi("u"),mi("&phi;",fontstyle = "normal"))`(r, theta, phi))*csgn(R+r*cos(theta))-csgn(r)*r*(diff(`#msub(mi("u"),mi("&theta;",fontstyle = "normal"))`(r, theta, phi), phi)))/(r*(R+r*cos(theta))), (2) = -(csgn(1, R+r*cos(theta))*(R+r*cos(theta))*`#msub(mi("u"),mi("&phi;",fontstyle = "normal"))`(r, theta, phi)+((R+r*cos(theta))*(diff(`#msub(mi("u"),mi("&phi;",fontstyle = "normal"))`(r, theta, phi), r))+cos(theta)*`#msub(mi("u"),mi("&phi;",fontstyle = "normal"))`(r, theta, phi))*csgn(R+r*cos(theta))-(diff(`#msub(mi("u"),mi("r"))`(r, theta, phi), phi)))*csgn(R+r*cos(theta))/(2*r*cos(theta)+2*R), (3) = (1/2)*(-(diff(`#msub(mi("u"),mi("r"))`(r, theta, phi), theta))+csgn(r)*r*(diff(`#msub(mi("u"),mi("&theta;",fontstyle = "normal"))`(r, theta, phi), r))+`#msub(mi("u"),mi("&theta;",fontstyle = "normal"))`(r, theta, phi)*(r*csgn(1, r)+csgn(r)))*csgn(r)/r})

(1)

``


 

Download

 

maple does not work at all

it displays this error

Error, (in StringTools:-FormatMessage) unknown option MAPLE
 

car_2som_opp := proc (U::list, V::list)  #construction d'un carré connaissant 2 sommets opposés 
local dist, eqCerU, eqCerV, r, sol, X, Y; 
dist := proc (M, N) sqrt(add((M[i]-N[i])^2, i = 1 .. 2)) end proc;
r := dist(U, V)/sqrt(2); 
eqCerU := (x-U[1])^2+(y-U[2])^2 = r^2; 
eqCerV := (x-V[1])^2+(y-V[2])^2 = r^2;
sol := solve([eqCerU, eqCerV], [x, y],allsolutions,explicit);  
map(allvalues,sol): 
X := [subs(op(sol[1]), x), subs(op(sol[1]), y)]; 
Y := [subs(op(sol[2]), x), subs(op(sol[2]), y)]; 
display(plot([U, X, V, Y, U],scaling = constrained, axes = none)) 
end proc:

car_2som_opp([-5,6],[7,-3]);"error
 

Need help solving this problem with a maple proc using the Crank–Nicolson method for the differential part and any other quadrature  for the integral part and thank you so much in advance any ideas or thoughts would be helpful

How I can convert attached maple code into series form to the latex format.

I want to convert to this format series form) not for certain M and N.

 

 

Hi. what is reason maple unable to integral in answer of dsolve? when I try to use dsolve for solve my equation in answer of maple there are expressions of integral that isnot calculated

&int;(-625 R^2 ro (-2 cp3 x1 lambdaopt+sin((pi t)/10) cp2+(16 cp2)/5) (-cos((pi t)/10)+cos(2 pi)+((-t/5+4) x1+(8 t)/25-32/5) pi) lambdaopt pi b11 (e)^(-((16+5 sin((pi t)/10)) cp1)/(10 lambdaopt x1))+625 (-R^2 k11 (cos((pi t)/10))^3+k11 R^2 (cos(2 pi)-((t-20) (x1-8/5) pi)/5) (cos((pi t)/10))^2+((32 sin((pi t)/10) R^2 k11)/5+(281 R^2 k11)/25+4 lambdaopt^2 (a11 x1+a13 x3)) cos((pi t)/10)-(32 k11 R^2 (cos(2 pi)-((t-20) (x1-8/5) pi)/5) sin((pi t)/10))/5+(281 (x1-8/5) (R^2 k11+(100 lambdaopt^2 (a11 x1+a13 x3))/281) pi t)/125) x1^2)&DifferentialD;t

Hi,

I'm trying to solve an expression with only one unknown variable, but for some reason solve and fsolve are unable to return a solution.

The expression I'm trying to solve is:

p := k*T*m__hhw*(ln(1+exp(E__fv-E__hh0)/(k*T))+ln(1+exp(E__fv-E__hh1)/(k*T)))/(Pi*h__bar^2*L__z) = 0.3e25

where all variables are predefined and am trying to solve for E__fv. However, when I use solve, I get the error: "Warning, solutions may have been lost", and when I try to use fsolve, it simply returns the expression as an answer and I am unable to find the numerical value for E__fv. Any helps or tips are appreciated.

If it helps, the defined variable values are:

k = 1.3806E-23;

T = 300;

h__bar = 1.05456E-34

L__z = 6E-9

m__hhw = 3.862216E-31;

E__hh0 = 3.012136E-21

E__hh1 = 1.185628E-20

 

A similar expression in the previous line was able to solve correctly and return a numerical value, so I'm not sure why solve/fsolve can't solve this one.

Similar solved expression:

solve(0.3e25 = m__cw*k*T*ln(1+exp(E__fc-E__c0)/(k*T))/(Pi*h__bar^2*L__z), E__fc);
                          -42.88488490
 

I have this polynomial equation: (x-2)^2*(x-3)+epsilon =0, I want to draw a bifurcation diagram in the (epsilon , x) plane.

 

How to implement this in maple 2018?

 

Thanks!

 

5 6 7 8 9 10 11 Last Page 7 of 44