(* irkdg.wl -- Unfinished Matrix Derivative Computation *)

n=3
Y=Array[y,n]
F=Function[{t,z},{z[[2]]*z[[3]]*Sin[t]-z[[1]]*z[[2]]*z[[3]],
  -z[[1]]*z[[3]]*Sin[t]+(1/20)*z[[1]]*z[[3]],
  z[[1]]^2*z[[2]]-(1/20)*z[[1]]*z[[2]]}]
K1=Array[k,n]
K2=Array[k,n,n+1]

(* Please check next week for the final version *)
