MaplePrimes Questions

How to get the inflection points for this function.?

fprime_expr:=x^sin(x)*(cos(x)*ln(x)+sin(x)/x);

i tried symbolic and numeric , but no answer 

for X:= solve( fprime_expr =0, x);

 

This worksheet plots the non-overlapped images of two spheres, but plotted procedure p fails to successfully pass the choice of color to the implicitplot3d command.

How can this be accompished?

Sphere_exclusions.mw

Hi there!

I am trying to express a given polynomial not in terms of the monomes, but in terms of some orthogonal polynomials p_i, i=1..n.

If I have let's say x^5+x^4+2x^2 then I first divide that term by p_5, then the remainder by p_4 and so on, until I have the wanted coefficients. My problem is, that Maple complains about my code:

Error, (in Basenwechsel) invalid input: op expects 1 or 2 arguments, but received 0.

NEUZMinus:= proc(Unten, Oben, f,G,Liste,n)::real;
  #Unten:= Untere Intervallgrenze; Oben:= Obere Intervallgrenze; g:= zu integrierende Funktion;
  #G:= Gewicht; n:= Hinzuzufügende Knoten;
 
  Basenwechsel:=proc(Dividend, m);
 
  print(Anfang,Dividend,p[m]);
  Koeffizient:=quo(Dividend, p[m],x);
  print(Koeffizient);
  Rest:=rem(Dividend, p[m],x);
  print(Rest);
  if m=0 then
    Basenwechsel:=[Koeffizient];
  else
    
    Basenwechsel:=[Koeffizient,op(Basenwechsel(Rest,m-1))];
   print(Basenwechsel);
  end if;
  print(Durchlauf)
  end proc;
p[-1]:=0;
p[0]:=1;
for i from 1 to max(n,numelems(Liste)) do
  p[i]:=x^i-add(int(x^i*p[j]*diff(G,x),x=Unten..Oben)*p[j]/int(p[j]^2*diff(G,x),x=Unten..Oben),j=0..i-1);
  print(p[i]);
end do;
print(Liste[1],numelems(Liste));
Hn:=mul(x-Liste[i],i=1..numelems(Liste));
print(Hn);
 Koeffizienten:=Basenwechsel(Hn,numelems(Liste));
print(Koeffizienten)
end proc
 

Funnily enough, when I add a print command before the op-function, it does work:

 

NEUZMinus:= proc(Unten, Oben, f,G,Liste,n)::real;
  #Unten:= Untere Intervallgrenze; Oben:= Obere Intervallgrenze; g:= zu integrierende Funktion;
  #G:= Gewicht; n:= Hinzuzufügende Knoten;
 
  Basenwechsel:=proc(Dividend, m);
 
  print(Anfang,Dividend,p[m]);
  Koeffizient:=quo(Dividend, p[m],x);
  print(Koeffizient);
  Rest:=rem(Dividend, p[m],x);
  print(Rest);
  if m=0 then
    Basenwechsel:=[Koeffizient];
  else
    print(Basenwechsel(Rest,m-1));
    Basenwechsel:=[Koeffizient,op(Basenwechsel(Rest,m-1))];
   print(Basenwechsel);
  end if;
  print(Durchlauf)
  end proc;
p[-1]:=0;
p[0]:=1;
for i from 1 to max(n,numelems(Liste)) do
  p[i]:=x^i-add(int(x^i*p[j]*diff(G,x),x=Unten..Oben)*p[j]/int(p[j]^2*diff(G,x),x=Unten..Oben),j=0..i-1);
  print(p[i]);
end do;
print(Liste[1],numelems(Liste));
Hn:=mul(x-Liste[i],i=1..numelems(Liste));
print(Hn);
 Koeffizienten:=Basenwechsel(Hn,numelems(Liste));
print(Koeffizienten)
end proc

Does anyone know how I can get rid of the print command without getting no result? And how does a simple print command change the outcome of an algorithm in the first place?

 

Thank you in advance,

Daniel Reksten

 

Why are units disappearing in this expression?

The answer should be 0 kN, not 0.

Need help with write calculations in maple and ploting

The array ( a 3 by 4 )  has elements:

1,1,1,1

2,4,8,16

3 ,9,27,81

How to do this a nested loop ? : how to reasoning this in Maple with do loops 

Hello,

While working for an assignment I had to use a piecewise function which gives a result based on a randomly generated value. The following is a much simpler version of what I've been working on, but it gives the same error.

A random value is uniformly distributed between 0 and 100. The piecewise returns a 1 if the value is between 0 and 50 and it returns 2 if the value is between 50 and 100. As far as I understand this function should only ever give 1 or 2, never something else. However when I loop this a few times Maple regularly returns a 0, which doesn't make sense to me. I've printed the values that return a 0 but none of these should break the piecewise. Can someone please explain to me what's going wrong here and how to fix it?

My code:

restart; randir := piecewise(0 <= r1() and r1() < 50, 1, 50 <= r1() and r1() < 100, 2);
for i to 1000 do r1 := rand(0. .. 100.0); if randir = 0 then print(fail[i], r1()) end if end do;
 

I've included a failure check to see when and at what values it returns a 0, and as you can see it happens very often.

Hi there.

I noticed that in Maple 2020 default font size of document (Times New Roman 12) looking bigger than in previous versions. How can I make the view the same as in previous versions?

Thank you

What does the error message, "Error, (in monspt) invalid input: ArrayTools:-Reverse expects its 1st argument, A, to be of type rtable, but received [[0, 4], [1, 3], [2, 2], [3, 1], [4, 0], [0, 3], [1, 2], [2, 1], [3, 0], [0, 2], [1, 1], [2, 0], [0, 1], [1, 0], [0, 0]]
" exactly mean, and how to find the required rtable type, and how to make it work? Thanks a lot!

The plot is not showing two lines perpendicular ?

Can anyone help with ploting the movement of the body along the segment AB, BC this task?

restart;
alpha:=Pi/6: V[A]:=1: l:=3: f:=0.2: d:=2.5:
# The movement of the body along the segment AB
de:=m*diff(x1(t),t,t)=G*sin(alpha)-F[fr]:
G:=m*g: F[fr]:=m*g*f*cos(alpha): g:=9.81:
ics:=x1(0)=0, D(x1)(0)=V[A]:
de:=evalf(de/m);
B:=rhs(de);
sol:=dsolve({de,ics}, x1(t)); # The solution of de
fsolve(rhs(sol)=l);
tau:=evalf[3](%[2]);
V['B']:=B*tau+V[A];
 

diff(diff(x1(t), t), t) = 3.205858157

 

3.205858157

 

x1(t) = (3205858157/2000000000)*t^2+t

 

-1.715094465, 1.091236545

 

1.09

 

4.494385391

(1)

# The movement of the body along the curve BC
sys:=m*diff(x(t),t,t)=0, m*diff(y(t),t,t)=m*g;
ics_sys:=x(0)=0, y(0)=0, D(x)(0)=V['B']*cos(alpha), D(y)(0)=V['B']*sin(alpha):
Sol:=evalf(dsolve({sys, ics_sys})); # The solution of sys
T:=evalf(solve(eval(x(t),Sol)=d));
h:=evalf(eval(eval(y(t),Sol),t=T));
T:=evalf[2](T);
h:=evalf[3](h);   
 

m*(diff(diff(x(t), t), t)) = 0, m*(diff(diff(y(t), t), t)) = 9.81*m

 

{x(t) = 3.892251925*t, y(t) = 4.905000000*t^2+2.247192696*t}

 

.6423016927

 

3.466940605

 

.64

 

3.47

(2)

 


 

Download sol_Maple.mw

Hi there.

How supress this warning in Maple 2020 when I use "e" just like a variable not an exponent:

Warning, if e is meant to be the exponential e, use command/symbol completion or palettes to enter this special symbol, or use the exp function

?

Thank you

Hi there.

I found bugs in covariant differentiation of vectors that defined in another coordinate system. I attached three documents - two of them with bugs and one with normal behavior.

In cov_diff_bug_1.mw I defined Killing vector in Minkowski spacetime. Then I transformed metric and this vector to spherical coordinate system and calculated covariant derivative of vector and obtained symmetric tensor matrix. But for Killing vector it should be antisymmetric.

In cov_diff_bug_2.mw I set metric in spherical coordinate system by hand instead of transformation and obtained zero covariant derivative.

And in cov_diff_no_bug.mw I set metric and vector by hand without any transformation and obtained normal results - antisymmetric covariant derivative and zero sum of two covariant derivatives with rearranged indexes.

So in Maple 2020 coordinate transformation don't doing well in combining with covariant differentiation. But in previous versions of Maple this behavior could be opposite)

I'm asking developers to consider this issues.

Thank you.

cov_diff_bug_1.mw

cov_diff_bug_2.mw

cov_diff_no_bug.mw

 

Dears,
Someone to know install modules as pyautogui? 
ImportModule("pyautogui")
In my software gives the following message "Error, (in Python:-ImportModule) ModuleNotFoundError: No module named 'pyautogui' ['Traceback (most recent call last):\n', '  File "<string>", line 1, in <module>\n', "ModuleNotFoundError: No module named 'pyautogui'\n"]"
Or someone know about automation keyboard in Maplesoft?
Kind regards

 Dear Maple community, 
I am struggling with the numerical solution of the complex (linear) partial differential equation to obtain results. For example, I see such as this statement "Error, (in pdsolve/numeric/plot) unable to compute solution for t<HFloat(0.0): signum is not differentiable at complex values" pops up. I did my homework but I didn't get anything. I will be happy if someone could give me some suggestions or advice. Thank you in advance. 

First 431 432 433 434 435 436 437 Last Page 433 of 2308