Items tagged with plot plot Tagged Items Feed

I need a help about running the following codes..

 

restart:
d:=1:
s:=0.8:
c1:=sqrt((s)/Pi/r^3)/2 * (-sin(r*(u-1)^2/4/(s)+Pi/4)+sin(r*(u+1)^2/4/(s))/sqrt(2)):


c22:=-((u+1)/(4*Pi*r*sqrt(y-1)))*(exp(-r*y*(u+1)^2/(8*s)))*(sin(r*y*(u+1)^2/8/(s*(y-1))+Pi/4)/2+sin(r*y*(u-1)^2/8/(s*(y-1))+Pi/4)):
c2:=Int(c22, y= 1..100):


c33:=(sin(r*y*(u+1)^2/8/(s*(y-1))+Pi/4)+sin(r*y*(u-1)^2/8/(s*(y-1))+Pi/4))*(((d*s)/(Pi*(r^2)*y*sqrt(y-1))*exp(-r*y*(u+1)^2/8/s))-((d^2*sqrt((2*s^3)/(Pi*(y-1)*y^3*r^5)))*exp(d*(u+1)+(2*d^2*s)/(y*r))*erfc((r*y*(u+1)+4*d*s)/(2*sqrt(2*s*r*y))))):
c3:=Int(c33, y= 1..100):


g1:=c1+c2+c3:
g2:= unapply(Int(g1, r= 1..infinity, epsilon=1e-4,digits=7),u):
plot(g2, 0..10, numpoints=100);


 

Thank you!

VIM.3.mwHi,How do I make a graph and error table for exact and by using this method i.e VIM 

I was solving thos differential equiation using dsolve:

diff(P(t),t)=b*P(t)+m*t

dsolve({diff(P(t),t)=b*P(t)+m*t,P(t0)=P0});

where t0=2000 and P0=100 and m=9/50

 

Then, I found b setting at t=2001, P=2*P0=200, I applied fsolve to do this and the result was:

b=-1.5979...

 

To this point, Maple computes correctly. The problem starts when I tried to plot the result in a single 2D-graph (P vs t) using the computed value for b, the expression involved is:

P(t)=.1126460951*t-0.7049523744e-1-1.090463096*10^1390*exp(-1.597924898*t)


I wrote

plot( .1126460951*t-0.7049523744e-1-1.090463096*10^1390*exp(-1.597924898*t) , t=2000..2001);

The result was a blank plot, only the axes appeared even if I used the y=A..B option in the plot command. Is there a way to get the desired plot without having to evaluate the function at several points and then using this set of points to generate the plot: plot({[a,b],[c,d],[f,g]...})? (that worked, of course).

Thanks and regards.

P.S. I have used Maple for over 15 years now and I've plotted expressions even more complicated than this one, I found this very strange. I guess the 10^1390 has something to do this time, but I'm not sure...

 

plot(t+1, t); plot(t+1, t = -1 .. 1);

Error, (in plot) incorrect first argument t+1
Error, (in plot) incorrect first argument t+1

my project supervisor has made it mandatory that i build a maple maplet of the bessel equation:

Hoe do you plot two vectors of data as a bar plot or historygram. I tried the statistics package but could not plot a bar plot that shows the proper relation of numbers in vector x to numbers in relation to vector y. Vector x contained years as data, and vector y contained for example crime data per year.

The equation tan(y) = 2*tan(x) defines y implicitly as a function of x.  Well, perphas "defines" is too strong a word, since there are multiple solutions for y.  However, if I am not mistaken, there exists a unique continuous solution y(x) that goes through the origin, that is, y(0)=0, and is defined for all x.

Question 1: How do we plot the graph of y(x)?

I have a roundabout solution as follows.  Differentiate the equation tan(y(x)) = 2*tan(x) with respect to x and arrive at a first order differential equation in y(x).  Solve the differential equation with the initial condition y(0)=0.  Surprisingly, Maple obtains an explicit solution:

which we can plot:

plot(rhs(%), x=0..2*Pi);

Question 2: Is there a neat way of getting that solution with algebra only, without appealing to differential equations?

 

Hello there..I tried to plot this program n somehow it didnt work.. I dont know how to fix it (I'm a beginner user)..Can anyone help me?? Please..

Thanks a lot

--------------------------------------------------

hmin := .2;

tclose := 150;

topen := 120;

B0 := 0; B1 := 800; dB := (B1-B0)*(1/100);

B := dB*i+B0;

xcoor := [seq(B, i = 0 .. 99)];

ycoor := [seq(fsolve(tclose*ln((1-(1-hmin)*exp((-B+APD)/topen))/hmin) = APD), i = 0 .. 99)];

with(plots);

P := plot(xcoor, ycoor);

display(P)

Hello everyone, 

In Maple8, I tried to plot this logistic map and an error occured (Error, (in Bifurcation) `plots` does not evaluate to a module).

What is wrong into this code?

Thank you

 

restart: with(plots):Warning, the name changecoords has been redefined

> Bifurcation := proc(initialpoint,xexpr,ra,rb,acc)
> local p1,hr,A,L1,i,j,phi:
> global r,L2:
> hr := unapply(xexpr,x);
> A := Vector(600):
> L1 := Vector(acc*500):
> for j from 1 to acc+1 do
> r := (ra + (j-1)*(rb-ra)/acc):
> A[1] := hr(initialpoint):
> for i from 2 to 500 do
> A[i] := evalf(hr(A[i-1])):
> end do:
> for i from 1 to 400 do
> L1[i+400*(j-1)] := [r,A[i+100]]:
> end do:
> end do:
> L2 := {seq(L1[i], i = 1..acc*400)}:
> p1 := plots:-pointplot(L2, 'symbol' = solidcircle, 'symbolsize' = 8, 'color' = blue):
> unassign('r'):
> return(p1):
> end proc:
> P1:= Bifurcation(1/2,r*x*(1-x),2.5,4,250):
>
Error, (in Bifurcation) `plots` does not evaluate to a module

 

Hello everyone!

 

Do you have any idea to solve and plot a 2-time logistic map:

x(t+ 1_\alpha)= r*x(t)*(1-x(t))  ,t=(t^1,t^2)  ?

 

Thank you

 

 

 

M := 10^2; plot(exp(M), t);

When we execute the above code we get graph in output naturally. But when execute the following code  , there is no graph in output, why?

M := 10^3; plot(exp(M), t);

Thanks in advance for any suggestion.

Dear Maple users

I am delighted that Maple has builtin commands to plot so many polyhedrons in 3D. Here I am talking about the polyhedraplot command in the plots package. I was however disappointed that the socalled Truncated Icosahedron is not supported (not present in the supported list ...). My first question is:

1. Why isn't it supported?

It seems more relevant than many of the other polyhedrons which are supported. It is a member of the Archimedian Solids (see https://en.wikipedia.org/wiki/Truncated_icosahedron). Besides it is the basic structure for soccer footballs. I found out that a TruncatedIcosahedron command is available in the geom3d package. This command is able to deliver data for the faces and more. With this command I succeded in writing a small program to actually plot this polyhedron in 3D:

-----------------------------------------------------

restart;
with(geom3d):
with(plots):
TruncatedIcosahedron(football,point(C,(0,0,0)),1):

PlotFootball:=proc(tr::float)
    
    local i::integer,
         FootballFace::Vector(datatype=float[8]),
         plotFace::Vector(datatype=float[8]);    
    
    for i from 1 to 32 do
        FootballFace[i]:=faces(football)[i];
        plotFace[i]:=polygonplot3d(Matrix(FootballFace[i]),axes=none,scaling=constrained,transparency=tr);
    end do;
    
    display(seq(plotFace[i],i=1..32)):
    
end proc:

PlotFootball(0.00);

-----------------------------------------------------

Since I am not really experienced in programming in Maple, here is my last question:

2. Can I simplify something in my code above?

 

Best wishes,

Erik

 

People im with error to show a plot, this is a wor for universty about mechanic materials, and i'm with this error... do not show the plot... before, i've a error because the variable x has values, and i've to unassign to have this.. without unassign, has error in domin of the plot... i dont know how to do this ! And ive to send the work today. 

Someone can help me about this?


Value of variables:
functions
plot algoritm...
plot algoritm

Thanks !

Hi

I can't understand difference between plots!

Please expailn it.

Thanks.

plot({sqrt(x+2*sqrt(x-1))+sqrt(x-2*sqrt(x-1)), sqrt(x-2*sqrt(x-1)), sqrt(x+2*sqrt(x-1))}, x = -3 .. 3)

 

Hi, I am trying to figure out how to find the maxima and minima of the function f(x)=(x^3-10*x^2-2*x+1)/(4*x^3+5*x+1) and its derivative.

a) Find the local maxima and minima of f. 

b) Find the vertical asymptote

 

I graphed the function and located the  local max to be (-1,1) and local min (1,-1) however, I can't figure out how to plot the derivative of the function. I tried to use the diff expression but it leaves the graph empty when i attempt to plot it.

Also, just curious how am i able to join text words and maple commands together? For example, I want to be able to say "The derivative of the original function is [maple command here]" ? I tried to switch from text to math however nothing worked.

 

Thanks in advance!

1 2 3 4 5 6 7 Last Page 1 of 163