archstevej

147 Reputation

4 Badges

17 years, 67 days

MaplePrimes Activity


These are questions asked by archstevej

Hi all,

I have a system of differential equations:

Eq1:=diff(R(t), t) = -0.22*W(t)-3000+0.3*R(t)+3488*exp(0.295e-3*t)+.744*(diff(W(t), t))

Eq2:=diff(W(t), t) = 0.304*W(t)-0.1916e-1*R(t)-2981*exp(0.295e-3*t)+2180-.36*H(t)+.6*(diff(R(t), t))

To solve, I could plug Eq1 into Eq2. However, maple will also solve this system without plugging Eq1 into Eq2. But, I see the answers are different. So, which method is correct?

Hi all,

I am assuming this is a pretty easy thing to do. I am solving a system of differential equations, and then want to export my results to maple. I basically want the solution to the system for 40 time steps (counting by 5).

By hand this is

Hi all,

I am solving a system of ODEs numerically, and then plotting the results over time. I would really like Maple to show me the equation behind the plotted graphs. Through my searches I have not been able to find the code for that. Any help would be greatly appreciated. For example, how do I see the equation for W(t)?

sol1a := dsolve([Eq1a, Eq2a, Eq3a, H(0) = 4879, W(0) = 105000, R(0) = 0], numeric, method = rosenbrock, output = listprocedure)

plots[odeplot](sol1a, [t, W(t)], t = 0 .. 40, color = blue, style = point);

Hi all,

I am doing some research where I solve three differential equations simultaneously.. The results are three optimal time paths, W(t), R(t), and H(t). I want to plug these three into an integral, and solve (to get the net present value of my initial problem). I have no problem solving the diff. equations and getting values, given initial conditions. I am getting an error in the integration. I know I am missing something in the way I set it up. Anyways, I am putting my steps here below. Thanks for taking a look!

Hi All, 

I am working on a numerical simulation of a system of ODEs that blows up when I do a numerical solution. One of my advisors says to multiply  the ODEs in my system by -1 to reverse the role of the eigen values. I only vaguely understand this. Would you be able to point me in the right direction to figure this out? When I do this (-1 in front of two of the ODEs), my results are much more stable. I will attach the ODE system here. The parameter values follow.

 Thanks for any input,

Steve

 > Restart;

1 2 3 Page 2 of 3