I have two coupled dimensionless differential equations, how can I run them to obtain a plot of amplitude against time

x1=y=yo*cos(w*t);

x2=dy/dt;

x3=q=qo*cos(w*t+phi);

x4=dq/dt;

x5=d^2y/dt^2=-w^2*yo*cos(w*t)

x6=d^2q/dt^2;

x5=-(2*zeta*delta+(gamma/mu))*x2-delta^2*x1+b*x3

x6=-e(x3^2-1)*x4-x3+A(x5)