MaplePrimes Questions

Hi,

I asked a question before, how to shift subscript of u(i,j):

P: u[i,j] -->  u[i+1,j]

https://www.mapleprimes.com/questions/223136-How-To-Define-An-Operator

How about an arbitrary expression?

P: u[i,j]^2 v[i,j]+u[i-1,j]  -->  u[i+1,j]^2 v[i,j]+u[i,j]

thanks.

 

See below, why the difference + what is the distinction between the 2?  Notice that they yield different results with the is command.

restart; P := a[0]+sum(4*(int(cos(2*Pi*(2*k-1)*x/T), x = 0 .. t))/T, k = 1 .. n); Q := a[0]+4*(int(sum(cos(2*Pi*(2*k-1)*x/T), k = 1 .. n), x = 0 .. t))/T; is(expand(P = Q)), is(combine(P = Q)); is(value(eval(P = Q, n = 10)))

false

(1)

restart; P := a[0]+Sum(4*(Int(cos(2*Pi*(2*k-1)*x/T), x = 0 .. t))/T, k = 1 .. m); Q := a[0]+Int(Sum(4*cos(2*Pi*(2*k-1)*x/T)/T, k = 1 .. m), x = 0 .. t); is(expand(P = Q)), is(combine(P = Q)); is(value(eval(P = Q, m = 10)))

true

(2)

sum(4*(int(cos(2*Pi*(2*k-1)*x/T), x = 0 .. t))/T, k = 1 .. n)

((1/2)*I)*(-ln(-(1+exp((2*I)*Pi*t/T))/(exp((2*I)*Pi*t/T)-1))+exp((2*I)*Pi*(1+2*n)*t/T)*LerchPhi(exp((4*I)*Pi*t/T), 1, 1/2+n)+ln((1+exp((2*I)*Pi*t/T))/(exp((2*I)*Pi*t/T)-1))-LerchPhi(exp(-(4*I)*Pi*t/T), 1, 1/2+n)*exp(-(2*I)*Pi*(1+2*n)*t/T))/Pi

(3)

a[0]+Sum(4*(Int(cos(2*Pi*(2*k-1)*x/T), x = 0 .. t))/T, k = 1 .. m)

a[0]+Sum(4*(Int(cos(2*Pi*(2*k-1)*x/T), x = 0 .. t))/T, k = 1 .. m)

(4)

``


 

Download sum_versus_Sum.mw

 

where i can change the value?

I'd like to find the harmonic mean of the lengths of the three altitudes of a triangle.  I believe this to be equal to the radius of the incenter of the triangle.  I haven't used the geometry package  much, but notice it has some interesting features like Harmonic mean. (Incidentally, it has a vertex, and vertices command, but these don't seem to be recognised in Maple 7 - not even in the Help section.)  I would normally tackle thissue  problem using coordinate geometry, and give a diagram of the triangle using the plots [display] command.  Id like some advice about the interactivity of commands in geometry and plots packages eg  can I find the coordinates of the base of an altitudemail using the geometry package?  There is an altitude command which gives the àltitude's equation.  For this the coding would be

:triangle(ABC, [point (A,0,0),point (B,,2,0), point (C,1,3)]:

altitude (hA1,A, ABC);

detail (hA1);

   I think I have answered my own question!  I see there is a DefinedAs command which would give me the coords.   Am I correct in thinking the geometry package is the best method of doing this.  Later I was hoping to display animation, to see how the Harmonic mean changes with varying different triangles.

Thanks for reading this and Happy New Year.

David  

Hallo,

I have a list of graphs made correctly with the package GraphTheory: g=[g[1],g[2],g[3]], where g[1], g[2] and g[3] are graph objects. GraphTheory[DrawGraph] inside a for-loop doesn't work. GraphTheory[DrawGraph](g) draws all graphs simultaneously. How can I draw the single graphs of the list successively when the number of members in g is a variable?

Thanks.

Hello everyone,

Compliment of the season. I received the error message while trying to solve some Pade approximant problem and try to find what could possibly be the error but could not.

Please, anyone with good informations on the possible solutions should help.

Thanks. The Link is below

 Error.mw

Hello.
Can I somehow print several matrices next to each other without getting the quotations marks printet too?
Attached what I currently have.


 

with(LinearAlgebra)

A := RandomMatrix(4); B := RandomMatrix(4); C := RandomMatrix(4); E := RandomMatrix(4)

print("A="*A, "   B="*B, "   C="*C, "   E="*E)

"A="*Matrix(%id = 18446746685123259198), "   B="*Matrix(%id = 18446746685123259438), "   C="*Matrix(%id = 18446746685123259678), "   E="*Matrix(%id = 18446746685123259918)

(1)

``


 

Download matrix_print_attempt.mw

 

 Hi,  

Is there a way of determining whether a function of two variables is continuous or not? Can we find the discontinuity points 

of such multivariable functions? (Considering the graphs of the function is not always helpful.)

Hi, I need to calculate some discrete variational problems, but it's a hard work just by hand.

Is there any toolboxes about the discrete variational computation? Or how to design a program to do that?

In particular, I need to realize the functions:

1. the linearized operator and the adjoint operator of some expressions.

       For example:

       the linearized operator of u(n,m+1)u(n+1,m) is u(n+1,m)Tm+u(n,m+1)Tn

       the adjoint operator is u(n+1,m-1)T^(-1)m+u(n-1,m+1)T^(-1)n,

       (The latter is just the discrete Euler operator.)

      Tm, Tn is the shift operator, T^(-1) is the inverse operator of T.

2. discrete Euler-Lagrange equation.

I am not familar with the maple, please help!

Thanks.

Hello

What is the fairest way to divide a shelf of books up between n people?

Partition_problem2.mw

Below is the output from test relation that seems contradictory.  All I did was swap the order of operation from summing an integral to taking the integral of the sum.  The summation and integration variables are independent so I would think the statement for S4 would yield a TRUE result.  Also, I got a FALSE return on a well known trig identity.

Am I missing something subtle or even obvious?

equivalence_discrepancy.mw

How do I introduce the first derivative of M when i<> k and aij when i <> j and aii in Maple?

 

 

How can i plot various 3 multiple graphs using Maple to Compare ADM, EXACT SOLUTION AND RKF45 using the table below:

 Time  ADM Solution                                   RKF45 Numerical Solution        EXACT Solution

               ADM Solution                     RKF45 Numerical Solution                EXACT Solution

 t           s               i              r                 s            i                r             s               i                   r

0.1     0.0769540597    0.1477783335     0.1393069312      0.0763361154   0.1477867626   0.1393072151   0.0763361266    0.1477867830     0.1393071934 

0.3     0.1093424148    0.1120401102     0.1509302274      0.1072278404   0.1121142033   0.1509369166   0.1072278479    0.1121142168     0.1509369024 

0.5     0.1392030568    0.0853083077     0.1558066181      0.1353291378   0.0855066806   0.1558355785   0.1353291558    0.08550671332   0.1558355439

0.7     0.1662374563    0.0652194693     0.1562235596      0.1604342222   0.0655908735   0.1562974878   0.1604342352    0.06559089617   0.1562974637 

0.9     0.1903821623    0.0500537619     0.1537887672      0.1825594352   0.0506356391   0.1539346707   0.1825594528    0.05063567192   0.1539346352

1.1     0.2117168860    0.0385555127     0.1496220815      0.2018541863   0.0393727175   0.1498707561   0.2018542024    0.03937274753   0.1498707234 

1.3     0.2303874000    0.0298096396     0.1444864012      0.2185409755   0.0308687065   0.1448804021   0.2185409880    0.03086872986   0.1448803765 

1.5     0.2465077820    0.0231661161     0.1388614678      0.2328759081   0.0244336214   0.1394893808   0.2328759200    0.02443364533   0.1394893543

 

 

 

Hello everyone.and complements

Please I am trying to obtain series expansion of the expression below in u but encounter difficulties particularly when b=0. I am very optimistic that when b=0 there will be a result not division by 0. Can I get help on the code?

Thank you in anticipation of your quick and positive responses and suggestions

# for k=2 CHEBY HYBRID WITH mu=(1-(1/2)*sqrt(2)))) AND v=(1+(1/2)*sqrt(2))))
restart:
omega:=u/h:
t:=(sum(a[j]*x^j,j=0..3)+a[4]*sin(omega*x)+a[5]*cos(omega*x)):
F:=diff(t,x):
G:=diff(t,x,x):
p1:=simplify(eval(t,x=q))=y[n]:
p2:=simplify(eval(t,x=q+(1-(1/2)*sqrt(2))*h))=y[n+mu]:
p3:=simplify(eval(t,x=q+h))=y[n+1]:
p4:=simplify(eval(t,x=q+(1+(1/2)*sqrt(2))*h))=y[n+v]:
p5:=simplify(eval(F,x=q+2*h))=f[n+2]:
p6:=simplify(eval(G,x=q+2*h))=g[n+2]:

vars:= seq(a[i],i=0..5):
Cc:=eval(<vars>, solve({p||(1..6)}, {vars})):
for i from 1 to 6 do
	a[i-1]:=Cc[i]:
end do:
Cf:=t:

K:=collect(combine(simplify(eval(Cf,x=q+2*h),size),trig),{y[n],y[n+mu],y[n+1],y[n+v],f[n+2], g[n+2]},factor):


Num := numer(K):
Den := denom(K):

N := 20:   # order of expansion
Num_N :=(convert(series(Num, u, N),polynom)):
Den_N := (convert(series(Den, u, N),polynom)):
b:=y[n+2]=(convert(series(Num_N/Den_N, u, N),polynom)):

eval(b,u=0); 

 

We have three vectors: position, velocity and accelleration

Then we calculate three unit vector of osculating plane and the center of curvature . Can you help me in creating a smart plot to show the idea

Es_1_1_2.mw

Thank's a lot and Regards

 

First 831 832 833 834 835 836 837 Last Page 833 of 2375