Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I am unable to compute integral, it displays the function itself.

burrpdf := alpha*k*((x-gamma)/beta)^(alpha-1)/beta(1+((x-gamma)/beta)^alpha)^(k+1)

 

mean_burr_temp := simplify(x*burrpdf)

 

int(mean_burr_temp, x = gamma .. t)

Hello Everybody!
Please help me to find a way to write the Taylor's series of some function in factiorial or Gamma notation.
Please see the attached file and suggest me something.
taylor_series_general.mw

I've been using Maple 2018 only a few days now, was mostly using Maple 2016 and never had any issues with returning an output (at least one that wasnt my fault), but with 2018 I've seen a few times where a simple task would return an output of "__SELECTION" and then my input. No idea what this is or why its happening and all I really need to know is how to prevent it from happening. If its something with my preferences or settings that needs to be addressed, thats fine, but otherwise this is getting to be a real pain. See attached 

 

input -->

restart; with(plots);
H := 11;
B := x; W := x;
t__E := (13100*(B*B)*H*(.23-.10)*1.2)/(1500*alpha);
A__t := 2*B*W+2*H*(B+W);
OF := alpha/A__t;
for alpha from .1 by .1 to 2 do XX[alpha] := solve({0 < x, OF < 0.6951871658e-2, t__E < 200}, x) end do

result-->

         XX[0.1] := {0.3222041888 < x, x < 1.155224359}
         XX[0.2] := {0.6354894776 < x, x < 1.633733956}
         XX[0.3] := {0.9405579047 < x, x < 2.000907284}
         XX[0.4] := {1.238023973 < x, x < 2.310448718}
         XX[0.5] := {1.528429210 < x, x < 2.583160196}
         XX[0.6] := {1.812253750 < x, x < 2.829710218}
         XX[0.7] := {2.089925662 < x, x < 3.056436363}
         XX[0.8] := {2.361828525 < x, x < 3.267467912}
         XX[0.9] := {2.628307644 < x, x < 3.465673077}
         XX[1.0] := {2.889675191 < x, x < 3.653140182}
         XX[1.1] := {3.146214498 < x, x < 3.831445747}
         XX[1.2] := {3.398183646 < x, x < 4.001814567}
         XX[1.3] := {3.645818516 < x, x < 4.165220660}
         XX[1.4] := {3.889335368 < x, x < 4.322453755}
         XX[1.5] := {4.128933055 < x, x < 4.474164703}
         XX[1.6] := {4.364794925 < x, x < 4.620897435}
         XX[1.7] := {4.597090458 < x, x < 4.763112053}
         XX[1.8] := {4.825976698 < x, x < 4.901201868}
 

question -->

anyone knows a simple way how I can plot the points of the lower boundaries and upper boundaries of x as two line graphs (without copy pasting all the values into separate lists and then using the listplot command)?

thanks

I have a function involving sinh(x)/cosh(x) to evaluate (see attached SKS77.mw and S77.pdf). I got different values depending on whether or not the function is algebraically simplified. What's the right way to evaluate the exponential function in this case? Thanks

Occasionally I use the Variables palette to inspect some variables after a run. In Maple2018 it seems every Vector or Array returns something like "Empty variable structure" or similar. I never saw this in prior versions of Maple.

A bug? Or am I missing something??

M.D.

I have four spindle tori governed by the following equations:

f1=x^4+2*x^2*y^2+2*x^2*z^2+y^4+2*y^2*z^2+z^4-2*x^3-2*x*y^2-2*x*z^2-(79/25)*x^2-(104/25)*y^2-(8/5)*z^2+(104/25)*x

f2=x^4+2*x^2*y^2+2*x^2*z^2+y^4+2*y^2*z^2+z^4-2*x^2*y-2*y^3-2*y*z^2-(104/25)*x^2-(79/25)*y^2-(8/5)*z^2+(104/25)*y

f3=x^4+2*x^2*y^2+2*x^2*z^2+y^4+2*y^2*z^2+z^4+2*x^3+2*x*y^2+2*x*z^2-(79/25)*x^2-(104/25)*y^2-(8/5)*z^2-(104/25)*x

f4=x^4+2*x^2*y^2+2*x^2*z^2+y^4+2*y^2*z^2+z^4+2*x^2*y+2*y^3+2*y*z^2-(104/25)*x^2-(79/25)*y^2-(8/5)*z^2-(104/25)*y

I could plot the surfaces using implicitplot3d and I can imagine the volume common to these surfaces but I could not visualize it. So, am looking for a way to plot the volume covered by the surfaces such that f1<0, f2<0, f3<0 and f4<0. I know that it's easy in case of a 2D filled plot but is there any way this could be done for the 3D case? Any mathematical advice as to how to characterize or calculate this volume would also be great.

 

 

Why Maple does not simplify to 1

beta^(1/2)*(1+beta^(1/2))/(beta+beta^(1/2))

I'm trying to use the continued fractions expansion to simplify a fraction, and I need to cut the expansion off for large values (for example when a denominator is greater than 20). I have a list of the outputs but I'm struggling to find the correct commands to use in order to do this.

p := ContinuedFraction((q-1)/2^t);
print(p);
cf := Term(p, 0..100);
print(cf);
R := Convergent( p, 100);
r := denom(R);
printf("The order is %d"\n, r);
 

This is the section of code, any help would be appreciated, thanks.

Dear Users!

Hoped everyone fine with everything! I want to define a square matrix P whose elements are

p[i,j]=<φ[i],ψψ[j]>

printlevel := 2; for i from 0 while i <= 2^k-1 do for j from 0 while j <= M-1 do varphi[i, j] := t^j end do end do;
printlevel := 2; for i from 0 while i <= 2^k-1 do for j from 0 while j <= M-1 do phi[M*i+j+1] := varphi[i, j] end do end do;
printlevel := 2; for i from 0 while i <= 2^k-1 do for j from 0 while j <= M-1 do psi[i, j] := 2^((1/2)*k)*sqrt(2*j+1)*(sum((-1)^(j+i1)*factorial(j+i1)*(2*t-i)^i1/(factorial(j-i1)*factorial(i1)^2), i1 = 0 .. j)) end do end do;
printlevel := 2; for i from 0 while i <= 2^k-1 do for j from 0 while j <= M-1 do `&psi;&psi;`[M*i+j+1] := psi[i, j] end do end do;
 

take k=2,M=3

There is a bug in inttrans:-hilbert:

restart;

inttrans:-hilbert(sin(a)*sin(t+b), t, s);
# should be:
sin(a)*cos(s+b);   expand(%);

sin(a)*cos(s)

 

sin(a)*cos(s+b)

 

sin(a)*cos(s)*cos(b)-sin(a)*sin(s)*sin(b)

(1)

########## correction ##############

`inttrans/expandc` := proc(expr, t)
local xpr, j, econst, op1, op2;
      xpr := expr;      
      for j in indets(xpr,specfunc(`+`,exp)) do
          econst := select(type,op(j),('freeof')(t));
          if 0 < nops(econst) and econst <> 0 then
              xpr := subs(j = ('exp')(econst)*combine(j/('exp')(econst),exp),xpr)
          end if
      end do;
      for j in indets(xpr,{('cos')(linear(t)), ('sin')(linear(t))}) do
          if type(op(j),`+`) then
              op1:=select(has, op(j),t); ##
              op2:=op(j)-op1;            ##
              #op1 := op(1,op(j));
              #op2 := op(2,op(j));
              if op(0,j) = sin then
                  xpr := subs(j = cos(op2)*sin(op1)+sin(op2)*cos(op1),xpr)
              else
                  xpr := subs(j = cos(op1)*cos(op2)-sin(op1)*sin(op2),xpr)
              end if
          end if
      end do;
      return xpr
end proc:

#######################################

inttrans:-hilbert(sin(a)*sin(t+b), t, s); expand(%);

-(1/2)*cos(a-b)*sin(s)+(1/2)*sin(a-b)*cos(s)+(1/2)*cos(a+b)*sin(s)+(1/2)*sin(a+b)*cos(s)

 

sin(a)*cos(s)*cos(b)-sin(a)*sin(s)*sin(b)

(2)

 


Download hilbert.mw

 

Hello

I was trying to introduce vector r_vec that has 3 components in x,y,z

I've attached my file, I've 2 questions here

first, why isn't the vector shown in as r_vec = () ei + () ej + () ek instead appears as a column vector

second, why doesn't it accept differentating

thank you
 

restart````

r := sqrt(VectorCalculus:-`+`(x^2, VectorCalculus:-`+`(y, VectorCalculus:-`-`(VectorCalculus:-`*`(z, 1/VectorCalculus:-`*`(tan, alpha))))^2))

(x^2+(y-z/(tan*alpha))^2)^(1/2)

(1)

 

theta := arctan((y-z)/x)

arctan((y-z)/x)

(2)

Z := VectorCalculus:-`*`(z, 1/VectorCalculus:-`*`(sin, alpha))

z/(sin*alpha)

(3)

with(VectorCalculus): 

r_vec := `<,>`(VectorCalculus:-`*`(r, `cos&theta;`), VectorCalculus:-`+`(VectorCalculus:-`*`(VectorCalculus:-`*`(r, sin), theta), VectorCalculus:-`*`(VectorCalculus:-`*`(Z, cos), alpha)), VectorCalculus:-`*`(VectorCalculus:-`*`(Z, sin), alpha))

Matrix(3, 1, {(1, 1) = sqrt(x^2+(y-z/(tan*alpha))^2)*`cos&theta;`, (2, 1) = sqrt(x^2+(y-z/(tan*alpha))^2)*sin*arctan((y-z)/x)+z*cos/sin, (3, 1) = z})

(4)

diff(r_vec, r)

Error, (in VectorCalculus:-diff) invalid input: diff received (x^2+(y-z/(tan*alpha))^2)^(1/2), which is not valid for its 2nd argument

 

``


 

Download tst1.mwtst1.mw

To demonstrate Maple 2018’s new Python connectivity, we wanted to integrate a large Python library. The result is the DeepLearning package - this offers an interface to a subset of the Tensorflow framework for machine learning.

I thought I’d share an application that demonstrates how the DeepLearning package can be used to recognize the numbers in images of handwritten digits.

The application employs a very small subset of the MNIST database of handwritten digits. Here’s a sample image for the digit 0.

This image can be represented as a matrix of pixel intensities.        

The application generates weights for each digit by training a two-layer neural network using multinomial logistic regression. When visualized, the weights for each digit might look like this.

Let’s say that we’re comparing an image of a handwritten digit to the weights for the digit 0. If a pixel with a high intensity lands in

  • an intensely red area, the evidence is high that the number in the image is 0
  • an intensely blue area, the evidence is low that the number in the image is 0

While this explanation is technically simplistic, the application offers more detail.

Get the application here

I have the following pertubation problem I want to use maple to expand for me.

 

We have epsilon := eps;

x(t,eps):= x_{-1}(t)/eps+x_0(t)+x_1(t)*eps

z(t,eps):=z_{-1}(t)/eps+z_0(t)+z_1(t)*eps

I want to expand a Taylor series of the following function upto some arbitray order of eps, i.e O(eps^3) or higher (depending on my mood :-)), around t=0, f(x(t,eps),z(t,eps),cos(t/eps),sin(t/eps)).

 

Anyone has any suggestion how to use maple 2017.3 to do this?

Thanks!

 

 

My code.mw

 

Footnote: Can we use variation of parameters method in this question?

First 824 825 826 827 828 829 830 Last Page 826 of 2216