Analytical Pipe flow matlab solution plotter (Bessel)

I need code written in matlab to plot pipe flow solutions (hagan poiseuille, startup flow, oscillating pressure gradient) . I have the equations you will need. I will provide the parameters, and you will need to implement them. I need this completed very quickly and efficiently. Once all is said and done, the plots should be time history plots that can be overlaid with experimental results.