The difference to a linear system of diffs eqs. is that two state variables are multiplied. Therefore, I think that this system falls under the nonlinear systems for which no general solution can be found analytically. Is this correct?
I could not find information on the net on systems of DEs of such a form. Can somebody send me some links about this kind of system with a multiplication of two state vars?

There probably isn't an analytic solution for the system, but there is some analysis you can do for this. One thing you can do is find the fixed points of the system: points (x*,y*,z*) such that (x',y',z') = (0,0,0) at these points.

0=a1*x*z
0=b1*x + b2*y
0=c1 + c2*z + c3*y*z

From eq 3, we see that we can't have z* = 0 (unless c1 = 0). The first equation then means that x* = 0, which in the second equation demands y* = 0. Then, the third equation needs z* = -c1/c2. So, if your system starts at the point (x,y,z) = (0,0,-c1/c2), it will stay there for all time. Then, you can consider the stability of this point - what happens when you perturb the system away from it? To do this, we linearize the system about the fixed point. We can write a linear system in the form

[tex] \mathbf{\dot{x}} = \mathsf{A}\mathbf{x}[/tex]

where [tex]\mathbf{x} = (x-x^\ast,y-y^\ast,z-z^\ast)^\mathcal{T}[/tex] and

Once you have this matrix, you find its eigenvalues and eigenvectors. The eigenvalues will tell you the stability of the point, and the eigenvectors will tell you along which direction trajectories will flow. For example, for a general eigenvalue [itex]\lambda = s + i \omega[/itex] with an eigenvector along the x direction (which may not be the case for this system), for s > 0 the system is unstable in the x direction (trajectories leave the fixed point), whereas for s < 0 trajectories will flow toward the fixed point in the x direction. For s = 0, the linear theory predicts closed orbits about the fixed point, but this could turn out to be incorrect due to our neglect of higher order terms.

There are many possible types of behaviours that could occur, especially with a 3D system. For example, the system could become chaotic, like the Lorenz system. The matter is further complicated by the fact that your parameters here are all arbitrary. The stability of the fixed point will generally depend on the values your parameters take.

For a good book about analysing dynamical systems like these, see "Nonlinear Dynamics and Chaos" by Steven Strogatz.