Mokusko

10 Reputation

2 Badges

10 years, 58 days

MaplePrimes Activity


These are replies submitted by Mokusko

@Carl Love 

thanks a lot for your answer, it's really helpful for me

@tomleslie 

Thank for your advice. I'm agreeing with all what you had wrote in your post. The second bullet is all true, and I see this problem, but I don't know how to solve it.

@ NULL

restart

alpha := .3;

.3

 

3.5

 

6

 

1.1

(1)

for i to 50 do IC := x(z)+y[i](z); ICon := x(0) = 1, y[i](0) = 0; for j to 50 do if i <> j then d1 := diff(x(z), z) = -G*x(z)*y[i](z)/IC-alpha*x(z); d2 := diff(y[i](z), z) = G*y[i](z)*y[j](z)/IC-alpha*y[i](z); dsys := {d1, d2}; F := dsolve({ICon, op(dsys)}, [x(z), y[i](z)], numeric) end if end do end do

Error, (in dsolve/numeric/process_input) unknown y[2] present in ODE system is not a specified dependent variable or evaluatable procedure

 

``


But I don't need to solve y[j](z) from this equations. I just need to find x(z), y[i](z), it is in dsolve. y[j](z) is helping to take into consideration an interaction between y[i](z) with different indecies.

If finding some analytic solution is impossible, then I need just to plot this solution like a graph.

Download 05.mw

Page 1 of 1