Error using / Matrix dimensions must agree for Vr

clc

clear all

close all

Vo= 10;

R= 1;

x = linspace(-4,4,1000);

y = linspace(-4,4,1000);

r = sqrt(x.^2 +y.^2);

theta = atan(y/x);

Vr = 1-((R^2)/(r.^2))*Vo*cos(theta);

% Vt = -(1- (R^2/r.^2)*Vo*sin(theta));

% Vx = (Vr*cos(theta) - Vt*sin(theta));

% Vy = (Vr*sin(theta) - Vt*cos(theta));

disp(Vr);

Answer 1

[Star Strider]

Vectorise!

Vr = 1-((R^2)./(r.^2)).*Vo.*cos(theta);
             ↑        ↑   ↑  ←  HERE

That should do what you want.

See Array vs. Matrix Operations for details.