MaplePrimes Questions

My Question

how may i know which is the theta prime or f double prime or theta only or f prime and theta double prime?

 i have attach double prime

tq maple
 

restart

with(student)

with(plots)

inf := 3

equ1 := (diff(f(eta), `$`(eta, 3)))/((1-`ϕ`)^2.5*(1-`ϕ`+`ϕ`*rho[s]/rho[f]))+(diff(f(eta), eta, eta))*f(eta)-(diff(f(eta), eta))^2+1+M*(1-(diff(f(eta), eta))) = 0

equ2 := (1+4/(3*N*k))*(diff(theta(eta), eta, eta))+Pr*(1-`ϕ`+`ϕ`*rho[s]*Cp[s]/(rho[f]*Cp[f]))*(diff(theta(eta), eta))*f(eta)/k+Br*(diff(f(eta), eta, eta))^2/(k*(1-`ϕ`)^2.5) = 0

Bcs := f(0) = 0, (D(f))(0) = `ε`, (D(f))(inf) = 1, theta(inf) = 0, theta(0) = 1

Pr := 6.2; Cp[s] := 385; Cp[f] := 4179; `ϕ` := .1; rho[f] := 997.1; rho[s] := 8933; k[s] := 400; k[f] := .613; k := (k[s]+2*k[f]-2*`ϕ`*(k[f]-k[s]))/(k[s]+2*k[f]+`ϕ`*(k[f]-k[s])); Br := .1; M := 1; N := 1

func := proc (v) options operator, arrow; rhs((dsolve({equ1, equ2, subs(`ε` = v, [Bcs])[]}, numeric))(0)[3])/(1-`ϕ`)^2.5 end proc; plot(func, -1 .. 1, title = typeset((diff(f(eta), eta, eta))*versus*`ε`/(1-'`ϕ`')^2.5), titlefont = [times, italic, 18])

 

func2 := proc (v) options operator, arrow; rhs((dsolve({equ1, equ2, subs(`ε` = v, [Bcs])[]}, numeric))(0)[5]) end proc; plot(func2, -1 .. 1, title = typeset((diff(theta(eta), eta))*versus*`ε`), titlefont = [times, italic, 18])

 

``


 

Download plot.mw
 

restart

with(student)

with(plots)

inf := 3

equ1 := (diff(f(eta), `$`(eta, 3)))/((1-`ϕ`)^2.5*(1-`ϕ`+`ϕ`*rho[s]/rho[f]))+(diff(f(eta), eta, eta))*f(eta)-(diff(f(eta), eta))^2+1+M*(1-(diff(f(eta), eta))) = 0

equ2 := (1+4/(3*N*k))*(diff(theta(eta), eta, eta))+Pr*(1-`ϕ`+`ϕ`*rho[s]*Cp[s]/(rho[f]*Cp[f]))*(diff(theta(eta), eta))*f(eta)/k+Br*(diff(f(eta), eta, eta))^2/(k*(1-`ϕ`)^2.5) = 0

Bcs := f(0) = 0, (D(f))(0) = `ε`, (D(f))(inf) = 1, theta(inf) = 0, theta(0) = 1

Pr := 6.2; Cp[s] := 385; Cp[f] := 4179; `ϕ` := .1; rho[f] := 997.1; rho[s] := 8933; k[s] := 400; k[f] := .613; k := (k[s]+2*k[f]-2*`ϕ`*(k[f]-k[s]))/(k[s]+2*k[f]+`ϕ`*(k[f]-k[s])); Br := .1; M := 1; N := 1

func := proc (v) options operator, arrow; rhs((dsolve({equ1, equ2, subs(`ε` = v, [Bcs])[]}, numeric))(0)[3])/(1-`ϕ`)^2.5 end proc; plot(func, -1 .. 1, title = typeset((diff(f(eta), eta, eta))*versus*`ε`/(1-'`ϕ`')^2.5), titlefont = [times, italic, 18])

 

func2 := proc (v) options operator, arrow; rhs((dsolve({equ1, equ2, subs(`ε` = v, [Bcs])[]}, numeric))(0)[5]) end proc; plot(func2, -1 .. 1, title = typeset((diff(theta(eta), eta))*versus*`ε`), titlefont = [times, italic, 18])

 

``


 

Download plot.mw

 

Maple gives me poor accuracy for simple arithmetic problems.

My calculator gives me 0.7 and 0.9 respectively.

I guess I could do evalf(convert(2.59/3.7,fraction)) = 0.7000000000

but this seems excessive.

Sup Brainiacs

I found this "solved" problem for backing rectangular items efficiently. The first example (without rotation of items).

https://www.researchgate.net/file.PostFileLoader.html?id=59938b635b49527571134c33&assetKey=AS%3A527790252490752%401502846410852

I think I have the constraints the same as the text, but I get different solutions. my questions are at the end of the worksheet. Thanks

I 2D_BP.mw

How does one look through the Maple 2017 help for a code structure. I have been looking for how to construct a While loop and I am not able to find any example that implements while loop nor the code structure for it.

Any help is appreciated.

I'm getting the message

Secure Connection Failed

when trying to get into this forum using the latest version of Firefox.
For writing this I used Internet Explorer.

The message continues with:

An error occurred during a connection to mapleprimes.com. You have received an invalid certificate. Please contact the server administrator or email correspondent and give them the following information: Your certificate contains the same serial number as another certificate issued by the certificate authority. Please get a new certificate containing a unique serial number. Error code: SEC_ERROR_REUSED_ISSUER_AND_SERIAL
    The page you are trying to view cannot be shown because the authenticity of the received data could not be verified.
    Please contact the website owners to inform them of this problem
.

 

Has anybody else experienced the same?
What to do? 

Please, anyone with usefull informations to hint me. Looking forward to your favorable response. 

Thanking you in advance

Hello,

how could I force Maple to solve this pde

pdsolve(y*(diff(f(x, y, z, w), x))-x*(diff(f(x, y, z, w), y))+w*(diff(f(x, y, z, w), z))-z*(diff(f(x, y, z, w), w)) = 0)

When I hit enter, it happens nothing.

Matrix([[xx3[1,2],xx3[1,3]],[xx3[2,2],xx3[2,3]]])
Matrix(2, 2, {(1, 1) = (1/6)*sqrt(3)+(1/2)*I, (1, 2) = (1/6)*sqrt(3)-(1/2)*I, (2, 1) = (1/6)*sqrt(3)-(1/2)*I, (2, 2) = (1/6)*sqrt(3)+(1/2)*I})
expect output be
but these example are wrong
((1/6)*sqrt(3)+(1/2)*I)*Matrix([[1,-I].[-I,1]])
but these example are wrong
MatrixMatrixMultiply(Matrix([[(1/6)*sqrt(3),(1/2)],[(1/6)*sqrt(3),(1/2)]]),Matrix([[1,-I],[1,I]]));
concept like this output

How do you evaluate a function expression at given variable values?

f := 2*sin(x)-x^2/10.00:

eval(f, x = 5) gives 2*sin(5)-2.500000000

evalf(eval(f, x = 5)) gives -4.417848549

So is there a way we can just numerically calculate f(x) where x = 5?

Below is MAPLE code to simplify a series.  MAPLE expresses the result in terms of functions which many people are not familiar with.  Is there a way to express the answer in terms of more conventional functions expecially if N is a positive integer?


 

Cn := ((-I)*(1/2))*(2*(I*Pi*n*tau-(2*I)*Pi*n)*cos(Pi*n*tau/T)-T*(2*I)*sin(Pi*n*tau/T)+(4*I)*Pi*n)/(Pi^2*n^2); S4 := a[0]+sum(Cn*sin(2*Pi*n*x/T), n = 1 .. k); a[0] := 0; T := 4; tau := 2; Cn; S5 := unapply(S4, k, x); T := simplify(S5(N, x))

convert(T, StandardFunctions);

(-polylog(2, exp(-((1/2)*I)*Pi*(x-1)))*N^2-exp(-((1/2)*I)*Pi*N*(x+1))*LerchPhi(exp(-((1/2)*I)*(x+1)*Pi), 2, N)*N^2+polylog(2, exp(((1/2)*I)*(x+1)*Pi))*N^2+exp(-((1/2)*I)*Pi*N*(x-1))*LerchPhi(exp(-((1/2)*I)*Pi*(x-1)), 2, N)*N^2+polylog(2, exp(-((1/2)*I)*(x+1)*Pi))*N^2-exp(((1/2)*I)*Pi*N*(x+1))*LerchPhi(exp(((1/2)*I)*(x+1)*Pi), 2, N)*N^2-polylog(2, exp(((1/2)*I)*Pi*(x-1)))*N^2+exp(((1/2)*I)*Pi*N*(x-1))*LerchPhi(exp(((1/2)*I)*Pi*(x-1)), 2, N)*N^2+exp(((1/2)*I)*Pi*N*(x+1))-exp(((1/2)*I)*Pi*N*(x-1))+exp(-((1/2)*I)*Pi*N*(x+1))-exp(-((1/2)*I)*Pi*N*(x-1))-I*exp(-((1/2)*I)*x*Pi*N)*LerchPhi(exp(-((1/2)*I)*x*Pi), 1, N)*N^2*Pi-I*ln(1-exp(-((1/2)*I)*x*Pi))*N^2*Pi+I*exp(((1/2)*I)*x*Pi*N)*LerchPhi(exp(((1/2)*I)*x*Pi), 1, N)*N^2*Pi+I*ln(1-exp(((1/2)*I)*x*Pi))*N^2*Pi-I*exp(((1/2)*I)*x*Pi*N)*N*Pi+I*exp(-((1/2)*I)*x*Pi*N)*N*Pi)/(N^2*Pi^2)

(1)

``


 

Download simplify.mw

How can I specify curve of function  when drawing set of functions in  one plot ,( curve1, curve2,...) in Legend

When attempting to numerically solve for a function using fsolv it is possible that the function has multiple roots.  So to focus on a particular region you specify a range such as:

xmax := fsolve(S, x = 0 .. 1/2)

Is it possible fsolve may not resolve the solution due to the fact that delta x is not small enough or does fsolve autonomously adjust delta x in order to find the solution?  If not, how do you manually dictate the delta x for the interval specified?

Hello,

Could you please help me with the following problem? I'm new to Maple and I need some help.

Find the real solutions of the following system(it's a nonlinear system):

x^2 + y^2 + z^2 = 4

x + y + z = 0

x*sin(y*z) = -1

Request help in solving the equation:

Eq_H_1 := C[11*m]*(int((t-(j-1)*T)*alpha*(-alpha*(K*T*beta+T*beta*j-K*T-T*beta-T*j-beta*t+T+t))^(-beta/(beta-1))*exp(-R[m]*t), t = (j-1)*T .. (K+j-1)*T))

 

Thanks

Hello,

Could you please help me with the following problem? I'm new to Maple and i need some help.

Solve the equation x^3 - a*x + 1 = 0 , in x. Determine the particular solution for a=1,2,... .Graphically represent the polynom that appears in the equation, in a case where the equation has a real root and in a case where the equation has 3 real roots.

Thank you !

First 896 897 898 899 900 901 902 Last Page 898 of 2433