Maple 15 Questions and Posts

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

How can I plot a vertical line say x = a in Maple 15 worksheet?

Thanks

Maple 15 does not allow me to do even small changes after executing even the smallest bit of command. Is there a way to overwrite this default setting and force Maple to undo? Do the newer versions of Maple provide this functionality and what is the purpose of undo getting deactivated after a command, at least it should be enabled for the input?

Greetings,

Maple 15 allows the following syntax

omega := sqrt(w0^2*(1+((z-zf)/z0)^2))

But it does not allow this one 

omega^2 := sqrt(w0^2*(1+((z-zf)/z0)^2))

Why is this so for functions and variables? Is there any way around this, I am really bugged by this issue. The Maple Math told me that the syntax was invalid so I had type the commands directly.

 

Is it possible in Maple 15 to solve an equation with a parameter for a given set of parameters? How can this be passed to the solve function, should I use some kind of list?  After obtaining the solution how can I assign the solutions to variables such as x1 for the first value of the parameter, x2 for the second value of the parameters and so on. Furthermore, is this possible with the fsolve command?

 

Thanks

Hello I am a Maple 15 user and I am using the command fsolve to solve for the intersection of two curves over a specified interval in x, namely from 0 to the lim defined in the Maple document. The specified interval contains asymptotes and when I specify the full interval only one of the three solutions is returned even if I can see that there are three distinct solutions by looking at the plot of RHS and LHS. Should I use another technique to find the solution or is my implementation of fsolve command wrong?

Thanks in advance


restart

with(ListTools):

n1 := 1:

n2 := 1.50:

n3 := 1.40:

lambda := 1.3:

k0 := 2*Pi/lambda:

d := 3:

x0 := k0*d:

arg1 := sqrt(x0^2*(n2^2-n1^2)):

arg2 := sqrt(x0^2*(n2^2-n3^2)):

lim := FindMinimalElement([arg1, arg2]):

sqr1 := sqrt(x0^2*(n2^2-n1^2)-x^2):

sqr2 := sqrt(x0^2*(n2^2-n3^2)-x^2):

LHS := tan(x):

RHS := (sqr1+sqr2)/(x*(1-sqr1*sqr2/x^2)):

plot([LHS, RHS], x = 0 .. lim, y = -6 .. 6)

 

fsolve(RHS = LHS, x = (1/2)*Pi .. 3*Pi*(1/2))

2.634254816

(1)

fsolve(RHS = LHS, x = 3*Pi*(1/2) .. 9*Pi*(1/4))

5.222527128

(2)

fsolve(RHS = LHS, x = 9*Pi*(1/4) .. lim)

7.598486053

(3)

``


Download HW4Q2.mw

test.mw

Hello. As you can see the minimum of set X is 2.

the time pair in set W whose difference is 2 corresponds to elements 54 and 55. How do i locate and print these automatically?

Digits:=4:

U := [1.010, 1.110, 1.210, 1.310, 1.410, 1.510, 1.610, 1.710, 1.810, 1.910, 2.020, 2.120, 2.220, 2.320, 2.420, 2.520, 2.620, 2.720, 2.820, 2.920, 3.030, 3.130, 3.230, 3.330, 3.430, 3.530, 3.630, 3.730, 3.830, 3.930, 4.040, 4.140, 4.240, 4.340, 4.440, 4.540, 4.640, 4.740, 4.840, 4.940, 5.050, 5.150, 5.250, 5.350, 5.450, 5.550, 5.650, 5.750, 5.850, 5.950, 6.060, 6.160, 6.260, 6.360, 6.460, 6.560, 6.660, 6.760, 6.860, 6.960, 7.070, 7.170, 7.270, 7.370, 7.470, 7.570, 7.670, 7.770, 7.870, 7.970, 8.080, 8.180, 8.280, 8.380, 8.480, 8.580, 8.680, 8.780, 8.880, 8.980, 9.090, 9.190, 9.290, 9.390, 9.490, 9.590, 9.690, 9.790, 9.890, 9.990, 10.01, 11.11, 12.21]:

W := select(proc (t) options operator, arrow; .59 >= frac(t) end proc, U);

[1.010, 1.110, 1.210, 1.310, 1.410, 1.510, 2.020, 2.120, 2.220, 2.320, 2.420, 2.520, 3.030, 3.130, 3.230, 3.330, 3.430, 3.530, 4.040, 4.140, 4.240, 4.340, 4.440, 4.540, 5.050, 5.150, 5.250, 5.350, 5.450, 5.550, 6.060, 6.160, 6.260, 6.360, 6.460, 6.560, 7.070, 7.170, 7.270, 7.370, 7.470, 7.570, 8.080, 8.180, 8.280, 8.380, 8.480, 8.580, 9.090, 9.190, 9.290, 9.390, 9.490, 9.590, 10.01, 11.11, 12.21]

(1)

X:=[seq(trunc(W[i+1])*60+frac(W[i+1])*100-(trunc(W[i])*60+frac(W[i])*100),i=1..nops(W)-1)];

[10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 2.00, 70.00, 70.00]

(2)

NULL

min(X),{W[54],W[55]};

2.00, {9.590, 10.01}

(3)

 



Download test.mw

  Given this equation:       (chi^2*(xi^2-1)^2+1-Omega^2)*(xi^2-Omega^2)-xi^2 = 0

   chi^2=1/1200;  I want to plot xi vs Omega. Can anyone help me out?

Hi there,

Is there a way in which i can solve the following optimal control problem numerically with Maple?

 dH/dt=λ-µ H-(1-u1)β H V+δ I,

dI/dt=(1-u1)β H V-σ I,

‎dV/dt=(1-u2)k I-γ‎V,

dλ1/dt=-1+λ1µ+β‎‎V(1-u1)(λ1-λ2),‎

dλ2/dt=λ1δ+λ2σ-λ3(1-u2)k‎,‎

‎dλ3/dt=β H(λ1-λ2)(1-u1)‎‎+λ3γ.

where

u1=βHV(λ1-λ2)/A1,

u2=-λ3kI/A2,

σ = α + δ,

and

H(0)=1.7*10^8,    I(0)=0,    V(0)=400.

 

λ=5*10^5,    µ=0.003,     β=4*10^(-10),    α=0.043,    δ= 0.2                  k=6.24,       γ‎=0.65.

A1=900, A2=1000.

 

Answers and advice are very appreciated. 

Thank you all for reading.

Aylin

Is it possible to evaluate a function at multiple points described by an array or something of that sort and have Maple return the evaluations as an array. I need approximations of a function at various values of its argument so it would be nice to do it with a single command.

Thanks

How can I find stepwise solution of the following integral?

Int(exp(-xi^2*b^2)*cos(xi^2*a*t)*cos(xi*x), xi = 0 .. infinity)

where xi greek letter.. 

I am Maple 15 student edition user and the packages is not loaded when I use the execute all command from the GUI. I have to manually press the enter key to have by package loaded before I use a function defined in the package. How can I work around this issue?

Thanks

P.S: You can see the Maple document attached to this post.

NumericalSimulation.mw

In Maple 15 is it possible to differentiate with respect to say x - y instead of x

 

Thanks

Greetings,

I am a Maple 15 Student Edition user on both Windows and Linux. I wonder if it is possible to see the actual command send to Maple when I click on a particular entry in the right click context menu. If it is possible to do so, how can I have it displayed in document/worksheet mode? Does it matter if I am inputting in 1D or 2D?

Thanks in advance

How can I plot a volume (many surfaces) in X,Y, Z axis where X,Y,Z are functions in 4 variables (a,b,c,d), and the domain for the 4 variables are 

-90<=a>=90, -10<=b>=10, -12<=c>=12, -90<=d>=0,

 

X := proc (a, b, c, d) options operator, arrow; 324.*cos(b)*sin(c)*cos(d)+324.*sin(b)*sin(d)+323.5*cos(b)*sin(c) end procX := proc (a, b, c, d) options operator, arrow; 324.*cos(b)*sin(c)*cos(d)+324.*sin(b)*sin(d)+323.5*cos(b)*sin(c) end proc

Y := proc (a, b, c, d) options operator, arrow; (324*1.*sin(a)*sin(b)*sin(c)+324*1.*cos(a)*cos(c))*cos(d)+(-1)*324.*sin(a)*cos(b)*sin(d)+323.5*sin(a)*sin(b)*sin(c)+323.5*cos(a)*cos(c)-100 end proc

Y := proc (a, b, c, d) options operator, arrow; (324*1.*sin(a)*sin(b)*sin(c)+324*1.*cos(a)*cos(c))*cos(d)+(-1)*324.*sin(a)*cos(b)*sin(d)+323.5*sin(a)*sin(b)*sin(c)+323.5*cos(a)*cos(c)-100 end proc

Z := proc (a, b, c, d) options operator, arrow; (324*cos(a)*sin(b)*sin(c)-324*sin(a)*cos(c))*cos(d)-324*cos(a)*cos(b)*sin(d)+323.5*cos(a)*sin(b)*sin(c)+(-1)*323.5*sin(a)*cos(c)+150 end proc

Z := proc (a, b, c, d) options operator, arrow; (324*cos(a)*sin(b)*sin(c)-324*sin(a)*cos(c))*cos(d)-324*cos(a)*cos(b)*sin(d)+323.5*cos(a)*sin(b)*sin(c)+(-1)*323.5*sin(a)*cos(c)+150 end proc

According to this site,"It is known that every even number can be written as a sum of at most six primes". 

http://www.theage.com.au/national/education/christians-goldbachs-magic-sum-20140903-3es2t.html

i wanted to test this using maple.

restart:
> PF := proc (a::integer)

> local cst,obj,res;
> cst := add(x[i], i = 1 .. numtheory:-pi(prevprime(a))) <= 6;
> obj := add(x[i]*ithprime(i), i = 1 .. numtheory:-pi(prevprime(a)))-a;
> res := Optimization:-LPSolve(obj, {cst ,obj>=0}, assume={nonnegative,integer}); end proc:
> PF(30);
[0, [x[1] = 0, x[2] = 0, x[3] = 6, x[4] = 0, x[5] = 0, x[6] = 0,x[7] = 0, x[8] = 0, x[9] = 0, x[10] = 0]]

the third prime is 5 and 6 of them make 30. as an aside, it would be nice to know how to get maple to output "30 = 6x5".

this is obviously pretty limited, because 30 can be written as the sum of two primes (7+23 and 11+19) [GOLDBACH], but using DS's GlobalSearch for all solutions takes a long time to compute. also I have to nominate the highest prime.

any suggestions?

First 11 12 13 14 15 16 17 Last Page 13 of 47