Hi Everyone!

I would like your help again.

Considering a Free-Pinned-Pinned-Free beam (page 88 in the pdf file). In case of a matrix 12x12 how could I find the coefficients (a1, a2, a3, a4, b1, b2, b3, b4, c1, c2, c3, c4) using MAPLE in order to plot the mode shapes of the Figure 3.22 (a) (page 70 in the pdf file)? in case of a matrix 16x16 and 20x20, the procedure is the same?

I tried to plot the mode shapes but I failed because I believe they should be similar to Figure 3.22 (a) (page 70 in the pdf file)

sys_GE := {-a1 + a3 = 0, -a2 + a4 = 0, 0.9341161484*a1 + 0.3569692163*a2 + 1.519943120*a3 + 1.819402948*a4 = 0, 0.6669014188*b1 - 0.7451459573*b2 + 5.530777989*b3 + 5.620454178*b4 = 0, 0.9938777922*c1 - 0.1104849953*c2 + 62.17096851*c3 + 62.17901032*c4 = 0, -0.9341161484*a1 - 0.3569692163*a2 + 1.519943120*a3 + 1.819402948*a4 + 0.9341161484*b1 + 0.3569692163*b2 - 1.519943120*b3 - 1.819402948*b4 = 0, -0.3569692163*a1 + 0.9341161484*a2 + 1.819402948*a3 + 1.519943120*a4 + 0.3569692163*b1 - 0.9341161484*b2 - 1.819402948*b3 - 1.519943120*b4 = 0, 0.3569692163*a1 - 0.9341161484*a2 + 1.819402948*a3 + 1.519943120*a4 - 0.3569692163*b1 + 0.9341161484*b2 - 1.819402948*b3 - 1.519943120*b4 = 0, -0.7451459573*b1 - 0.6669014188*b2 + 5.620454178*b3 + 5.530777989*b4 + 0.7451459573*c1 + 0.6669014188*c2 - 5.620454178*c3 - 5.530777989*c4 = 0, -0.6669014188*b1 + 0.7451459573*b2 + 5.530777989*b3 + 5.620454178*b4 + 0.6669014188*c1 - 0.7451459573*c2 - 5.530777989*c3 - 5.620454178*c4 = 0, 0.7451459573*b1 + 0.6669014188*b2 + 5.620454178*b3 + 5.530777989*b4 - 0.7451459573*c1 - 0.6669014188*c2 - 5.620454178*c3 - 5.530777989*c4 = 0}

solve(sys_GE, {a1, a2, a3, a4, b1, b2, b3, b4, c1, c2, c3, c4});

{a1 = 0.1014436637 c4, a2 = -0.1143870369 c4,

a3 = 0.1014436637 c4, a4 = -0.1143870369 c4,

b1 = 0.07095134140 c4, b2 = -0.1260395712 c4,

b3 = 0.1510591164 c4, b4 = -0.1737777468 c4,

c1 = 0.2102272829 c4, c2 = -0.2816561313 c4,

c3 = -1.003990622 c4, c4 = c4}

26_06_2020_Transcedental_equation_matrix_12x12_artigo_2.mw

Hello

I want to know if there is any option to convert the Maple code into MATLAB version to work across platforms. I have tried the code generation option in Maple and I could not understand how to use the single line code generated. I am working on an already written Maple code and conversion of code to use on both platforms is kind of important. Please advice me on the same.

Regards!

We have just run into a problem (loss of data), when 2 users have opened the same workbook.

Does Maple support Multiuser access to a workbook, and if not - why isn't there any lock mechanism in place to prevent that?

I've asked about this 5 years ago

https://www.mapleprimes.com/questions/203628-Maple-2015-No-Longer-Autoscroll-As

And according to the reply there at the time, this is supposed to have been fixed in Maple 2015.1:

**It appears that the scrolling issue has been fixed in the next version of Maple**

But ever since then, whenever I use Maple, the worksheet do not auto-scroll as the program runs and prints outout to the screen. I have to keep using the mouse to move the bar down to see the latest prints.

Not only that, the bar itself (on the right side) shrinks in size as more outpout is printed.

It becomes smaller and smaller with time, and becomes very hard to grab it with the mouse since it become so thin. Once all the output is deleted, the bar become long again.

So I think there is an option to enable autoscrolling somewhere.

But I am not able to find it. I looked at options, display, general and interface and see nothing there to turn it on.

Sometimes I get the feeling that the folks at Maplesoft do not use Maple interface the way we users use it, becuase this is such a basic interface problem and any one who uses Maple interface for more than one day must see these problems all the time.

So how come these are still not fixed? Here is a screen shot showing the bar on the right after long print out, and showing how it was before.

I use only worksheet mode, not document mode. On windows 10.

Hi, I want to use Maple to derive the eigenvectors of a 3x3 matrix A. I can obtian the 3 eigenvalues correctly by using eigenvalues(A), and the results are "u, -a+u, a+u". But after inputing eigenvectors(A), I can only get "[-a+u, 1, {r}], [a+u, 1, {r}], [u, 1, {r}]". All the 3 eigenvectors are written as {r}, I dont know why...and how can I tranform {r} into visiable results?

I would be very appreciate if anyone could give me some help!

it took a while to track this. When I run these 3 ode's in this sequence, Maple gives

**Error, (in Engine:-Dispatch) too many levels of recursion**

and it has to be this sequence. I also found I might have to run the code 2 or 3 times to get this error if it does not show first time.

Will show a movie below also.

Here is the code

restart;
ode1:=cos(y(x)) - x*sin(y(x))*diff(y(x), x) = sec(x)^2:
ode2:=y(x)*sin(x/y(x)) + x*cos(x/y(x)) - 1 + (x*sin(x/y(x)) - x^2*cos(x/y(x))/y(x))*diff(y(x), x) = 0:
ode3:=x/(x^2 + y(x)^2) + y(x)/x^2 + (y(x)/(x^2 + y(x)^2) - 1/x)*diff(y(x), x) = 0:
func:=y(x);
the_ic:={};
the_time_limit:=60;
advisor:=timelimit(the_time_limit,DEtools:-odeadvisor(ode1)):
solution_by_CAS := timelimit(the_time_limit,:-dsolve([ode1,y(0)=0],func,singsol=all)):
result_of_ode_test := timelimit(the_time_limit,:-odetest~([solution_by_CAS],ode1)):
advisor:=timelimit(the_time_limit,DEtools:-odeadvisor(ode2)):
solution_by_CAS := timelimit(the_time_limit,:-dsolve([ode2,y(0)=0],func,singsol=all)):
result_of_ode_test := timelimit(the_time_limit,:-odetest~([solution_by_CAS],ode2)):
advisor:=timelimit(the_time_limit,DEtools:-odeadvisor(ode3)):
solution_by_CAS := timelimit(the_time_limit,:-dsolve([ode3,y(1)=0],func,singsol=all)):
result_of_ode_test := timelimit(the_time_limit,:-odetest~([solution_by_CAS],ode3));

And here is a movie (this time it happend on second run) sometimes it happens on first run. so it is random.

Maple 2020.1

Why does it happen and is there a work around?

Noticed few errors generated on some ode's. I do not know if these are expected or known, I thought may be someone might want to check them.

restart;
ode1:=(a^2*x+y(x)*(x^2-y(x)^2))*diff(y(x),x)+x*(x^2-y(x)^2) = a^2*y(x);
ode2:=x^7*y(x)^2*diff(y(x),x)^3+(1-3*x^6*y(x)^3)*diff(y(x),x)^2+3*x^5*y(x)^4*diff(y(x),x)-x^4*y(x)^5 = 0;
ode3:=diff(y(x),x)^4+x*diff(y(x),x)-3*y(x) = 0;

sol:=timelimit(60,:-dsolve(ode1,y(x))):

Error, (in dsolve) numeric exception: division by zero

sol:=timelimit(60,:-dsolve(ode2,y(x))):
timelimit(60,:-odetest([sol][5],ode2))

Error, (in IntegrationTools:-Indefinite:-AlgebraicFunction) argument is not an algebraic

sol:=timelimit(60,:-dsolve(ode3,y(x))):
timelimit(60,:-odetest(sol[1],ode3));
timelimit(60,:-odetest(sol[2],ode3))

Error, (in odetest) received not an ODE w.r.t x(_T)

Error, (in odetest) received not an ODE w.r.t y(_T)

This last one is only due to how the solution is given. So I do not think it is a real problem.

Maple** 2020.1**, Physics **708**

windows 10

Maple Worksheet - Error

Failed to load the worksheet **/maplenet/convert/check.mw **.

Download check.mw

Dear All.

I hope we are all staying safe.

Please I want to solve Sine Gordon Equation using a numerical method I constructed (non-classical), I need to compare the result of the method with the exact solution to generate the errors. However, since the exact solution has two variables, x, and t, I really don't know how to accommodate the two in my coding.

Can someone be of help in this regard?

Thank you and kind regards

Download Discretization_of_Sine_Gordon_Equation.mw

Download Sine_Gordon_Implementation_Trial.mw

Hi there,

Can someone explain to me why is it happening? If smaller equal than coloured well, if smaller then just the equal line appears. Why?

with(plots);
dm1 := .23; pm12 := .7; dm12 := .23;
inequal(pm1*(dm12+am21)-pm12*(dm1+am21) <= 0, pm1 = 0 .. 1, am21 = 0 .. 1, color = "Nautical 1");

inequal(pm1*(dm12+am21)-pm12*(dm1+am21) < 0, pm1 = 0 .. 1, am21 = 0 .. 1, color = "Nautical 1")

Thank you!

Etude d'un cas particulier
a := 5: b := 7:
k := 9:
A := [a, 0]: B := [0, b]: #A et B fixes
P := [t, 0]: Q := [0, k/t]:#P et Q 2 points mobiles
cir := -a*x-b*y+x^2+y^2 = 0:
sol := solve(subs(y = 5, cir), x):
cen := [solve(diff(cir, x)), solve(diff(cir, y))]:
x0 := sol[1]: y0 := 5:
M := [x0, y0]:
R := sqrt(cen[1]^2+cen[2]^2):
beta := arctan(diff(solve(EQ(M, cen), y), x)):
Recherche des valeurs de t pour que les 2 droites soient perpendiculaires
eq := t^2*(y0-b)+t*(a*b-a*y0+b*x0-k)-x0*(a*b-k) = 0;
sol := solve(eq, t);
t := sol[1]; tp := sol[2];
P1 := [t, 0]; Q1 := [0, k/t];
PQ1 := simplify(x*(-a*b+b*t+k)+y*t*(t-a)-t*(-a*b+b*t+k)) = 0:#1ere tangente
PQ2 := simplify(x*(-a*b+b*tp+k)+y*tp*(tp-a)-tp*(-a*b+b*tp+k)) = 0:#2ième tangente
P2 := [tp, 0]; Q2 := [0, k/tp];
CIR := implicitplot(cir, x = -4 .. 8, y = -4 .. 12, color = red);
Fig := proc (alpha) local Dr1, DR1, Dr2, DR2, N, u0, v0, Po, t, tp, sol; global a, b, k, cen, R; u0 := cen[1]+R*cos(alpha); v0 := cen[2]+R*sin(alpha); N := [u0, v0]; sol := solve(t^2*(v0-b)+t*(b*u0-a*v0+a*b-k)-u0*(a*b-k) = 0, t); t := sol[1]; tp := sol[2]; Dr1 := simplify(x*(-a*b+b*t+k)+y*t*(t-a)-t*(-a*b+b*t+k)) = 0; DR1 := implicitplot(Dr1, x = -4 .. 8, y = -4 .. 12, color = brown); Dr2 := simplify(x*(-a*b+b*tp+k)+y*tp*(tp-a)-tp*(-a*b+b*tp+k)) = 0; DR2 := implicitplot(Dr2, x = -4 .. 8, y = -4 .. 12, color = pink); Po := pointplot([N[]], symbol = solidcircle, color = [black], symbolsize = 8); display([Po, DR1, DR2]) end proc;
DrPQ1 := implicitplot(PQ1, x = -4 .. 22, y = -4 .. 12, color = blue);
DrPQ2 := implicitplot(PQ2, x = -4 .. 22, y = -4 .. 12, color = blue);
Points := pointplot([A[], B[], M[], P1[], P2[], Q1[], Q2[], cen[]], symbol = solidcircle, color = [green], symbolsize = 10);
T := plots:-textplot([[A[], "A"], [B[], "B"], [M[], "M"], [P1[], "P1"], [P2[], "P2"], [Q1[], "Q1"], [Q2[], "Q2"], [cen[], "cen"]], font = [times, 10], align = {below, left});
n := 19;
display([seq(Fig(2*i*Pi/n), i = 0 .. n), Fig(beta), CIR, DrPQ1, DrPQ2, Points, T], scaling = constrained, size = [500, 500]);
I would find out the focus of the ellipse. Thank you.