Items tagged with plotting

Feed

I don't recall ever changing any settings in my copy of Maple 16. But for some reason, I tried to make a plot and it spit out line-printed output rather than a regular plot.

How do I change this back? I have no use for line-printed output.

lineprinted_fail.mw
 

plot([E1(x), E2(x), E3(x)], x = 0 .. 5, title = "Nitrogen Nucleus Second-Order Energy Shifts")


                                                                              
                  Nitrogen Nucleus Second-Order Energy Shifts                 
  |                                                                           
  +                                                                     HHHHH
6e-08                                                               HHHHH     
  +                                                            HHHHHH         
4e-08                                                     HHHHHH              
  +                                                  HHHHHH                   
  +                                           HHHHHHH                         
2e-08                                  HHHHHHHH                               
  +                           HHHHHHHHHH                                      
  +             HHHHHHHHHHHHHHH                                               
  ***************************************************************************
0 +              1    HHHHHHHHHH2HH           3              4              5
  +                               HHHHHHHHH                                   
-2e-08                                    HHHHHHHH                            
  +                                              HHHHHHH                      
-4e-08                                                 HHHHHH                 
  +                                                         HHHHHH            
  +                                                              HHHHH        
-6e-08                                                                HHHHH   
  +                                                                       HHH
                                                                              

 

NULL

``


 

Download lineprinted_fail.mw

 

Hello,
I have a second problem which is related to converting horizontal axis to degrees ( or in terms of Pi :  0 - Pi/2 ).

Thank you for your helps

Bengu


restart;
s1 := (1/3*(4*(1.66+10^(-1))^2-x^2))^(3/2)/(x^2*(1+.66+10^(-1)))-sin(2*a);
s2 := (1/3*(4*(1.66+10^(-2))^2-x^2))^(3/2)/(x^2*(1+.66+10^(-2)))-sin(2*a);
s3 := (1/3*(4*(1.66+10^(-3))^2-x^2))^(3/2)/(x^2*(1+.66+10^(-3)))-sin(2*a);
k1 := series(s1, x);
                    (-2)                                 
       4.769067180 x     + (-0.5773502688 - 1. sin(2. a))

                           2    / 4\
          + 0.01164914507 x  + O\x /
k2 := series(s2, x);
                    (-2)                                 
       4.293792441 x     + (-0.5773502692 - 1. sin(2. a))

                           2    / 4\
          + 0.01293857500 x  + O\x /
k3 := series(s3, x);
                    (-2)                                 
       4.247636750 x     + (-0.5773502693 - 1. sin(2. a))

                           2    / 4\
          + 0.01307916821 x  + O\x /
p1 := convert(k1, polynom);
  4.769067180                                                2
  ----------- - 0.5773502688 - 1. sin(2. a) + 0.01164914507 x 
       2                                                      
      x                                                       
p2 := convert(k2, polynom);
  4.293792441                                                2
  ----------- - 0.5773502692 - 1. sin(2. a) + 0.01293857500 x 
       2                                                      
      x                                                       
p3 := convert(k3, polynom);
  4.247636750                                                2
  ----------- - 0.5773502693 - 1. sin(2. a) + 0.01307916821 x 
       2                                                      
      x                                                       
with(plots, implicitplot);
implicitplot([p1, p2, p3], a = 0 .. (1/2)*Pi, x = 0 .. 3, style = point, symbol = [box, cross, circle], labeldirections = [horizontal, vertical], labels = ["θ (radians)", "Switching  field (T)"], color = [blue, red, black], symbolsize = 13);

Dear sir, 

 

I request your kind help with regard to the above cited subject.

 

 

With warm regards.

 

Mr.M.Anand

Associate Professor in Mathematics.

Aurora's Technological & Research Institute, 

Hyderabad-INDIA

What coding will display a spherical cap of an arbitrary sphere (an arbitrarily located centre and arbitrary radius) where the cap covers an arbitrary solid angle and the radius to the cap's centre is any arbitrary radius of the sphere?

Assume the sphere is defined by its radius and its azimuth and polar angles.

hello experts,

I was using maple for a physical  problem,

and things turned very complicated with a equation with bessel function in it,

like this BesselI(1, (0.9067480359e-2+0.9067480359e-2*I)*sqrt(f))andBesselI(0., (0.9067480359e-2+0.9067480359e-2*I)*sqrt(f)),

which include complex,

the whole equation is as followed:

the variable is f and RV,,dependent variable is RV.

how am I supposed to plot RV when f=100..4000?

it is certain RV has real part and imaginary part,maybe i need a 3Dplot?

please let me know if you have any idea.

best regards,


 

-7.873519774*10^18*RV^4+(2.676513624*10^12-3.842712573*10^15*(-1)^(3/4)*BesselI(1, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))/(Pi*sqrt(f)*BesselI(0, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))*(4+623.8617593*(-1)^(3/4)*BesselI(1, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))/(BesselI(0, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))*sqrt(f)))))*RV^2+80864.83845+1.440831316*10^9*(-1)^(3/4)*BesselI(1, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))/(Pi*sqrt(f)*BesselI(0, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))*(4+623.8617593*(-1)^(3/4)*BesselI(1, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))/(BesselI(0, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))*sqrt(f)))) = 0
 

want to plot the relationship between RV and f,but how?

 

complexplot(RV, f = 100 .. 4000, labels = ["Re", "Im"])

complexplot(RV, f = 100 .. 4000, labels = ["Re", "Im"])

(1)

NULL


 

Download bessel_in_equation.mw

 

work.mw

I'm having problems plotting these three functions on the same graph. I can get two of them, but as soon as I try to add the third, it throws a couple errors no matter what I do. It doesn't help that the docs are confusing. The first two functions are approximating diff(y(x),x) = -2xy(x)+1 with the Euler method (step size of 0.05 and 0.1) that I implemented using a for loop. I gather the data into an array and then plot those. That works. But once I try to set y(x) equal to a Taylor approximation, I can plot that graph but the values for the other functions are strangely absent. I've attached my work. What is wrong with my syntax and how can I plot all three functions at once?

I want to plot a function f(x,y)=x^3-2x^2y. But I only want to plot values where $x^2+y^2 <= 1$. So far I specified the plotting range in the standard way, but I really need to specify it like I described. I want to use plot3d.

I am trying to plot some graphs for a differential equations class. I need to plot two equations (soln1a, soln1b1) in one graph, with the two curves overlayed on each other. I am trying to plot them on a graph named "gr1c". The error I get says "Error, (in plots/multiple) empty plot". Clicking on the error results in a webpage saying that "There is no help page available for this error".

When I plot each equation separately, there are no errors. The graphing command that works for soln1a is 

gr1a:=plot(rhs(soln1a),dom1,color=blue);

 

and the graphing command that works for soln1b1 is

gr1b1:=plot(rhs(soln1b_1),dom1,color=purple);
 

My code is as follows. I would greatly appreciate any help ASAP to remove the error. Thank you!

 

with(plots):with(DEtools):

K:=9;
deG:=diff(theta(t),t,t) + mu*diff(theta(t),t)+K*sin(theta(t))= 0;
deL:=diff(theta(t),t,t) + mu*diff(theta(t),t)+K*theta(t)= 0;
Iv:=theta(0)=0.75, D(theta)(0)=2.0;
dom1:=t=0..10;
soln1a:=dsolve({eval(deL,mu=0),Iv});

soln1b_1:=dsolve({eval(deL,mu=1),Iv});

gr1c:=multiple(plot[soln1a,soln1b1],dom1,color=[blue,purple]);
 

dear all,

I have a problem when I try to find points on boundary of the surface.

the surface ploted from a matrix as follow:

 

please help me to generate a bondary curve of this surface.

Thanks

 

 

 

I have made a program to plot y(t) vs x(t). Can you please modify to plot y vs x data from an excel file.

Given x in A column and y in B column in 5 rows.

Example data: [[1,1],[2,4],[3,9],[4,16],[5,26]]

StLine_VRK.msim

Thanks.

Ramakrishnan V

I want to set the position of the legend inside the figure window for example say top right in the figure window. Maple chose the legend to be at right, left, top, bottom outside the figure window. How can I place the legend of my own wish?

thanks,

Hi, everybody.

I have a problem when I try to plot3d a sub-surface from a surface as follows:


 

S := proc (u, v) options operator, arrow; Matrix([[40*u], [80*v], [10*u^2*v+20*u*v+15]]) end proc;

proc (u, v) options operator, arrow; Matrix([[40*u], [80*v], [10*u^2*v+20*u*v+15]]) end proc

(1)

S(u, v)

Matrix(3, 1, {(1, 1) = 40*u, (2, 1) = 80*v, (3, 1) = 10*u^2*v+20*u*v+15})

(2)

 

p:= proc(u,v) if u<v then S(u,v) else S(u,v)+10 end if end proc:

h:= proc(u) 2*u  end proc:

plot3d(p, 0 .. 1, 0 .. h)

Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

 

 

``


 

Download sub-surface.mw

sub-surface.mw

 

Please help me!

 

Thanks and have a nice day.

 

 

I have final project to make a media for learning mathematic using maple. But I'm so confused to make a net of cuboid, anybody can help me? please :D

Dear All,

Please, could anyone point me out how to draw the beautiful Apollonian Gasket using Maple?

Thank you.

Very kind wishes,

Wang Zhe

We have a new video about plotting with procedures, specifically on how to avoid an error related to evaluation that many people make. 

The worksheet I used in the video is available here:  PlottingWithProcedures.mw

By the way, this mistake doesn't occur just in the context of plotting. Users of Optimization commands (and other commands that allow functions to be expressed as either an expression in a variable or as a procedure) frequently run into this problem too.

 

1 2 Page 1 of 2