Question: how to use dsolve for arbitrary parameter then plot3d

Hi experts in the field of solving complex system of ODE

I have this sytem of var. are  (n(t),u(t),z(t)) with delta and t are (real)

then want to plot 3d n(t) as a function of( t, delta)

restart:assume(delta,real):
phi:=0:lambda:=0.1:N:=5:M:=sqrt(N*(N+1))*exp(I*phi):omegap:=10:
var:={n(t),u(t),z(t)}:
dsys:={diff(n(t),t)=-2*(n(t)-N)+(u(t)-M)*exp(-2*I*omegap*t/lambda)+((z(t)-conjugate(M))*exp(2*I*omegap*t/lambda)),diff(u(t),t)=-2*(1-I*delta)*u(t)+2*(n(t)-N)*exp(2*I*omegap*t/lambda)+2*M,diff((z(t),t))=-2*(1+I*delta)*z(t)+2*(n(t)-N)*exp(-2*I*omegap*t/lambda)+2*conjugate(M)}:
res1:=dsolve(dsys union {n(0)=0,u(0)=0,z(0)=0},numeric,output=listprocedure);

P3:=plots3d(n(t),t=0..1,delta=-10..10,axes=boxed,tickmarks = [3, 2], color = black, thickness = 1, linestyle = solid, titlefont = [Helvetica, roman, 18], labeldirections = [horizontal, vertical], labelfont = [Helvetica, roman, 24]);