Matlab diary
rkf.m
comprkf.m
plot of solution