% testCentripetal;.m   test centripetal.m
% Stephanie Countess for CISC106, 10 am, 10/01/2007

% Test formula to calculate the centripetal force of an object around a
% circle

%%%%%%%%%%%%%%%%%%%%%%%%
% define the tolerance %
%%%%%%%%%%%%%%%%%%%%%%%%

% tolerance is "how far away from the expected value" the result
% is allowed to be.  Since the calculation involves a square root,
% not all inputs will produce results that can be represented exactly.

tolerance = 0.001;

%%%%%%%%%%%%%%%%%
% Run the tests %
%%%%%%%%%%%%%%%%%

% Test 1


actual = centripetal(2,5,10);
expected = 5;

diff = abs(expected-actual);

if (diff < tolerance)
  disp('test 1 passed');
else
  disp('test 1 failed');
  expected
  actual
end
  

% Test 2

actual = centripetal(1,2,1);
expected = 4;

diff = abs(expected-actual);

if (diff < tolerance)
  disp('test 2 passed');
else
  % signal that test failed, and print both expected and actual
  disp('test 2 failed');
  expected
  actual
end
  

% end testCentripetal.m