Maple 2018 Questions and Posts

These are Posts and Questions associated with the product, Maple 2018

Does there exist a Maple command that on its own calculates the residual standard error of two regression lists? 

Given the following functions, graph them and identify relative and absolute extrema (if any).

f(x)=3x^3-2x^2+5x-7     [-3,6]

How I can do ?

Thank you.

 

Substitution of . 5,6,7) into Eqs. 1–(4), gives the new equation as functions of the generalized coordinates,
u_m,n(t);  v_m,n ( t), and w_m,n ( t). These expressions are then inserted in the Lagrange equations (see Eq. 8)) a set of N second-order coupled ordinary differential equations with both quadratic   and cubic nonlinearities.

In Eq (8) q are generalized coordinate such as uvw  and q = {`u__m,n`(t), `v__m,n`(t), `w__m,n`(t)}^T.

\where the elements of the vector,q_i are the time-dependent generalized coordinates.

L_Maple
 

U = (1/2)*(int(int(int(E*(`∂`(u(x, y, t))/`∂`(x)+(1/2)*(`∂`(w(x, y, t))/`∂`(x))^2+`∂`(w(x, y, t))/`∂`(x)*(`∂`(w__0(x, y, t))/`∂`(x))-z*(diff(w(x, y, t), x, x))+v(x, y, t)*(`∂`(v(x, y, t))/`∂`(y)+(1/2)*(`∂`(w(x, y, t))/`∂`(y))^2+`∂`(w(x, y, t))/`∂`(y)*(`∂`(w__0(x, y, t))/`∂`(y))-z*(diff(w(x, y, t), y, y))))*(`∂`(u(x, y, t))/`∂`(x)+(1/2)*(`∂`(w(x, y, t))/`∂`(x))^2+`∂`(w(x, y, t))/`∂`(x)*(`∂`(w__0(x, y, t))/`∂`(x))-z*(diff(w(x, y, t), x, x)))/(-nu^2+1)+E*(`∂`(nu(x, y, t))/`∂`(y)+(1/2)*(`∂`(w(x, y, t))/`∂`(y))^2+`∂`(w(x, y, t))/`∂`(y)*(`∂`(w__0(x, y, t))/`∂`(y))-z*(diff(w(x, y, t), y, y))+v(x, y, t)*(`∂`(u(x, y, t))/`∂`(x)+(1/2)*(`∂`(w(x, y, t))/`∂`(x))^2+`∂`(w(x, y, t))/`∂`(x)*(`∂`(w__0(x, y, t))/`∂`(x))-z*(diff(w(x, y, t), x, x))))*(`∂`(v(x, y, t))/`∂`(y)+(1/2)*(`∂`(w(x, y, t))/`∂`(y))^2+`∂`(w(x, y, t))/`∂`(y)*(`∂`(w__0(x, y, t))/`∂`(y))-z*(diff(w(x, y, t), y, y)))/(-nu^2+1)+E*(`∂`(u(x, y, t))/`∂`(y)+`∂`(v(x, y, t))/`∂`(x)+`∂`(w(x, y, t))/`∂`(x)*(`∂`(w(x, y, t))/`∂`(y))+`∂`(w__0(x, y, t))*`∂`(w(x, y, t))/(`∂`(x)*`∂`(y))+`∂`(w__0(x, y, t))*`∂`(w(x, y, t))/(`∂`(x)*`∂`(y))-2*z*(diff(w(x, y, t), x, y)))^2/(2*(1+nu))+E*l^2*(diff(w(x, y, t), x, y))^2/(1+nu)+E*l^2*(diff(w(x, y, t), x, y))^2/(1+nu)+E*l^2*(diff(w(x, y, t), y, y)-(diff(w(x, y, t), x, x)))^2/(2*(1+nu))+E*l^2*(diff(v(x, y, t), y, y)-(diff(u(x, y, t), x, x)))^2/(8*(1+nu))+E*l^2*(diff(v(x, y, t), x, y)-(diff(u(x, y, t), y, y)))^2/(8*(1+nu)), z = -(1/2)*h .. (1/2)*h), y = 0 .. b), x = 0 .. a))

U = (1/2)*(int(int((1/12)*(-E*(-v(x, y, t)*(diff(diff(w(x, y, t), y), y))-(diff(diff(w(x, y, t), x), x)))*(diff(diff(w(x, y, t), x), x))/(-nu^2+1)-E*(-v(x, y, t)*(diff(diff(w(x, y, t), x), x))-(diff(diff(w(x, y, t), y), y)))*(diff(diff(w(x, y, t), y), y))/(-nu^2+1)+4*E*(diff(diff(w(x, y, t), x), y))^2/(2+2*nu))*h^3+E*(`∂`(u(x, y, t))/`∂`(x)+(1/2)*`∂`(w(x, y, t))^2/`∂`(x)^2+`∂`(w(x, y, t))*`∂`(w__0(x, y, t))/`∂`(x)^2+v(x, y, t)*(`∂`(v(x, y, t))/`∂`(y)+(1/2)*`∂`(w(x, y, t))^2/`∂`(y)^2+`∂`(w(x, y, t))*`∂`(w__0(x, y, t))/`∂`(y)^2))*(`∂`(u(x, y, t))/`∂`(x)+(1/2)*`∂`(w(x, y, t))^2/`∂`(x)^2+`∂`(w(x, y, t))*`∂`(w__0(x, y, t))/`∂`(x)^2)*h/(-nu^2+1)+E*(`∂`(nu(x, y, t))/`∂`(y)+(1/2)*`∂`(w(x, y, t))^2/`∂`(y)^2+`∂`(w(x, y, t))*`∂`(w__0(x, y, t))/`∂`(y)^2+v(x, y, t)*(`∂`(u(x, y, t))/`∂`(x)+(1/2)*`∂`(w(x, y, t))^2/`∂`(x)^2+`∂`(w(x, y, t))*`∂`(w__0(x, y, t))/`∂`(x)^2))*(`∂`(v(x, y, t))/`∂`(y)+(1/2)*`∂`(w(x, y, t))^2/`∂`(y)^2+`∂`(w(x, y, t))*`∂`(w__0(x, y, t))/`∂`(y)^2)*h/(-nu^2+1)+E*(`∂`(u(x, y, t))/`∂`(y)+`∂`(v(x, y, t))/`∂`(x)+`∂`(w(x, y, t))^2/(`∂`(x)*`∂`(y))+2*`∂`(w__0(x, y, t))*`∂`(w(x, y, t))/(`∂`(x)*`∂`(y)))^2*h/(2+2*nu)+2*E*l^2*(diff(diff(w(x, y, t), x), y))^2*h/(1+nu)+E*l^2*(diff(diff(w(x, y, t), y), y)-(diff(diff(w(x, y, t), x), x)))^2*h/(2+2*nu)+E*l^2*(diff(diff(v(x, y, t), y), y)-(diff(diff(u(x, y, t), x), x)))^2*h/(8+8*nu)+E*l^2*(diff(diff(v(x, y, t), x), y)-(diff(diff(u(x, y, t), y), y)))^2*h/(8+8*nu), y = 0 .. b), x = 0 .. a))

(1)

T = rho*h*(int(int((`∂`(u(x, y, t))/`∂`(t))^2+(`∂`(v(x, y, t))/`∂`(t))^2+(`∂`(w(x, y, t))/`∂`(t))^2, y = 0 .. b), x = 0 .. a))

T = rho*h*(int(int(`∂`(u(x, y, t))^2/`∂`(t)^2+`∂`(v(x, y, t))^2/`∂`(t)^2+`∂`(w(x, y, t))^2/`∂`(t)^2, y = 0 .. b), x = 0 .. a))

(2)

F = (1/2)*c*(int(int((`∂`(u(x, y, t))/`∂`(t))^2+(`∂`(v(x, y, t))/`∂`(t))^2+(`∂`(w(x, y, t))/`∂`(t))^2, y = 0 .. b), x = 0 .. a))

F = (1/2)*c*(int(int(`∂`(u(x, y, t))^2/`∂`(t)^2+`∂`(v(x, y, t))^2/`∂`(t)^2+`∂`(w(x, y, t))^2/`∂`(t)^2, y = 0 .. b), x = 0 .. a))

(3)

W = int(int(w(x, y, t)*f__1(x, y, t)*cos(omega*t), y = 0 .. b), x = 0 .. a)

W = int(int(w(x, y, z)*f__1(x, y, z)*cos(omega*t), y = 0 .. b), x = 0 .. a)

(4)

u(x, y, t) = sum(sum(`u__m,n`(t)*sin(m*Pi*x/a)*sin(n*Pi*y/b), n = 1 .. N), m = 1 .. M)

u(x, y, t) = -(1/4)*(cos(Pi*y*N/b)*cos(Pi*x/a)*sin(Pi*y/b)*sin((M+1)*Pi*x/a)-cos(Pi*y*N/b)*cos((M+1)*Pi*x/a)*sin(Pi*x/a)*sin(Pi*y/b)+cos(Pi*x/a)*sin(Pi*y*N/b)*cos(Pi*y/b)*sin((M+1)*Pi*x/a)-cos((M+1)*Pi*x/a)*sin(Pi*x/a)*sin(Pi*y*N/b)*cos(Pi*y/b)-cos(Pi*y*N/b)*sin(Pi*y/b)*sin((M+1)*Pi*x/a)-cos(Pi*x/a)*sin(Pi*y*N/b)*sin((M+1)*Pi*x/a)-cos(Pi*x/a)*sin(Pi*y/b)*sin((M+1)*Pi*x/a)+cos((M+1)*Pi*x/a)*sin(Pi*x/a)*sin(Pi*y*N/b)+sin(Pi*x/a)*sin(Pi*y/b)*cos((M+1)*Pi*x/a)-sin(Pi*y*N/b)*cos(Pi*y/b)*sin((M+1)*Pi*x/a)+sin(Pi*y*N/b)*sin((M+1)*Pi*x/a)+sin(Pi*y/b)*sin((M+1)*Pi*x/a))*`u__m,n`(t)/((cos(Pi*x/a)-1)*(cos(Pi*y/b)-1))+(1/4)*(-cos(Pi*y*N/b)*sin(Pi*y/b)*sin(Pi*x/a)-sin(Pi*y*N/b)*cos(Pi*y/b)*sin(Pi*x/a)+sin(Pi*y*N/b)*sin(Pi*x/a)+sin(Pi*y/b)*sin(Pi*x/a))*`u__m,n`(t)/((cos(Pi*x/a)-1)*(cos(Pi*y/b)-1))

(5)

v(x, y, t) = sum(sum(`v__m,n`(t)*sin(m*Pi*x/a)*sin(n*Pi*y/b), n = 1 .. N), m = 1 .. M)

v(x, y, t) = -(1/4)*(cos(Pi*y*N/b)*cos(Pi*x/a)*sin(Pi*y/b)*sin((M+1)*Pi*x/a)-cos(Pi*y*N/b)*cos((M+1)*Pi*x/a)*sin(Pi*x/a)*sin(Pi*y/b)+cos(Pi*x/a)*sin(Pi*y*N/b)*cos(Pi*y/b)*sin((M+1)*Pi*x/a)-cos((M+1)*Pi*x/a)*sin(Pi*x/a)*sin(Pi*y*N/b)*cos(Pi*y/b)-cos(Pi*y*N/b)*sin(Pi*y/b)*sin((M+1)*Pi*x/a)-cos(Pi*x/a)*sin(Pi*y*N/b)*sin((M+1)*Pi*x/a)-cos(Pi*x/a)*sin(Pi*y/b)*sin((M+1)*Pi*x/a)+cos((M+1)*Pi*x/a)*sin(Pi*x/a)*sin(Pi*y*N/b)+sin(Pi*x/a)*sin(Pi*y/b)*cos((M+1)*Pi*x/a)-sin(Pi*y*N/b)*cos(Pi*y/b)*sin((M+1)*Pi*x/a)+sin(Pi*y*N/b)*sin((M+1)*Pi*x/a)+sin(Pi*y/b)*sin((M+1)*Pi*x/a))*`v__m,n`(t)/((cos(Pi*x/a)-1)*(cos(Pi*y/b)-1))+(1/4)*(-cos(Pi*y*N/b)*sin(Pi*y/b)*sin(Pi*x/a)-sin(Pi*y*N/b)*cos(Pi*y/b)*sin(Pi*x/a)+sin(Pi*y*N/b)*sin(Pi*x/a)+sin(Pi*y/b)*sin(Pi*x/a))*`v__m,n`(t)/((cos(Pi*x/a)-1)*(cos(Pi*y/b)-1))

(6)

w(x, y, t) = sum(sum(`w__m,n`(t)*sin(m*Pi*x/a)*sin(n*Pi*y/b), n = 1 .. N), m = 1 .. M)

w(x, y, t) = -(1/4)*(cos(Pi*y*N/b)*cos(Pi*x/a)*sin(Pi*y/b)*sin((M+1)*Pi*x/a)-cos(Pi*y*N/b)*cos((M+1)*Pi*x/a)*sin(Pi*x/a)*sin(Pi*y/b)+cos(Pi*x/a)*sin(Pi*y*N/b)*cos(Pi*y/b)*sin((M+1)*Pi*x/a)-cos((M+1)*Pi*x/a)*sin(Pi*x/a)*sin(Pi*y*N/b)*cos(Pi*y/b)-cos(Pi*y*N/b)*sin(Pi*y/b)*sin((M+1)*Pi*x/a)-cos(Pi*x/a)*sin(Pi*y*N/b)*sin((M+1)*Pi*x/a)-cos(Pi*x/a)*sin(Pi*y/b)*sin((M+1)*Pi*x/a)+cos((M+1)*Pi*x/a)*sin(Pi*x/a)*sin(Pi*y*N/b)+sin(Pi*x/a)*sin(Pi*y/b)*cos((M+1)*Pi*x/a)-sin(Pi*y*N/b)*cos(Pi*y/b)*sin((M+1)*Pi*x/a)+sin(Pi*y*N/b)*sin((M+1)*Pi*x/a)+sin(Pi*y/b)*sin((M+1)*Pi*x/a))*`w__m,n`(t)/((cos(Pi*x/a)-1)*(cos(Pi*y/b)-1))+(1/4)*(-cos(Pi*y*N/b)*sin(Pi*y/b)*sin(Pi*x/a)-sin(Pi*y*N/b)*cos(Pi*y/b)*sin(Pi*x/a)+sin(Pi*y*N/b)*sin(Pi*x/a)+sin(Pi*y/b)*sin(Pi*x/a))*`w__m,n`(t)/((cos(Pi*x/a)-1)*(cos(Pi*y/b)-1))

(7)

diff(`∂`(T(x, y, t))/`∂`(`#mscripts(mi("q"),mi("j"),none(),none(),mo("."),none(),none())`), t)-`∂`(T(x, y, t))/`∂`(`#mscripts(mi("q"),mi("j"),none(),none(),mo("."),none(),none())`)+`∂`(U(x, y, t))/`∂`(`#mscripts(mi("q"),mi("j"),none(),none(),mo("."),none(),none())`)+`∂`(U(x, y, t))/`∂`(`#mscripts(mi("q"),mi("j"),none(),none(),mo("."),none(),none())`)+`∂`(F(x, y, t))/`∂`(`#mscripts(mi("q"),mi("j"),none(),none(),mo("."),none(),none())`) = `∂`(W(x, y, t))/`∂`(`#mscripts(mi("q"),mi("j"),none(),none(),mo("."),none(),none())`), j = 1, () .. (), N

(D(`∂`))(T(x, y, t))*(diff(T(x, y, t), t))/`∂`(`#mscripts(mi("q"),mi("j"),none(),none(),mo("."),none(),none())`)-`∂`(T(x, y, t))/`∂`(`#mscripts(mi("q"),mi("j"),none(),none(),mo("."),none(),none())`)+2*`∂`(U(x, y, t))/`∂`(`#mscripts(mi("q"),mi("j"),none(),none(),mo("."),none(),none())`)+`∂`(F(x, y, t))/`∂`(`#mscripts(mi("q"),mi("j"),none(),none(),mo("."),none(),none())`) = `∂`(W(x, y, t))/`∂`(`#mscripts(mi("q"),mi("j"),none(),none(),mo("."),none(),none())`), j = 1, () .. (), N

(8)

NULL


 

Download L_Maple

 

 

Hi, 

I was able to determine a cubic spline fit, F(v), to x1 and y1. Now I have vector x2 which I would like to use F(v) to calculate y2 as another Vector[row]. I am having trouble accomplishing this task. Any help is greatly appreciated. Thanks.
 

restart

 x1 := Vector[row]([0.8e-1, .28, .48, .68, .88, 1, 1.2, 1.4, 1.6, 1.8, 2, 2.2, 2.4, 2.6, 2.8, 3, 3.2, 3.4, 3.6, 3.8, 4, 4.2]);

 y1 := Vector[row]([-10.081, -10.054, -10.018, -9.982, -9.939, -9.911, -9.861, -9.8, -9.734, -9.659, -9.601, -9.509, -9.4, -9.293, -9.183, -9.057, -8.931, -8.806, -8.676, -8.542, -8.405, -8.265]);

 

m := ArrayTools[Dimensions](x1);

maxx := rhs(m[1]);

 

F := proc (v) options operator, arrow; CurveFitting:-Spline(x1, y1, v, degree = 3) end proc;

 

x2 := Vector[row]([seq(log10(2*10^x1[k]), k = 1 .. maxx)])

 

y2:=?

 

Pts1 := plot(x1, y1, style = point, symbol = diamond, gridlines = true, color = red);

plt_sp := plot(F(v), v = x1[1] .. x1[maxx], color = blue);

plots:-display(Pts1, plt_sp)``

"# How to calculate Vector y2 using spline fit F with x2"? "    x1:=Vector[row]([0.08,0.28,0.48,0.68,0.88,1,1.2,1.4,1.6,1.8,2,2.2,2.4,2.6,2.8,3,3.2,3.4,3.6,3.8,4,4.2]):    y1:=Vector[row]([-10.081,-10.054,-10.018,-9.982,-9.939,-9.911,-9.861,-9.8,-9.734,-9.659,-9.601,-9.509,-9.4,-9.293,-9.183,-9.057,-8.931,-8.806,-8.676,-8.542,-8.405,-8.265]):    m:=ArrayTools[Dimensions](x1):  maxx:=rhs(m[1]):      F:=v->CurveFitting:-Spline(x1,y1, v,degree=3):    x2:=Vector[row]([seq(log10(2*10^(x1[k])),k=1..maxx)]):                   #` PLOT RESULTS`   Pts1:=plot(x1,y1,style=point,symbol = diamond, gridlines=true, color = red):       plt_sp:=plot(F(v),v=x1[1]..x1[maxx],color = blue):     plots:-display(Pts1,plt_sp);     "

 

``

``


 

Download splfit.mw

In Maple 2018, I was playing around with some sums of infinite series, and I came across a result that made me wonder if Maple was perhaps using some other definition or understanding of the sum of a series in its calculation. Take a look at the screenshot linked below:

https://ibb.co/hMdkQHn

That first series is most certainly divergent since the limit as n approaches infinity of n^2/(n+1) is not equal to 0. And just to confirm my own sanity, I even checked some of the partial sums of the series, which sure enough are diverging. And yet for the infinite sum, Maple is giving this finite result.

I even checked a more familiar alternating series, the alternating harmonic series, which Maple does correctly calculate to be ln(2).

What am I missing here? Is Maple using a different definition for the sum of the series than the limit of the partial sums as n approaches infinity? Or is there a mistake with how I've written something that I'm not noticing?


 

"I am trying to solve these two algebraic equations for A and B "

 

eq1 := -(64*A*B^2*k[1]+28*B^2*k[2]+112*A*k[1]+105*k[2])*A/(105*B)-(32*A^2*B^2*k[1]+28*A*B^2*k[2]+56*A^2*k[1]+105*A*k[2]+210*a)/(105*B) = 0

-(1/105)*(64*A*B^2*k[1]+28*B^2*k[2]+112*A*k[1]+105*k[2])*A/B-(1/105)*(32*A^2*B^2*k[1]+28*A*B^2*k[2]+56*A^2*k[1]+105*A*k[2]+210*a)/B = 0

(1)

eq2 := -(64*A^2*B*k[1]+56*A*B*k[2])*A/(105*B)+(32*A^2*B^2*k[1]+28*A*B^2*k[2]+56*A^2*k[1]+105*A*k[2]+210*a)*A/(105*B^2) = 0

-(1/105)*(64*A^2*B*k[1]+56*A*B*k[2])*A/B+(1/105)*(32*A^2*B^2*k[1]+28*A*B^2*k[2]+56*A^2*k[1]+105*A*k[2]+210*a)*A/B^2 = 0

(2)

solve({eq1, eq2}, {A, B})

{A = RootOf(192*k[1]^2*_Z^3+440*k[1]*k[2]*_Z^2+(480*a*k[1]+210*k[2]^2)*_Z+315*a*k[2]), B = (1/2)*RootOf(1536*k[1]^2*RootOf(192*k[1]^2*_Z^3+440*k[1]*k[2]*_Z^2+(480*a*k[1]+210*k[2]^2)*_Z+315*a*k[2])^2+2512*k[1]*k[2]*RootOf(192*k[1]^2*_Z^3+440*k[1]*k[2]*_Z^2+(480*a*k[1]+210*k[2]^2)*_Z+315*a*k[2])+21*_Z^2+3840*a*k[1]+105*k[2]^2)/k[2]}

(3)

``

How can I find A and B explicitly

 

 

``


 

Download solve_equattion.mw

Using the definition of a derivative as a limit i.e., lim h->0  f(x+h)-f(x)/h  .Find the derivative of the following functions:

a) f(x)=3x^3-2x^2+5x-7

 

how do i solve and type in maple? please help

Hi, I've created an array from a do loop, where the number in the array is the number of degrees for which I've calculated the answer, I now wish to graph the angle against the value in the array. i.e. the value a[50]=156.21 should graph to the point (50,156.21). Ideally as well I'd like for it to be joined up by a smooth curve of best fit, but I'll take what I can get, I can't seem to figure out how to plot it at all at the moment.

See code below:

restart;

for i from 50 by 5 to 85 do
ThetaBn := (1/180)*i*Pi;
s := cos(2*ThetaBn)*x+(2*sin(ThetaBn)*sin(ThetaBn))*sin(x);
a[i] := 180.0*fsolve(s = 0, x, 1 .. 6)*(1/Pi)
end do;

 

Thanks

I would like to understand why the 'solve' command is unable to find any solution to the system of equation { x^2=2, x^3=2*sqrt(2) }  (which obviously has a unique solution x=sqrt(2) ). The syntax that I used is
> solve({x^2 = 2, x^3 = sqrt(2)^3}, [x]);
and the output is the empty list.

I suspect that this is related to the presence of the algebraic number sqrt(2). Does anyone have a better understanding ?

I am using Maple version 2018.0, build ID 1298750.

Thank you.

I discovered incidentally that the command  Matrix(3, 3, `-`) (the number 3 is purely illustrative) returned the same result than the command Matrix(3, 3, (i,j) -> i-j).
In the same way `+` realizes (i,j) -> i+j), `*` realizes (i,j) -> i*j), ...

More surprisingly `.` realizes (i,j) -> i*j while I'm in worksheet mode, with "old" maple input style, and that the command 2.3 does not answer 6 but concatenates 2 and 3.

Is this a known behaviour or an undocumented feature?

 

 

 

 

 

Hello,

I am attempting to check solutions to a system of ODEs using odetest. However, odetest only appears to partially substitute the provided solution. Furthermore, it appears to be related to the alphabetical order of the functions.

For instance, here I have two functions, phiL and phiM, that satisfy Laplace's equation and are coupled through the boundary conditions, BCs.

laplace := {-phiL(z) + diff(phiL(z),z$2)=0, -phiM(z) + diff(phiM(z),z$2)=0}:
BCs := {phiL(d1)=0,phiM(-d1)=0,phiL(0)=phiM(0), D(phiL)(0)-D(phiM)(0)=-n}:
sol := {
phiM(z) = n/2/coth(d1)*(cosh(z)+coth(d1)*sinh(z)),
phiL(z) = n/2/coth(d1)*(cosh(z)-coth(d1)*sinh(z))
}:

odetest(sol, laplace union BCs,{phiL(z),phiM(z)});

This returns

{0, 1/2*(2*phiL(0)*coth(d1)-n)/coth(d1), D(phiM)(0)-1/2*n}

Here, phiL(0) and phiM(0) are unevaluated even though the provided solutions are valid there.

Furthermore, while renaming phiL to an alphabetically earlier name (eg, phiJ) causes the corresponding change in the output. However, renaming it to something alphabetically after phiM (eg, phiN) causes the terms in the output to switch. That is, changing phiL to phiN in the above code results in

{0, 1/2*(2*phiM(0)*coth(d1)-n)/coth(d1), D(phiN)(0)+1/2*n}

Therefore, it seems to be related to the way Maple internally stores the list of variables.

Is this a bug? Or is there something I'm missing?

Thanks!

My question has two steps:

STEP 1:  The multiplication  of is defined as follows

 

if n<>l, then

.

if n=l and m<=s,

Question 1: I wrote a code for calculating the multiplication  of. Is it right?

The code for Step 1  

restart;

multiply:=proc(n,m,l,s) local g,a: 
a:=unapply(doublefactorial(2*j-1)/factorial(j),j):
g:=unapply((a(m-j)*a(j)*a(s-j)/a(m+s-j))*(2*m+2*s-4*j+1)/(2*m+2*s-2*j+1),j):

if n<>l then 0 else
sqrt((2*m+1)*(2*s+1))*2^(K/2-1).add((g(j)/sqrt(m+s-2*j+1/2))*phi[n, m+s-2*j],j=0..m) 
end if
end proc:
 
n:=2:
l:=2:
m:=1:
s:=1:
multiply(n,m,l,s);

when I compared the results which I got and the results which is given in the book as follows, I think it is right.

Step 2:

We know that the outer product matrix is calculated as follows 

  

We found the elements of the outer product matrix in Step 1. 

Question 2 : I want to write the elements which are derived in step 1 to the outer product matrix in step 2. In here, the outer product matrix is NxN matrix. N=(M+1).2^(K-1) where K, M are any integers.

Hello

Could you please help me to solve the following problem

I want to catch the optimized values for df, tf, hw, and tw using considered constraint.

Thank you very much.

Dear Community,

 

I'm trying to import a simple Excel (xlsx) file into my Maple worksheet.

1st problem: Maple reads it only, if I put it into my root directory. As soon as I try to read it from the directory I'm working in, and provide full file name (Path + file name)  it fails. Why?

2nd problem: The real surprise comes, when I manage to import it. Its a 150 x 2 matrix, but as soon as I try to retrieve a value, I get an error message "Error, bad index into Matrix". Why is this? Should be very simple. Files attached.

 

Tx in advance

best regards

Andras

ImportTest.mw

RC_Rate_Schedule.xlsx


Suggestions for a Maple program entering the coordinates of three points and giving the equation of the circle passing through these three points, non-aligned?
Thank you.

 

 

First 6 7 8 9 10 11 12 Last Page 8 of 39