Advanced Search

Parham2016

70 Reputation

6 Badges

9 years, 212 days
Contact Parham2016

MaplePrimes Activity


These are questions asked by Parham2016

Contrudiction between plot and NextZero...

Parham2016 70 Maple
plot heun
November 10 2015
0 3

Wonderful.mw

Hi all

 

I don't know what is its meaning???

 

The plot and the solved Zero-roots aren't coincident with each other...

 

Zero-roots: 7.532332868 Wheras from the plot we noticed that there's a root near 14...

 

And why this Code can't find other Zero-roots?

 

Thanks a lot

How to get data from plot ?...

Parham2016 70 Maple
November 09 2015
1 1

Plot.mw

Hi all


As you know  in Matlab many Maple Functions didn't define. One of that is   HeunT  function. Now I want to plot that in Matlab, I'm surfing in the net, Unfortunately I can't find it...!

Please help me to get all the data points from my Plot in Maple so that I could plot those in Matlab.

Analysis of beam with two supports...

Parham2016 70 Maple
beam mechanics
November 02 2015
2 0

Hi all

I want to analysis this beam  (a simple and a tip rotational spring) in Maple to find the first four Frequencies...

Please help me

Plot Heun in Maple...

Parham2016 70 Maple
heun
October 24 2015
1 1

Heun.mw

plot3d([r, theta, -3.3203*HeunT(.4995036958*3^(2/3)/([1.1233, 6.8291, 12.2089, 17.4216, 32.9018, 53.2699, 64.3282]*[1])^(4/3), 0, (.3138423830*[1.1233, 6.8291, 12.2089, 17.4216, 32.9018, 53.2699, 64.3282])*[1]*3^(1/3), -.3258398511*3^(2/3)*r)*exp(0.5743565187e-1*r*(2.710463448*r^2+3))+4.9407*exp(-(.1148713037*(1.355231724*r^2+3/2))*r)*HeunT(.4277706929*3^(2/3), 0, .3525391488*3^(1/3), .3258398511*3^(2/3)*r)], r = 0 .. 1, theta = 0 .. 2*Pi, coords = cylindrical)

 

Download Heun.mw

Hi all

I want to know if there's any plot for the function:

 

Zero Roots of equation...

Parham2016 70 Maple
rootfinding
October 21 2015
2 10

mysolve.mw

``

k := -81.58765950*HeunT(2.405631288*10^5*W^(2/3)*3^(2/3), 0, 980.9170330*W^(1/3)*3^(1/3), -0.5544992207e-1*3^(2/3)*W^(1/6))*exp(81.58842672*sqrt(W))*sqrt(W)-81.58765950*exp(-81.58842670*sqrt(W))*HeunT(2.405631288*10^5*W^(2/3)*3^(2/3), 0, 980.9170330*W^(1/3)*3^(1/3), 0.5544992207e-1*3^(2/3)*W^(1/6))*sqrt(W)

-81.58765950*HeunT(240563.1288*W^(2/3)*3^(2/3), 0, 980.9170330*W^(1/3)*3^(1/3), -0.5544992207e-1*3^(2/3)*W^(1/6))*exp(81.58842672*W^(1/2))*W^(1/2)-81.58765950*exp(-81.58842670*W^(1/2))*HeunT(240563.1288*W^(2/3)*3^(2/3), 0, 980.9170330*W^(1/3)*3^(1/3), 0.5544992207e-1*3^(2/3)*W^(1/6))*W^(1/2)

 (1)

plot(k, W = 0 .. 8, color = ["OrangeRed"], labels = ["W", "Determinant"]);

  

``

``

Digits := 10:

K := proc (W) options operator, arrow; (-1)*81.58765950*HeunT(2.405631288*10^5*W^(2/3)*3^(2/3), 0, 980.9170330*(W^(1/3))*3^(1/3), (-1)*0.5544992207e-1*(3^(2/3))*W^(1/6))*exp(81.58842672*sqrt(W))*sqrt(W)+(-1)*81.58765950*exp((-1)*81.58842670*sqrt(W))*HeunT(2.405631288*10^5*W^(2/3)*3^(2/3), 0, 980.9170330*(W^(1/3))*3^(1/3), 0.5544992207e-1*(3^(2/3))*W^(1/6))*sqrt(W) end proc; G := table([0 = 0]); for k to 4 while G[k-1] < 10000 do G[k] := RootFinding:-NextZero(K, G[k-1], maxdistance = 10000, guardDigits = 3) end do

Error, cannot determine if this expression is true or false: FAIL < 10000

  

``

NULL

 

Download mysolve.mw

 

 

 

 

4 5 6 7 8 9 10 Page 6 of 10