MaplePrimes Questions

Hi,

Due to an unexpected maintenance operation, I had to uninstall Maple 18 from my Windows 7-64 bit PC.

Later on I installed it without any problems but, to my surprise now I can't configure it properly. This is:

* Enabling Maple Text as Input (classical input method)
* Removing numbers from equations
* Setting Maple language to English (it took Spanish by default because of Windows)
* Hiding left panel

As usual, this is performed under Tools/Options/Interface and so on...

After I modify my preferred settings, apply globally and close-open again, the program is again in its original form.

How can I do in this case? I have installed it several times. It is also worth noting that the maintenance that I performed was related to deep Windows registry modifications. 

Thanks and regards.

 

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

Would you teach me how to collect and isolate specific terms in an equation?

For example, this is a given equation:

v1(t) - R1*(-i3(t)/n13 - i2(t)/n12) - L1*diff(-i3(t)/n13 - i2(t)/n12, t) = n12*(v2(t) - R2*i2(t) - L2*diff(i2(t), t))

From this equation, I wonder if there is a way to collect and isolate the derivative terms at the RHS, like this:

(- L1/n12+(n12 * L2))*diff(i2(t), t) + (L1/n13)*diff(i3(t), t) = -v1(t) + R1*(-i3(t)/n13 - i2(t)/n12) + n12*v2(t) - n12*R2*i2(t)

Thank you in advance!

 

Hi 

Is there a way to solve it numericall

evalf(int(exp((5+I*6)*x)*sin(1+erf(x)),x=-to..to));

evalf(int(exp((5+I*6)*x)*cos(1+erf(x)),x=-to..to));

 

Would you tell me how to achieve this?

> eqn1 := e1*i1 + e2*i2 + e3*i3 = 0
Then, 

solve( { eqn1 } , { i1 } , explicit ) ;
                     /       e2 i2 + e3 i3\ 
                    { i1 = - ------------- }
                     \            e1      / 
Then, I want to replace 'e2/e1' with this expression:

> eqn2 := n12 = e1/e2;

In other words, I want to express the following:

                       /       e2 i2 + e3 i3\ 
                    { i1 = - ------------- }
                     \            e1      / 

as this:

                      /          -i2     - e3 i3\ 
                    { i1 =  ------   + ------ }
                     \            n12        e1   / 

 

Thank you in advance!

 

Does maple has surrogate optimization built in function or some codes. Thank you

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


 

Dear all,

How can I decide that the product a*n appears on the list [1-a*n*xy, x+y]? I used from the following function but the result is false.

I am looking forward to hearing from you.

Sincerely yours

i have a problem with a nonlinear equation containing Heaviside step function. fsolve does not return anything and SolveEquation command also returns incorrect answer. How can i solve this equation? any method or help is appreciated.


 

 

restart:

with(DirectSearch):

F:=x*(Heaviside(x-(S0))/(1-DD)+Heaviside(-x+S0)/(1-HH*DD));

x*(Heaviside(x-S0)/(1-DD)+Heaviside(-x+S0)/(-DD*HH+1))

(1)

EQ:=F-E*Y;

-E*Y+x*(Heaviside(x-S0)/(1-DD)+Heaviside(-x+S0)/(-DD*HH+1))

(2)

HH:=0.2;
DD:=0.579350262599427;
S0:=1e8;
E:=70e9;
Y:=0.001668993263529;

.2

 

.579350262599427

 

0.1e9

 

0.70e11

 

0.1668993263529e-2

(3)

EQ

-116829528.5+x*(2.377274752*Heaviside(x-0.1e9)+1.131055455*Heaviside(-x+0.1e9))

(4)

fsolve(EQ)

fsolve(-116829528.5+x*(2.377274752*Heaviside(x-0.1e9)+1.131055455*Heaviside(-x+0.1e9)), x)

(5)

SolveEquations(EQ)

Warning, complex or non-numeric value encountered; trying to find a feasible point

 

[1.3868049384289111*10^13, Vector(1, {(1) = -3723983.000}), [x = 9.999999999999999*10^7], 231]

(6)

eval(EQ,x = 1.00000000000000*10^8)

Float(undefined)

(7)

 


 

Download fsolve.mw

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

 

I am new to Maple. I have been plotting 3-dimensional graphs, and would like to plot points on the graph, but I don't see how to do that.

Thanks.

Hi all

for a given two parameter in the attached code ( G and PA) to compute the coefficient Cr

I would like to display the text to describe the corresponding situation ( maybe we one add a table )

 

Many tahnks for your help

 

 

 

testprint.mw

I want to write a code for 3D animation of solutions of PDEs.  Specially, I wonder 3D animation of wave equations etc.

Example: I can write a code for 2d animation but not write 3d animation.

a:=-2:
b:=2:
alpha:=0.1:
beta:=1:
pde:=diff(uu(x,t),t)=alpha*diff(uu(x,t),x,x)+beta*uu(x,t)*(1-uu(x,t));
IC:=uu(x,0)=(sech(10*x))^2;
BC:=uu(a,t)=0,uu(b,t)=0;  

s:=pdsolve(pde,{IC,BC},numeric ,timestep=1/100, spacestep=1/100);
s :- animate( uu(x,t) ,t=0..tbas, frames=40, labels=["x", "u(x,t)"], labelfont=[TIMES,ROMAN,14]); 
g3:=s :- plot3d(uu(x,t), x=-xbas..xbas,t=0..tbas, shading=zhue, axes=boxed, labels=["x","t","u(x,t)"], labelfont=[TIMES,ROMAN,16]);

 But how to write a code for 3d animation?

Thanks.

 

plot([v[Iterations-4], v[Iterations-3], v[Iterations-2], v[Iterations-1], v[Iterations-1]], t = 1 .. 2.9, color = [red, green, blue, black, purple]); g11 := plot([v[Iterations-4], v[Iterations-3], v[Iterations-2], v[Iterations-1], v[Iterations-1]], t = 1 .. 2.9, color = [red, green, blue, black, purple], numpoints = 200);

with(plots); p11 := display({g11}, axes = boxed, thickness = 2); dat1 := `~`[plottools:-getdata]([p11]); for i while i <= 5 do A[i] := dat1[i, 3]; Y[i] := A[i][1 .. -1, 2] end do; X := A[1][1 .. -1, 1]; MM1 := `<|>`(X, Y[j2]$j2 = 1 .. 5);

        [[                                                        
        [[                                                        
dat1 := [["curve", [1. .. 2.9, 1.25843241253069330813854745573 .. 
        [[                                                        

                                    [ 200 x 2 Matrix       ]]  
                                    [ Data Type: sfloat    ]]  
  3.64143345197220537035927195897], [ Storage: rectangular ]], 
                                    [ Order: Fortran_order ]]  

  [                                                              
  [                                                              
  ["curve", [1. .. 2.9, 1.26033290400709038909782882918 .. 6.355\
  [                                                              

                               [ 200 x 2 Matrix       ]]  [
                               [ Data Type: sfloat    ]]  [
  77353302842195250192992449], [ Storage: rectangular ]], [
                               [ Order: Fortran_order ]]  [

  "curve", [1. .. 2.9, 1.26039280181499972068892216922 .. 9.3496\

                              [ 200 x 2 Matrix       ]]  [
                              [ Data Type: sfloat    ]]  [
  4246576526181629314962273], [ Storage: rectangular ]], [
                              [ Order: Fortran_order ]]  [

  "curve", [1. .. 2.9, 1.26039397590711135163918382958 .. 10.112\

                              [ 200 x 2 Matrix       ]]  [
                              [ Data Type: sfloat    ]]  [
  7233044121210212769952674], [ Storage: rectangular ]], [
                              [ Order: Fortran_order ]]  [

  "curve", [1. .. 2.9, 1.26039397590711135163918382958 .. 10.112\

                              [ 200 x 2 Matrix       ]]]
                              [ Data Type: sfloat    ]]]
  7233044121210212769952674], [ Storage: rectangular ]]]
                              [ Order: Fortran_order ]]]
Error, 1st Matrix dimension size (1..1) is too small to fit source size (1..200)
Error, number of indices exceeds rank
ExportMatrix("C:/Users/DELL/Dropbox/Common Data (UH)/Hamid made model/For Hamid Thesis/P4 Picard for ODEs/Tecplot//P5 case 1.dat", MM1)
Error, invalid input: ExportMatrix expects its 2nd argument, M, to be of type {Matrix, list(Matrix)}, but received Vector[row](6, {(1) = X, (2) = Y[1], (3) = Y[2], (4) = Y[3], (5) = Y[4], (6) = Y[5]})
 

First 666 667 668 669 670 671 672 Last Page 668 of 2425