%Stokes Law - Rise speed of bubbles
%Author: Tom D Pering - University of Sheffield
%You are free to use and alter with acknowledgement
function [u] = stokes (rho, rs, dSlug, Eta)
%rho = denisty of the magma (or liquid)
%rs = radius of the bubble
%dSlug = density of the slug
%Eta = the viscosity of the magma
g=9.81;
u=(2*(rho-dSlug)*g*(rs.^2))/(9*Eta); %Eta is magma viscosity
clear g
end