Two Polarizing Plates

%% Two Polarizing Plates
% Determines the intensity of a wave passing through polarizer A and B ans
% then interchanged

% Inputs and Variables
I0 = input('The average intensity of the polarized wave\n');
theta1 = input('Angle of polarizer A\n')*pi/180;
theta2 = input('Angle of polarizer B\n')*pi/180;

% Outputs and Equations
I1 = I0*(cos(theta1))^2
I2 = I1*(cos(theta2-theta1))^2
I1prime = I0*(cos(theta2))^2;
If = I1prime*(cos(theta1-theta2))^2
fprintf(['The average intensity after...\npolarizer A: %.4e W/m^2\n'...
         'polarizer B: %.4e W/m^2\n'...
         'the interchange: %.4e W/m^2\n'], I1, I2, If)

%{
The trigonomic functions are in radians
%}