Question: Plotting functions which have been numeric solved

Hello,

I have a little problem. I have a system of ODEs, I can solve (with numeric solve) it but I can't plot it...
Any idea to help ?

Thank you !!

Here is my code:

with(DEtools);
beta := 1; lambda := 5; nu := 1; Delta := 1;


sys := {C(0) = 0, In(0) = .2, diff(C(u), u) = beta*In(u)*s(u)-C(u)/nu, diff(In(u), u) = C(u)/nu-In(u)/lambda+Delta*(diff(diff(In(u), u), u)), diff(s(u), u) = -beta*In(u)*s(u), s(0) = .8, (D(In))(0) = .2}

dsys := dsolve(sys, numeric, [In(u), s(u), C(u)])

plot(dsys[1], u = 0 .. .5, axes = boxed)

And I have this error:

Warning, expecting only range variable u in expression dsys[1] to be plotted but found name dsys[1]


PS: Maple solves the system without any problem:

 dsys(.5);

        [                                              
        [u = 0.5, In(u) = HFloat(0.33305297276372425),
        [                                              

           d                                       
          --- In(u) = HFloat(0.33950703755414946),
           du                                      

          s(u) = HFloat(0.7022184797811994),

                                           ]
          C(u) = HFloat(0.0781621551198558)]
                                           ]