%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% doit.m
% This will do lots of stuff.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fprintf(1,'\n\nStarting...\n');

clear;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 0 Example
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%fprintf(1,'\n\nExample...\n');
%n = 400; h = 10 / n;
%y = [1 0 0.001 0];
%tacoma([0 400], y, h, 3)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 1 Zero Initial Angle
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fprintf(1,'\n\nZero Initial Angle...\n');

%n = 400; h = 10 / n;
%y = [0 0 0 0];
%tacoma([0 200], y, h, 3)

fprintf(1,'\n\nResult: Nothing happens to the default system.\n');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 2a Replace trapstep With rk4step
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fprintf(1,'\n\nRk4step...\n');

%n = 400; h = 10 / n;
%y = [1 0 0.001 0];
%rk4tacoma([0 200], y, h, 3, 11, 3, 0)

fprintf(1,'\n\nResult: It still works using rk4step.\n');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 2b Find a Range for "a" (F)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fprintf(1,'\n\nFinding a Range for "a"...\n');

%n = 400; h = 10 / n;
%y = [1 0 0.01 0];
%F = 9.45;
%W = 3;
%while(F <= 9.45)
%    F
%    rk4tacoma([0 1000], y, h, 3, F, W, 1)
%    F = F + 0.5;
%end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 2c Angle When Not Excited
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fprintf(1,'\n\nAngle When Not Excited...\n');

fprintf(1,'\n\nWhen not excited the angle goes to zero.\n');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 3 Plot (A, w)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fprintf(1,'\n\nPlot(A, w)...\n');

%n = 400; h = 10 / n;
%y = [1 0 0.01 0];
%F = 15.6;
%W = 2;
%while(F <= 15.6)
%    F
%    rk4tacoma([0 1000], y, h, 3, F, W, 1)
%    F = F + 0.5;
%end

%n = 400; h = 10 / n;
%y = [1 0 0.01 0];
%F = 38;
%W = 4;
%while(F <= 99)
%    F
%    rk4tacoma([0 2000], y, h, 3, F, W, 1)
%    F = F + 0.1;
%end

%n = 400; h = 10 / n;
%y = [1 0 0.01 0];
%F = 83.8;
%W = 5;
%while(F <= 83.8)
%    F
%    rk4tacoma([0 500], y, h, 3, F, W, 1)
%    F = F + 0.1;
%end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 4 Increasing Damping
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fprintf(1,'\n\nIncreasing Damping...\n');

fprintf(1,'\n\nFinished!\n');