Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Help page says "minus~" is an element-wise operators in Maple
Why this command doesn't return [{1}, {3}] but generates an error ?
[{1, 2}, {2, 3}] minus~ {2};
Error, dimension bounds must be the same for all container objects in an elementwise operation

OK, you can answer I could use  map(u -> u minus {2}, [{1, 2}, {2, 3}]), but it wouldn't explain the reason of the (my?) error.

Thanks in advance

 

 

 

duaxisplot([x=1,x=2,x=3], ....);

just want to plot multiple vertical line with a time series at the same time

How could one tell Maple to simplify  log(x)+log(y) to log(x*y)?

As Maple own apps says  https://www.maplesoft.com/support/help/maple/view.aspx?path=MathApps%2FProductRuleForLogarithms

And known identity  https://en.wikipedia.org/wiki/List_of_logarithmic_identities

 

This is what I tried

restart;
simplify( ln(x) + ln(y), symbolic);
simplify( ln(x) + ln(y), ln);
simplify( ln(x) + ln(y), size);

Maple knows this identity, becuase

simplify( ln(x*y)- (ln(x)+ln(y)) , symbolic)
       0

Note, this identity has no assumptions on it. It is defined as is.

Maple help page says

ln(x y) ==> ln(x) + ln(y) provided     0 < x  and   signum(y) is unknown

so, is Wikipedia and all the other references I saw wrong for not showing this assumption on x>0 whey they state the product rule for logs?

I want to go the other way

 ln(x) + ln(y) ==> ln(x y) 

So this releation should really be written as

ln(x) + ln(y) <==> ln(x y) 

 

How could I substitute the (t1*e) into t2 at the index of e?  I am trying many ways to do that, but it doesn't work. please help me.
 

restart:

alias(epsilon = e, omega = w, omega[0] = w0, t[1] = t1, t[2] = t2, alpha[1] = a1, alpha[2] = a2, alpha[3] = a3, beta[1] = b1, beta[2] = b2, F[0] = f0); e := proc (t1, t2) options operator, arrow; e end proc; w0 := proc (t1, t2) options operator, arrow; w0 end proc; a := proc (t1, t2) options operator, arrow; a end proc; f := proc (t1, t2) options operator, arrow; f end proc; mu := proc (t1, t2) options operator, arrow; mu end proc

 

 

(1/8)*alpha[2](t[1], t[2])*b^3*exp((3*I)*omega*t[1])+(1/8)*alpha[2](t[1], t[2])*b^3*exp(-(3*I)*omega*t[1])+(1/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp((3*I)*t[1])+(1/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp(-(3*I)*t[1])+3*alpha[2](t[1], t[2])*b^3*exp(-I*omega*t[1])*(1/8)+3*alpha[2](t[1], t[2])*b^3*exp(I*omega*t[1])*(1/8)-I*alpha[1](t[1], t[2])*A(t[2])*exp(-I*t[1])*(1/2)+3*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(-I*omega*t[1])*(1/4)-I*(diff(A(t[2]), t[2]))*exp(-I*t[1])+3*alpha[2](t[1], t[2])*A(t[2])*exp(-I*t[1])*b^2*(1/4)-I*alpha[1](t[1], t[2])*b*omega*exp(-I*omega*t[1])*(1/2)+I*alpha[1](t[1], t[2])*b*omega*exp(I*omega*t[1])*(1/2)+3*alpha[2](t[1], t[2])*A(t[2])^3*exp(-I*t[1])*(1/8)+I*(diff(A(t[2]), t[2]))*exp(I*t[1])+3*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(-I*t[1]-(2*I)*omega*t[1])*(1/8)+3*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(-I*t[1]+(2*I)*omega*t[1])*(1/8)+q[1](t[1], t[2])+3*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(-(2*I)*t[1]-I*omega*t[1])*(1/8)+3*alpha[2](t[1], t[2])*A(t[2])^3*exp(I*t[1])*(1/8)+3*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(-(2*I)*t[1]+I*omega*t[1])*(1/8)+3*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(I*t[1]-(2*I)*omega*t[1])*(1/8)+3*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(I*t[1]+(2*I)*omega*t[1])*(1/8)+diff(q[1](t[1], t[2]), t[1], t[1])+3*alpha[2](t[1], t[2])*A(t[2])^2*b*exp((2*I)*t[1]-I*omega*t[1])*(1/8)+3*alpha[2](t[1], t[2])*A(t[2])^2*b*exp((2*I)*t[1]+I*omega*t[1])*(1/8)+I*alpha[1](t[1], t[2])*A(t[2])*exp(I*t[1])*(1/2)+3*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(I*omega*t[1])*(1/4)+3*alpha[2](t[1], t[2])*A(t[2])*exp(I*t[1])*b^2*(1/4);

(3/4)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(-I*omega*t[1])+(3/4)*alpha[2](t[1], t[2])*A(t[2])*exp(-I*t[1])*b^2+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(-I*t[1]-(2*I)*omega*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(-I*t[1]+(2*I)*omega*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(-(2*I)*t[1]-I*omega*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(-(2*I)*t[1]+I*omega*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(I*t[1]-(2*I)*omega*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(I*t[1]+(2*I)*omega*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp((2*I)*t[1]-I*omega*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp((2*I)*t[1]+I*omega*t[1])+(3/4)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(I*omega*t[1])+(3/4)*alpha[2](t[1], t[2])*A(t[2])*exp(I*t[1])*b^2-((1/2)*I)*alpha[1](t[1], t[2])*A(t[2])*exp(-I*t[1])+((1/2)*I)*alpha[1](t[1], t[2])*A(t[2])*exp(I*t[1])+I*(diff(A(t[2]), t[2]))*exp(I*t[1])+(1/8)*alpha[2](t[1], t[2])*b^3*exp((3*I)*omega*t[1])+(1/8)*alpha[2](t[1], t[2])*b^3*exp(-(3*I)*omega*t[1])+(1/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp((3*I)*t[1])+(1/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp(-(3*I)*t[1])+(3/8)*alpha[2](t[1], t[2])*b^3*exp(-I*omega*t[1])+(3/8)*alpha[2](t[1], t[2])*b^3*exp(I*omega*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp(-I*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp(I*t[1])-I*(diff(A(t[2]), t[2]))*exp(-I*t[1])+((1/2)*I)*alpha[1](t[1], t[2])*b*omega*exp(I*omega*t[1])-((1/2)*I)*alpha[1](t[1], t[2])*b*omega*exp(-I*omega*t[1])+diff(diff(q[1](t[1], t[2]), t[1]), t[1])+q[1](t[1], t[2])

(1)

 

 

eq:=combine(subs(w=3+e*sigma,%),exp);

(3/4)*alpha[2](t[1], t[2])*A(t[2])*exp(-I*t[1])*b^2+(3/4)*alpha[2](t[1], t[2])*A(t[2])*exp(I*t[1])*b^2-((1/2)*I)*alpha[1](t[1], t[2])*A(t[2])*exp(-I*t[1])+((1/2)*I)*alpha[1](t[1], t[2])*A(t[2])*exp(I*t[1])+(3/4)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(-I*(epsilon*sigma+3)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(-I*t[1]-(2*I)*(epsilon*sigma+3)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(-I*t[1]+(2*I)*(epsilon*sigma+3)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(-(2*I)*t[1]-I*(epsilon*sigma+3)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(-(2*I)*t[1]+I*(epsilon*sigma+3)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(I*t[1]-(2*I)*(epsilon*sigma+3)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(I*t[1]+(2*I)*(epsilon*sigma+3)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp((2*I)*t[1]-I*(epsilon*sigma+3)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp((2*I)*t[1]+I*(epsilon*sigma+3)*t[1])+(3/4)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(I*(epsilon*sigma+3)*t[1])+I*(diff(A(t[2]), t[2]))*exp(I*t[1])+(1/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp((3*I)*t[1])+(1/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp(-(3*I)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp(-I*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp(I*t[1])-I*(diff(A(t[2]), t[2]))*exp(-I*t[1])+(1/8)*alpha[2](t[1], t[2])*b^3*exp((3*I)*(epsilon*sigma+3)*t[1])+(3/8)*alpha[2](t[1], t[2])*b^3*exp(I*(epsilon*sigma+3)*t[1])+(1/8)*alpha[2](t[1], t[2])*b^3*exp(-(3*I)*(epsilon*sigma+3)*t[1])+(3/8)*alpha[2](t[1], t[2])*b^3*exp(-I*(epsilon*sigma+3)*t[1])+((1/2)*I)*alpha[1](t[1], t[2])*b*(epsilon*sigma+3)*exp(I*(epsilon*sigma+3)*t[1])-((1/2)*I)*alpha[1](t[1], t[2])*b*(epsilon*sigma+3)*exp(-I*(epsilon*sigma+3)*t[1])+diff(diff(q[1](t[1], t[2]), t[1]), t[1])+q[1](t[1], t[2])

(2)

expand(%);

(3/4)*alpha[2](t[1], t[2])*A(t[2])*exp(I*t[1])*b^2+((1/2)*I)*alpha[1](t[1], t[2])*A(t[2])*exp(I*t[1])+I*(diff(A(t[2]), t[2]))*exp(I*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp(I*t[1])+(1/8)*alpha[2](t[1], t[2])*A(t[2])^3*(exp(I*t[1]))^3+(1/8)*alpha[2](t[1], t[2])*A(t[2])^3/(exp(I*t[1]))^3+(3/8)*alpha[2](t[1], t[2])*A(t[2])^3/exp(I*t[1])-I*(diff(A(t[2]), t[2]))/exp(I*t[1])+diff(diff(q[1](t[1], t[2]), t[1]), t[1])+q[1](t[1], t[2])+(1/8)*alpha[2](t[1], t[2])*b^3*(exp(I*epsilon*sigma*t[1]))^3*(exp(I*t[1]))^9+(3/8)*alpha[2](t[1], t[2])*b^3*exp(I*epsilon*sigma*t[1])*(exp(I*t[1]))^3+(1/8)*alpha[2](t[1], t[2])*b^3/((exp(I*epsilon*sigma*t[1]))^3*(exp(I*t[1]))^9)+(3/8)*alpha[2](t[1], t[2])*b^3/(exp(I*epsilon*sigma*t[1])*(exp(I*t[1]))^3)+(3/4)*alpha[2](t[1], t[2])*A(t[2])*b^2/exp(I*t[1])-((1/2)*I)*alpha[1](t[1], t[2])*A(t[2])/exp(I*t[1])+(3/4)*alpha[2](t[1], t[2])*A(t[2])^2*b/(exp(I*epsilon*sigma*t[1])*(exp(I*t[1]))^3)+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2/((exp(I*epsilon*sigma*t[1]))^2*(exp(I*t[1]))^7)+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*(exp(I*epsilon*sigma*t[1]))^2*(exp(I*t[1]))^5+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b/(exp(I*epsilon*sigma*t[1])*(exp(I*t[1]))^5)+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(I*epsilon*sigma*t[1])*exp(I*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2/((exp(I*epsilon*sigma*t[1]))^2*(exp(I*t[1]))^5)+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*(exp(I*epsilon*sigma*t[1]))^2*(exp(I*t[1]))^7+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b/(exp(I*epsilon*sigma*t[1])*exp(I*t[1]))+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(I*epsilon*sigma*t[1])*(exp(I*t[1]))^5+(3/4)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(I*epsilon*sigma*t[1])*(exp(I*t[1]))^3+((3/2)*I)*alpha[1](t[1], t[2])*b*exp(I*epsilon*sigma*t[1])*(exp(I*t[1]))^3-((3/2)*I)*alpha[1](t[1], t[2])*b/(exp(I*epsilon*sigma*t[1])*(exp(I*t[1]))^3)+((1/2)*I)*alpha[1](t[1], t[2])*b*exp(I*epsilon*sigma*t[1])*(exp(I*t[1]))^3*epsilon*sigma-((1/2)*I)*alpha[1](t[1], t[2])*b*epsilon*sigma/(exp(I*epsilon*sigma*t[1])*(exp(I*t[1]))^3)

(3)

eq:=combine(subs({e^(t1*e)=e^t2},%),exp);

(3/4)*alpha[2](t[1], t[2])*A(t[2])*exp(-I*t[1])*b^2+(3/4)*alpha[2](t[1], t[2])*A(t[2])*exp(I*t[1])*b^2-((1/2)*I)*alpha[1](t[1], t[2])*A(t[2])*exp(-I*t[1])+((1/2)*I)*alpha[1](t[1], t[2])*A(t[2])*exp(I*t[1])+(1/8)*alpha[2](t[1], t[2])*b^3*exp((3*I)*t[1]*epsilon*sigma+(9*I)*t[1])+(3/8)*alpha[2](t[1], t[2])*b^3*exp(I*epsilon*sigma*t[1]+(3*I)*t[1])+(1/8)*alpha[2](t[1], t[2])*b^3*exp(-(3*I)*t[1]*epsilon*sigma-(9*I)*t[1])+(3/8)*alpha[2](t[1], t[2])*b^3*exp(-I*epsilon*sigma*t[1]-(3*I)*t[1])+I*(diff(A(t[2]), t[2]))*exp(I*t[1])+(1/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp((3*I)*t[1])+(1/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp(-(3*I)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp(-I*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^3*exp(I*t[1])-I*(diff(A(t[2]), t[2]))*exp(-I*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(-(2*I)*epsilon*sigma*t[1]-(7*I)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*exp((2*I)*epsilon*sigma*t[1]+(5*I)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(-I*epsilon*sigma*t[1]-(5*I)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(I*epsilon*sigma*t[1]+I*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*exp(-(2*I)*epsilon*sigma*t[1]-(5*I)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])*b^2*exp((2*I)*epsilon*sigma*t[1]+(7*I)*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(-I*epsilon*sigma*t[1]-I*t[1])+(3/8)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(I*epsilon*sigma*t[1]+(5*I)*t[1])+diff(diff(q[1](t[1], t[2]), t[1]), t[1])+q[1](t[1], t[2])+(3/4)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(-I*epsilon*sigma*t[1]-(3*I)*t[1])+(3/4)*alpha[2](t[1], t[2])*A(t[2])^2*b*exp(I*epsilon*sigma*t[1]+(3*I)*t[1])+((3/2)*I)*alpha[1](t[1], t[2])*b*exp(I*epsilon*sigma*t[1]+(3*I)*t[1])-((3/2)*I)*alpha[1](t[1], t[2])*b*exp(-I*epsilon*sigma*t[1]-(3*I)*t[1])+((1/2)*I)*alpha[1](t[1], t[2])*b*epsilon*sigma*exp(I*epsilon*sigma*t[1]+(3*I)*t[1])-((1/2)*I)*alpha[1](t[1], t[2])*b*epsilon*sigma*exp(-I*epsilon*sigma*t[1]-(3*I)*t[1])

(4)

How could I substitute the (t1*e) into t2 at the index of e?  I am trying many ways to do that, but it doesn't work. please help me.   (t1*e)=t2


 

Download question.mw

Hi all

Hope the best for all. 

I have a single product of two sum, where the second sum admits one variable, namely j, from first one.

how can I it via maple?

thanks for any help.

I have a function g(x,t), which is the numerical result of a PDE.

How to calculate the function f (t), defined by:

f:=t->int(g(x,t),x=a..b)

I tried to calculate that way, but that did not work.

Thanks for the help.

 f(a)(x) = a*x mod 256  for all x in Z[256]

how can write this loop that has written in matlab to maple??
i=1;
for e=T1P1:0.0399:T1P;
    E(i)=e;
    gama1(i,:)=a-atan(e/Rb1);
     i=i+1;
end
 

HELLO

I AM A STUDENT AND WAS LOOKING FOR USING MAPLE2018.

I REQUESTED FOR A MAPLE EVALUATION. MY QUESTIONS ARE:

FOR HOW LONG CAN I USE THIS EVALUATION VERSION?

SECONDLY, WHAT SOFTWARE FACILITIES I WOULD NOT GET WITH THIS EVALUATION VERSION? WHAT ARE THE DIFFERENCES BETWEEN THE PRICED ENTERPRISE VERSION AND MAPLESOFT EVALUATION VERSION?

THIRDLY, WHAT ARE THE DIFFERENCES BETWEEN MAPLESOFT STUDENT VERSION AND THE ULTIMATE PRICED ENTERPRISED VERSION?

PLEASE REVERT BACK TO ME AT THE EARLIEST.

Hi,

this has been asked before but none of the other answers worked for me. I want to generate code from expressions that were simplified using some assumptions. I the assumptions remain on the variables Maple replaces variable names wich I don't understand.

Here is my script where I tried to remove the assumptions based on the answers of previous questions. Unfortunately, even though assumptions are removed in the end, Code generation still replaces the variables.

Thank you in advance!

Here is my minimal exaple Code:

 

restart:

 

# test removing assumptions

funa:=x+y

x+y

(1)

 Maple_proc := codegen[makeproc](funa,(ListTools[Flatten]([eval(alist,1)])));

proc (alist) x+y end proc

(2)

# make some assumptions

assume(x,real);assume(y,real);

x,y

x, y

(3)

#place variables in a list

alist:=[x,y];

[x, y]

(4)

# make cool interferences with assumptions...

# clear assumptions

# try 1 from https://www.mapleprimes.com/questions/207601-Remove-Assumptions-

nms:=convert(indets(alist,name),list);
nmsS:=convert~(nms,string);
L:=StringTools:-Substitute~(nmsS,"~",""); #Removing "~"
L1:=parse~(L);
S:=nms=~L1;
Expr:=subs(S,expr);

[x, y]

 

["x~", "y~"]

 

["x", "y"]

 

[x, y]

 

[x = x, y = y]

 

expr

(5)

hasassumptions(x)

true

(6)

# try 2

nops(alist)

2

(7)

for i from 1 to nops(alist) do
parse(cat(StringTools[Substitute]~(alist[i],"~",""),":='",StringTools[Substitute]~(alist[i],"~",""),"'"));
end do;

'x'

 

'y'

(8)

# check if assumtions are there

hasassumptions(x)

true

(9)

#try 3 https://www.mapleprimes.com/questions/39555-Unassume

for u in alist do
`property/object`[u]:=evaln(`property/object`[u]);
`property/OrigName`[u]:=evaln(`property/OrigName`[u]);
end do;

`property/object`[x]

 

`property/OrigName`[x]

 

`property/object`[y]

 

`property/OrigName`[y]

(10)

hasassumptions(x)

false

(11)

hasassumptions(y)

false

(12)

alist

[x, y]

(13)

funb:=x+y

x+y

(14)

 Maple_proc := codegen[makeproc](funb,(ListTools[Flatten]([eval(alist,1)])));

proc (x, y) x+y end proc

(15)

 Maple_proc := codegen[makeproc](funa,(ListTools[Flatten]([eval(alist,1)])));

proc (x, y) x+y end proc

(16)

 

 CodeGeneration[C](Maple_proc, optimize, declare=[x::float,y::float], defaulttype=float, deducetypes=false, coercetypes=false, output="fun.c");

Warning, the following variable name replacements were made: x~ -> cg, y~ -> cg1

 

 

>

 

 

 

 

 

 

 

 


 

Download testunassuming.mw

please I need the code to calculate the convergence order of a multipoint iterative method to solve nonlinear equation systems.

 

my email is: yjmsigma@gmail.com

 

is there a way to change the color of just the labels on the axis ticks only? 

I wanted to do this for DEplot, but I started with normal plot since it is easier. I could change the colors of the axis, but this changes not only the tick labels but the whole axis itself:

restart;
plot(sin(x),x=-Pi..Pi,axis=[tickmarks=[color=red]], color=blue);

So the axis itself changed color as well. It is OK for the ticks as well as the ticks labels (the numbers of the ticks). But prefer the axis itself not to be red.

Once there is a way to do the above, I wanted to apply to DEtools:-DEplot. Why? Becuase it is hard to read the Axis tick labels on DEplot. I thought that changing the color of the tick labels will it easy. But DEtools:-DEplot does not even accept axis command:

ode:= diff(y(x),x)=2*x:
DEtools:-DEplot(ode,y(x),x = -2 .. 2,y = -2 .. 2, [[0.1,0]],
               labels=["",""],
               linecolour = red,
               color = blue,
               'arrows' = 'medium',
               axesfont=['Times', 'bold', 12]
               );

 

You can see it is hard to read the axis tick labels above, as arrows go over some of them,. If I make them say RED, I think it will make it more clear.

Now when adding the axis options, it gives error

DEtools:-DEplot(ode,y(x),x = -2 .. 2,y = -2 .. 2, [[0.1,0]],
               linecolour = red,
               color = blue,
               arrows = MEDIUM,
               axesfont=['Times', 'bold', 10],
               axis=[tickmarks=[color=red]]
               );

Error is generated, axis is not supported command.

Any idea how to make tick labels red for DEtools:-DEplot? It will be OK if I can get both the axis and the ticks RED, at least it will be better than it is now.

One workaround way I currently use is this:

ode:= diff(y(x),x)=2*x:
p1:=DEtools:-DEplot(ode,y(x),x = -2 .. 2,y = -2 .. 2, [[0.1,0]],
               labels=["",""],
               linecolour = red,
               color = blue,
               'arrows' = 'medium'                    
               );
p2:=plot(0,x=-2..2,y=-2..2,axis=[tickmarks=[color=red]]);
plots:-display([p2,p1]);

 

But this seems like a hack to me. There should be a way to directly set tickmarks color in DEplot?

Help me extracting the data from the Maple figure. I was used the 'implicitplod3D' command to draw the figure.  I want to save the parameters of sigma, A  and f0 in Exce
 

af_eq := (3*alpha[2](t[1], t[2])*A^3*(1/8)-A*sigma)^2+(1/4)*alpha[1](t[1], t[2])^2*A^2 = (1/4)*f0^2

(0.4871250000e-1*A^3-A*sigma)^2+0.2500000000e-2*A^2 = (1/4)*f0^2

(1)

with(plots, implicitplot,implicitplot3d):

 alpha[1]:=0.1;alpha[2]:=0.1299;

.1

 

.1299

(2)

##3D

implicitplot3d(af_eq, sigma = -1.5 .. 1.5,f0 =0..0.5, A = 0 .. 5,numpoints = 20000,style=surfacecontour,contours = 15);

 

 

 


 

Download extract_data_from_the_figure.mw

l format. then, I can use the data draw a picture in other software. Many thanks.

First 656 657 658 659 660 661 662 Last Page 658 of 2097