The paper is devoted to the results of the EAS neutron component investigations by means of the PRISMA-32 array. The array consists of 32 en-detectors and enables to record delayed thermal neutrons accompanying showers. For registration of thermal neutrons, the scintillator based on $^{6}Li$ isotope as a target is used in the detectors. Some results of the processing of data accumulated over a long period of time are presented: the lateral distribution function of neutrons in EAS and preliminary results on EAS neutron multiplicity spectrum and distribution of showers in e/n ratio.

We report on a measurement of thermal neutrons, generated by the hadronic component of extensive air showers (EAS), by means of a small array of EN-detectors developed for the PRISMA project (PRImary Spectrum Measurement Array), novel devices based on a compound alloy of ZnS(Ag) and $^{6}$LiF. This array has been operated within the ARGO-YBJ experiment at the high altitude Cosmic Ray Observatory in Yangbajing (Tibet, 4300 m a.s.l.). Due to the tight correlation between the air shower hadrons and thermal neutrons, this technique can be envisaged as a simple way to estimate the number of high energy hadrons in EAS. Coincident events generated by primary cosmic rays of energies greater than 100 TeV have been selected and analyzed. The EN-detectors have been used to record simultaneously thermal neutrons and the air shower electromagnetic component. The density distributions of both components and the total number of thermal neutrons have been measured. The correlation of these data with the measurements carried out by ARGO-YBJ confirms the excellent performance of the EN-detector.

A novel type of EAS array (PRISMA-32) has been constructed on the base of NEVOD-DECOR experiment (MEPhI,Moscow) and is now taking data. It consists of 32 specially designed scintillator en-detectors able to measure two main EAS components: hadrons (n) and electrons (e). First results on thermal neutron lateral as well as temporal distributions are presented. Obtained exponential neutron lateral distributions are consistent with that expected for normal hadron production with exponential transverse momentum distribution. As there are no other experimental data on thermal neutron distributions and so, to compare results with other measurements, we additionally obtained electron lateral distribution function (using the same detectors) and compared it with NKG - function. Recorded neutron temporal distributions are very close to that obtained with data of our previous prototypes.

Preliminary results of one year anisotropy measurement in the energy range 10^{13} -10^{14} eV as a function of energy are presented. The results are compared for two methods of data analysis: the standard one with meteo correction approach in use and another one so-called "East minus West" method. Amplitudes and phases of anisotropy for three median energies E = 25 TeV, E = 75 TeV and E = 120 TeV are reported. Brief consideration of amplitude-phase dependence of anisotropy on energy is expounded.

When using a recently developed method of Doppler-Zeeman mapping (Vasilchenko et al., 1996) for analysis of a real star and real observational data, we are confronted with limitations due to the model simplifications and unavoidable errors in observed spectra. We discuss the errors introduced by probable inaccurracies of the mathematical model: analytical fit of the local Stokes parameters, influence of magneto-optical effect, ignorance of the true atmosphere model to compute local Stokes profiles, non-uniform surface brightness. The magnetic field configuration is found in the form of arbitrarily shifted dipole and sum of dipole and quadrupole, along with the distribution of Si, Ti, Cr and Fe over the surface of the star. Lines of different elements lead to the same magnetic field configuration, which is reliably determined for the part of the stellar surface which faces the observer. This allows to compare the magnetic field and chemical maps of the surface of HD 215441. A large-scale ring structure with the magnetic pole at its center is clearly seen on the abundance maps. Si, Cr and Ti are highly deficient where the magnetic field lines are vertical (near the magnetic pole) while Fe is highly overabundant there.

An efficient method has been developed for solving the inverse problem of Doppler-Zeeman mapping of magnetic, chemically peculiar stars. A regularized iteration method is used to simultaneously solve the integral equations for the Stokes I, V, Q and U parameters. The validity of analytical fits to the local profiles of the Stokes parameters is substantiated. The algorithm had been tested on models and makes it possible to obtain simultaneously, from the observed Stokes profiles, a map of the distribution of a chemical element and the parameters of an arbitrarily shifted magnetic dipole.