%% regression for .0031'' orifice
fitEqn = 'exp(-x/tau)';
fitOps = fitoptions('Method','NonlinearLeastSquares',...
    'Lower',[0],...
    'Upper',[inf],...
    'StartPoint',[1]);
fitType = fittype(fitEqn,'options',fitOps)
[fitInfo0031,GoF0031] = fit(raw0031(:,1),raw0031(:,2)/raw0031(1,2),fitType)
figure(7);
plot(fitInfo0031,raw0031(:,1),raw0031(:,2)/raw0031(1,2));
%% regression for .0062'' orifice
fitEqn = 'exp(-x/tau)';
fitOps = fitoptions('Method','NonlinearLeastSquares',...
    'Lower',[0],...
    'Upper',[inf],...
    'StartPoint',[1]);
fitType = fittype(fitEqn,'options',fitOps)
[fitInfo0062,GoF0062] = fit(raw0062(:,1),raw0062(:,2)/raw0062(1,2),fitType)
figure(7);
plot(fitInfo0062,raw0062(:,1),raw0062(:,2)/raw0062(1,2));
%% regression for .0128'' orifice
fitEqn = 'exp(-x/tau)';
fitOps = fitoptions('Method','NonlinearLeastSquares',...
    'Lower',[0],...
    'Upper',[inf],...
    'StartPoint',[1]);
fitType = fittype(fitEqn,'options',fitOps)
[fitInfo0128,GoF0128] = fit(raw0128(:,1),raw0128(:,2)/raw0128(1,2),fitType)
figure(7);
plot(fitInfo0128,raw0128(:,1),raw0128(:,2)/raw0128(1,2));
%% regression for .0240'' orifice
fitEqn = 'exp(-x/tau)';
fitOps = fitoptions('Method','NonlinearLeastSquares',...
    'Lower',[0],...
    'Upper',[inf],...
    'StartPoint',[1]);
fitType = fittype(fitEqn,'options',fitOps)
[fitInfo024,GoF024] = fit(raw024(:,1),raw024(:,2)/raw024(1,2),fitType)
figure(7);
plot(fitInfo024,raw024(:,1),raw024(:,2)/raw024(1,2));