CODE:
clear; clc;
text= load('V.dat');
x=text(:,1);
r=text(:,2);
vx=text(:,3);
dvx=text(:,5);
dvr=text(:,6);
v11=text(:,8);
v21=text(:,9);
%%
Fx = scatteredInterpolant(x(,r(,dvx();
Fr = scatteredInterpolant(x(,r(,dvr();
%%
xx=linspace(min(x(),max(x(),500);
rr=linspace(min(r(),max(r(),500);
[xgg,rgg]=meshgrid(xx,rr);
xstream = Fx(xgg,rgg);
ystream = Fr(xgg,rgg);
%%
scrsz = get(0,'ScreenSize');
figure1 = figure('Position',[0.05*scrsz(3) 0.05*scrsz(4) 0.9*scrsz(3) 0.9*scrsz(4)]);
[xs,rs] = meshgrid(x,r);
[dvxs,dvrs] = meshgrid(dvx,dvr);
quiver(x,r,dvx,dvr,'r');
numstream=400;
strx=randi([2,99],numstream,1);
stry=randi([0,12],numstream,1);
strx=[strx,strx];
stry=[stry,-stry];
h=streamline(xgg,rgg,xstream,ystream,strx,stry);
set(h,'LineWidth',0.5,'Color','k')
axis equal
axis tight
box on