I am having some trouble with solving a somewhat heavy differential equation system, which is consisted of 10 variables and other abstract parameters as follows:
A[t] == a0 + a1 * L[t]
B[t] == c0 + c1 * L[t]
M[t] == b0 + b1 * H[t]
Q[t] == (A[t] - 1) * B[t]
R[t] == x * q * A[t] * B[t] + f * M[t]
F'[t] == d * (Q[t] - R[t])
G[t] == v * F[t]
H[t] == m0 + m1 * L[t] - m2 * G[t]
L[t] == p0 + p1 * J[t] - p2 * F[t]
J'[t] == n * (y - L[t])
I used DSolve as follows:
DSolve[{A[t] == a0 + a1 * L[t], B[t] == c0 + c1 * L[t], M[t] == b0 + b1 * H[t], Q[t] == (A[t] - 1) * B[t],
R[t] == x * q * A[t] * B[t] + f * M[t], F'[t] == d * (Q[t] - R[t]), G[t] == v * F[t], H[t] == m0 + m1 * L[t] - m2 * G[t],
L[t] == p0 + p1 * J[t] - p2 * F[t], J'[t] == n * (y - L[t])}, {A[t], B[t], M[t], Q[t], R[t], F[t], G[t], H[t], L[t], J[t]}, t]
Mathematica can't solve this and just return the same code. Am I missing something?
