t = 0:.01:1;
l_value = l(1);
v1 = velocity( r, angular_velocity, t, l_value );
max_flow1 = max(v)*piston_area;
disp(['For l = ' num2str(l_value) ' the peak flow rate is ' num2str(max_flow1) ' in^3/sec' ]);
l_value = l(2);
v2 = velocity( r, angular_velocity, t, l_value );
max_flow1 = max(v)*piston_area;
disp(['For l = ' num2str(l_value) ' the peak flow rate is ' num2str(max_flow1) ' in^3/sec' ]);
l_value = l(3);
v3 = velocity( r, angular_velocity, t, l_value );
max_flow1 = max(v)*piston_area;
disp(['For l = ' num2str(l_value) ' the peak flow rate is ' num2str(max_flow1) ' in^3/sec' ]);
l_value = l(4);
v4 = velocity( r, angular_velocity, t, l_value );
max_flow1 = max(v)*piston_area;
disp(['For l = ' num2str(l_value) ' the peak flow rate is ' num2str(max_flow1) ' in^3/sec' ]);
l_value = l(5);
v5 = velocity( r, angular_velocity, t, l_value );
max_flow1 = max(v)*piston_area;
disp(['For l = ' num2str(l_value) ' the peak flow rate is ' num2str(max_flow1) ' in^3/sec' ]);
plot(t,v1*piston_area, 'r:' );
plot(t,v2*piston_area, 'g:');
plot(t,v3*piston_area, 'g:');
plot(t,v4*piston_area, 'g:');
plot(t,v5*piston_area, 'g:');