I tried to solve the following equation to find an explicit expression for x in terms of n. But the answer is a long relation in terms of RootOf and _Z. 

(1/4)*((-x^2+1)*(-4*Pi*(x^2-1)+n*(x^4-2*x^2+5))*cosh(n*x)*sinh(Pi*x)+sinh(n*x)*((-x^2+1)*(-4*n*(x^2-1)+Pi*(x^4-2*x^2+5))*cosh(Pi*x)-2*x*(x^4-2*x^2-3)*sinh(Pi*x)))/(x^2-1)^2=0

Any comment is welcomed.

F(0) := a; F(1) := b; F(2) := c; F(3) := d

for k from 0 to 1 do F(k+4) := -(N[1]*G(k)+Re*(sum(F(k-m)*(m+1)*(m+2)*(m+3)*F(m+3), m = 0 .. k))-Re*(sum((k-m+1)*F(k-m+1)*(m+1)*(m+2)*F(m+2), m = 0 .. k)))/((1+N[1])*(k+1)*(k+2)*(k+3)*(k+4)) end do

How to plot a graph for this equation with different values of N_1 and Re number

Dear Users!

Hoped everyone fine here. I have three main questions regarding the maple code given bellow:

restart; with(LinearAlgebra); with(plots);

alpha := 1; beta := 1; theta := 1/2;

UU := sinh(x)*sinh(y)*sinh(z)*exp(-1.*t);

NN := 3; L := 0; R := 1; T := 1; N := NN; Mx := NN; My := NN; Mz := NN; `Δx` := (R-L)/Mx; `Δy` := (R-L)/My; `Δz` := (R-L)/Mz; `Δt` := (R-L)/N;

kappa[1] := 1; kappa[2] := 2/x^2; kappa[3] := 1/x^2; kappa[X] := x^2+y^2+z^2+1; kappa[Y] := x^2+y^2+z^2+1; kappa[Z] := x^2+y^2+z^2+1; kappa[4] := 0; NL := 3;

ics := [seq(seq(seq([u[i, j, k, 0] = eval(UU, [x = i*`Δx`, y = j*`Δy`, z = k*`Δz`, t = 0]), u[i, j, k, -1] = eval(u[i, j, k, 1]-2*`Δt`*(eval(diff(UU, t), t = 0)), [x = i*`Δx`, y = j*`Δy`, z = k*`Δz`, t = 0])][], i = 0 .. Mx), j = 0 .. My), k = 0 .. Mz)];

bcs := [seq(seq(seq([u[0, j, k, n] = eval(UU, [x = 0, y = j*`Δy`, z = k*`Δz`, t = n*`Δt`]), u[Mx, j, k, n] = eval(UU, [x = L, y = j*`Δy`, z = k*`Δz`, t = n*`Δt`])][], j = 0 .. My), k = 0 .. Mz), n = 1 .. N), seq(seq(seq([u[i, 0, k, n] = eval(UU, [x = i*`Δx`, y = 0, z = k*`Δz`, t = n*`Δt`]), u[i, My, k, n] = eval(UU, [x = i*`Δx`, y = L, z = k*`Δz`, t = n*`Δt`])][], i = 1 .. Mx-1), k = 0 .. Mz), n = 1 .. N), seq(seq(seq([u[i, j, 0, n] = eval(UU, [x = i*`Δx`, y = j*`Δy`, z = 0, t = n*`Δt`]), u[i, j, Mz, n] = eval(UU, [x = i*`Δx`, y = j*`Δy`, z = L, t = n*`Δt`])][], i = 1 .. Mx-1), j = 1 .. My-1), n = 1 .. N)];
Sol := {u[1, 1, 1, 1] = 0.2366497936e-1, u[1, 1, 1, 2] = 0.7589975856e-2, u[1, 1, 1, 3] = 0.6029906475e-3, u[1, 1, 2, 1] = 0.3778786317e-1, u[1, 1, 2, 2] = 0.7126415819e-2, u[1, 1, 2, 3] = -0.1197885714e-2, u[1, 2, 1, 1] = 0.3778786315e-1, u[1, 2, 1, 2] = 0.7126415820e-2, u[1, 2, 1, 3] = -0.1197885718e-2, u[1, 2, 2, 1] = 0.6038763054e-1, u[1, 2, 2, 2] = 0.4264591907e-2, u[1, 2, 2, 3] = -0.3509477851e-2, u[2, 1, 1, 1] = 0.3171958616e-1, u[2, 1, 1, 2] = -0.1327161715e-1, u[2, 1, 1, 3] = -0.4628647419e-2, u[2, 1, 2, 1] = 0.4979852397e-1, u[2, 1, 2, 2] = -0.3060811899e-1, u[2, 1, 2, 3] = -0.344914876e-4, u[2, 2, 1, 1] = 0.4979852397e-1, u[2, 2, 1, 2] = -0.3060811898e-1, u[2, 2, 1, 3] = -0.3449150010e-4, u[2, 2, 2, 1] = 0.7882396741e-1, u[2, 2, 2, 2] = -0.6192340018e-1, u[2, 2, 2, 3] = 0.1156615222e-1}

Using set of points given in ics, bcs and Sol

1. I want to contruct a vector at any time level (by fixing fourth suffix like u[i,j,k,n]) for i = 0..Mx,j=0..My,k=0..Mz and then find its L2 and L[infinity] norms.

2. Next I want contruct a vector by fixing two suffixes like u[i,j,k,n]) for i = 0..Mx,j=0..My and plot a surface in 3D

3. Finally I want to construct a vector by fixing three suffixes like u[i,j,k,n]) for i = 0..Mx, and plot a curve in 2D.

I'm waiting for your positive respone. I shall be very thankfull to you in advance.

Special request to:
@acer @Carl Love @Kitonum @Preben Alsholm

Hi everyone, I have a mechanical problem, I solved it using hybrid method (differential transformation method for time and finite difference method for space) I need some help please.

How do i pipe maple output to the clipboard and vice versa (assigning clipboard content to a variable declared in a maple worksheet)?

1. The member() function doesn't find the value in a table if that value is associated with the key (aka index)  '0'
2. Doesn't seem to matter if the table is created explicitly (ie by calling table()) or implicitly (ie by simple indexed assignment)
3. member() doesn't seem to have the same issue with a zero-based Array()
4. I'm pretty sure that this is a bug, but I'd like some opinions before I report it
5. This behaviour has been around for a while: same thing happens all the way back to Maple 18. (I can't check anything earlier)

Check the ouput of member(10,t1) in the attached

Download memberProp.mw

Recently I discovered you can get the wolfram language and Mathematica for free if you buy a Raspberry Pi. 

Will there be something similar with Maple?

I'm guessing since Mathematica has secured a contract with Raspberry Pi, that would exclude competitors, ie Maple

Hi everybody 

l I tried to install Maple 2018.2 on the golden master of Catalina and it didn't work : the installation process ended after the entering of the password to authorize the installation. In fact Maple 2018.2 still contains 32 bit elements. Is there a solution ? Thank you

David

Second terms coefficient = 2/3*x when (....)^2

9/16*x^2 when (...)^3

64/125*x^3 when (...)^4

625/1296*x^4

Which Bernstein has these result?

Which limit function of x^n = (n^(n-1))/((n+1)^(n-1))*x^(n-1)?

how to find and calculate?

Hi,

    This is not an answer. I can't ask question, since I can't add any 'tags'! This problem appeared long time ago, but I hadn't solved it until now. So I ask here.

    I have a question about the function coeffs(). I don't know why it cann't get coefficients for me.

question.mw

Thanks.

Variable disappear after read file and then pass parameter to create thread

Hello people in mapleprimes,
I have a question.

I use maple2019 with mac os 10.14.6.

$ sw_vers
ProductName:	Mac OS X
ProductVersion:	10.14.6
BuildVersion:	18G95

With maple2019, errors appears.

> kernelopts(version);
        Maple 2019.1, APPLE UNIVERSAL OSX, Jun 6 2019, Build ID 1403154

> assume(a>-1,b>0);
> additionally(a<=1);
> about(a);
Originally a, renamed a~:
  is assumed to be: FAIL

> assume(tau<1,tau>0,s<1,s>0):
> a_e1:=tau*s*(1+tau)<tau*s+tau+s-1:
> b_e2:=expand(lhs(a_e1)-rhs(a_e1))<0:
>  b_e3:=collect(b_e2,s,factor):
> solve(b_e2,s) assuming tau<1;
Error, (in assuming) when calling 'property/ConvertProperty'. Received: 'FAIL
is an invalid property'

On the other hand, with maple2018, they do not.

> kernelopts(version);
       Maple 2018.2, APPLE UNIVERSAL OSX, Nov 16 2018, Build ID 1362973

> assume(a>-1,b>0);
> additionally(a<=1);
> about(a);
Originally a, renamed a~:
  is assumed to be: RealRange(Open(-1),1)

> assume(tau<1,tau>0,s<1,s>0):
> a_e1:=tau*s*(1+tau)<tau*s+tau+s-1:
> b_e2:=expand(lhs(a_e1)-rhs(a_e1))<0:
> b_e3:=collect(b_e2,s,factor):
> solve(b_e2,s) assuming tau<1;
                                    1
                               [{-------- < s~}]
                                 1 + tau~

Are errors due to some bugs in maple2019 when being used with mac os 10.14.6?

Take care.

taro

Hello,

I'm working in MapleSim to model a mechanism called Schatz which has a six-Revolut, known whit TURBULA name.
I modeled the mechanism in CREO to have a CAD part to use in maplesim.
everything is okay when I run the simulation without a speed constant. as shown in the picture above.
when I add a speed source, MapleSim finds an error (Internal error occurred during the simulation).
Can someone help me please?

Hi everybody!

Quick question for you guys. I have access to Maple 2019 on my main computer and Maple 2018 on a remote Desktop (but I don't want to use because it's super slow and bad (normal rdp stuff)). The problem I have is for the same document, same line, same instruction, Maple 2019 on my computer is failling to do it correctly (unable to do the matrix correctly) and in the Remote Maple 2018 it's doing it just fine. I really don't understand. If someone can give me a little help it would be highly apreciated!

Maple 2018 (what it's supposed to do):

Maple 2019 (what it is not suppose to do):