Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

help me! 

 

I have a problem with the system, looking forward to everyone's help!

expand((x-c)^2+(y-d)^2-R^2) = 0; algsubs(-R^2+c^2+d^2 = f, %); P := proc (x, y) options operator, arrow; -2*x*c-2*y*d+x^2+y^2+f = 0 end proc; 2 2 P := (x, y) -> -2 x c - 2 y d + x + y + f = 0 P(a*cos(theta), b*sin(theta)); G := unapply(%, theta); #usage des formules d'Euler simplify(expand(4*(exp(I*theta))^2*subs(cos(theta) = (exp(I*theta)+exp(-I*theta))*(1/2), sin(theta) = (exp(I*theta)-exp(-I*theta))/(2*I), G(theta)))); poly := sort(subs(exp(I*theta) = X, exp((2*I)*theta) = X^2, exp((3*I)*theta) = X^3, exp((4*I)*theta) = X^4, %)); coeff(lhs(poly), X^4)/tcoeff(lhs(poly)); # exp(I*theta1),exp(I*theta2),exp(I*theta3),exp(I*theta4) sont les racines de ce polynôme unitaire : exp(I*theta1)*exp(I*theta2)*exp(I*theta3)*exp(I*theta4) =1 exp(I*(θ1+θ2+θ3+θ4)=1 d'où θ1+θ2+θ3+θ4 ≡ 2*Pi

I am trying to construct a procedure which makes a boxplot from a list using the statistics package. 

I have come to this here. However it shows both the x and y - axis. Is there any way to make it show a horizontal boxplot like in my code. But only displays the x- axis ? Meaning the values from the plot? Plus any to make it not stick to x-axis like my code does? I tried the command scale, but it doesn't work in statistics package 

    BoxP:=proc(xd::list) 
        uses Statistics:
        return plots[display](BoxPlot(xd,deciles = false,orientation=horizontal),size=[600,400],axes=normal);
    end proc;

Assume 

matrix([[1,0],[0,1]]) * matrix[[y],[y']]
= a system 

 

how to get back matrix([[1,0],[0,1]]) from this system in general method if the system is very large?

Is there a function to convert recurrence equation to a system of equations?

 

I think we will write nested for loops or nested sequences. But I can' t achieve.

Could you help me?

EDITED AND EXTENDED:

In fact, the original question is as follows:

and

where it is a standard inner product on L^2 and u(x,t) is a function. For example; u(x,t) =x*t etc.

(The screenshots is taken from a book)

I want to write a code for finding Matrix U whose elements are u_ij.

Hello!

How can I make MAPLE to create the solution of the following system?

Program of Transfer matrix method for solving the vibrational behavior of a cracked beam. For more details and comprehension check the paper entitled ''A transfer matrix method for free vibration analysis and crack identification of stepped beams with multiple edge cracks and different boundary conditions ''
 

restart; with(LinearAlgebra)NULLNULL

W11 := A[1]*cos(nu*x)+A[2]*sin(nu*x)+A[3]*cosh(nu*x)+A[4]*sinh(nu*x);

theta11 := -A[1]*nu*sin(nu*x)+A[2]*nu*cos(nu*x)+A[3]*nu*sinh(nu*x)+A[4]*nu*cosh(nu*x)  NULL

M11 := EI*(-A[1]*nu^2*cos(nu*x)-A[2]*nu^2*sin(nu*x)+A[3]*nu^2*cosh(nu*x)+A[4]*nu^2*sinh(nu*x))
S11 := EI*(A[1]*nu^3*sin(nu*x)-A[2]*nu^3*cos(nu*x)+A[3]*nu^3*sinh(nu*x)+A[4]*nu^3*cosh(nu*x))
 

MD11 := subs(A[1] = 1, A[2] = 0, A[3] = 0, A[4] = 0, W11); MD12 := subs(A[1] = 0, A[2] = 1, A[3] = 0, A[4] = 0, W11); MD13 := subs(A[1] = 0, A[2] = 0, A[3] = 1, A[4] = 0, W11); MD14 := subs(A[1] = 0, A[2] = 0, A[3] = 0, A[4] = 1, W11)
NULL

MD21 := subs(A[1] = 1, A[2] = 0, A[3] = 0, A[4] = 0, theta11); MD22 := subs(A[1] = 0, A[2] = 1, A[3] = 0, A[4] = 0, theta11); MD23 := subs(A[1] = 0, A[2] = 0, A[3] = 1, A[4] = 0, theta11); MD24 := subs(A[1] = 0, A[2] = 0, A[3] = 0, A[4] = 1, theta11)
 

MD31 := subs(A[1] = 1, A[2] = 0, A[3] = 0, A[4] = 0, M11); MD32 := subs(A[1] = 0, A[2] = 1, A[3] = 0, A[4] = 0, M11); MD33 := subs(A[1] = 0, A[2] = 0, A[3] = 1, A[4] = 0, M11); MD34 := subs(A[1] = 0, A[2] = 0, A[3] = 0, A[4] = 1, M11)
 

MD41 := subs(A[1] = 1, A[2] = 0, A[3] = 0, A[4] = 0, S11); MD42 := subs(A[1] = 0, A[2] = 1, A[3] = 0, A[4] = 0, S11); MD43 := subs(A[1] = 0, A[2] = 0, A[3] = 1, A[4] = 0, S11); MD44 := subs(A[1] = 0, A[2] = 0, A[3] = 0, A[4] = 1, S11)
 

 

TM := Matrix(4, 4, [[MD11, MD12, MD13, MD14], [MD21, MD22, MD23, MD24], [MD31, MD32, MD33, MD34], [MD41, MD42, MD43, MD44]])

C := Matrix(4, 4, [[1, 0, 0, 0], [0, 1, c44, 0], [0, 0, 1, 0], [0, 0, 0, 1]])NULL

NULL    TM2 := subs(x = 0, TM); TM3 := subs(x = L, TM)

with(MTM)

TM4 := inv(TM)NULLNULL 

TM5 := inv(TM2)NULL

Y11 := MatrixMatrixMultiply(TM3, TM4)

    Y22 := MatrixMatrixMultiply(C, TM)   

Y33 := MatrixMatrixMultiply(Y11, Y22)

Y44 := MatrixMatrixMultiply(Y33, TM5)NULL

BB11 := Y44[3, 3]

BB12 := Y44[3, 4]

BB21 := Y44[4, 3]

BB22 := Y44[4, 4] 

BB := Matrix(2, 2, [[BB11, BB12], [BB21, BB22]]) 

NULL

R11 := det(BB) 

NULL

L := .18; L1 := 1; h := 0.5e-2; b := 0.2e-1; rho := 957.5; area = b.h; m := rho*h*b; EI := 0.2682e10*b*h^3mu := ((m.(omega^2))*L^4/EI)^(1/4); x := .5; c44 := .1; c11 := 0NULLNULL

plot(R11, omega = 1 .. 100)

 

 

 

``

NULL

NULL


 

Download transfer.mw

Hi, 

I'm currently solving an equations where the boundary conditions is at infinity. I'm trying to solve it by using dsolve but i can't seem to find a solutions. Here is my equations:

ode1 := diff(f(eta), eta$3)+(diff(f(eta), eta$2))*f(eta)-(diff(f(eta), eta))^2-M . (diff(f(eta), eta))-A . (diff(f(eta), eta)+(1/2)*(eta . (diff(f(eta), eta$2)))) = 0;

ode2 := diff(theta(eta), eta$2)+Pr*(f(eta) . (diff(theta(eta), eta))-(diff(f(eta), eta)) . theta(eta)-A . (theta(eta)+1/2 . eta . (diff(theta(eta), eta)))) = 0;

and my boundary conditions are:

bcs := f(0) = 0, (D(f))(0) = 1, ((D@@2)(f))(inf) = 0, theta(0) = 1, theta(inf) = 0;

The value of Pr=7, M=1 and A=[0,1,2,4]

I really need your help, please. 

Thank you :)

The solution to the logistic map .The solution now oscillates but doesn't appear to show any discernible pattern. The value of Xn seems to "jump around". This  called chaotic.

For convenience, I am looking to extract a sequence of numbers that is generated by a simple procedure.

The attached shows such an example.

In this (simple) case, the output I require is [1,4,9,16,25].

Can anyone suggest a way to do this?

Thank you all ...

MaplePrimes_Example.mw

E-_Original_RSA.pdf

 

If I want to calculate and know the exact time duration of running step 8, do the current algorithm, 

st := time();

elapsed_time := (time() - st)*sec

is good enough? precise enough? or accurate enough? 

Any other recommendations? 

https://math.stackexchange.com/questions/301068/how-do-you-find-a-corresponding-recurrence-relation-for-some-random-algorithm/301709

according to this link, how to parse or walk through the algorithm in maple to generate recurrence relation formula?

How to remove all punctuation marks in a string read from a text file 

and then store in a list of words according to the ordering of the original text?

Ahead of the upcoming deadline 14 jan 2020 for switching from Windows 7 to Windows 10,  I am investigating whether Maple 2019 will run on an HP Stream laptop: Celeron CPU N3060 @ 1.6GHz, RAM 2GB, 64 bit.running Windows 10.

My guess is that the above spec is too low to successfully run Maple 2019.  But if anyone has had success running Maple on a similar spec device please let me know.

Melvin

First 494 495 496 497 498 499 500 Last Page 496 of 2097