% 30-second demo of survey.m, showing first 10 modes.
% CCK 2016-Mar-10

k1 = 1
k2 = 5
L  = 2
safety = 40
outfile = 'surveytest.mat'

figure(1)
start_time = datestr(now)
[tension, k] = survey(k1, k2, L, safety, outfile);
finish_time = datestr(now)

figure(2)
hold off
plot(k,log10(tension),'r')

