Kitonum

14009 Reputation

24 Badges

11 years, 190 days

MaplePrimes Activity


These are replies submitted by Kitonum

These 2 expressions are not equal:

restart;
A:=Ve*k21*ke*x1+Ve*k21*x1^2+k12*ke^2*x1[t]+2*k12*ke*x1*x1[t]+k12*x1^2*x1[t]+k21*ke^2*x1[t]+2*k21*ke*x1*x1[t]+k21*x1^2*x1[t]+Ve*ke*x1[t]+ke^2*x1[t,t]+2*ke*x1*x1[t,t]+x1^2*x1[t,t]:
B:=x1[t,t]*(x1+k2)^2+(k12+k21)*x1[t]*(x1+ke)^2+Ve*x1[t]*ke+(k21*Ve)*x1*(x1+ke):

expand(A-B);

Output:    -k2^2*x1[t, t]-2*k2*x1*x1[t, t]+ke^2*x1[t, t]+2*ke*x1*x1[t, t]

@acer  We can replace the last 2 lines of your code with one line:

subs(1-isolate(e_nsp2,gamma), isolate(e_nsp1,gamma), set_3);

Maple sequentially performs these 2 substitutions.

@lhasif  Thank you. The surface will look more expressive if you use the following code:

with(Student:-Calculus1):
f := x->(x-1)^2:
g := x->x+1:
A:=SurfaceOfRevolution(f(x), x = 0 .. 3, distancefromaxis = 4, surfaceoptions=[color=khaki, style=patch], output=plot):
B:=SurfaceOfRevolution(g(x), x = 0 .. 3, distancefromaxis = 4, surfaceoptions=[color=khaki, style=patch], output=plot):
plots:-display(A,B, caption=``, lightmodel=light3, scaling=constrained);

                      

For even greater expressiveness, we can make a cutout on this surface in order to clearly see the curve to be rotated and the axis of rotation. But in this case, we must specify the surface and these curves by explicit equations:

f := x->(x-1)^2;
g := x->x+1;
A:=plot3d([[x,(4-f(x))*cos(phi),(4-f(x))*sin(phi)+4],[x,(4-g(x))*cos(phi),(4-g(x))*sin(phi)+4]],x=0..3,phi=-Pi/2..5*Pi/6, style=surface, color=khaki):
B:=plots:-spacecurve([[t,0,f(t)],[t,0,g(t)]], t=0..3, color=red, thickness=3): # The line to be rotated
C:=plots:-spacecurve([t,0,4], t=-1..4, linestyle=3, color=blue, thickness=2): # Axis of rotation
plots:-display(A,B,C, axes=normal, scaling=constrained, view=[-1..4,-4..4,0..8], orientation=[-125,75]);

                 

 

@Teep You can do it similarly:


 

``

restart; with(plots)

X := Vector([1, 2, 3, 4, 5, 6, 7, 8, 9, 10])

Y1 := Vector([1, 3, 6, 2, 9, 4, 2, 7, 4, 9])

Y2 := Vector([1, 2, 1, 4, 11, 10, 7, 8, 9, 12])

display(plot([], x = 1 .. 10), seq(`$`(Statistics:-ColumnGraph([Y1[1 .. i], Y2[1 .. i]]), 5), i = 1 .. 10), insequence)

 

``


 

Download MaplePrimes_Example_animChart.mw

 

@Teep    I don’t understand what you call an animated column chart. Write the code for a non-animated chart, I'll try to animate it.

@Carl Love  OP wrote  "...n is natural and k=0..n"

Of course, your interpretation is also possible, but only the OP  can clarify what was meant.

@Carl Love  Of course, it is always useful to have an alternative solution. In my opinion, it’s easier to interpret your line 

B:= eval(L, t= solve(L.v = C.v)); 

(without mentioning the plane) how "the vector  L - C  is perpendicular to the vector  v " and rewrite it as

B:= eval(L, t= solve((L-C).v)); 

@tomleslie  I corrected my answer.

@Zeineb To get a solution, we can take arbitrary vectors, but  n  must be specific. This solution can be easily written as a procedure in which  n  will be one of the parameters. But when using this procedure, we must indicate the specific value of .

@JTB2222   Similarly.

@ActiveUser 

You wrote a lot of monotonous code, in each line of which you assign a certain expression to  g , and then do some substitution. If you think Maple is calculating something incorrectly, select only this line from the entire code and ask a specific question.

@Zen 

evalindets(g*a*(a*g+g*b+u^2*(Df-Dr))*b/(Df*Dr*(a+b)*u),`+`, t->simplify(t, {a+b=L}));

                             g*a*(L*g+u^2*(Df-Dr))*b/(u*L*Dr*Df)

@Christian Wolinski  You can have several versions of Maple at the same time in one system. On a fragment of my desktop there are labels of 4 versions (Maple 2017 in two options - Classic and Standard):

@Christian Wolinski  I do not know if this is possible. But you can always tweak Standard worksheet to make it look exactly like Classic worksheet.

@ActiveUser  Probably all your errors are related to the incorrect use of the  subs  command to substitute into the derivative that has not yet been calculated. Use the more "clever" the  eval  command for this.

See example:

A:=subs(x=1, diff(m(x),x));
B:=eval(diff(m(x),x), x=1);
m:=x->x^2;
A;
B;

                                      

 

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