close all
clear all

g = 9.81;
rho = 1025;  %kg/m^3
d = 5;  %Sphere diameter
V = pi*d^3/12;  %m^3
m = rho*V;

f = .05:.01:.55;
w = 2*pi*f;
k = w.^2/g;
kd = k*d;

[a,b] = sphere_heave_radiation_coeffs(d,w,g,rho);



plot(w.^2*d/g,a/(rho*V),w.^2*d/g,b./(rho*V*w));