MaplePrimes Questions

I have code 
with(plots);

implicitplot3d((r*sin(theta)*cos(phi))^2+(r*sin(theta)*sin(phi))^2+(r*cos(theta))^2=36,r=-3..3, phi=0..2*Pi,theta=0..Pi,coords=spherical);

How can I fix the code to display the picture

How to plot f(x)=x in the interval -2<x<2 with a period 4?

The plot should look like this

Hello,
How do i generate 4d plot of f(x,y,z)=1 is it possible made it in 3d?

So I am trying to find the structure of a material from x-ray diffraction data. To accomplish this I have to solve a system of nonlinear equations. Since there is some unaccuracy in the measurements I wanted to use a numerical tool, so I am trying to solve using fsolve. To check that I am doing this correctly I am trying to solve a set of similar equations where I know the solution. My code is shown below.

eqs := {(2^2/a^2+2^2/c^2-2*(2*2)*cos(x)/(a*c))/sin(x)^2 = 1/2.2393^2, (2^2/a^2+4^2/c^2+4*(2*2)*cos(x)/(a*c))/sin(x)^2 = 1/1.5968^2, (1^2/a^2+sin(x)^2/b^2+2^2/c^2+(2*2)*cos(x)/(a*c))/sin(x)^2 = 1/2.7896^2, 1/sin(x)^2*(2^2*sin(x)^2/b^2) = 1/2.8650^2}

fsolve(eqs)

For this system I know that a = 5.44 b = 5.73 c = 7.89 x=pi/2, but fsolve only returns the input.

Am I doing something wrong or should I use a different function? I am quite new to Maple, so have some patience with me.

Thanks in advance.

Hello! Prompt please as in Maple can determine the distribution function of the resulting histogram distribution? I know about cdf function, but how to act in relation to the histogram do not know.

Histogram:

restart;
with(stats);
with(stats[statplots]);
data2 := [30, 30.5, 31, 31.5, 32, 32.5, 32.6, 33, 33.1, 33.3, 33.6, 34, 35, 36];
histogram(data2, area = count);

In other words, I need smoothing the histogram, get graph and get on it to obtain an analytical expression.

Aslam o Alikum!

Hope everyone going fine with sound health I want to make the graph of following system of ODEs

Eq1 := diff(F(eta), eta, eta, eta, eta)-M*(eta*(diff(F(eta), eta, eta, eta))+3*(diff(F(eta), eta, eta))+(diff(F(eta), eta))*(diff(F(eta), eta, eta))-F(eta)*(diff(F(eta), eta, eta, eta)))-Ha^2*(diff(F(eta), eta, eta));
Eq2 := diff(G(eta), eta, eta)+Pr*M*(F(eta)*(diff(G(eta), eta))-eta*(diff(G(eta), eta)))+Pr*Ec*(diff(F(eta), eta, eta))^2+Nb*(diff(H(eta), eta))+diff(G(eta), eta)+Nt*(diff(G(eta), eta))^2;
Eq3 := diff(H(eta), eta, eta)+M*Sc*(F(eta)*(diff(H(eta), eta))-eta*(diff(H(eta), eta)))+Nt*(diff(F(eta), eta, eta))/Nb;
 

using the BC's

IC1 := F(0) = 0, ((D@@2)(F))(0) = 0, (D(G))(0) = 0, (D(H))(0) = 0;
     F(0) = 0, @@(D, 2)(F)(0) = 0, D(G)(0) = 0, D(H)(0) = 0
IC2 := F(1) = 1, (D(F))(1) = 0, G(1) = 1, H(1) = 1;
 

for Ec := .1; Nt := .1; Nb := .1; Sc := .5; Pr := 10; M := .5; and different values of Ha like Ha:=0, Ha:=2, Ha:=4, Ha:=6 and Ha:=8 in one coordinate. 

I am waiting your positive response.


Hello all, first time posting on here so if I make any mistakes please let me know. So I have created a function called "g" and plotted it on a graph. After looking at the two graphs I have plotted (one with discont=false and one with discont=true), I can see 3 discontinuities are present. However, attempting to solve for these discontinuities by setting the denominator of "g" = 0 only returns 1 of the points. I'm not sure why this is happening and how would I go by solving for the other 2 points? Thank you.

f := exp(x-188.5)-4*x^2+1507.964*x-142122.3+cos(x)^2;

exp(x-188.5)-4*x^2+1507.964*x-142122.3+cos(x)^2

(1)

g := (-f+x^2-376.991*x+35530.6)/(f+x^2-376.991*x+35530.6);

(-exp(x-188.5)+5*x^2-1884.955*x+177652.9-cos(x)^2)/(exp(x-188.5)-3*x^2+1130.973*x-106591.7+cos(x)^2)

(2)

``

plot(g, x = 186 .. 194, y = -25 .. 25,discont=false)

 

plot(g, x = 186 .. 194, y = -25 .. 25,discont=true)

 

vasymp := fsolve( denom(g)=0);

192.2026966

(3)

``

Dear all,

My script has multiple animations, in which the range of the variable is larger then I wish to see. Therefore, I have used the view=[x1..x2,y1..y2] command, which has different values for every plot. Plotted seperately, this gives me exactly what I want. However, if I combine the plots using display(A,B,C) they loose their individual views. Instead, maple chooses the view of plot A, and plots B an C in the range defined by view A. 

Is there a way to combine plots and still maintain individual views, within the largest view?

Kind regards,
Bastiaan Overdorp

Am trying to valid a research work done by kuiken(1968)

Kuiken_(1968).pdf

where we have this two eauations:

restart;
Digits := 35;
with(ODETools);
with(student);
with(plots);
inf := 4;
equ1 := diff(f[0](eta), `$`(eta, 3))+theta[0](eta);
equ2 := diff(theta[0](eta), `$`(eta, 2))+3*f[0](eta)*(diff(theta[0](eta), eta));
Bcs1 := f[0](0) = 0, (D(f[0]))(0) = 0, theta[0](0) = 1, theta[0](inf) = 0, (D(D(f[0])))(inf) = 0;
S1 := dsolve({Bcs1, equ1, equ2}, {f[0](eta), theta[0](eta)}, type = numeric, method = bvp[midrich]);
proc(x_bvp)  ...  end;
S1(0);
[                            d                   
[eta = 0., f[0](eta) = 0., ----- f[0](eta) = 0., 
[                           deta                 

    d   /  d            \                                          
  ----- |----- f[0](eta)| = 0.82449782146165697398999365896678734, 
   deta \ deta          /                                          

  theta[0](eta) = 1.0000000000000000000000000000000000, 

    d                                                         ]
  ----- theta[0](eta) = -0.71098574970825563256340736114251047]
   deta                                                       ]
S1(inf);
[                                                            
[eta = 4., f[0](eta) = 1.7815670728545914261072119522795076, 
[                                                            

    d                                                      
  ----- f[0](eta) = 0.51061876174095320088291844433043562, 
   deta                                                    

    d   /  d            \                           
  ----- |----- f[0](eta)| = 0., theta[0](eta) = 0., 
   deta \ deta          /                           

    d                                                             
  ----- theta[0](eta) = -0.000054818176138173095945902421930470836
   deta                                                           

  ]
  ]
  ]
 

 

Pls, I need to find the function of the limit of f[0](eta) at eta tend to infinity. checked equation 45 of the attached document and for the two equation pls checked equation 36 and 37 for the ODE equation solved above.

Kuiken_solution for equation 36 and 37.pdf

Is it possible to animate plot like this parametically:
animate(plot, [[cos(t), sin(t), t = 0 .. A]], A = 0 .. 2*Pi, scaling = constrained, frames = 50)

however given the x and y components as solutions of an implicit equation.

I know I could run RootOf. But it seems that there is a part missing due to a branch cut :-/

What do people find a good screen size is? 

Personally I find two monitors would be ideal for using Maple and especially Maplesim.  I'm going to be adding a second monitor, probably a 20".  Maybe 24" is better?  I think 20" would be enough.  Just wondering how many monitors people out there work with and their monitor sizes.

After changing the time and date on my computer, I get a message saying “Invalid or missing license file.” whenever I try to open Maple. I have tried activating it with my license code. It then gives me a message saying “Activation succesful! Please restart Maple for the new license file to take effect.”. After opening it again, I will once again get the original message. I have also tried reinstalling it twice. One time with the file my school provides, and with a file from Maples website. I get the same error message no matter what I do. 

I haven't been able to access Maple for several weeks now, so I'm getting a bit behind in math class.

How should I go about fixing this?

Hi all,
I have a question about editing strings in Maple2016.
When I make a variable "DataName" which is as below in Maple worksheet,
DataName:=Array(["Data1","Data2","Data3"]):

I want to make the variable "DataSet" whose elements are as below.
DataSet:=[1="Data1",2="Data2",3="Data3"]

Of course I can type the elements of "DataSet", but if the number of DataName increases, it requires a lot of time to do that, so I want to make the "DataSet" by using general command so that I can obtain "DataSet" automatically.
I tried to use the command "cat", "||", "map", and so on... but I could not find how to do that.
Is there anyone who can solve this problem? Please let me know how to solve this problem. I appreciate any idea you may have.

Always thanks for the nice answer

hi..

is correct this answer for differential equations??

i think order of result should be in (10^6 or 10^9 or higher) range....

please check it

thanks

hpp.mw
 

restart

L := 100*10^(-9):

Eq1 := {-(1017/1600000000000000000000000000000000000000000)*(diff(w(x), x, x, x, x, x, x))+(26169/40000000000000000000000000)*(diff(w(x), x, x, x, x))-0.8325000000e-4*omega^2+1.560937500*10^(-21)*omega^2*(diff(w(x), x, x)), w(0) = 0, w(1/10000000) = 0, (D(w))(0) = 0, (D(w))(1/10000000) = 0, ((D@@2)(w))(0) = 0, ((D@@2)(w))(1/10000000) = 0}:

sys := subs(omega^2 = omega2, Eq1);

{-(1017/1600000000000000000000000000000000000000000)*(diff(diff(diff(diff(diff(diff(w(x), x), x), x), x), x), x))+(26169/40000000000000000000000000)*(diff(diff(diff(diff(w(x), x), x), x), x))-0.8325000000e-4*omega2+0.1560937500e-20*omega2*(diff(diff(w(x), x), x)), w(0) = 0, w(1/10000000) = 0, (D(w))(0) = 0, (D(w))(1/10000000) = 0, ((D@@2)(w))(0) = 0, ((D@@2)(w))(1/10000000) = 0}

 

{-(5085/8)*(diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y))+6542250*(diff(diff(diff(diff(g1(y), y), y), y), y))-0.8325000000e-4*omega2+0.1560937500e-6*omega2*(diff(diff(g1(y), y), y)), 10000000*(D(g1))(0) = 0, 10000000*(D(g1))(1) = 0, 100000000000000*((D@@2)(g1))(0) = 0, 100000000000000*((D@@2)(g1))(1) = 0, g1(0) = 0, g1(1) = 0}

 

{-(5085/8)*(diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y))+6542250*(diff(diff(diff(diff(g1(y), y), y), y), y))-0.8325000000e-4*omega2+0.1560937500e-6*omega2*(diff(diff(g1(y), y), y))}, {10000000*(D(g1))(0) = 0, 10000000*(D(g1))(1) = 0, 100000000000000*((D@@2)(g1))(0) = 0, 100000000000000*((D@@2)(g1))(1) = 0, g1(0) = 0, g1(1) = 0}

 

{g1(0) = 0, g1(1) = 0, (D(g1))(0) = 0, (D(g1))(1) = 0, ((D@@2)(g1))(0) = 0, ((D@@2)(g1))(1) = 0}

 

{diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y), diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), diff(diff(diff(diff(g1(y), y), y), y), y), diff(diff(diff(g1(y), y), y), y), diff(diff(g1(y), y), y), diff(g1(y), y)}

 

{-(1017/1600000000000000000000000000000000000000000)*(diff(diff(diff(diff(diff(diff(w(x), x), x), x), x), x), x))+(26169/40000000000000000000000000)*(diff(diff(diff(diff(w(x), x), x), x), x))-0.8325000000e-4*omega2+0.1560937500e-20*omega2*(diff(diff(w(x), x), x))}

 

{-(5085/8)*(diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y))+6542250*(diff(diff(diff(diff(g1(y), y), y), y), y))-0.8325000000e-4*omega2+0.1560937500e-6*omega2*(diff(diff(g1(y), y), y))}

 

{diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y) = 10292.62537*(diff(diff(diff(diff(g1(y), y), y), y), y))-0.1309734513e-6*omega2+0.2455752212e-9*omega2*(diff(diff(g1(y), y), y))}

 

{diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y) = 10292.62537*(diff(diff(diff(diff(g1(y), y), y), y), y))-0.1309734513e13*omega3+2455752212.*omega3*(diff(diff(g1(y), y), y))}

 

{diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y) = 10292.62537*(diff(diff(diff(diff(g1(y), y), y), y), y))-0.1309734513e13*omega3+2455752212.*omega3*(diff(diff(g1(y), y), y)), g1(0) = 0, g1(1) = 0, (D(g1))(0) = 0, (D(g1))(1) = 0, ((D@@2)(g1))(0) = 0, ((D@@2)(g1))(1) = 0}

 

{((D@@3)(g1))(0), ((D@@3)(g1))(1), ((D@@4)(g1))(0), ((D@@4)(g1))(1), ((D@@5)(g1))(0), ((D@@5)(g1))(1)}

 

((D@@3)(g1))(0)

 

((D@@3)(g1))(1)

 

((D@@4)(g1))(0)

 

((D@@4)(g1))(1)

 

((D@@5)(g1))(0)

 

((D@@5)(g1))(1)

 

((D@@5)(g1))(1), ((D@@4)(g1))(0), ((D@@5)(g1))(0), ((D@@3)(g1))(1), ((D@@4)(g1))(1), ((D@@3)(g1))(0)

 

HFloat(-8.852947665097804e-24), HFloat(-8.991820290300328e-22), HFloat(8.852947665097804e-24), HFloat(-9.672787782157173e-20), HFloat(-8.991820290300328e-22), HFloat(9.672787782157165e-20)

(1)

sqrt(8.85294766509780*10^(-21)*10^19);

.2975390338

(2)

NULL


 

Download hpp.mw

 

Dear all

I have n, p two integers greater than one

I would like to minimize the following real number by a fixed positive constant if its possible

(n/(n+p))^(n/p)-(n/(n+p))^((n+p)/p)

Many thanks

First 1054 1055 1056 1057 1058 1059 1060 Last Page 1056 of 2428