Maple 17 Questions and Posts

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

Can Maple 17.0 work on Windows 10 Home Edition? or it can be only Windows 10 Pro Edition

Trying to plot the function 7*x^2+22*xy+7*y^2+14*xz*sqrt(3)+14*yz*sqrt(3)-5*z^2 = 180, I tried using implicitplot3d to plot it, with ranges I'm quite certain should contain the surface but for whatever reason all I get is a blank plot with no graphics. Here's my input:

 

Loading plots

with(plots);
[animate, animate3d, animatecurve, arrow, changecoords,

complexplot, complexplot3d, conformal, conformal3d,

contourplot, contourplot3d, coordplot, coordplot3d,

densityplot, display, dualaxisplot, fieldplot, fieldplot3d,

gradplot, gradplot3d, implicitplot, implicitplot3d, inequal,

interactive, interactiveparams, intersectplot, listcontplot,

listcontplot3d, listdensityplot, listplot, listplot3d,

loglogplot, logplot, matrixplot, multiple, odeplot, pareto,

plotcompare, pointplot, pointplot3d, polarplot, polygonplot,

polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus,

semilogplot, setcolors, setoptions, setoptions3d, spacecurve,

sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]
implicitplot3d(7*x^2+22*xy+7*y^2+14*xz*sqrt(3)+14*yz*sqrt(3)-5*z^2 = 180, x = -50 .. 50, y = -50 .. 50, z = -50 .. 50, axes = normal);

 

Not sure what I'm doing wrong exactly. 

As far as I can tell from the help pages, Maple 17, which I am using, can perform only one-dimensional Fourier transformations. Has that changed in latter versions?

I ask because I would like to find the 4D Fourier transformation of Heaviside(t^2 - x^2 - y^2 - z^2), the argument being the Minkowski line element. I have made quite some attempts with pen and paper, but the results are not stable. One calculational strategy of mine seems to depend on the order of integration: alongside some other part, a 4D Dirac delta function either appears or not. Another strategy produces an altogther differently looking expression.

It would therefore be nice to have some computational guidance.

The following is the log file after running a code file,

Extrema_Network := table([])

Extrema_Network[Han] := [20.0, 385.61]

Extrema_Network[Liv] := [20.0, 385.61]

Extrema_Network[Vir] := [20.0, 385.61]

values := [[2.8274333874308139146163, 2 Pi]]

theta_step := 0

phi_step := 0

counter_theta := 0

counter_phi := 0


distance eff distance_eff
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
Im in has par
Im in has par
the time offset for Han was .71848e-2*cos(lat_begin)*cos(long_begin)+.12746e-1*cos(lat_begin)*sin(long_begin)-.15394e-1*sin(lat_begin)
the time offset for Liv was .24690e-3*cos(lat_begin)*cos(long_begin)+.18302e-1*cos(lat_begin)*sin(long_begin)-.10809e-1*sin(lat_begin)
the time offset for Vir was -.15117e-1*cos(lat_begin)*cos(long_begin)-.28029e-2*cos(lat_begin)*sin(long_begin)-.14657e-1*sin(lat_begin)
the angles with the coordinate axis are [long_begin-3.1416, 1.5708+lat_begin]
Using logChirp

common factor in Signal() was .43329e-21

Using logChirp

common factor in Signal() was .43329e-21

Using logChirp

common factor in Signal() was .43329e-21

snr_network was (187.00-.48032e-3*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^2*sin(long_begin)-378.50*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)*sin(long_begin)-247.94*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^3*cos(long_begin)^3-56.997*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^3*cos(long_begin)-31.130*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)*cos(long_begin)^3+32.512*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^3*sin(long_begin)+196.54*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)*cos(long_begin)+.65598e-3*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^2*cos(long_begin)-.13119e-2*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^2*cos(long_begin)^3+64.585*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^3*cos(long_begin)^2*sin(long_begin)-61.911*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)*cos(long_begin)^2*sin(long_begin)+110.68*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^3*sin(long_begin)^2*cos(long_begin)-375.04*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)*sin(long_begin)^2*cos(long_begin)+.85041e-3*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^2*cos(long_begin)^2*sin(long_begin)+.10495e-3*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^2*sin(long_begin)^2*cos(long_begin)-.49247e-3*sin(lat_begin)^3+260.07*cos(long_begin)^4+.16559e-3*cos(long_begin)^4*sin(lat_begin)+.16559e-3*sin(lat_begin)^3*cos(long_begin)^4+.90222e-3*sin(lat_begin)^3*cos(long_begin)^2-.12334e-2*sin(lat_begin)*cos(long_begin)^2-268.44*cos(long_begin)^2+102.34*sin(long_begin)*cos(long_begin)-590.30*sin(lat_begin)^2*cos(long_begin)^2-312.01*sin(lat_begin)^4*cos(long_begin)^2-102.19*sin(lat_begin)^4*sin(long_begin)^2+260.07*sin(lat_begin)^4*cos(long_begin)^4+650.59*sin(lat_begin)^2*cos(long_begin)^4+114.40*cos(long_begin)^3*sin(long_begin)+33.292*sin(long_begin)^2*cos(long_begin)^2-.11326e-2*cos(long_begin)^3*sin(lat_begin)*sin(long_begin)-.85371e-4*sin(long_begin)^2*cos(long_begin)^2*sin(lat_begin)-.11326e-2*sin(lat_begin)^3*cos(long_begin)^3*sin(long_begin)+.49010e-3*sin(lat_begin)^3*sin(long_begin)*cos(long_begin)-.85371e-4*sin(lat_begin)^3*sin(long_begin)^2*cos(long_begin)^2-394.84*sin(lat_begin)^2*sin(long_begin)*cos(long_begin)+287.80*sin(lat_begin)^4*sin(long_begin)*cos(long_begin)+114.40*sin(lat_begin)^4*cos(long_begin)^3*sin(long_begin)-219.38*sin(lat_begin)^2*cos(long_begin)^3*sin(long_begin)+33.292*sin(lat_begin)^4*sin(long_begin)^2*cos(long_begin)^2+1085.5*sin(lat_begin)^2*sin(long_begin)^2*cos(long_begin)^2+.64246e-3*sin(lat_begin)*sin(long_begin)*cos(long_begin)+157.13*sin(long_begin)^2-54.941*sin(lat_begin)^2*sin(long_begin)^2+78.538*sin(lat_begin)^4+27.150*sin(lat_begin)^2+.57530e-3*sin(lat_begin))^(1/2)
low_fisher for detector Han was


Error, (in fprintf) number expected for floating point format
13

32

31

17

hou := 0

mini := 0

seci := 0

memory used=26.3MB, alloc=32.3MB, time=0.47




Not quite sure why I am getting this error. Also this code was written in maple 14 and I am using maple 17. Could that also explain why I am getting all these errors? Is it possible to downgrade maple 17 to maple 14? Sorry if my question above on fprintf is a bit too cryptic as it is research work. I will upload the code file in a reply.

After running Maple in a shell file, I come up with this error that I do not understand on my Mac,

gap_long := 0.117647058823529 Pi

gap_lat := 0.0588235294117647 Pi

lat_begin := 0.441176470588235 Pi

long_begin := -Pi

lat_begin_0 := 0.441176470588235 Pi

long_begin_0 := -Pi

long_max := 0.882352941176471 Pi

lat_max := -0.441176470588235 Pi

33

Warning, `parameter` is implicitly declared local to procedure `set_par_eff`

distance eff distance_eff
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
Im in has par
Im in has par
Error, invalid input: eval expects its 2nd argument, eqns, to be of type
{integer, equation, set(equation)}, but received par_eff_post
13

32

31

17

hou := 0

mini := 0

seci := 0

memory used=4.0MB, alloc=32.3MB, time=0.23



If needed, I can attach more files if my question is still a bit too cryptic. Please let me know asap as this is urgent. Thank you so much,
-Z

Dear All,

I am going to solve the following systems of ODEs but get the error: Newton iteration is not converging.
Could you please share your idea with me. In the case of AA=-0.2,0,0.2,0.4,...; I could get the solution.
Thank you in advance.


restart;
with(plots);
Pr := 2; Le := 2; nn := 2; Nb := .1; Nt := .1; QQ := .1; SS := .1; BB := .1; CC := .1; Ec := .1; MM := .2;AA:=-0.4;

Eq1 := diff(f(eta), `$`(eta, 3))+f(eta).(diff(f(eta), `$`(eta, 2)))-2.*nn/(nn+1).((diff(f(eta), eta))^2)-MM.(diff(f(eta), eta)) = 0; Eq2 := 1/Pr.(diff(theta(eta), `$`(eta, 2)))+f(eta).(diff(theta(eta), eta))-4.*nn/(nn+1).(diff(f(eta), eta)).theta(eta)+Nb.(diff(theta(eta), eta)).(diff(h(eta), eta))+Nt.((diff(theta(eta), eta))^2)+Ec.((diff(f(eta), `$`(eta, 2)))^2)-QQ.theta(eta) = 0;
Eq3 := diff(h(eta), `$`(eta, 2))+Le.f(eta).(diff(h(eta), eta))+Nt/Nb.(diff(theta(eta), `$`(eta, 2))) = 0;

bcs := f(0) = SS, (D(f))(0) = 1+AA.((D@@2)(f))(0), theta(0) = 1+BB.(D(theta))(0), phi(0) = 1+CC.(D(phi))(0), (D(f))(etainf) = 0, theta(etainf) = 0, phi(etainf) = 0

Error, (in dsolve/numeric/ComputeSolution) Newton iteration is not converging

I have a system of equations e.g.

A^2+B*A+C=0

where A,B,C are Matrices and I want to solve for A.

Sure I can write every equations in brakets [..=0], but isn'T it possible to just use the matrix notation?

Hello,

 

I tried to plot the problem presented below:

restart; with(plots); C := setcolors(); with(LinearAlgebra);

formula1 := 2.6*BodyWeight*abs(sin(4*Pi*t));
2.6 BodyWeight |sin(4 Pi t)|
BodyWeight := 80*9.81;
plot(formula1, t = 0 .. 2);


eq2 := formula1-SpringConstant*y(t)-DampConstant*(diff(y(t), t)) = Mass*(diff(y(t), `$`(t, 2)));
2040.480 |sin(4 Pi t)| - SpringConstant y(t)

/ d \ / d / d \\
- DampConstant |--- y(t)| = Mass |--- |--- y(t)||
\ dt / \ dt \ dt //
DampConstant := 50;
50
Mass := .200;
Springt := 200;
200
SpringConstant := Youngsmodulus*Surface/DeltaLength;
DeltaLength := 0.2e-1-y(t);
Surface := .15;
Youngsmodulus := 6.5*10^6/(t+1)+6.5*10^6;
plot(Youngsmodulus, t = 0 .. 10000);

eq2;
2040.480 |sin(4 Pi t)|

/ 6 \
|6.5000000 10 6|
0.15 |------------- + 6.5000000 10 | y(t)
\ t + 1 / / d \
- ----------------------------------------- - 50 |--- y(t)| =
0.02 - y(t) \ dt /

/ d / d \\
0.200 |--- |--- y(t)||
\ dt \ dt //

incs := y(0) = 0, (D(y))(0) = 0;
eq4 := dsolve({eq2, incs}, y(t), type = numeric, method = lsode[backfull], maxfun = 0);
proc(x_lsode) ... end;

plots:-odeplot(eq4, [t, y(t)], 0 .. 5);

 When I try to plot it beyond t=5, Maple gives the following error:

Warning, could not obtain numerical solution at all points, plot may be incomplete

Does anyone know how to plot it even further?

 

 

vz := 2*(-eta^2+1);

D_im := .22;

r0 := 1;

pde := diff(vz*Y(eta, z), z)-D_im*((diff(eta*(diff(Y(eta, z), eta)), eta))/eta+diff(Y(eta, z), `$`(z, 2)))/r0 = 0;

pde := expand(%);

ibc := [Y(1, z) = 0, (D[1](Y))(0, z) = 0, Y(eta, 0) = 1, (D[2](Y))(eta, 0) = 0];

sol := pdsolve(pde, ibc, numeric, time = z, range = 0 .. 1);

pds := sol:-value(z = 0, output = listprocedure);

sol:-plot(z = 0.1e-3, numpoints = 50, color = blue, view = 0 .. 1)

So I was trying to solve this conservation equation for the radial coordinate eta and the z coordinate being treated as time. The flow is in z direction. Now unfortunately it is diverging. Not sure why though. What am I doing wrong?

Im trying to draw a shpere but it always saying: 

Error, (in plot3d) unexpected option: z = -2 .. 2


this is the equation: x^2+y^2+z^2-4=0

i'm writing this way

plot3d(x^2+y^2+z^2-2^2, x = -2 .. 2, y = -2 .. 2, z = -2 .. 2)


what should I do? this is my first time with this software

 

best from Brazil,
Nina

Hi all,

Let $A$ be a 0-1 square matrix of order $n$. I want to obtain a matrix $D$ from powers $A$, $A^2$, $\dots$, $A^{n-1}$, where the ($i,j$)-element of $D$ is the smallest $k$ for which the ($i,j$)-element of $A^k$ is nonzero. I only consider the non-diagonal elements of $D$.For example, if the matrix $A$ is Matrix( [[0, 1, 0, 0, 0], [1, 0, 1, 0, 0], [0, 1, 0, 1, 0], [0, 0, 1, 0, 1], [0, 0, 0, 1, 0]]), then $D$ shoule be Matrix([[0, 1, 2, 3, 4], [1, 0, 1, 2, 3], [2, 1, 0, 1, 2], [3, 2, 1, 0, 1], [4, 3, 2, 1, 0]]). However, I cannot obtain this result.

The code is as follows.

p := proc ()

   local A, B, D, m, n, k, r;

   A := Matrix([[0, 1, 0, 0, 0], [1, 0, 1, 0, 0], [0, 1, 0, 1, 0], [0, 0, 1, 0, 1], [0, 0, 0, 1, 0]]);

   r := LinearAlgebra:-RowDimension(A);

   D := A;

   for k from 2 by 1 to r-1 do

       B := A^k;

       for m from 1 by 1 to r do

           for n from 1 by 1 to r do

               if m <> n and B[m, n] <> 0 and D[m, n] = 0 then

                    D[m, n] := k

                end if;

           end do;

        end do;

   end do;

   D;

end proc;

By executing this procedure, I obtain D=Matrix([[0,1,2,3,3],[1,0,1,2,3],[2,1,0,1,2],[3,2,1,0,1],[3,3,2,1,0]]), which is not I want.

Thanks.

When the Physics package is loaded, Maple returns the following error in connection with the simple statement 2*D[1]:

Error, (in TypeTools/ac_var_local) unable to determine the anticommutative character of D[1]

Why that? Although to no avail, for D[1] there is no problem. Note that the above error results even without having set up any anticommutative prefixes or the likes in Physics Setup.

PS: I am using Maple 17.

This is the coding till i do dhe decryption process. 

Do(plaintext=GetProperty("message",value));
Do(plaintext=convert(GetProperty("message",value),name));
Do(plaintextInt = convert(plaintext, bytes));
Do(plaintextBin = `~`[convert](plaintextInt, binary));
Do(plaintextBin2 = map2(nprintf, "%07d", plaintextBin));
Do(n0 = plaintextBin2[]);
Do(length1 = length(n0));
Do(plaintextCode = cat("", plaintextBin2[]));
Do(length2 = length(plaintextCode));
Do(z = convert(plaintextCode, decimal, binary));
Do(z1 = z+1);
Do(z2 = z1+%sk1);
Do(z3 = convert(z2, base, 2));
Do(b = cat("", z3[]));
Do(z4 = length(b));
Do(z5 = [Bits:-GetBits(-z2, -1 .. 0, bits = z4)]);
Do(z6 = cat("", z5[]));
Do(z7 = convert(z6, decimal, binary));
Do(%C = `mod`(Power(z7, %txte), %txtN));
Do(%C1 = `mod`(Power(%sk1, %txte), %txtN));

Do(%m = `mod`(Power(%C, %d), %N));

Do(%sk2=`mod`(Power(%C1,%d),%N));

Then nw i need to decrypt back to the original message with the coding:

Do(z8 = [Bits:-GetBits(-%m,-1 .. 0, bits = z4)]);
Do(c = cat("", z8[]));
Do(z9 = convert(c, decimal, binary));
Do(z10 = z9-sk2);
Do(z11 = z10-1);
Do(z12 = [Bits:-GetBits(z11, -1 .. 0, bits = length2)]);
Do(d = cat("", z12[]));
Do(plaintextBin2 = [StringTools:-LengthSplit(d, length1)]);
Do(plaintextInt2 = `~`[convert](plaintextBin2, decimal, binary));
Do(%message1 = convert(plaintextInt, bytes));

when i execute the program it shows the error

so how should I solve this as although i think that it should be problem of parsing the number z4 in the sentence that i highlighed, but whenever i correct it it still can't work.Thus anyone who know please help.Thanks.

 

restart:with(plots):
eq:=(diff(f(eta),eta$2))-a*f(eta)+b*(1+diff(f(eta),eta)^2)^(-1/2)=0;
bc:=f(1)=0,D(f)(0)=0;
ans := dsolve(eq);

I am trying to find the root of an equation. The problem is, I keep getting the error

"Error, (in fsolve) Can't handle expressions with typed procedures"

and 

Warning, solutions may have been lost


whenever I try to solve it. Anyone have any ideas? My worksheet is here:  1.mw

First 19 20 21 22 23 24 25 Last Page 21 of 60