MaplePrimes Questions

Meanwhile, thank you so much for everything.
I know I'm asking a lot but if you have time, you can help me do this?

Building a system of interactive components that, taken a function, two points 'a' and 'b' values ​​and an integer n, the calculations point between a and b in which the function assumes the minimum value by using the following procedure:

• It divides the values ​​between a and b into n equal parts (these will distance the one with the other (b-n)/2);
• calculates the function in each of these points;
• located between these values,  what is the minimum (in case of a tie, take the one closest to a)

I think i have to create a vector for each  part and  prehaps with a fcycle for, calculate the function, finelly i'll use minimize with all function.
Do you think is the correct procedure? If yes, how can I do it?

yVal := 0.01

xVal := 0.01

p8 := plot([fdiff(('rhs')((pds:-value(f(x, y)))(x, yVal)[3]), [y, y], x = z)], z = 0 .. 20, color = [red])

I cant seem to plot the graph for yVal, but if I were to switch it around to (xVal,y), the graph works perfectly. How should I input the code such that I can obtain the f''(x,0) graph. 

Any help will be greatly appreciated :) thanks

fyp3.mw

Dear All, 

I work on the period doubling phenomenon by integrating the dynamical system at some particular parameter values. My main commands are as follows

restart:
with(DEtools): with(plots): with(plottools):
a:=-1: b:=-3: c:=3: d:=1: omega:=1: v1:=1: f:=-4: epsilon:=0.01: 
sys:=diff(u1(t),t)=v1*u1(t)-(omega+u2(t)^2)*u2(t)+u1(t)*(a*(u1(t)^2+u2(t)^2)+b*z(t)^2),
     diff(u2(t),t)=(omega+u1(t)^2)*u1(t)+v1*u2(t)+u2(t)*(a*(u1(t)^2+u2(t)^2)+b*z(t)^2),
     diff(z(t),t)=z(t)*(-v1+c*(u1(t)^2+u2(t)^2)+z(t)^2)+epsilon*z(t)*(v2+f*z(t)^4):

solC:=dsolve({eval(sys, v2=2.0500014987),
              u1(0)=0.6,u2(0)=0.6,z(0)=0.1},
              type=numeric, method=rkf45, maxfun=10^7,
              range=350..750):
p1:=odeplot(solC, [sqrt(u1(t)^2+u2(t)^2),z(t)],
               refine=2,
               color=burgundy,
               thickness=1);
 

The resulting curves are quite blurring and furry. I have tried to increase the "refine" or the "numpoints" values, but with minimal improvement. I would like to know if there a way to obtain clean and sharp integral curves?  

I would be very grateful if someone could help me! 

Thank you in advance.

With kind wishes,

Wang Zhe

Hi, I try Aproximate integral with simpsons method, but i want value in function the a with number solution:

 

with(Student[Calculus1]);
ApproximateInt(a*tan(x)-2*x, x = -100 .. 100, method = simpson, output = plot, partition = 500);

This error.

 

How we can construct conservation laws of fractional differential equation with Maple?

Hi I try this integral:

m2 := int(exp(-(1/2)*z^2)*((exp(B*J*sqrt(q)*z))^2-1)/(sqrt(Pi)*sqrt(2)*((exp(B*J*sqrt(q)*z))^2+1)), z = -infinity .. infinity)

But not resolve.

How can i do?

Regards.

Basically what I'm trying to do is use Newton's method to find the root of f(x)=sqrt(x)+ln(0.1x) starting at x0=5 within a tolerance of 0.001.

 

f := proc (x) options operator, arrow; sqrt(x)+ln(.1*x) end proc

> xk := 1.0;
print(`output redirected...`); # input placeholder
                              1.0
> for k to 6 do xk1 := xk-f(xk)/(D(f))(xk); xk := xk1 end do;
 

 

When I do that, Maple barfs out pages of nonsense when I'm looking for it to give numerical values. It seemed to work fine with a different function, so maybe that's part of the problem? But I need to do it with the function I asked about.

Respected member!
Please help me to find the solution of attached problem,  I am a new user so please forgive any mistakes.maple.mwmaple.mw

Hello,

I have been trying to solve a simple nonlinear equation. Im interested in the solution per say rather than the plot but when I browsed about the commands to use, this came up. I tried it in my case and it is giving me the following errors:

ode.mw
 

``

restart;

``

with(plots);

[animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, densityplot, display, dualaxisplot, fieldplot, fieldplot3d, gradplot, gradplot3d, implicitplot, implicitplot3d, inequal, interactive, interactiveparams, intersectplot, listcontplot, listcontplot3d, listdensityplot, listplot, listplot3d, loglogplot, logplot, matrixplot, multiple, odeplot, pareto, plotcompare, pointplot, pointplot3d, polarplot, polygonplot, polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus, semilogplot, setcolors, setoptions, setoptions3d, shadebetween, spacecurve, sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]

(1)

eq5:=C*sqrt(y(x)*((diff(y(x),x))^2+1))-y(x)=0;

C*(y(x)*((diff(y(x), x))^2+1))^(1/2)-y(x) = 0

(2)

C:=1;

1

(3)

bcs:=y(-1)=1, y(1)=1;

y(-1) = 1, y(1) = 1

(4)

dsys:={eq5,bcs};

{(y(x)*((diff(y(x), x))^2+1))^(1/2)-y(x) = 0, y(-1) = 1, y(1) = 1}

(5)

dsol:=dsolve(dsys, numeric); odeplot(dsol,[x,y(x)],0..1,color=red,axes=box);

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

 

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

 

``


 

Download ode.mw

 

The question explores the family of differential equations dy/dx = sqrt(1+(a*x)+(2*y)) for various values of the parameter a.

This figure shows the tangent field in the case a=1.

sketch a tangent field in the case a=-2.

need the following diagram on maple:

 

hello,

Given f(x)=2x^2-4x, how do I define and plot the functions f(x), 2f(x), 4f(x), -f(x) and -4f(x) in maple 15?

Thank You.

 

How wired icon of motor to icon of sources?thanks

Can anyone solve,

eqn := diff(x(t), t, t)+omega^2*x(t)-epsilon*mu*(diff(x(t), t))+epsilon*alpha*x(t)^2*(diff(x(t), t))+epsilon*Zeta*cos(Omega*t)*x(t)=0;

This equation using MMS (Method of Multiple Scales) in maple.

Pls use Polar transformation function as (1/2)*a(t)*exp(-I*beta(t))

Please reply.

Hi, I try expand in serie this function:

mtaylor(tanh(B*J*sqrt(q)*z+B*Jo*m), [q, m], 9)

 

But this is error.

 

How can i do?

Best Regards

Compare seq(-GAMMA(k+1, x), k = 0 .. -2) and seq(-GAMMA(k+1, x), k = -2 .. 0)

sum(-GAMMA(k+1, x), k = -2 .. 0) = -1/x*Ei(2,x)-Ei(1,x)-exp(-x) may be what you expect

First 886 887 888 889 890 891 892 Last Page 888 of 2308