AbstractA method by which the numerical difficulties that one encounters while attempting to solve the differential equations which describe a diffusion flame can be circumvented is the quasi‐linearization technique. Application of this method in conjunction with difference approximations transforms the equations into tri‐diagonal matrices which can be inverted by the Thomas method. Since this is an iterative method initial profiles are generated by employing Burke‐Schumann kinetics. This method is found to be superior to the marching techiques that are frequently used in the analyses of diffusion f