How to solve two differential equation simultaneously contain complex numbers?

how to solve these two differential equation to get the value of q,p:

dq/dz=-j*b1*q-j*k*p dp/dz=-j*b2*p-j*k*q where j notation of complex numbers (root -1) , b1 ,b2 and k constants

and q(0)=0,p(0)=1 and z=[0 10]

i try to solve these equations analytically and get the solution in the attached image but how to get the same solution using matlab so where is the error . where A,B,D are constants related to b1,b2,k