please use matlab to code this program

Doublet and Uniform Flow (flow past a cylinder) Psi = Psi_d + Psi- = Psi_1 + Psi_2 = - A sin theta/r + Uy = - A sin theta/r + Ur sin theta Psi = U (r - A/Ur) sin theta Psi = Ur(1 - a^2/r^2)sin theta a = Squareroot A/U phi = phi_d + phi_u = phi_1 + phi_2 = a A cos theta/r - U x = - A cos theta/r = Ur cos theta phi = - U (r + A/Ur) cos theta = - Ur(1 + a^2/r^2) cos theta You have to plot stream lines for a flow past a cylinder for two different values