gkokovidis

2355 Reputation

13 Badges

20 years, 132 days
Draeger Medical Systems, Inc.

 

 

Regards,
Georgios Kokovidis
Dräger Medical

MaplePrimes Activity


These are answers submitted by gkokovidis

Can you upload your worksheet that causes this error.  Use the Green Up-arrow  in the menu.  Maple version?  

Try opening the attached and run it.
 

Download Intsin(x).mw

Regards,

Georgios

 

Round-off error. Try increasing Digits to 20 or higher.

 

Digits:=20;

at the top of your worksheet.

Regards,

Georgios

You have a system of two 2nd degree equations, with two unknowns.  Your solutions look correct and match the intersection of the plot, as in the attachment.

Download implicitplotSol.mw
 

Regards,

Georgios

 

restart;
with(plots):
with(plottools):
k := disk([0, 0], 3, color = white):

k1 := polygon([[0, 0], [3, 4], [3, 1]]):
display([k, k1]);

Regards,

Georgios

I am guessing as to what is desired, but here is my interpretation of the question.

Assuming the imported data looks something like this constructed example below:

restart:

A:=Matrix(3, shape = identity);

Matrix(3, 3, {(1, 1) = 1, (1, 2) = 0, (1, 3) = 0, (2, 1) = 0, (2, 2) = 1, (2, 3) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 1})

(1)

subs(0=``,A);

"[[[1,,],[,1,],[,,1]]]"

(2)

 

Download NULL.mw

Regards,

Georgios

Looks like an issue with the ResponsePlot command.  Here is a workaround.

restart:

with(inttrans):
H1:=1.000000000*10^10/(s + 100000.)^2;
H2:=1/(1 + 6.319637500*10^(-13)*s^2 + 1.1240*10^(-6)*s);
input:=laplace(Heaviside(t),t,s):
yStep1:=invlaplace(input*H1,s,t);
yStep2:=invlaplace(input*H2,s,t);
plot({yStep1,yStep2},t= 0..0.00005,
     title=`step response`,
     labels=[`time`,`y(t)`],gridlines=true);

0.1000000000e11/(s+100000.)^2

 

1/(1+0.6319637500e-12*s^2+0.1124000000e-5*s)

 

1.-1.*exp(-100000.*t)*(1.+100000.*t)

 

1.-0.7908395474e-6*exp(-889291.5140*t)*(1263927.379*sin(889679.6304*t)+1264479.*cos(889679.6304*t))

 

 

 

Download stepfunc.mw

Regards,

Georgios

Your X and Y are the same, so plotting one against the othter does not make much sense unless you are looking for a straight line (which you are not).

The attached plots one of your data sets.  Is this what you want?

Plot.mw
You can also use the command dataplot, like this:

dataplot(1..110,X);

which will yield the same plot as in the attached file.


Regards,

Georgios

 

Which version of Maple are you using, and what commands have you entered that give you a parabola?

restart:
with(plots,implicitplot):
implicitplot(x^2+y^2=1,x=-1..1,y=-1..1);

will draw a circle.  Depending on your version of Maple, you might have to use the scaling=constrained option at the end, like this:

implicitplot(x^2+y^2=1,x=-1..1,y=-1..1,scaling=constrained);

Can you post your code so we have an idea of what is causing you problems.

Regards,

Georgios

Maple does show an answer.  odeplot does not show your plot because you are missing with(plots) at the beginning of your worksheet.

If you do not want to load the plots package, call odeplot like this:

plots:-odeplot(sol_num, [t, phi1(t), epsilon(t)], t = 0 .. 0.5)

The error is due to how you are using odetest.  Try commenting it out for now, and rerun you worksheet.

Regards,

Georgios

 

From the Maplesoft site:  https://www.maplesoft.com/applications/view.aspx?SID=3473&view=html

Use this as a starting point.

 

//Georgios

Adding to what was answered already, with an example that you can modify for your own puposes.

restart:with(Student[Calculus1]):

eq := sqrt(3)*exp(-t/2)*sin(sqrt(3)*t/2)/3;
 

(1/3)*3^(1/2)*exp(-(1/2)*t)*sin((1/2)*3^(1/2)*t)

                                                                                                                                                           (1)

plot(eq,t=0..20,y=0.3..-0.3);

 

ans:=Roots(eq,0..20);

[0, (2/3)*Pi*3^(1/2), (4/3)*Pi*3^(1/2), 2*Pi*3^(1/2), (8/3)*Pi*3^(1/2), (10/3)*Pi*3^(1/2)]

                                                                                                     (2)

NumberOfZeroCrossingsOverGivenRange:=nops(ans);

6

                                                                                                                              (3)

 

Download ZeroCrossing.mw

Upload your code, so we can look at it.  

If you are creating a code region in Maple, and want to have multiple lines, then you must use the "Shift-Enter" key combination to go to the next line, and continue to do so until the end of your program.

You can use an external editor, save your file, and read it into Maple.

Another option is to use an external editor, and then cut and paste the code, and then execute.

Regards,

Georgios

Below is my interpretation of what the OP is trying to do.  I'm not sure if the word Vector is relevant here.

The attached creates a sequence of points, like acer shows, spaced at 0.5.  Then there is a function mapped to these

points ( 5*x^2+3).  Then x and y are plotted against each other.   

restart:

x:=[seq(0..20,0.5)];

[0, .5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5, 7.0, 7.5, 8.0, 8.5, 9.0, 9.5, 10.0, 10.5, 11.0, 11.5, 12.0, 12.5, 13.0, 13.5, 14.0, 14.5, 15.0, 15.5, 16.0, 16.5, 17.0, 17.5, 18.0, 18.5, 19.0, 19.5, 20.0]

(1)

y:=[seq(5*i^2+3,i=x)];

[3, 4.25, 8.00, 14.25, 23.00, 34.25, 48.00, 64.25, 83.00, 104.25, 128.00, 154.25, 183.00, 214.25, 248.00, 284.25, 323.00, 364.25, 408.00, 454.25, 503.00, 554.25, 608.00, 664.25, 723.00, 784.25, 848.00, 914.25, 983.00, 1054.25, 1128.00, 1204.25, 1283.00, 1364.25, 1448.00, 1534.25, 1623.00, 1714.25, 1808.00, 1904.25, 2003.00]

(2)

xylist := zip( (x,y) -> [x,y], x, y);

[[0, 3], [.5, 4.25], [1.0, 8.00], [1.5, 14.25], [2.0, 23.00], [2.5, 34.25], [3.0, 48.00], [3.5, 64.25], [4.0, 83.00], [4.5, 104.25], [5.0, 128.00], [5.5, 154.25], [6.0, 183.00], [6.5, 214.25], [7.0, 248.00], [7.5, 284.25], [8.0, 323.00], [8.5, 364.25], [9.0, 408.00], [9.5, 454.25], [10.0, 503.00], [10.5, 554.25], [11.0, 608.00], [11.5, 664.25], [12.0, 723.00], [12.5, 784.25], [13.0, 848.00], [13.5, 914.25], [14.0, 983.00], [14.5, 1054.25], [15.0, 1128.00], [15.5, 1204.25], [16.0, 1283.00], [16.5, 1364.25], [17.0, 1448.00], [17.5, 1534.25], [18.0, 1623.00], [18.5, 1714.25], [19.0, 1808.00], [19.5, 1904.25], [20.0, 2003.00]]

(3)

plot(xylist);

 

 

Download plotFuncData.mw

@digerdiga The procedure by Preben Alsholm, from a previous posting here, does what you are looking for.

Regards,

Georgios

restart: with(plots):

f:=(x)->-(2.50000000000000000000000000000*(-2.424176758*10^30*x^5+5.982448001*10^29*x^3+(5.336360731*10^33*I)*exp((4.0405213687327619451756171962*I)/x)+(3.144879905*10^98*I)*exp((48.6318841588794957405477169504*I)/x)+(8.489643736*10^35*I)*exp((4.04052136873276194517561719622*I)/x)))*x^6/((1.009157478*10^101*I)*exp((48.6318841588794957405477169504*I)/x)*x^6-2.760182453*10^32*x^11-(7.651751413*10^37*I)*exp((4.0405213687327619451756171962*I)/x)*x^4+1.887831524*10^32*x^9-(8.567811466*10^36*I)*exp((4.04052136873276194517561719622*I)/x)*x^4-(1.207009672*10^100*I)*exp((48.6318841588794957405477169504*I)/x)*x^4+(6.538572988*10^29*I)*x^6+(1.246166933*10^31*I)*x^10-1.837651455*10^31*x^7+(1.504191760*10^38*I)*exp((4.04052136873276194517561719622*I)/x)*x^6+1.574891043*10^34*exp((4.0405213687327619451756171962*I)/x)*x^5-(7.612867062*10^30*I)*x^8+4.556143298*10^99*exp((48.6318841588794957405477169504*I)/x)*x^5-(2.488678360*10^105*I)*exp((52.6724055276122576857233341466*I)/x)+6.791123646*10^36*exp((4.04052136873276194517561719622*I)/x)*x^5-3.454612059*10^36*exp((4.0405213687327619451756171962*I)/x)*x^3-4.710342238*10^98*exp((48.6318841588794957405477169504*I)/x)*x^3+(3.488285966*10^35*I)*exp((4.0405213687327619451756171962*I)/x)*x^6+5.512268188*10^106*exp((52.6724055276122576857233341466*I)/x)*x);

complexplot(f,0.2..0.4);

Regards,

Georgios

 

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