Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

How can I see 8 1/2 x 11" page margins to write a simple letter size report?

Hello everyone,

first, I'd like to mention that I am relatively new to Maple and am therefore thankful for any advice you might have!

I am trying to integrate the term (k_1^2 * r) from a to infinity, see the picture below as well as the attached file. Maple seems to have some issues with that. However, if I break the integral down into more manageable parts it suddenly works! Why is that? How can I get Maple to solve this immeadiately?  I suspect the culprit lies in the term that contains (-Ei(-B*r)*r^(-1)) where Ei is the exponential integral as defined in Maple. The resulting  hypergeometric function seems suspicious.

The problem is that I have to evaluate 21 integrals of this type (k_x*k_y*r) and breaking them down manually becomes pretty cumbersome, especially as the number of terms in the expanded expressions increases. Is there a way to automate this procedure? I guess I would need to extract individual terms and automatically plug them into the integral expression. That should the last resort, however.

The specific problem (everything included for context, weird stuff happens after equation 15):

Maple_Problem.mw

As for the variables: E, t, and R are real positive numbers. a and B are already assumed as real and positive. A_0, C_0, A_2, and C_2 are real numbers (could be negative, I do not know yet since they must be determined later on). a_0 is definitely real, but it may be negative. r is the polar coordinate, so it is also real and positive, but adding this assumption did not yield a better result.

 

Thank you for your help!

 

How to use random variables for different parameters in Maple to find min and max value, like i have R1 min, nom and max value,

I have ratio of R1/R2, temperature coefficient for R1 for min and max value?

I've trying to compose operators in Maple and I've hit a snag. Here's a simplified example:

dx := f -> diff(f, x);
dy := f -> diff(f, y);
eq1 := x * dy
eq2 := y * dx

eq1(eq2)

The output of the final statement is x(y*dx)dx. It took me a while to figure out that Maple is treating it like so:

(x * dy)(y * dx)
(x(y * dx)) * (dy(y * dx))
(x(y * dx)) * (dx)
x(y * dx) * dx

The problem is I don't want it applying anything to x, the expression should reduce to x * dx. Same as if x was replaced by a constant, for example:

(4 * dy)(y * dx) = 4 * dx

In the context of my full code, Maple already knows x is not "appliable" e.g. is(x::appliable) returns false. How do I tell Maple to treat x as not-a-function?

Hello everyone,

How would I do something like this in Maple:

i.e. do the product (or sum) from k=1 to n, skipping the value i. This one has me really scratching my head...

Thanks for your help!


 

``

lambda := .3:

omega := lambda+mu+xi:

alpha := 2*sqrt(lambda*mu)

.9165151390

(1)

``

B[1] := BesselI(k-1, alpha*(u-y))

BesselI(k-1, .9165151390-.9165151390*y)

(2)

B[2] := BesselI(k+1, alpha*(u-y))

BesselI(k+1, .9165151390-.9165151390*y)

(3)

``

F := evalf(Int(sum((B[1]-B[2])*exp(-omega*(u-y)), k = 1 .. infinity), y = 0 .. u))

Int(sum((BesselI(k-1., .9165151390-.9165151390*y)-1.*BesselI(k+1., .9165151390-.9165151390*y))*exp(-1.200000000+1.200000000*y), k = 1 .. infinity), y = 0. .. 1.)

(4)

``

``

``


 

Download int.mw

when i write :

I want answer me   '2'

but it answer on this way:

what i should do for it answer just '2' and ynderstand dont need conjugate 'x and lambda'

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

First 543 544 545 546 547 548 549 Last Page 545 of 2217