The innermost stable circular orbit equation of a test particle is obtained for an approximate Kerr-like spacetime with quadrupole moment. We derived the effective potential for the radial coordinate by the Euler-Lagrange method. This equation can be employed to measure the mass quadrupole by observational means, because from this equation a quadratic polynomial for the quadrupole moment can be found. As expected, the limiting cases of this equation are found to be the known cases of Kerr and Schwarzschild.

A new rotation version of the Curzon-Chazy metric is found. This new metric was obtained by means of a perturbation method, in order to include slow rotation. The solution is then proved to fulfill the Einstein field equations using a REDUCE program. Furthermore, the applications of this new solution are discussed.

A metric representing a slowly rotating object with quadrupole moment is obtained using a perturbation method to include rotation into the weak limit of the Erez-Rosen metric. This metric is intended to tackle relativistic astrometry and gravitational lensing problems in which a quadrupole moment has to be taken into account.

A new Kerr-like metric with quadrupole moment is obtained by means of perturbing the Kerr spacetime. The form of this new metric is simple as the Kerr metric. By comparison with the exterior Hartle-Thorne metric, it is shown that it could be matched to an interior solution. This approximate metric may represent the spacetime of a real astrophysical object with any Kerr rotation parameter a and slightly deformed.

A program to generate codes in Fortran and C of the full Magnetohydrodynamic equations is shown. The program used the free computer algebra system software REDUCE. This software has a package called EXCALC, which is an exterior calculus program. The advantage of this program is that it can be modified to include another complex metric or spacetime. The output of this program is modified by means of a LINUX script which creates a new REDUCE program to manipulate the MHD equations to obtain a code that can be used as a seed for a MHD code for numerical applications. As an example, we present part of output of our programs for Cartesian coordinates and how to do the discretization.

Computer simulations of plasmas are relevant nowadays, because it helps us understand physical processes taking place in the sun and other stellar objects. We developed a program called PCell which is intended for displaying the evolution of the magnetic field in a 2D convective plasma cell with perfect conducting walls for different stationary plasma velocity fields. Applications of this program are presented. This software works interactively with the mouse and the users can create their own movies in MPEG format. The programs were written in Fortran and C. There are two versions of the program (GNUPLOT and OpenGL). GNUPLOT and OpenGL are used to display the simulation.

Gravitational lenses are presently playing an important role in astrophysics. By means of these lenses the parameters of the deflector such as its mass, ellipticity, etc. and Hubble's constant can be determined. Using C, Xforms, Mesa and Imlib a computer program to visualize this lens effect has been developed. This program has been applied to generate sequences of images of a source object and its corresponding images. It has also been used to visually test different models of gravitational lenses.

A REDUCE code for the Newman-Janis algorithm is described. This algorithm is intended to include rotation into nonrotating solutions of the Einstein field equations with spherically symmetry or perturbed spherically symmetry and has been successfully applied to many spacetimes. The applicability of the code is restricted to metrics containing potentials of the form 1/r.

A new Kerr-like metric with quadrupole moment is obtained by means of perturbing the Kerr spacetime. By comparison with the exterior Hartle-Thorne metric, it is showed that it could be matched to an interior solution. This metric may represent the spacetime of an astrophysical object.

N-body simulations predict that dark matter haloes are described by specific density profiles on both galactic- and cluster-sized scales. Weak gravitational lensing through the measurements of their first and second order properties, shear and flexion, is a powerful observational tool for investigating the true shape of these profiles. One of the three-parameter density profiles recently favoured in the description of dark matter haloes is the Einasto profile. We present exact expressions for the shear and the first and second flexions of Einasto dark matter haloes derived using a Mellin-transform formalism in terms of the Fox H and Meijer G functions, that are valid for general values of the Einasto index. The resulting expressions can be written as series expansions that permit us to investigate the asymptotic behaviour of these quantities. Moreover, we compare the shear and flexion of the Einasto profile with those of different mass profiles including the singular isothermal sphere, the Navarro-Frenk-White profile, and the S\'ersic profile. We investigate the concentration and index dependences of the Einasto profile, finding that the shear and second flexion could be used to determine the halo concentration, whilst for the Einasto index the shear and first and second flexions may be employed. We also provide simplified expressions for the weak lensing properties and other lensing quantities in terms of the generalized hypergeometric function.

We present a simple model for the non-thermal emission from the historical supernova remnant SN 1006. We constrain the synchrotron parameters of the model with archival radio and hard X-ray data. Our stationary emission model includes two populations of electrons, which is justified by multi-frequency images of the object. From the set of parameters that predict the correct synchrotron flux we select those which are able to account, either partly or entirely, for the gamma-ray emission of the source as seen by HESS. We use the results from this model as well as the latest constraints imposed by the Fermi observatory and conclude that the TeV emission cannot be accounted for by neutral pion decay produced by high-energy cosmic rays with a single "soft" power-law distribution (i.e., with a particle index greater than 2 or so).

Recent high-resolution N-body CDM simulations indicate that nonsingular three-parameter models such as the Einasto profile perform better than the singular two-parameter models, e.g. the Navarro, Frenk and White, in fitting a wide range of dark matter haloes. While many of the basic properties of the Einasto profile have been discussed in previous studies, a number of analytical properties are still not investigated. In particular, a general analytical formula for the surface density, an important quantity that defines the lensing properties of a dark matter halo, is still lacking to date. To this aim, we used a Mellin integral transform formalism to derive a closed expression for the Einasto surface density and related properties in terms of the Fox H and Meijer G functions, which can be written as series expansions. This enables arbitrary-precision calculations of the surface density and the lensing properties of realistic dark matter halo models. Furthermore, we compared the S\'ersic and Einasto surface mass densities and found differences between them, which implies that the lensing properties for both profiles differ.

In recent high resolution N-body CDM simulations, it has been had found that nonsingular three-parameter models, e.g. the Einasto profile has a better performance better than the singular two-parameter models, e.g. the Navarro, Frenk and White in the fitting of a wide range of dark matter halos. A problem with this profile is that the surface mass density is non-analytical for general values of the Einasto index. Therefore, its other lensing properties have the same problem. We obtain an exact analytical expression for the surface mass density of the Einasto profile in terms of the Fox H-function for all values of the Einasto index. With the idea of facilitate the use of the Einasto profile in lensing studies, we calculate the surface mass density, deflection angle, lens equation, deflection potential, magnification, shear and critical curves of the Einasto profile in terms of the Meijer G-function for all rational values of the Einasto index. The Meijer G-function have been implemented in several commercial and open-source computer algebra systems, thus the use of the lensing properties of the Einasto profile in strong and weak lensing studies is straighforward. We also compare the S\'ersic and Einasto surface mass densities profiles and found differences between them. This implies that the lensing properties are not equal for both profiles.

The gravitational lens equation for the tilted Plummer family of models can be reduced to a one-dimensional non-linear equation. For certain values of the slope of the radial profile it can be reduced to a polynomial form. Both forms are advantageous to find the roots, i.e. the images for a given model. The critical curve equations can also be reduced to a non-linear or polynomial form, and therefore it is useful to find the caustics. This lens model family has ample use in gravitational lens theory, and can produce up to five images.