MaplePrimes Questions

matrixplot(A, style=surfacewireframe,heights=histogram);

when A has zero values they are graphed, is there any way to not have them show? Choose any matrix for A I just want to see non-zero elements.

 

Also when looking down on the 3d graph I'd like the tops of the histogram cells to be colored as if I'm looking at a 2d plot. When I do it now I can't tell any difference between any of the cells. This looks to be a lighting issue I suppose. There is some variation but it is very difficult to tell. Some light settings work better than others. Z-Hue seems to be the best but still a little washed out.

I've always found that the way maple handles rotating graphs to be very unintuitive. It is very difficult for me to get the exact view I want by rotating the graph.

 

Usually i want to view for 1 2d projection then rotate to view the other. In most programs one can do this easier and there usually is some modifier key. In maple this is not the case and how it handles rotation is bizzaire.

For example, if I set the orientation to be x axis and start to rotate by dragging the mouse horizontally the 3d graph starts to rotate all 3 axis. It seems to depend on where the mouse is exactly at in the graph. That is, if I drag only in the horizontal but the mouse starts at the top then it is different than when it starts at the bottom. This is not intuitive it it requires me to know exactly how maple does things and thing work backwards...

Many programs have a modifier key that will snap one direction/angle to not change. This is idea. rather than have all 3 angles changing arbitrarily(it is nice when one just wants to perview the graph in 3D but not good when one wants to get a specific angle).

 

Is there any easy way to get what I want? I know I can set the orientation manually. That is not what I'm after. I'm after changing the rotation algorithm that maple uses or to add a modifier key to restrict rotation to only 2 angles.

 

 

Draw the surface area of the paraboloid z=12-x^2-y^2 which is inside the cone z=sqrt(x^2+y^2)?

Hello,

I was working one of the problems for my course in structural dynamics and came across the following function that needs to be plotted. How can we do it in maple 2018.

 

i have a acceleration of record earthquake(in excel ) and i get fourier and power spectra from another software (seismosignal) .Now i want to get it in maple but i try a lot but i couldn't.
i very Eager to learn it

rr.xlsx

I compared the numerical value of an integral expression in Gradshteyn and Rhyzhik (8.431.5) for the BesselI(v,z) function to the BesselI(v,z) function in Maple 2019.2 and got different answers.  The GR expression is the difference of two integrals and it looks like Maple only evaluates the first integral, which is correct only if the order v is an integer.  Please see the attached.

BesselI.mw

At first, I appology my poor English and I am a non-proessional on mathematics and programing. I am an old retired engineer/researcher.

I found strange phenomena with maple.
Could someone please verify this result?
This is sound like a stupid joke.
This queation is related to factorization of composite number N=p*q
p, q are large prime number respectively.

Degits:=200;
1.    Let p is large prime such as 50 degits
2.    Let R is small integer such as 3
3.    Let q is nearest prime R*p
4.    Let N=p*q
5.    Let f=frac(N/(q+s))

j:=1;
m:=2^b*2^c;
r:=m/2^b;
b:=8; c:=0;
F:=[ ];
for s from -m to m by r do
f:= frac(N/(q+s));
      if f=0 OR 1.0 then f:=100000; # zero dev exeption;
      elif f <= 0.5 then f:=log10(f);
      elif f > 0.5 AND f < 1.0 then f:=-log10(1.0-f);
      else print(“error”) break;
      end if:
F:=[F[ ], [s, f]];  # maple plot format
j::=j+1;
end do:

6.    Draw a Graph s vs F

You wil find a oscillating graph with step at s=0;

7. b:=8 and c:=75; and try same caluclation.

You will find a graph with dip at s=0;
Caution! Resolution and Scanwidth can not go together.

8.    Let qx=q+10^25;
9.    Nx:=p*qx;
10.    Let ps:=sqrt(N/3.0) and let qs:=3.0*sqrt(N/3.0);
11.    fs:=frac(Nx/(qx+s));
Caution!  Nx not equal Ns:=p*qs
12.     Draw a graph s vs Fs with Nx, p, qs
You will find a similar graph dip point shifted from s=0
13. Verify qx=ps-dip point s

My questions are
1.    Is it a kaind of factorization of Nx?
2.    This phenomena hava posibirity to make a vulneravility for RSA crypto?
3.    Power spectrum of F( c=0) indicates R=q/p; with an peak
4.    Is it an rediscovery of “Wheel”?  Are thea any papers similar to this phenomena

I appreciate to read this questions.

Environment: Windows 10 for workstation and Maple 2019(64bit) CPU: Intel Xeon.

15 Dec 2019  some miss spelling corrected

                       A graph added   frac(N(q+s)) vs s b=8 c=75 

 

Dear maple user ,while solving the system of pdes i am getting the errors and unable to plot the curves. please help me to  rectify the errors and  plot the curves.

Thanks in advance

 

restart;
with(PDETools):
G:=1300:beta:=0.0075:epsilon:=4.9*(10)^(-6):Pa:=100:alpha:=24:Pv:=97:
sys1:= {G*(beta^2*u(r,z)*diff(u(r,z),z)+v(r,z)*diff(u(r,z),r)) = -diff(p(r,z),z) +beta^2*diff(u(r,z), z$2)+diff(u(r,z),r$2)+(1/r)*diff(u(r,z),r),G*(beta^2*u(r,z)*diff(v(r,z),z)+v(r,z)*diff(v(r,z),r)) = -diff(p(r,z),r) +beta^2*diff(v(r,z), z$2)+diff(v(r,z),r$2)+(1/r)*diff(v(r,z),r)-(v(r,z)/r^2),beta^2*diff(u(r,z),z)+diff(v(r,z),r)+(v(r,z)/r)=0  }:

IBC:=(D[1](u))(0, z) = 0,v(0,z)=0,phi*(D[1](u))(1, z)+u(1,z) = 0,v(1,z)=epsilon*(p(1,z)+(Pa/alpha)-1),p(r,-1)=0,p(1,z)=(Pv-Pa)/alpha:
sol := pdsolve(sys1, IBC, numeric):
r:=0:

p1 := sol:-plot(u(r, z), phi = 0, numpoints = 100, z = -1 .. 1, color = ["Blue"], legend = ["u(r,z)"]):

p2 := sol:-plot(u(r, z), phi = 0.15, numpoints = 100, z = -1. 1, color = ["red"], legend = ["u(r,z)"]):

p3 := sol:-plot(u(r, z), phi = 0.4, numpoints = 100, z = -1 .. 1, color = ["green"], legend = ["u(r,z)"]):

plots:-display({p1, p2,p3});

 

 

"a=4 x+1/diff(f(y),y)";
b := 3*(4*x + 1)/(diff(f(y), y)*(3*x + 1));
subs((4*x + 1)/diff(f(y), y) = a, b);
 

I am trying to figure out which computational method is used in function Issimilar, (which determines similarity of matrices) in Maple 2018 and whether or not its a rational one.

how can i use curve fitting tools for get m and n in this equation  for excel data

equation

i want to get m and n?

New_Microsoft_Excel_Worksheet.xlsx

 

How would one code a real time clock (digital's the easiest - of course, or analog) within maple?

Hi,

Why extracting the features of the events doesn't return the good results when executed within a loop (see the pink test) ?

Thanks in advance.

PS: please, spare me  replies of the type "you can solve this equation formally"

 

restart;

interface(version);

`Standard Worksheet Interface, Maple 2015.2, Mac OS X, December 21 2015 Build ID 1097895`

(1)

sys := { diff(x(t), t) = 1, x(0) = 0 }:
evs := [ [x(t)-0.1, none],  [x(t)-0.3, none], [x(t)-0.5, none] ]:
sol := dsolve(sys, numeric, events=evs):

plots:-odeplot(sol, [t, x(t)], t=0..0.5, gridlines=true);

 

# times that fired the events

sol(1): # initialization

sol(eventfired=[1]);
sol(eventfired=[2]);
sol(eventfired=[3]);
 

[HFloat(0.1)]

 

[HFloat(0.3)]

 

[HFloat(0.49999999999999994)]

(2)

# Same times computed  within a loop

for i from 1 to 3 do
  te := op(sol(eventfired=[i]));
end do;

HFloat(0.1)

 

HFloat(0.3)

 

HFloat(0.49999999999999994)

(3)

# Values of x(t) computed  within a loop
#
# Why are calues for events 2 and 3 wrong ?
for i from 1 to 3 do
  te := op(sol(eventfired=[i]));
# xe := sol(te);             # this doesn't return the correct result,
# xe := subs(sol(te), x(t)); # this doesn't work neither
  xe := eval(x(t), sol(te)); # this doesn't work neither
end do;

HFloat(0.1)

 

HFloat(0.1)

 

HFloat(0.0)

 

HFloat(0.0)

 

HFloat(0.0)

 

HFloat(0.0)

(4)

 


 

Download Incomprehensible.mw

 

Dear experts

I am using Maple to solve a complex equation. My idea is to separate real and imaginary parts and then solve a set of the equation when both real and imaginary parts are zero. the following are the equation and the way I made real and imaginary parts;

(K*( Q*sinh(K)*cosh(Q)-K*cosh(K)*sinh(Q))*(1+s*K^2)    +p*(-4*K^2*Q*(K^2+Q^2)        +Q*(Q^4+2*K^2*Q^2+5*K^4)*cosh(K)*cosh(Q)        -K*(Q^4+6*K^2*Q^2+K^4)*sinh(K)*sinh(Q)))/(K^2*Q*cosh(Q))

eq:= (K*( Q*sinh(K)*cosh(Q)-K*cosh(K)*sinh(Q))*(1+s*K^2)+p*(-4*K^2*Q*(K^2+Q^2)+Q*(Q^4+2*K^2*Q^2+5*K^4)*cosh(K)*cosh(Q)-K*(Q^4+6*K^2*Q^2+K^4)*sinh(K)*sinh(Q)))/(K^2*Q*cosh(Q)):

so the K and the Q are both complex variables and p and s are constant.

p := 0.1019367992e-3, s := 7.135575943      K:=Kr+I*Kim    Q:= sqrt(K^2-I*h^2*2*Pi/1.0e-6)

K:=Kr+I*Kim;

Q:= sqrt(K^2-I*h^2*2*Pi/1e-6);

therefore the real and imaginary parts of the equation are 

A:=evalc(Re(eq)):   B:=evalc(Im(eq)):

finally, I tried to solve it as following

sys:={eval(A,[p=nu^2/g/h^3,s=sigma/rho/g/h^2])=0,eval(B,[p=nu^2/g/h^3,s=sigma/rho/g/h^2])=0}:

sol2:=(fsolve(sys,{Kr=0..1,Kim=4..5}));

sys:={eval(A,[p=nu^2/g/h^3,s=sigma/rho/g/h^2])=0,eval(B,[p=nu^2/g/h^3,s=sigma/rho/g/h^2])=0}:

sol2:=fsolve(sys,{Kr=0..5,Kim=0..5},maxsols=5);

 

the problem is that Maple can not solve it and returns the command. I  know that there is solutions. How can I solve this equation?

 

the maple file is attached.mapleprime.mw

==============================================================

I guess Maple use Newton method to solve equation or system of equations. Is there an alternative? I mean what are the possible methods?

Hello users. I have a question on my work.

I'm trying to construct the equation and plot it. And I got 2 errors(warning).

Please help me how to solve this problem check the image below and attached file.

 

 

Question_plot_the_curvature.mw

First 608 609 610 611 612 613 614 Last Page 610 of 2427