MaplePrimes Questions

Hi, I'm trying to plot the following functions in 3 dimensions using the following:

plot3d(kappa*sech(kappa*(-k^2*t+x)), x = -10 .. 10, t = -10 .. 10)

plot3d((2*(-K^2+k^2))*(K*cosh(k*(-k^2*t+x))+k*cosh(K*(-K^2*t+x)))/((K+k)^2*cosh(k*(-k^2*t+x)-K*(-K^2*t+x))+(-K+k)^2*cosh(k*(-k^2*t+x)+K*(-K^2*t+x))+4*K*k), x = -10 .. 10, t = -10 .. 10)

plot3d(2*(diff(arctan(lambda*sin(mu*((-3*lambda^2+mu^2)*t+x))*sech(lambda*((-lambda^2+3*mu^2)*t+x))/mu), x)), x = -10 .. 10, t = -10 .. 10);

Each produces an empty plot with a message similar to the following:
Warning, expecting only range variables [x, t] in expression 2*(lambda*cos(mu*((-3*lambda^2+mu^2)*t+x))*sech(lambda*((-lambda^2+3*mu^2)*t+x))-lambda^2/mu*sin(mu*((-3*lambda^2+mu^2)*t+x))*sech(lambda*((-lambda^2+3*mu^2)*t+x))*tanh(lambda*((-lambda^2+3*mu^2)*t+x)))/(lambda^2/mu^2*sin(mu*((-3*lambda^2+mu^2)*t+x))^2*sech(lambda*((-lambda^2+3*mu^2)*t+x))^2+1) to be plotted but found names [lambda, mu]

in each case kappa, k, K, lambda, mu are real constants 

 

Thanks for any help with this
 

 

Is there any facility to apply Finite Volume Method to Partial idifferential equation on MAPLE?
Any comand?

Any Code?

I know abut semilogpot - it gives a semilog scale on the horizontal axis.

Is there any way to get semilog scale on the vertical axis?

Season's greeting

Modify the procedure that implements the secant and tangent routines in such a way that instead of the number of iterations in the beginning give the given accuracy E, with which the approximation is to be determined.

As a result, the procedure should return the last approximation along with the iteration number.

===

===

 

I need help.


How can we remove 0=0 from the above by single comand if it lie in any position from the set?

 

How to eliminate the I in the numerator?

Hello,

Suppose we have a set of quadratic equations of the form:

 

U_11 * a * q + U_12 * b * q + U_13 * c * q + U_14 * d * q + U_15 * a * w + U_16 * b * w + U_17 * c * w + U_18 * d * w + .. = C_i

.

.

.

 

Where terms in uppercase means constant while lowercase letters correspond to unknowns.

 

Now, I want to make a change of variable, so instead of having non-linear term `a * q` we would have `s_ij`, meaning i-th equation and j-th unknown after the change was done.

 

I'm trying to do so, because in my case I will be left with a linear system with a small number of unknowns(relative to N) and N equations, and so I could then solve it easily.

 

Any help is appreciated, thank you!

I have a worksheet that creates many dynamic GIF images.  They use up a huge large amount of RAM, causing MAPLE to grind to a halt.

Instead, on its creation, I would like to export each image file to my external file system, and then delete it from my MAPLE worksheet.  I wish then to use commands in MAPLE to start up and view the image files using IrfanView outside MAPLE.  This would, significantly reduce my RAM usage. (I.e. I wish to export the iGIF mage viewing to IrfranView.)

Can anyone help please?  (I have tried to follow MAPLE documentation, but have been unable find a working solution.)

I am running under Windows 7.

MRB

NULL

 

Typesetting:-delayDotProduct(with, DEtools, true):

with(IntegrationTools):

z := 'z';

z

(1)

whattype(z)

symbol

(2)

``

eq := proc (z) options operator, arrow; evalf(Int(-(2*I)*exp((-1)*.5*I*t)/(sqrt(t^2+1)*exp(sqrt(t^2+1))), t = -500 .. z)) end proc

proc (z) options operator, arrow; evalf(Int(-(2*I)*exp((-1)*.5*I*t)/(sqrt(t^2+1)*exp(sqrt(t^2+1))), t = -500 .. z)) end proc

(3)

f(500)

f(500)

(4)

a := 2*t^2+5*t

2*t^2+5*t

(5)

s := t^2+t+1

t^2+t+1

(6)

d := t^2+t-9

t^2+t-9

(7)

``

 

 

sys := {diff(x(t), t) = eq(z)*y(t)+a*x(t), diff(y(t), t) = s*x(t)+d*y(t)}

{diff(x(t), t) = (Int(-(2.*I)*exp(-(.5*I)*t)/((t^2+1.)^(1/2)*exp((t^2+1.)^(1/2))), t = -500. .. z))*y(t)+(2*t^2+5*t)*x(t), diff(y(t), t) = (t^2+t+1)*x(t)+(t^2+t-9)*y(t)}

(8)

IC_1 := {x(-1) = 0, y(-1) = 1}

{x(-1) = 0, y(-1) = 1}

(9)

dsol3 := dsolve(`union`(sys, IC_1), numeric, parameters = [z], method = rkf45, range = -500 .. 500, output = procedurelist)

"dsol3:=proc(x_rkf45) ... end proc"

(10)

dsol3(parameters)

[z = undefined]

(11)

``

dsol3(parameters = [500])

Error, (in evalf/int) invalid arguments

 

dsol3(500);

Error, (in dsol3) parameters must be initialized before solution can be computed

 

``

``

``

``

``

``

``

``

``

``

``

``

``

``

``

``

``

``

Download eslam_().mw

Hi. could anyone help on this code? fsovle doesn't solve it;
12.mw

Dear Users,

I have difficulty in finding numerical integration of a function f(r,t) which is a function of position r and time t. Function f(r,t) consists 100 terms (for example : BesselJ(0, 151.5793716314014*r)+BesselJ(0, 151.5793716314014*r)*r^2+......100 terms). For a particular time t=t1, f(r,t1) is calculated and then integrated as follows:

I am using evalf(Int(f(r,t1),r=0..1)

Maple takes a lot of time  to evaluate it as it is integrating it in one shot!  Is there a way to

a) pick the terms individually and integrate it

b) then sum these individual terms up together

c) How reliable is evalf(int(f(r,t1),r=0..1)) is? Is evalf (Int()..)  the best way to evaluate integration?

thanks.

Greetings, any one help  me to remove the error  in plots display,here is the codes 

i want 2 rows and 2 columns graphs with 9 lines in each graphs. 

i am attaching the maple questions also  like this i want 

 https://www.mapleprimes.com/questions/224443-How-Do-I-Get-A-Legend-Only-Three-Curves--In-Maple

https://www.mapleprimes.com/questions/225078-How-To-Remove-Error-In-Plot-Display

 

 

restart;
 h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
 K1:=((4/h(z)^4)-(sin(alpha)/Gamma2)-h(z)^2+Nb*h(z)^4):
 lambda:=unapply(Int(K1,z=0..1), Gamma2):
 L1:=0.2:
 
 alpha:=Pi/6:
 with(plots):
 display
  ( Vector[row]
    ([seq [ seq
        ( plot
          ( [ seq
              ( eval(lambda(Gamma2), Nb=j),
                j in [0.1,0.2,0.3]
              )
           d1 in [0.1,0.2,0.5])] ],
            delta2=0.02..0.1,
            legend=[Nb=0.1,Nb=0.2,Nb=0.3],
            labels=[typeset(`δ1`), typeset(conjugate(`Δp`))],
            title=typeset("Effect of ", ''alpha'', " when ", Gamma,"2=", Gamma2)
          ),
          Gamma2 in [10,20,30,40]
        )
      ]
    )
  );

I need Maple code of Bisection method.

I need to know how to compute and manupulate a symbolic equations in vectorial forms in maple. For istance, I need to compute the derivate of the following expression:

r1=T21r2

where ris a vector of three components, representing the position in the reference frame (RF) 1, T21 is the rotation matrix of RF2 w.r.t  RF1 and r2 is the position vector in RF2. In maple I need to derivate (in time) the expression twice to obtain the acceleration.

d(r1)=d(T2)r2 + T21d(r2)

I don't know which command is requested or how should I declare the variables. 

After this I need to substitute the expression:

d(T2) =(T2wo )x

where 'x' is the vectorial product and wo is the angular velocity of thr RF. Of course after the substitution I will need to derivate again the expression.

Is there a way to procede in maple? All this is needed to build a mathematical model of a body in space.

Thank you

EDIT: If the question is not clear, please let me know. This is my first contact with maple and perhaps I am not even understanding what I can do with this software

I have written the following commands in Maple 18.

implicit_func := x^3+y^3 = 9*x*y;
c := 2;
s := evalf(solve(subs(x = c, implicit_func)));
m1 := evalf(eval(implicitdiff(implicit_func, y, x), {x = c, y = s[1]}));
                          0.8000000000
m2 := evalf(eval(implicitdiff(implicit_func, y, x), {x = c, y = s[2]}));
                          -1.257321410
m3 := evalf(eval(implicitdiff(implicit_func, y, x), {x = c, y = s[3]}));
                          0.4573214099

How can I graph all "implicit_func, y - s[1] = m1*(x-c), y - s[2] = m2*(x-c), y - s[3] = m3*(x-c), and the points (c,s[1]), (c,s[2]), and (c,s[3])" in a plane? (Each of them in a different color)     

First 739 740 741 742 743 744 745 Last Page 741 of 2427