Use quiver () function to plot vector field (xy plane) of O2 flux in lab5.

What I have currently:

rt= 50; % in micrometers

mu= 3.75*10.^8; %in mPa

delta_p= 1000; %delta p mmHg/um

l= 200; %length of the vessel (0-200 um)

v= ((rt.^2)/4.*u).*delta_p.*(1-(r/rt).^2); %velocity of vessel

x=linspace(0,200,20); % x points for the meshgrid

y=linspace(-rt,rt,-200); %y point for the meshgrid

%graph selection

[X, Y]=meshgrid (x,y);

vrx=blood(Y,rt,mu,delta_p);

U=Vrx;

V=U;

V(:,:)=0;

%vrx=v;

figure (1)

quiver(X,Y,u,v);

%part b graphing

%velocity of blood function

function [Vrx]=blood(r,rt,mu,delta_p)

Vrx=((rt^2)./(4*mu))*delta_p*(1-r/rt.^2);

end

What I had for Lab 5:

function[j]=lab05(x0,y0)

syms x y;

%muscle cell radius (R) - um - 1*10^-6m

R=25*(10^(-6));

%diffusion coeff(D) - cm^2*s^-1 - (1*10^-2m)^2*s^-1

D=(10^(-5))*(10^(-4));

%oxygen concentration outside cell (co2)0 - mol O2*mL^-1

co2=2.22*(10^(-7));

%oxygen consumption rate (sigma)- mole*mL^-1*s^-1

sigma=1.65*10^-8;

%oxygen concentration equation(phi)

phi=((sigma*(R^2))/4*D)*((((x^2)+(y^2))/(R^2))-1)+co2;

rtest= (x0^2)+(y0^2);

if rtest < (R)^2

disp('Will Work')

else

rtest > (R)^2; disp('Wont Work')

end

gradphi=gradient(phi,[x,y]);

%diffusing oxygen flux equation(j)

j=-D*gradphi;

jr0=subs(j,{x,y},{x0,y0});

divergence(j,[x,y])-phi;

disp('Flux of diffusing oxygen:') end