- /
-
Developing Flow
on 10 Nov 2023
- 4
- 9
- 0
- 0
- 281
drawframe(1);
Write your drawframe function below
function drawframe(f)
u = 30;
v = 0;
clf
hold on
xlim([-5 5])
ylim([-5 5])
u = f;
Q = 40;
for psi = 1:20:200
xVals = [];
yVals = [];
for theta = 0:0.001:pi
y = (psi - Q./(2.*pi).*theta)./u;
r = y./sin(theta);
x = r.*cos(theta);
xVals = [xVals,x];
yVals = [yVals,y];
end
plot(xVals,yVals,'Color','b')
plot(xVals,-yVals,'Color','b')
end
end