Kitonum

18914 Reputation

26 Badges

14 years, 37 days

MaplePrimes Activity


These are replies submitted by Kitonum

@Ahmed111 

D[2](q)  is the derivative with respect to the second variable of the function q .

Using the  $  operator, we can easily prevent the animation from coming so quickly (in the first example of tomleslie's answer, each frame is multiplied 7 times):

restart:
  with(Fractals:-LSystem):
  with(LSystemExamples):
  with(plots):
  display([seq(PlotExample(DragonCurve,j)$7, j=1..12)], insequence=true);

@Will_iii  I don't understand why you are typing  b:=solve(eq,b) = 1. Enough

eq := 5*a^2 + 2*b - 7 = 0;   
a:=1;                                            
b:=solve(eq,b);

Maple now knows that  b  is  1 .

This is probably a problem with your version of Maple. In Maple 2018.2 everything works as expected:

restart:

with(plots):

equ1 := BesselJ(sqrt(17)/2, 10*sqrt(t)*sqrt(2))/t^(1/4) + BesselY(sqrt(17)/2, 10*sqrt(t)*sqrt(2))/t^(1/4):

equ2 := BesselJ(sqrt(17)/2, 10*sqrt(t)*sqrt(2))/t^(1/4) + 5*BesselY(sqrt(17)/2, 10*sqrt(t)*sqrt(2))/t^(1/4):

equ3 := BesselJ(sqrt(17)/2, 10*sqrt(t))/t^(1/4) + 5*BesselY(sqrt(17)/2, 10*sqrt(t))/t^(1/4):

tmax   := 30:
colors := ["Red", "Violet", "Blue"]:

p1 := plot([equ1, equ2, equ3], t = 0 .. tmax, labels = [t, T[2](t)], tickmarks = [0, 0], labelfont = [TIMES, ITALIC, 12], axes = boxed, color = colors):

ymin := min(op~(1, op~(2, op~(2, [plottools:-getdata(p1)])))):
ymax := max(op~(2, op~(2, op~(2, [plottools:-getdata(p1)])))):
dy   := 2*ymax:

legend1 := typeset(C[3] = 1, ` , `, C[4] = 1, ` , `, Omega^2 = 50):
legend2 := typeset(C[3] = 1, ` , `, C[4] = 5, ` , `, Omega^2 = 50):
legend3 := typeset(C[3] = 1, ` , `, C[4] = 5, ` , `, Omega^2 = 25):

p2 := seq(textplot([tmax-2, ymax-k*dy/20, legend||k], align=left), k=1..3):

p3 := seq(plot([[tmax-2, ymax-k*dy/20], [tmax-1, ymax-k*dy/20]], color=colors[k]), k=1..3):
display(p1, p2, p3, view=[default, -ymax..ymax], size=[800, 500])

 

 

Download Legend_Inside_new.mw

@zenterix  I understood your question how to write the same in Maple syntax. It seems that MaplePrimes editor does not support writing formulas in the Latex language. But there are resources where you can do this, for example http://primat.org/mathred/mathred.html  . And here you can paste the corresponding screenshot.

@Will_iii  I'm wondering why Maple help is missing in Maple Flow, in which Maple syntax is used.

@Lali_miani
 
The first example:

f:=n->2*(n-1)/(2*n-1):
seq(f(n), n=1..10);


and so on.

@Carl Love  Thank you. I fixed the output in my code. When I posted it yesterday, I corrected something but forgot to replace the output. Now it's all right.

@south The  eval  command does roughly the same as the  subs  one. You can replace  eval(f(x), %)  by the subs(%, f(x))  in the code below .   

@etian2  You forgot to call the Physics[Vectors] subpackage:

restart;
with(Physics[Vectors]):
a_-b_;

                                           

 

@jalal  Should be:

A := animate(plot, [[T(a), [[a, f(a)]]], x = -20 .. a, style = [line, point], color = red, thickness = 3, symbolsize = 16], a = -20 .. 20, frames = 90, background = G, view = -20 .. 20):
 
B := animate(plot, [[g(a), [[a, f(a)]]], x = -20 .. a, style = [line, point], color = red, thickness = 3, symbolsize = 16], a = -20 .. 20, frames = 90, background = G, view = -20 .. 20):
 
display(A, B);

 

@acer  1. I did not set a goal to calculate this integral of a complex-valued function over the entire interval,, but simply wanted to explain to OP why Maple does not directly calculate this integral. I also wanted to show that for this calculation it is not necessary to call the VectorCalculus package, but we can simply use the int  command (or evalf(Int) for numerical calculation).

2. As for the integrand, I confused path integral with line integral (omitting the arc differential), as already pointed out by Carl.

@acer Thanks for the helpful critique. Below is the corrected version with your example:

restart:
CartProd:=proc(L)
local n, l, p;
option remember;
n:=nargs;
if n=1 then return args else
if n=2 then return [seq(seq([p,l], l=args[2]), p=args[1])] else
[seq(seq([op(p),l], l=args[n]), p=thisproc(args[1..n-1]))] fi; fi;
end proc:

LL:=seq([seq(x[i,j],j=1..3)],i=1..2):
CartProd(LL);

               

 

@mmcdara  Yes, sure. It's very easy to fix. 

1 2 3 4 5 6 7 Last Page 1 of 124