%% Function name
%       TandTinv

%% Revised:
%       21 January 2014

%% Author
%       Michael Bornick, Trey Moore, & Autar Kaw
%       Section: All
%       Semester: Fall 2012

%% Purpose
%       Given the angle of the ply in degrees output the transformation 
%       matrix and inverse transformation matrix

%% Usage
%       function [T,Tinv] = TandTinv(angle)
% Input Variable
%       angle=angle of ply in degrees
% Output Variables
%       [T]=transformation matrix
%       [Tinv]=inverse transformation matrix
% Keyword
%       transformation matrix
%       transformation of strains
%       transformation of stresses

%% License Agreement
%       http://www.eng.usf.edu/~kaw/OCW/composites/license/limiteduse.pdf

%% Code
function [T,Tinv] = TandTinv(angle)
% Sine of the angle of the lamina
s=sind(angle);
% Cosine of the angle of the lamina
c=cosd(angle);
% Transformation matrix
[T]=[c^2, s^2, 2*s*c; s^2, c^2, -2*s*c; -s*c, s*c, c^2-s^2;];
% Inverse transformation Matrix
[Tinv]=inv(T);
end