@tomleslie 

Dear Tom,

I need your help. I have a delay differential equation to solve and extract the value of the solution y(t) at a selected point of the independent variable t.  I am uploading a small sample code. 

Thanks.
 

Download delay-differential-equation.mw
 

Download delay-differential-equation.mw

Hi everyone: 

I want to obtain r1(t), V(t) , q(t) in terms of U0 and plot V(t) in terms of U0, how? 

eq1:= diff(r[1](t), t, t)+(.3293064114+209.6419478*U[0])*(diff(r[1](t), t))+569.4324330*r[1](t)-0.3123434112e-2*V(t) = -1.547206836*U[0]^2*q(t)
eq2:= 2.03*10^(-8)*(diff(V(t), t))+4.065040650*10^(-11)*V(t)+0.3123434112e-2*(diff(r[1](t), t)) = 0
eq3 := diff(q(t), t, t)+1047.197552*U[0]*(q(t)^2-1)*(diff(q(t), t))+1.096622713*10^6*U[0]^2*q(t) = -2822.855019*(diff(r[1](t), t, t))
ics:=r1(0)=0,D(r1)(0)=0,V(0)=0,q(0)=0,D(q)(0)=0;

Tnx...

Hi, everytimes I enter anythings, it turn out with Typesetting:-mparsed( bla1bla2, bla1bla2;"_noterminate")

example:

For an Array A, say, and some positive integer n, say, Maple interpretes A^n as raising each entry separately to the same power n. Without the Physics package loaded, A^n can also be written as A . A . ... . A (n times). But with the Physics package loaded, this equality is broken (at least in Maple 2017): If A is a 2D square Array, A . A all of a sudden is no longer equal to A^2, but rather to convert(A,Matrix)^2, i.e., to the square of the Array considered as a Matrix. The presence of the dot operator seems to make the Physics enviroment convert A to a Matrix. This seems to me to be a bug.

During running my ws I faced with memory error as below, where as my system have enough memory (120GB)

Warning, Run: unable to set assignto result due to error:  Maple was unable to allocate enough memory to complete this computation.  Please see ?alloc

Maple's help suggests :Software limits are imposed by the -T command-line argument, the datalimit argument to kernelopts and system imposed user limits (for example shell limits).
  But I could not understand how to increase software limit.

how to fix that?

find the fourier cosine and sine transform. Use your identites and integrals to eliminate all exponents.

piecewise function is {cosh(t) for 0 less or equal to t, t is less than 1,

                                   0 for 1is less than or equal to t.}

alpha+{6*RK[1]*alpha+2+(96/5)*R^2*K^2*alpha^2-(1/6)*R*alpha+64*R^3*K^3*alpha^3}*beta+{(24/5)*RK+44*R^2*K^2*alpha^3}*beta^2=0

i am currently using Maple 18, i have a problem on inserting 12 row by 12 column matrix and above is seem to be impossible. please can help or direct me on how to insert 12 row by 12 column matrix in maple. because my maple 18 seem to stop in 10 row by 10 column matrix.  thanks

Hi everybody.....

I want to obtain the roots of below equation, How? 

eq:=1+cos(beta)*cosh(beta)+.5720823799*beta*(cos(beta)*sinh(beta)-cosh(beta)*sin(beta))-0.1285578382e-1*beta^2*sin(beta)*sinh(beta)-0.9629800618e-4*beta^3*(cos(beta)*sinh(beta)+cosh(beta)*sin(beta))+0.1377259814e-4*beta^4*(1-cos(beta)*cosh(beta)) = 0

Tnx...

Dear users,

I have an issue with finding real part of a complex variable function. In calculating the real part I see two arguments and the plot is not smooth. How to get real part correct. The worksheet is attached.
 

Download Toya_complexPlot2.mw

Hi all, how to get name file  in directory?

Example. I have  and I want call then return

["equals", "function", "polynomial", "viete"]

or

["equals.txt", "function.txt", "polynomial.txt", "viete.txt"]

Thank you very much.

Dear authors,
How to solve this ode problem.

Download link ode.mw

In this problem the boundary condition is

Note: F=g in our problem.

eta approaches N.

Thank you.

I am trying to put a number of related 2-d plots into a 3-d frame so I can see them stacked up in the third dimension (which follows a parameter) and rotate things around.

The way I once did this successfully was to create the 2-d plots and then use plottools:-transform to move the individual plots in the third dimension, like so:

plt:=plot(something);

tr:=plottools:-transform((x,y) -> [x,2,y]); # the "2" gets changed for the other plots (not shown here).

plots:-display(tr(plt));

The only effect I can get is that the GUI gets confused and I have to close and reload the sheet to get it back again. I have a (complicated) sheet where this actually works, but I am not able to make it work even in the small example I am posting below.

Any hint of where I am going off trail is appreciated. Incidentally, this problem is what led to the corrupted sheet I had maybe a week ago.

Thanks,

Mac Dude.

display3d.mw

I unprotect the GAMMA, but still receives error:

Error, attempting to assign to `GAMMA` which is protected.  Try declaring `local GAMMA`; see ?protect for details.

Download soal.mw

What is the problem?

Hi, I'm investigating an equation of motion, and I'm attempting to plot t against different values of thetabn. I have acquired a list of my data points (in the real domain) of the form [thetabn1, t1], [thetabn2, t2] ... etc. But I am struggling to figure out how to plot it. Using dataplot gives me a graph but it is incorrect. I would like a scatterplot, ideally once I have refined my range of thetabn I would like to be able to join it up with a line to create a nice looking plot. Also any improvements to my method would be appreciated. v[i], thetavn and omegac are set variables and s[n] is my function which I am setting equal to zero, and solving for t.

with(RealDomain);
v[i] := 145000;
thetavn := (1/6)*Pi;
omegac := .1;
s[n] := v[i]*cos(thetavn)*(cos(2*thetabn)*tan(thetabn)+sin(2*thetabn)*sin(omegac*t)/omegac);
for thetabn from evalf((1/100)*Pi) by evalf((1/100)*Pi) to evalf(Pi) do assign('result', [op(result), [thetabn, solve(v[i]*cos(thetavn)*(cos(2*thetabn)*tan(thetabn)+sin(2*thetabn)*sin(omegac*t)/omegac) = 0, t)]]) end do;
dataplot(result);