Naby

48 Reputation

2 Badges

14 years, 43 days

MaplePrimes Activity


These are questions asked by Naby

I am trying to solve the differential equation below using the following code. deq := diff(c(x), x, x) = c(x)*exp(gamma*beta*(1-c(x))/(1+beta*(1-c(x)))) ic := (D(c))(1) = 0, c(0) = lambda dsol := dsolve({deq, ic}, numeric, output=array([0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0]), continuation = lambda)
I have a system of two ODE's. One is of the form dT/dt=a*(dc/dt)-b, where a and b are constants. The other is of the form dc/dt=-(c^2)*d*exp(j*(1/T-1/k)), where d,j, and k are constants. I need to find the maximum value for T. I can solve dT/dt via the separation of variables method to get and expression T=f(c,t). I then substituted this expression into dc/dt=f(c(t),t) and solved it numerically in Maple. How do I substitute the values of c(t) obtained in Maple into T=f(c,t) and find the maximum? Or is there a better method?
Page 1 of 1