## how to get nonempty plot?...

Asked by:

I would like to plot h over the real interval from 0 to 2*Pi, but

produces an empty plot...  How can I plot successfully?

## Is it possible to write a caption of a plot in lat...

Asked by:

Hello every body

I have a plot that contains two curve. I need to chracterize the curves with `min(D_{T})` and `max(D_{E})`. Is it possible to write legend of plot such that when include it in latex file,  T and E be presented in indices?

## Multiple plot in one figure...

Asked by:

Hello every body

I have three vectors `X`, `Y` and `Z`. I need to plot `X` vs `Z` and `Y` vs `Z` in one plot. How can I do this?

Thanks in advance

## How do I solve and plot a coupled differential eq...

Asked by:

,Hello friends

.I want to plot a coupled differential equation numerically  but I dont know whats wrong with the loops

numeric4.mw  .pdfnumeric4.pdf

## Evaluating expressions inside an Explore command...

Asked by:

Dear Maple users

Let's say I want to plot a line containing a parameter, for example (a+2)*x+6, where a is the parameter. When using the following Explore command:

Explore(plot((a+2)*x+6, x = -5 .. 5, view = -10 .. 10), parameters = [a = -4.0 .. 3], initialvalues = [a = 1])

I will be able to change the value of a by using the slider. The actual value of a is currently displayed as well. Now, However, I want to have the current value of the slope of the line displayed as I move the slider. How is this done? When using a print command inside the Explore command the value of a+2 is not displayed, rather "a+2". Can I use some kind of component?

I hope someone can help.

NB! I constructed this rather doll example, just to make it simple.

Regards, Erik

## how to plot this?...

Asked by:

ode1a := diff(y1(t), t) = round(rhs(odeparm1[1][1]))*y1(t)+round(rhs(odeparm1[1][2]))*y2(t)+round(rhs(odeparm1[1][3]))*y3(t);
ode2a := diff(y2(t), t) = round(rhs(odeparm1[1][4]))*y1(t)+round(rhs(odeparm1[1][5]))*y2(t)+round(rhs(odeparm1[1][6]))*y3(t);
ode3a := diff(y3(t), t) = round(rhs(odeparm1[1][7]))*y1(t)+round(rhs(odeparm1[1][8]))*y2(t)+round(rhs(odeparm1[1][9]))*y3(t);
try
ode1a := diff(y1(t), t) = rhs(odeparm1[1][1])*y1(t)+rhs(odeparm1[1][2])*y2(t)+rhs(odeparm1[1][3])*y3(t);
ode2a := diff(y2(t), t) = rhs(odeparm1[1][4])*y1(t)+rhs(odeparm1[1][5])*y2(t)+rhs(odeparm1[1][6])*y3(t);
ode3a := diff(y3(t), t) = rhs(odeparm1[1][7])*y1(t)+rhs(odeparm1[1][8])*y2(t)+rhs(odeparm1[1][9])*y3(t);
sys := DiffEquation([ode1a, ode2a, ode3a], inputvariable = [y1(t)], outputvariable = [y2(t), y3(t)]);
sysz := ToDiscrete(sys, ts); in_t := Sine(1, 1, 0, 0);
sol := Simulate(sys, [in_t]);
try
p1 := plots[odeplot](sol, [[t, y2(t)]], t = 0 .. t_sim, numpoints = 200, color = red);
print("succeed 1 2", i)
catch:
print("error draw at ", i)
end try;
try
p1 := plots[odeplot](sol, [[t, y3(t)]], t = 0 .. t_sim, numpoints = 200, color = red);
print("succeed 1 3", i)
catch:
print("error draw at ", i)
end try
catch: print("error at ", i);
print(lastexception);
print(ode1a);
print(ode2a);
print(ode3a);
end try;
try
ode1a := diff(y1(t), t) = rhs(odeparm1[1][1])*y1(t)+rhs(odeparm1[1][2])*y2(t)+rhs(odeparm1[1][3])*y3(t);
ode2a := diff(y2(t), t) = rhs(odeparm1[1][4])*y1(t)+rhs(odeparm1[1][5])*y2(t)+rhs(odeparm1[1][6])*y3(t);
ode3a := diff(y3(t), t) = rhs(odeparm1[1][7])*y1(t)+rhs(odeparm1[1][8])*y2(t)+rhs(odeparm1[1][9])*y3(t);
sys := DiffEquation([ode1a, ode2a, ode3a], inputvariable = [y2(t)], outputvariable = [y1(t), y3(t)]);
sysz := ToDiscrete(sys, ts);
in_t := Sine(1, 1, 0, 0);
sol := Simulate(sys, [in_t]);
try
p1 := plots[odeplot](sol, [[t, y1(t)]], t = 0 .. t_sim, numpoints = 200, color = red);
print("succeed 2 1", i)
catch:
print("error draw at ", i)
end try;
try
p1 := plots[odeplot](sol, [[t, y3(t)]], t = 0 .. t_sim, numpoints = 200, color = red);
print("succeed 2 3", i)
catch:
print("error draw at ", i)
end try
catch:
print("error at ", i);
print(lastexception);
print(ode1a);
print(ode2a);
print(ode3a)
end try;
try
ode1a := diff(y1(t), t) = rhs(odeparm1[1][1])*y1(t)+rhs(odeparm1[1][2])*y2(t)+rhs(odeparm1[1][3])*y3(t);
ode2a := diff(y2(t), t) = rhs(odeparm1[1][4])*y1(t)+rhs(odeparm1[1][5])*y2(t)+rhs(odeparm1[1][6])*y3(t);
ode3a := diff(y3(t), t) = rhs(odeparm1[1][7])*y1(t)+rhs(odeparm1[1][8])*y2(t)+rhs(odeparm1[1][9])*y3(t);
sys := DiffEquation([ode1a, ode2a, ode3a], inputvariable = [y3(t)], outputvariable = [y1(t), y2(t)]);
sysz := ToDiscrete(sys, ts);
in_t := Sine(1, 1, 0, 0);
sol := Simulate(sys, [in_t]);
try
p1 := plots[odeplot](sol, [[t, y1(t)]], t = 0 .. t_sim, numpoints = 200, color = red);
print("succeed 3 1", i)
catch:
print("error draw at ", i)
end try;
try
p1 := plots[odeplot](sol, [[t, y2(t)]], t = 0 .. t_sim, numpoints = 200, color = red);
print("succeed 3 2", i)
catch:
print("error draw at ", i)
end try
catch:
print("error at ", i);
print(lastexception);
print(ode1a);
print(ode2a);
print(ode3a)
end try

diff(y1(t), t) = 1.052936200*10^5*y1(t)+70106.19000*y2(t)+35169.00000*y3(t)
diff(y2(t), t) = 70106.19000*y1(t)+71031.61000*y2(t)+35511.00000*y3(t)
diff(y3(t), t) = 35169.00000*y1(t)+35511.00000*y2(t)+36100.00000*y3(t)
"the DEs contain functions with undefined values (probably caused by a discontinuity in the input that was differentiated). As a result, the numerical solution cannot be calculated. The DE system is: %1\"",[(&DifferentialD;)/(&DifferentialD;t) y1(t)=1.052936200 10^5 y1(t)+70106.19000 y2(t)+35169.00000 ({[[0,t<0],[sin(t),otherwise]]),(&DifferentialD;)/(&DifferentialD;t) y2(t)=70106.19000 y1(t)+71031.61000 y2(t)+35511.00000 ({[[0,t<0],[sin(t),otherwise]]),{[[0,t<0],[undefined,t=0],[cos(t),0<t]]=35169.00000 y1(t)+35511.00000 y2(t)+36100.00000 ({[[0,t<0],[sin(t),otherwise]]),y2(0)=0,y1(0)=0]

it has error when plot

## unable to evaluate the function to numeric values ...

Asked by:
ygraph1 := -.736312023696564122*exp(2.26140104440167664*10^5*tt)-.591826613918776445*exp(28994.5376895644186*tt)+.328002839648234568*o*exp(13767.7178702679158*tt);
ygraph2 := -.591859486202007235*exp(2.26140104440167664*10^5*tt)+.328381376616263988*exp(28994.5376895644186*tt)-.736116852194203974*o*exp(13767.7178702679158*tt);
ygraph3 := -.327943520064913564*exp(2.26140104440167664*10^5*tt)+.736143281263262450*exp(28994.5376895644186*tt)+.592069351595225779*o*exp(13767.7178702679158*tt);
FunctionAdvisor(branch_points, ygraph1);
plot(ygraph1, tt=-5..5);
plot(ygraph2, tt=-5..5);
plot(ygraph3, tt=-5..5);
Warning, unable to evaluate the function to numeric values in the region;

how to plot this system?

## how to deplot this case?...

Asked by:
ode1a := diff(y1(tt), tt) = 1.342398800*10^5*y1(tt)+89591.20000*y2(tt)+44647.44000*y3(tt);
ode2a := diff(y2(tt), tt) = 89591.20000*y1(tt)+89803.24000*y2(tt)+44901.60000*y3(tt);
ode3a := diff(y3(tt), tt) = 44647.44000*y1(tt)+44901.60000*y2(tt)+44859.24000*y3(tt);

DEplot3d({ode1a,ode2a,ode3a}, {y1(tt), y2(tt), y3(tt)}, tt=0..10,[[y1(0) = 0, y2(0) = 0, y3(0) = 0]],scene=[tt,y1(tt),y2(tt)]);
DEplot3d({ode1a,ode2a,ode3a}, {y1(tt), y2(tt), y3(tt)}, tt=0..10,[[y1(0) = 0, y2(0) = 0, y3(0) = 0]],scene=[tt,y1(tt),y3(tt)]);
DEplot3d({ode1a,ode2a,ode3a}, {y1(tt), y2(tt), y3(tt)}, tt=0..10,[[y1(0) = 0, y2(0) = 0, y3(0) = 0]],scene=[tt,y2(tt),y3(tt)]);

can it plot 3 functions ?
and why it return a straight line 3d graph

is there some interesting graph from this system?

## Error, (in plot) unexpected options...

Asked by:

Dear all,

I'm trying to plot dispersion curves in Maple but I'm having some trouble. The code is attached as a file also.

f1  is my main function and I want to plot Vx with regard to f as all the rest are known variables, but I'm not able to do so.

I hope that someone can tell me why do I have this error.

Thanks in advance!

restart;
f1 := (C33*Rp*kzp+C13*kx)*(Rm*kx+kzm)*sin(kzp*h)*cos(kzm*h)-(C33*Rm*kzm-C13*kx)*(Rp*kx+kzp)*sin(kzm*h)*cos(kzp*h)=0;
Rp := (-C11*kx^2-C55*kzp^2+omega^2*rho)/((C55+C13)*kx*kzp);
Rm := (-C11*kx^2-C55*kzm^2+omega^2*rho)/((C55+C13)*kx*kzm);
kzp := sqrt(((-M+sqrt(M^2-4*N))*(1/2))*kx^2);
kzm := sqrt(((-M-sqrt(M^2-4*N))*(1/2))*kx^2);
M := (C11*C33/rho^2-2*C55*C13/rho^2+C13^2-omega^2*(C33+C55)/(rho*kx^2))*rho^2/(C33*C55);
N := (omega^2/kx^2-C11/rho)*(omega^2/kx^2-C55/rho)*rho^2/(C33*C55);
C11 := 0.435e10;
C13 := 0.259e10;
C55 := 0.112e10;
C33 := 0.108e11;
rho := 923;
h := 0.7e-2*(1/2);
kx := omega/Vx;
omega := 2*Pi*f;
f1;
plot(f1, f = 10 .. 0.100e6);
Error, (in plot) unexpected options: [(.4633081900*(-0.1717311166e12*f^2/Vx^2-0.2210791386e11*(...

Asim_dispers.mw

## plot/ zoom and change the axises...

Asked by:

hi

I have an equation and I want to plot this in maple.

for example, this is my equation:

w=-4.918231636*z^3+2.486861786*z^2+.8573483099*z+2.341623514

1) I want to plot z/H (H is a number) versus w so I need to switch the axises.

2) I want to plot specific area (by using commands) for example vertical axis from -0.5..0.5 and horizontal from -2..3.

i don not want to use plot property because it should be aligned for every run and it is not permanent.

thank you very much

## is it possible to evaluate a floating number and p...

Asked by:

unknown := sum(exp(product(sum(kk, kk=1..jj),jj=1..y))*1/y!*x^y, y=0..infinity):
evalf(subs(x=1, unknown));
plot(unknown, x=-3..3, numpoints = 5);

## how to plot this function?...

Asked by:

g2 := arctanh((exp(2*y)+sqrt((exp(2*y))^2+exp(2*y)))/exp(2*y)-1)-1;
singular(g2);
FunctionAdvisor(definition, g2);
plot(g2, y=-5..5);

## How to print units in a table-structured plot ?...

Asked by:

Hi everybody,

I want to plot many curves in the form of a table.

For instance, let    p := plot(x, x=0..1, labels=[x*Unit('m'), y*Unit('m')]):

The command   plots:-display(p)  returns a plot  with labels x [[m]] and y [[m]], just as expected.

But the commands (remember I have several plots)
plots:-display(<p>);                              # and the generalization plots:-display(<p | p>)
DocumentTools:-Tabulate([p]);           # and the generalizationDocumentTools:-Tabulate([p,p])
both provide labels  x Units:-Unit('m') and  y Units:-Unit('m')

Is there a way to fix this and obtain plots withe the expected labels ?

Thank to all of you for yout help

## plot an indefinite integral ...

Asked by:

hi every one, i want to plot an indefinite integral  , it is some what complex and maple can not compute the answer, ( but numeric integration can be computed) , but we want to plot the output, what should we do ? tnx for help in advance

corrected.mw

## Changing coordinates for matrix plot...

Asked by:

I'm plotting matrices using the proc in the attached file.

How do I move the rows in the plot closer to each other?
The first picture show the rows plotted in the values [13,14...19].
I want them to be plotted in values like [13,13.5,14,14.5...].

Second picture shows how I would like it to be. Created it in paint, from a different matrix.Matrix_plot.mw

 1 2 3 4 5 6 7 Last Page 1 of 175
﻿