Items tagged with differential_equation differential_equation Tagged Items Feed

Consider the problem of a hard-hit baseball. The air-friction drag on a baseball is approximately given by the following formula

and subsequent differential equations : 

F:=-((C_d)*rho*Pi*(r)^2*v*v)/2;
m:=0.145;
v0:=65;
g:=9.81;
v_x:=diff(x(t),t);
v_y:=diff(y(t),t);
d2v_x:=-((C_d)*rho*Pi*(r^2)*(v_x)*sqrt((v_x)^2 +(v_y)^2))/(2*m);
d2v_y:=-((C_d)*rho*Pi*(r^2)*(v_y)*sqrt((v_x)^2 +(v_y)^2))/(2*m)-g;

where

C[d] is the drag coefficient (about 0.35 for a baseball)

•
rho[air] is the density of air (about 1.2 kg/
3m
)
•
r is the radius of the ball (about 0.037 m)

v is the vector velocity of the ball

Then if given that : 

Power hitters in baseball say they would much rather play in Coors Field in Denver than in sea-level stadiums because it is so much easier to hit home runs. The air pressure in Denver is about 10% lower than it is at sea level. The field dimensions at Coors Field are:

Left Field - 347 feet (106 m)
Left-Center - 390 feet (119 m)
Center Field - 415 feet (126 m)
Right-Center - 375 feet (114 m)
Right Field - 350 feet (107 m)

 1. Overlay two plots: one at sea level and one in Denver to show why power hitters prefer Coors field.

2. Find the initial magnitude of velocity, v0

needed to hit a home run to Right-Center, where v_x(0)=v0/sqrt(2) and v_y(0)=v0/sqrt(2)

I don't quite understand how to use the field dimensions for both 1 and 2 and am pretty clueless as to how to approach this question using the ordinary differential equations mentioned above.

 

 

Hey guys. I've tried to find an answer for this, but have struggled since our learning book is in danish, so the used terms may not be technically correct, so sorry for poor phrasing.

 

Anyway, how would you solve this problem in maple?

 

Find a solution for the differential equation:

d4y/dt4 - 16y = u' + u

With the effect*  u(t) = e3t + 3et

 

 

I've gotten this result (By hand calculation)

y(t) = 4/65*e3t - 6/15*et

 

 

Thanks for the help. It's my first post, so let me know if I should do something different next time :)

-Alex

 

Bonus question:

How do you calculate the transfer function in maple:

H(s) = (s + 1)/(s4-16)

 

*Don't know if 'effect' is the correct term. 

I am attempting to plot an initial value problem in Maple 18.  I have my equation defined, as well as a general solution and two particular solutions at y(0)=3/4 and y(0)=1/2.  To graph, I entered the command

DEplot(de,y(x),x=-3..3,{[0,1/2],[0,3/4]},dirgrid=[12,12],color=black,linecolor=blue,thickness=2);

but instead of returning a graph, the software gave me the error message

Error, (in DEtools/DEplot/CheckDE) extra unknowns found: sinx

The Maple support site lists this as an unknown error, and as a new user, I'm not sure what to do.  What does this mean?

Good day everyone, how can one check for the congence of numerical solution of this ODE in maple? See it here 1.mw

Best regards.

Hello guys

I have a coupled linear differentional equation which are in the 4th order. they are shown in the below:

P:=phi(x):
Q:=psi(x):

eq1:=a11*diff(P,x,x,x,x)+a22*diff(P,x,x)+a33*P+a44*diff(Q,x,x)+a55*Q:
eq2:=a44*diff(P,x,x)+a55*P+a66*diff(Q,x,x)+a77*Q:

eq1:=0:
eq2:=0:

The boundary values for this coupled equation are:
phi(a)=sigma1,phi(-a)=sigma1,diff(P,x)(a)=0,diff(P,x)(-a)=0,psi(a)=sigma2,psi(-a)=sigma2

Now consider:

a11:=6.36463*10^(-10):
a22:=-1.22734*10^(-9):
a33:=3.48604*10^(-10):
a44:=2.94881*10^(-11):
a55:=-5.24135*10^(-11):
a66:=-1.03829*10^(-9):
a77:=4.86344*10^(-10):
when I use dsolve for deriving a good answer in this equation, there are six real roots .How can I solve it with these boundary condition?

I need to extract phi(x) and psi(x) from this coupled equation.

Thanks

 

Hi there,
I have the following set of equation which I want to solve using Maple's dsolve command:

d[V(t)*C(t)]/dt = G - K *C(t)
dV(t)/dt = alpha - beta

where V is the volume, C is the concentration, and t is the time variable. G, K, alpha and beta are constant parameters of the problem.


The solution for V(t) is easy to find operating the second equation:
V(t) = V_0 + (alpha - beta)*t

But solving for C(t) is a bit harder.

I would like Maple to solve the system, but the result I get does not really make sense to me. My attempt is this: Maple_artificial_kidney.mw

Any thoughts about how to introduce the equations successfully?

Thank you,
jon

Hallo. I'm solving a initial value problem for system of 7 ODE:

dsn := dsolve({expand(maineq[1, 1]), expand(maineq[1, 2]), expand(maineq[1, 3]), expand(maineq[1, 4]), expand(maineq[1, 5]), expand(maineq[1, 6]), expand(maineq[1, 7]), T(0) = .5, u(0) = u0, Y[1](0) = .8, Y[2](0) = .2, Y[3](0) = 0, Y[4](0) = 0, Y[5](0) = 0}, numeric, method = lsode[backfull])

 

Is there easy way how to plot result?

 

 

 

ds(t)/dt = a*s(t)*(1 - s(t) - m(t)) - b*s(t) 

dm(t)/dt = c*s(t) - d*m(t)

 

need to find steady state of this system ( finding this simultaneously) in maple 

 

How can you do it? 

hello , to solve a differential euqations , maple gave me this error. please help me.

 

TG.mw

hello dear freinds

im new comer in maple.

i want to find  particular solution of an ode by following code:

ode := diff(u[1](t), t, t)+u[1](t) = -(1/4)*a^3*cos(3*beta[0]+3*t)-(3/4)*a^3*cos(beta[0]+t)

m := combine(convert(particularsol(ode), trig))

but maple solution is : m := u[1](t) = (81/32)*a^3*cos(-3*beta[0]+t)-(81/16)*a^3*cos(3*beta[0]+t)-(3/8)*a^3*t*sin(beta[0]+t)+(3/16)*a^3*cos(-beta[0]+t)-(27/16)*a^3*cos(beta[0]+t)+(1/32)*a^3*cos(3*beta[0]+3*t)

but  particular solution is :

u[1](t) = -(3/8)*a^3*t*sin(beta[0]+t)+(1/32)*a^3*cos(3*beta[0]+3*t)

is there any idear for finding the solution?

thanks in advance

Hey, how is can i see all the steps in maple? I would specially like to know it for differential equations.

For example we could use this one:

dl := 3*(diff(y(t), t, t))+6*(diff(y(t), t))+4*y(t) = 0 

Hello,

 

could you help me solve this error ? I don't understand what it means.

 


> eq3:=diff(x(t),t,t)+Gamma*diff(x(t),t)+omega[0]^2*(x(t)-(diff(x(t),t,t)+Gamma*diff(x(t),t)+omega[0]^2*x(t)+omega[0]^2*X[0])/omega[0]^2) = -omega[0]^2*X[0]:
> dsolve(eq3);
Warning, it is required that the numerator of the given ODE depends on the highest derivative. Returning NULL.

 

Thanks.

Dear all,

I am trying to find the intial velocity of a ball that is shot under an angle while only the start and end coordinates are given. The air resistance should also be taken into account. 

In order to do that I have build the following Maple sheet:

Assignment_question_1.mw 

I have used two differential equations that both include the variables v0 and t, and then try to solve them. Only I receive an answer in the form of RootOf, which I cannot remove with for example allvalues. 

I have been working for quite a long time on this but I am not coming any further, so is there anyone who can find what I am doing wrong/what I should be doing else? Or maybe my whole approach is not right?

Even a small step in the right direction would be appreciated a lot!

Thanks in advance,

Elise

Please, how do i compute,solve and derive d' Alembert wave formula with maple since wave=How do i derive d'Alembert formula using maple and also how to solve any wave problems eg IC u(x,0)=1/2e^-x^2 and ut(x,0)-e^-x^2 when c=4

1 2 3 4 5 6 Page 1 of 6