%Script to illustrate use of subplot to present results
clear all; clc; close all;

%Generate junk values of shear force etc.
x = linspace(0,2,21);
sforce = x; 
bm = x.^2;
torque = x.^3;

figure(1); 
clf;
suptitle('Main Shaft Loads');
subplot(3,1,1)
plot(x, sforce, '-k');
xlabel('x (m)');  ylabel('Shear (kN)');

subplot(3,1,2)
plot(x, bm, '-r');
xlabel('x (m)');  ylabel('Moment (kN*m)');

subplot(3,1,3)
plot(x, torque, '-b');
xlabel('x (m)');  ylabel('Torque (kN*m)');