Please show the stream function plots or provide it in MATLAB, thank you so much!

The flow field has the following polar coordinate velocities:

$V_r=(1-\frac{R^2}{r^2})V_{}cos()$

$V_{}=-(1+\frac{R^2}{r^2})V_{}sin()$

given that the stream functions is defined as follows in polar coordinates:

$V_r=\frac{1}{r}\frac{\text{d}\text{e}\text{l}}{\text{d}\text{e}\text{l}}$

$V_{}=-\frac{\text{d}\text{e}\text{l}}{\text{d}\text{e}\text{l}\text{r}}$

a$)$ plot the streamlines for this flow for values of $r=0\mathrm{.}5R,R,2R$

b$)$ calculate the velocity along the streamline $r=R$

c$)$ calculate the pressure coefficient along the streamline $r=R$

d$)$ prove whether the flow is rotational or irrotational.