La solución de sistemas lineales es similar a la de una sola de las ecuaciones.
En primer lugar vamos a bajar el grado de las ecuaciones mediante el establecimiento $u:=y',v:=z'$ y volver a escribir el sistema en la forma canónica:
$$u'(t)-pv(t)=0\\v'(t)+pu(t)=pq.$$
Encontrar la solución general del sistema homogéneo
$$u'(t)-pv(t)=0\\v'(t)+pu(t)=0.$$
Usted puede hacer que al intentar una solución exponencial como $u=u_0e^{\lambda t},v=v_0e^{\lambda t}$. Usted obtener:
$$\lambda u_0e^{\lambda t}-pv_0e^{\lambda t}=0\\\lambda v_0e^{\lambda t}+pu_0e^{\lambda t}=0,$$
es decir, después de la simplificación
$$\lambda u_0-pv_0=0\\\lambda v_0+pu_0=0.$$
Usted debe reconocer un Eigenproblem para un $2\times2$ matriz (dos elementos que faltan). Los Valores propios son $\pm ip$, y los correspondientes vectores propios se $(1,i), (1,-i)$.
La combinación, la solución general puede escribirse
$$u(t)=C_0\cos(pt)+C_1\sin(pt)\\v(t)=C_1\cos(pt)-C_0\sin(pt).$$
Ahora mirando el no-sistema homogéneo, debemos encontrar una constante de la solución particular:
$$u(t)=u_0\\v(t)=v_0.$$
Conectar el sistema obtenemos
$$-pv_0=0\\pu_0=pq,$$
y por la identificación,
$$u_0=q,v_0=0.$$
Ahora, para obtener el $y$$z$, combinar e integrar una vez en $t$.
Este es un resumen rápido de un enfoque general para resolver este tipo de sistemas, que pueden estar basados en el álgebra matricial técnicas.