clear all; close all;
%-----User input---------------
ID = 1.5; %cylinder ID in inches
D = .92; %distance in inches
P = 900; %pressure in psi
fpin = 'Gysin_900psi.txt'
%--------------------------------

%calc work
Cyl_area = pi() * (ID/2)^2; %in^2
F = P * Cyl_area; %pounds
W = F * D/12; %ft-lbs
W = 1.35582 * W; %Joules
fprintf('Mechanical work = %f J\n',W);



%calc regen J, Joules = Volts * Amps * Seconds
fid = fopen(fpin,'r');

tline = fgetl(fid);
A = sscanf(tline,'%s %s %f %f');
E = A
while ischar(tline)
    %disp(tline)
    tline = fgetl(fid);
end
fclose(fid);