MaplePrimes Questions

Could be helpful to see again those worksheets about programming in Maple.

It was back in 2000-2004 that have seen these workssheets on the Maple website

 I looked at Maple website, but nothing to find anymore for those programming worksheets?

conj := conjugate; d := a*x+b*y-c = 0; z := x+I*y; evalc(z+conj(z)); evalc(z-conj(z)); d := expand((1/2)*a*(z+conj(z))+b*(z-conj(z))/(2*I))-c; is(d = z(a-I*b)+conj(z)*(a+I*b)-2*c); varpi = a+I*b; is(d = z*conj(varpi)+conj(z)*varpi-2*c); How to perform calculations correctly ? Thank you.

Unfortunately, I can't solve the problem in MAPLE that appears in the last line of the picture (0=1 instead of a0=1) unless I write a0 seperately and add the partial sum starting from k=1 to it.

I assume the problem is with 00 in the power series but since MAPLE defines 00:=1 I'm not sure about it. Anyone knows a better solution or can tell me if I made a mistake somewhere?


 

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

First 556 557 558 559 560 561 562 Last Page 558 of 2433