1. A method of designing an integrated computational imaging system for producing images having an extended depth of field comprising steps of: formulating an imaging system having an aberration component dominated by a third order spherical aberration for producing a blurred intermediate image; and determining an amount of the third order spherical aberration that reduces variations among impulse responses over the extended depth of field; calculating an impulse response that departs from an arithmetic average of the impulse responses over the extended depth of field to unevenly weight the impulse responses over the extended depth of field; and incorporating the calculated impulse response in a digital processing algorithm for deblurring the intermediate image produced by the multifocal imaging subsystem and for calculating a recovered image having an extended depth of field.