Items tagged with plot plot Tagged Items Feed

Hello,

I am trying to plot 3d scattered data as a 2d scatterplot with the third dimension being displayed as colour or a line that separates the value of the two functions (for example, the higher value and lower value at that point). I attached the file. I have two functions and several parameters that are constant. Here, I want to use P[A] on the x-axis, and P[B] on the y-axis. The z dimension is the profit functions, either Profit[1] or Profit [2]. I want to show them in colors like the highest one is red and the lower one is green, or something similar.  P[A] and P[B] are probability values and P[B] is between 0 to 0.50 with 0.01 increments, while P[A] is always less than P[B]. The file is here: example_file.mw

Example;

P[A] := 0 ..P[B]-0.01

P[B] := 0 ..0.50

Anyone please help me to code this problem.Thanks,

YC

Hello all,

I am pretty new to Maple, but I am trying to understand something. I was able to find a numerical solution to a PDE in maple (with some community help, thanks guys!). I am trying to manipulate this data but am struggling with it. I thought if I could take this data to matlab it would be pretty easy for me to manipulate and do what I want. 

So my question is: How do I export my numerical solution (pds module) to matlab. Just taking the data is okay. I know you can evaluate the data at some points. 

 

I see there is "Matlab" command that converts code, it doesn't seem to like pds as an input though.

I also see an export matrix command. I guess it could be possible to create a matrix of plot data and convert it this way?

 

I had to do a change of coordinate system to solve the PDE because of boundary conditions. I'm trying to transform this data back to my regular x,y coordinate system to see if it matches some other simulations. 

Thanks in advance! And here is my file.


restart;
with(Physics):
Setup(mathematicalnotation = true):

V(Z, f);

If function ( varphi) is defined, use this one.

PDE1 := subs[inplace](x = xi*varphi(t), PDE1);
If function ( varphi) is defined, use this one.


PDE1 := subs[inplace](t = 2*vartheta*(1/omega), PDE1);


 

 

 

 

 

 

Thesis_Pde2_attempt.mw

how to plot ln(sqrt(x-12)/(-x^(2)+15x)

Whenever I try to plot I get 

"Warning, unable to evaluate the function to numeric values in the region; complex values were detected".

HELP PLEASE

The scale of vertical axis is logarithmic and function in each interval must be linear. When I use the command " mode=log " the shape of lines are deformed into curves, however when I use the command " curve ", the shape of lines remains straight. I can use right-click for changing mode to logarithmic scale but is it possible to use label font for vertical axis as well as size, resolution and other options in the command "curve" ?

Thanks alot

 

 

restart; with(plottools):

SN := curve([[3, (4.566256-4.544647)/(4.544647)], [4, (4.544933-4.544647)/(4.544647)], [5, (4.544653-4.544647)/(4.544647)], [6, (4.544649-4.544647)/(4.544647)], [7, 0]]): PLOT(SN);

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

for m to 5 do

x[m] := m+2 end do:

y[1, 1] := (4.566256-4.544647)/(4.544647):y[2, 1] := (4.544933-4.544647)/(4.544647): y[3, 1] := (4.544653-4.544647)/(4.544647): y[4, 1] := (4.544649-4.544647)/(4.544647): y[5, 1] := 0:

for j to 1 do for i to 4 do

L[i, j] := (y[i+1, j]-y[i, j])*(x-x[i])+y[i, j] end do:

PW[j] := piecewise(`and`(x >= x[1], x < x[2]), L[1, 1], `and`(x >= x[2], x < x[3]), L[2, 1], `and`(x >= x[3], x < x[4]), L[3, 1], `and`(x >= x[4], x < x[5]), L[4, 1]) end do:

plot(PW[1], x = x[1] .. x[5], color = black, axes = framed, font = [Times, 13], size = [650, 550], resolution = 1200, legendstyle = [location = top], legend = ['gamma' = -.4], linestyle = solid, thickness = 2, labels = [Number*of*Basis, '(P-P[min])/P[min]'], labelfont = [Times, 13], axis[2] = [mode = log]);

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

I have the following code:

restart;
PDE := diff(u(x, t), t) = diff(u(x, t), x, x)-sin(x+t)+cos(x+t);
IBC:= D[1](u)(0,t)=-sin(t),
D[1](u)(1,t)=-(u(1, t))^4+(cos(1+t))^4-sin(1+t),
u(x,0)=cos(x);

pds := pdsolve( PDE, [IBC], numeric, time = t, range = 0 .. 1,

spacestep = 0.1e-1, timestep = 0.1e-1,
errorest=true
)

 

And I want to plot the difference |pds(t,x) - cos(x+t)| in maple for x=1 and t=0..5

 

I thought to use the following piece of commands but I get an error:

P:=unapply(pds, t,x);
Q:=abs(P(t,x)-cos(x+t));
Q:-plot(x=1, t=0..1);

I get an error that Q isn't a module.

I thought that unapply is used for this case, can you help me with this simple task?

 

Thanks in advance.

 

hi all

 

i have tow equations and i want plot q vs.A , that A change beetwin 0..phi. vn and a are definite and positive.

plz plz help me

 

(diff(f,y))^2-(f[y]^4)/4+q^2/(f[y]^2)+(1+2*q^2)*(f[y]^2)/2=piecewise(y<-a,1/4+2*q^2,-a<y<a,vn*(f[y]^2)-vn+1/4+2*q^2,y>a,1/4+2*q^2);

phiB[infinity]=q*int(f,y=-infinity..infinity):
A:=phiB[infinity]:

Good Morning,

This is the script that I have used to to an integration, as the command ''int'' didn't work because the coputation was really heavy.


p_old := 0;

g := 0; #index

E_tot_net_Point := 0;

for g to 2 do

 h := g+1;

 loopstart := SP_time[g];

 loopend := SP_time[h];

  for tt from loopstart by 0.01 to loopend do

   t := tt;

   P_tot_net_Fkt := t -> C_motore[g]*omega_motore[g]

if t≠SP_time[g] and t ≠ SP_time[h] then

 P_tot_net_Point := P_tot_net_Fkt(t);

 P_tot_net_Point_next := P_tot_net_Fkt(t+0.001);

 plt_P_tot_net_plt[t] := pointplot([t, P_tot_net_Point]);

 E_tot_net_Point := E_tot_net_Point+(P_tot_net_Point*0.01);

 plt_E_tot_net_plt[t] := pointplot([t, E_tot_net_Point]);

 p_old := P_tot_net_Point;

end if;      

  end do:

end do:

plots:-display(entries(plt_E_tot_net_plt, 'nolist'), title = "Energy", color = "Red", symbol = point, symbolsize = 20, gridlines = true, size = [1000, 500], labels = ['[s]', '[Watt*s]'], font = [axes, 12])


I have plotted the function as a series of points over the orizontal axis (t).

 

Now I would like to calculate the last value of the curve (which is itself a point) as a function of the orizontal axes t. How can I do it?

Something like E_tot_net_Point = f(t)

 

Thank you for the help

How I can sketch the helix with parametric equations x=2cost  y=sint  z=t  and the line with parameric equation x=-2t  y=1  z=(pi/2)+t   on a three-dimensional coordinate system?

How I can graph parabolic cylinder y=x^2 and elipsoid x^2+4y^2+4z^2=16 on a three-dimensional coordinate system?

Hello

I have an SEIR model.

Equation 5 is for disease death but I would like to plot the cumulative numbers of disease death which will be the integral of Equation 5. I added the integral inside odeplot but it is not working. Any idea  about  how to compute the integral ?

Maple code is attached

Thank you

code.mw

Hi,

 

I was wondering how to go about plotting a Fourier Tranform in Maple.

My assignment is to plot a simple harmonic equation as a Fourier transform, depicting amplitude against fequency.

I've been given: x'' + w^2 x = 0

And want to obtain both the f(x) = a0 sin(wt) + b0cos(wt) form, and a graph of the the amplitude (c^2 =a0^2 + b0^2) against frequency.

I know how to do this on paper but not in Maple, so any help with line commands and layout would be very much appreciated.

 

Thanks

Hi,

  I would need to get access to the values of euler angles defining the orientation of a 3D plot.

I generate the 3D plot with a given orientation, then the user is allowed to rotate it with mouse orbit. How can I get the values of the new euler angles to use them in further computations?

In the plot/structure I cannot find any data concerned with orientation...

 

Thank you.

 

Marco

Hello,

I would like to plot a vertical line each time my function is null and is increasing.

Here an extract of my code:

v:=unapply(H*sin(w*t),t);
L:=0.080;
H:=0.020;
Vf:=0.3;
w:=10;
fsolve(v(t)=0,t=0.5..2);
zeros:=Roots(v(t),t=0.5..2);

vdot:=unapply(diff(v(t),t),t);
map(vdot,zeros);
LignesVerticales=implicitplot(x=To complete, colour=yellow,linestyle=3, thickness=2):

With Roots function, I could obtain the zero-crossings.

With the derative of the functions, I could obtain a vector which gives me when the derivative is positive or not.

I would like now to obtain the list of the zero-crossing values which are increasing but I have difficulties to obtain it. 

By list manipulation, it should be easy to do. 

May you help to plot obtain this list so that I can plot the desired vertical lines?

Thank you for your help

Hello
I need to plot this region :

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