Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

 

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

I am trying to do a small procedure which converts an input of degrees into the sin function into its corresponding unit-circle value. I can get this to work. 

Sinus:= proc(g::unit)   
return evalf(sin(g*Pi/180));
end proc; 

 

but is there a more elegant way to do it? Using maybe Units package ? Thanks in advance :)

Hello!

How can I make MAPLE to put out the numerical solution of the following system?
(a =0.12, c = 47.04)  Neither solve nor fsolve does the job!

Thank for your help!

>restart:
>S:=t->c*exp(-a*t)+18;

                      S := t -> c*exp(-a*t)+18

> sys:={S(2)=55,S(8)=36};

        sys := {c*exp(-8*a)+18 = 36, c*exp(-2*a)+18 = 55}

> solve(sys);

         {a = -1/2*ln(18/37*RootOf(18*_Z^3-37)^2), c = 37*RootOf(18*_Z^3-37)}

> fsolve(sys);

         fsolve({c*exp(-8*a)+18 = 36, c*exp(-2*a)+18 = 55},{a, c})

 

I have a differential equation involving several functions of the following form:

diff(h,z) = iAf + iBg,

where h, f and g are functions of the Cartesian coordinates x, y and R and the third coordinate corresponds to z = R for some fixed constant value R.  The derivative is then with respect to the coordinate z and A and B are constants, with i the usual imaginary unit.  Is there some way this equation could be solved explicitly with Maple?

I have a system of nonlinear differential equations with initial conditions:


 

restart

de0 := {(1-p)*(diff(f(x), x, x))+p*(beta*[5]*f(x)*(diff(f(x), x))^2+beta[11]*f(x)^7+beta[6]*(diff(f(x), x))^3+beta[7]*f(x)^5+f(x)+diff(f(x), x, x)+beta[3]*f(x)^3+beta[1]*f(x)+beta[2]*(diff(f(x), x))+mu[1]*(diff(f(x), x))+(3/4)*rho*f(x)^3+beta[9]*f(x)^2*(diff(f(x), x))^3+beta[12]*f(x)^4*(diff(f(x), x))^3+beta[13]*f(x)*(diff(f(x), x))^6+beta[10]*f(x)*(diff(f(x), x))^4+beta[8]*f(x)^4*(diff(f(x), x))+beta[4]*f(x)^2*(diff(f(x), x))-(1/4)*rho*(g(x)^3-3*f(x)^2*g(x))*sin(4*Omega*x)-alpha[2]*f(x)^3*cos(2*Omega*x)-(1/4)*rho*(3*f(x)*g(x)^2-f(x)^3)*cos(4*Omega*x)-alpha[1]*g(x)*sin(2*Omega*x)-(1/2)*alpha[2]*(f(x)^3+3*f(x)^2*g(x))*sin(2*Omega*x)+(3/4)*rho*f(x)*g(x)^2-alpha[1]*f(x)*cos(2*Omega*x)-F*Omega^2*cos(Omega*x+theta)), (1-p)*(diff(g(x), x, x))+p*(diff(g(x), x, x)+beta[11]*g(x)^7+beta[7]*g(x)^5+beta[6]*(diff(g(x), x))^3+beta[1]*g(x)+beta[2]*(diff(g(x), x))+beta[3]*g(x)^3+mu[2]*(diff(g(x), x))+(3/4)*rho*g(x)^3+g(x)+beta[9]*g(x)^2*(diff(g(x), x))^3+beta[13]*g(x)*(diff(g(x), x))^6+beta[8]*g(x)^4*(diff(g(x), x))+beta[10]*g(x)*(diff(g(x), x))^4+beta[4]*g(x)^2*(diff(g(x), x))+beta[12]*g(x)^4*(diff(g(x), x))^3+beta[5]*g(x)*(diff(g(x), x))^2-(1/4)*rho*(3*f(x)*g(x)^2-f(x)^3)*sin(4*Omega*x)+(3/4)*rho*g(x)*f(x)^2-F*Omega^2*sin(Omega*x+theta)-(1/2)*alpha[2]*(f(x)^3+3*f(x)*g(x)^2)*sin(2*Omega*x)-alpha[2]*g(x)^3*cos(2*Omega*x)-(1/4)*rho*(3*f(x)^2*g(x)-g(x)^3)*cos(4*Omega*x)-alpha[1]*f(x)*sin(2*Omega*x)+alpha[1]*g(x)*cos(2*Omega*x))}:

ibvc0 := {f(0), g(0), (D(f))(0), (D(g))(0)}:

n := 3:

F := unapply(add(b[k](x)*p^k, k = 0 .. n), x):

de := map(series, eval(de0, {f = F, g = G}), p = 0, n+1):

for k from 0 to n do if k = 0 then ibvc := expand(eval[recurse](ibvc0, {f = F, g = G, p = 0})) else ibvc := {b[k](0)-1, (D(b[k]))(0)-1, (D(c[k]))(0)-1, ((D@@2)(b[k]))(0)-1, ((D@@2)(c[k]))(0)-1, c[k](0) = 1} end if; sys := `union`(simplify(map(coeff, de, p, k)), ibvc); soln := dsolve(sys); b[k] := unapply(eval(b[k](x), soln), x); c[k] := unapply(eval(c[k](x), soln), x) end do; 'F(x)' = F(x); 'G(x)' = G(x)

Error, (in dsolve) invalid input: `PDEtools/NumerDenom` expects its 1st argument, ee, to be of type algebraic, but received [0]

 

F(x) = b[1](x)*p+b[2](x)*p^2+b[3](x)*p^3

 

G(x) = c[1](x)*p+c[2](x)*p^2+c[3](x)*p^3

(1)

``


How would I deal with this error message

Download 6_11_2019_nonlinear_d_e.mw
 

>  

 

 

 

 

 

Basically it spits out the subset of values for which a division by zero error will occur for the function you specify on  range you specify for each of it's arguments, but I get an ambigous error when ever exponentiation features in the function I specify, which of course dramatically reduces the application of the calculator. Division,addition,substraction and multiplication are currently the only available arithmetic operators availble for the function window that I know the error will not occur.

If some one can help it is much appriciated

 

DIVISION_BY_ZERO_CALCULATOR.mw

First 615 616 617 618 619 620 621 Last Page 617 of 2218