# Plot a stream function?

GROUPS:
 Hi, I have problem to plot this equation \[Psi](y,z)=( k[y]' z^2 (k[y]-z))/(8k[y] ); where k[y]'=Sqrt[(1-k[y]/p)(1-k[y] p Cos[\[Alpha]]- k[y] Log[k[y]/p])] and I am expected to get this kind of graphI get k[y] from this f[k_?NumericQ]:=p*NIntegrate[1/Sqrt[(1-t)(1-t p^2 Cos[\[Alpha]])-3 M p t Log[t]],{t,k/p,1}] ky=FindRoot[y==f[k],{k,0}]; Table[{y,ky[[1,2]]},{y,0,a1+\[Delta],\[Delta]}]//N where p=1.57516;a1=1.7824; Can anyone help me with this problem? Thanks in advance!