Index: src/sans/models/media/model_functions.rst
===================================================================
--- src/sans/models/media/model_functions.rst	(revision 6386cd87d17d00e1716dde9a230015a6fb255f4f)
+++ src/sans/models/media/model_functions.rst	(revision 7f42aadf1091b713f50ced43ccd0876d4c5a6fcd)
@@ -389,5 +389,5 @@
 NIST (Kline, 2006). Figure 1 shows a comparison of the output of our model and the output of the NIST software.
 
-.. image:: img/image005.JPG
+.. image:: img/image005.jpg
 
 Figure 1: Comparison of the DANSE scattering intensity for a sphere with the output of the NIST SANS analysis software.
@@ -436,5 +436,5 @@
 ==============  ========  =============
 
-.. image:: img/image009.JPG
+.. image:: img/image009.jpg
 
 *Figure. 1D plot using the default values above (w/200 data point).*
@@ -504,5 +504,5 @@
 ==============  ========  =============
 
-.. image:: img/image012.JPG
+.. image:: img/image012.jpg
 
 *Figure. 1D plot using the default values (w/200 data point).*
@@ -525,5 +525,5 @@
 The structure is:
 
-.. image:: img/raspberry_pic.JPG
+.. image:: img/raspberry_pic.jpg
 
 where *Ro* = the radius of the large sphere, *Rp* = the radius of the smaller sphere on the surface, |delta| = the
@@ -561,5 +561,5 @@
 ==============  ========  =============
 
-.. image:: img/raspberry_plot.JPG
+.. image:: img/raspberry_plot.jpg
 
 *Figure. 1D plot using the values of /2000 data points.*
@@ -625,5 +625,5 @@
 NIST (Kline, 2006). Figure 1 shows a comparison of the output of our model and the output of the NIST software.
 
-.. image:: img/image014.JPG
+.. image:: img/image014.jpg
 
 Figure 1: Comparison of the SasView scattering intensity for a core-shell sphere with the output of the NIST SANS
@@ -686,5 +686,5 @@
 *qmax* = 0.7 -1 and the above default values.
 
-.. image:: img/image015.JPG
+.. image:: img/image015.jpg
 
 *Figure: 1D plot using the default values (w/200 data point).*
@@ -692,5 +692,5 @@
 The scattering length density profile for the default sld values (w/ 4 shells).
 
-.. image:: img/image016.JPG
+.. image:: img/image016.jpg
 
 *Figure: SLD profile against the radius of the sphere for default SLDs.*
@@ -721,5 +721,5 @@
 The *I* :sub:`0` is calculated in the following way (King, 2002)
 
-.. image:: img/secondmeq1.JPG
+.. image:: img/secondmeq1.jpg
 
 where *scale* is a scale factor, *poly* is the sld of the polymer (or surfactant) layer, *solv* is the sld of the
@@ -748,5 +748,5 @@
 ==============  ========  =============
 
-.. image:: img/secongm_fig1.JPG
+.. image:: img/secongm_fig1.jpg
 
 REFERENCE
@@ -764,5 +764,5 @@
 solvent and the shells are interleaved with layers of solvent. For *N* = 1, this returns the VesicleModel (above).
 
-.. image:: img/image020.JPG
+.. image:: img/image020.jpg
 
 The 2D scattering intensity is the same as 1D, regardless of the orientation of the *q* vector which is defined as
@@ -791,5 +791,5 @@
 is the number of shells.
 
-.. image:: img/image021.JPG
+.. image:: img/image021.jpg
 
 *Figure. 1D plot using the default values (w/200 data point).*
@@ -818,19 +818,19 @@
 The 1D scattering intensity is calculated in the following way
 
-.. image:: img/image022.GIF
-
-.. image:: img/image023.GIF
+.. image:: img/image022.gif
+
+.. image:: img/image023.gif
 
 where, for a spherically symmetric particle with a particle density |rho|\ *(r)*
 
-.. image:: img/image024.GIF
+.. image:: img/image024.gif
 
 so that
 
-.. image:: img/image025.GIF
-
-.. image:: img/image026.GIF
-
-.. image:: img/image027.GIF
+.. image:: img/image025.gif
+
+.. image:: img/image026.gif
+
+.. image:: img/image027.gif
 
 Here we assumed that the SLDs of the core and solvent are constant against *r*.
@@ -838,5 +838,5 @@
 Now lets consider the SLD of a shell, *r*\ :sub:`shelli`, defined by
 
-.. image:: img/image028.GIF
+.. image:: img/image028.gif
 
 An example of a possible SLD profile is shown below where *sld_in_shelli* (|rho|\ :sub:`in`\ ) and
@@ -846,5 +846,5 @@
 For \| *A* \| > 0,
 
-.. image:: img/image029.GIF
+.. image:: img/image029.gif
 
 For *A* ~ 0 (eg., *A* = -0.0001), this function converges to that of the linear SLD profile (ie,
@@ -852,11 +852,11 @@
 so this case is equivalent to
 
-.. image:: img/image030.GIF
-
-.. image:: img/image031.GIF
-
-.. image:: img/image032.GIF
-
-.. image:: img/image033.GIF
+.. image:: img/image030.gif
+
+.. image:: img/image031.gif
+
+.. image:: img/image032.gif
+
+.. image:: img/image033.gif
 
 For *A* = 0, the exponential function has no dependence on the radius (so that *sld_out_shell* (|rho|\ :sub:`out`) is
@@ -864,28 +864,28 @@
 factor contributed by the shells is
 
-.. image:: img/image034.GIF
-
-.. image:: img/image035.GIF
+.. image:: img/image034.gif
+
+.. image:: img/image035.gif
 
 In the equation
 
-.. image:: img/image036.GIF
+.. image:: img/image036.gif
 
 Finally, the form factor can be calculated by
 
-.. image:: img/image037.GIF
+.. image:: img/image037.gif
 
 where
 
-.. image:: img/image038.GIF
+.. image:: img/image038.gif
 
 and
 
-.. image:: img/image039.GIF
+.. image:: img/image039.gif
 
 The 2D scattering intensity is the same as *P(q)* above, regardless of the orientation of the *q* vector which is
 defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 NB: The outer most radius is used as the effective radius for *S(Q)* when *P(Q)* \* *S(Q)* is applied.
@@ -909,9 +909,9 @@
 NB: *rad_core* represents the core radius (*R1*) and *thick_shell1* (*R2* - *R1*) is the thickness of the shell1, etc.
 
-.. image:: img/image041.JPG
+.. image:: img/image041.jpg
 
 *Figure. 1D plot using the default values (w/400 point).*
 
-.. image:: img/image042.JPG
+.. image:: img/image042.jpg
 
 *Figure. SLD profile from the default values.*
@@ -945,5 +945,5 @@
 and a shell thickness, *t*.
 
-.. image:: img/image018.JPG
+.. image:: img/image018.jpg
 
 The 2D scattering intensity is the same as *P(q)* above, regardless of the orientation of the *q* vector which is
@@ -970,5 +970,5 @@
 NB: *radius* represents the core radius (*R1*) and the *thickness* (*R2* - *R1*) is the shell thickness.
 
-.. image:: img/image019.JPG
+.. image:: img/image019.jpg
 
 *Figure. 1D plot using the default values (w/200 data point).*
@@ -1001,25 +1001,25 @@
 The 1D scattering intensity is calculated in the following way:
 
-.. image:: img/image022.GIF
-
-.. image:: img/image043.GIF
+.. image:: img/image022.gif
+
+.. image:: img/image043.gif
 
 where, for a spherically symmetric particle with a particle density |rho|\ *(r)*
 
-.. image:: img/image024.GIF
+.. image:: img/image024.gif
 
 so that
 
-.. image:: img/image044.GIF
-
-.. image:: img/image045.GIF
-
-.. image:: img/image046.GIF
-
-.. image:: img/image047.GIF
-
-.. image:: img/image048.GIF
-
-.. image:: img/image027.GIF
+.. image:: img/image044.gif
+
+.. image:: img/image045.gif
+
+.. image:: img/image046.gif
+
+.. image:: img/image047.gif
+
+.. image:: img/image048.gif
+
+.. image:: img/image027.gif
 
 Here we assumed that the SLDs of the core and solvent are constant against *r*. The SLD at the interface between
@@ -1028,13 +1028,13 @@
 1) Exp
 
-.. image:: img/image049.GIF
+.. image:: img/image049.gif
 
 2) Power-Law
 
-.. image:: img/image050.GIF
+.. image:: img/image050.gif
 
 3) Erf
 
-.. image:: img/image051.GIF
+.. image:: img/image051.gif
 
 The functions are normalized so that they vary between 0 and 1, and they are constrained such that the SLD is
@@ -1044,9 +1044,9 @@
 to the form factor *P(q)*
 
-.. image:: img/image052.GIF
-
-.. image:: img/image053.GIF
-
-.. image:: img/image054.GIF
+.. image:: img/image052.gif
+
+.. image:: img/image053.gif
+
+.. image:: img/image054.gif
 
 where we assume that |rho|\ :sub:`inter_i`\ *(r)* can be approximately linear within a sub-layer *j*.
@@ -1054,22 +1054,22 @@
 In the equation
 
-.. image:: img/image055.GIF
+.. image:: img/image055.gif
 
 Finally, the form factor can be calculated by
 
-.. image:: img/image037.GIF
+.. image:: img/image037.gif
 
 where
 
-.. image:: img/image038.GIF
+.. image:: img/image038.gif
 
 and
 
-.. image:: img/image056.GIF
+.. image:: img/image056.gif
 
 The 2D scattering intensity is the same as *P(q)* above, regardless of the orientation of the *q* vector which is
 defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 NB: The outer most radius is used as the effective radius for *S(Q)* when *P(Q)* \* *S(Q)* is applied.
@@ -1098,9 +1098,9 @@
 NB: *rad_core0* represents the core radius (*R1*).
 
-.. image:: img/image057.JPG
+.. image:: img/image057.jpg
 
 *Figure. 1D plot using the default values (w/400 point).*
 
-.. image:: img/image058.JPG
+.. image:: img/image058.jpg
 
 *Figure. SLD profile from the default values.*
@@ -1249,9 +1249,9 @@
 and |phi|. Those angles are defined in Figure 1.
 
-.. image:: img/image061.JPG
+.. image:: img/image061.jpg
 
 *Figure 1. Definition of the angles for oriented cylinders.*
 
-.. image:: img/image062.JPG
+.. image:: img/image062.jpg
 
 *Figure 2. Examples of the angles for oriented pp against the detector plane.*
@@ -1286,5 +1286,5 @@
 NIST (Kline, 2006). Figure 3 shows a comparison of the 1D output of our model and the output of the NIST software.
 
-.. image:: img/image065.JPG
+.. image:: img/image065.jpg
 
 *Figure 3: Comparison of the SasView scattering intensity for a cylinder with the output of the NIST SANS analysis*
@@ -1301,5 +1301,5 @@
 distribution *p(*\ |theta|,\ |phi|\ *)* = 1.0. Figure 4 shows the result of such a cross-check.
 
-.. image:: img/image066.JPG
+.. image:: img/image066.jpg
 
 *Figure 4: Comparison of the intensity for uniformly distributed cylinders calculated from our 2D model and the*
@@ -1351,5 +1351,5 @@
 ==============  ========  =============
 
-.. image:: img/image074.JPG
+.. image:: img/image074.jpg
 
 *Figure. 1D plot using the default values (w/1000 data point).*
@@ -1358,9 +1358,9 @@
 (Kline, 2006).
 
-.. image:: img/image061.JPG
+.. image:: img/image061.jpg
 
 *Figure. Definition of the angles for the oriented HollowCylinderModel.*
 
-.. image:: img/image062.JPG
+.. image:: img/image062.jpg
 
 *Figure. Examples of the angles for oriented pp against the detector plane.*
@@ -1388,5 +1388,5 @@
 The Capped Cylinder geometry is defined as
 
-.. image:: img/image112.JPG
+.. image:: img/image112.jpg
 
 where *r* is the radius of the cylinder. All other parameters are as defined in the diagram. Since the end cap radius
@@ -1397,9 +1397,9 @@
 The scattered intensity *I(q)* is calculated as
 
-.. image:: img/image113.JPG
+.. image:: img/image113.jpg
 
 where the amplitude *A(q)* is given as
 
-.. image:: img/image114.JPG
+.. image:: img/image114.jpg
 
 The < > brackets denote an average of the structure over all orientations. <\ *A*\ :sup:`2`\ *(q)*> is then the form
@@ -1409,9 +1409,9 @@
 The volume of the Capped Cylinder is (with *h* as a positive value here)
 
-.. image:: img/image115.JPG
+.. image:: img/image115.jpg
 
 and its radius-of-gyration
 
-.. image:: img/image116.JPG
+.. image:: img/image116.jpg
 
 **The requirement that** *R* >= *r* **is not enforced in the model! It is up to you to restrict this during analysis.**
@@ -1432,5 +1432,5 @@
 ==============  ========  =============
 
-.. image:: img/image117.JPG
+.. image:: img/image117.jpg
 
 *Figure. 1D plot using the default values (w/256 data point).*
@@ -1439,9 +1439,9 @@
 |theta| = 45 deg and |phi| =0 deg with default values for other parameters
 
-.. image:: img/image118.JPG
+.. image:: img/image118.jpg
 
 *Figure. 2D plot (w/(256X265) data points).*
 
-.. image:: img/image061.JPG
+.. image:: img/image061.jpg
 
 *Figure. Definition of the angles for oriented 2D cylinders.*
@@ -1485,5 +1485,5 @@
 the outer shell is given by *L+2t*. *J1* is the first order Bessel function.
 
-.. image:: img/image069.JPG
+.. image:: img/image069.jpg
 
 To provide easy access to the orientation of the core-shell cylinder, we define the axis of the cylinder using two
@@ -1520,5 +1520,5 @@
 NIST (Kline, 2006). Figure 1 shows a comparison of the 1D output of our model and the output of the NIST software.
 
-.. image:: img/image070.JPG
+.. image:: img/image070.jpg
 
 *Figure 1: Comparison of the SasView scattering intensity for a core-shell cylinder with the output of the NIST SANS*
@@ -1531,5 +1531,5 @@
 2D output using a uniform distribution *p(*\ |theta|,\ |phi|\ *)* = 1.0. Figure 2 shows the result of such a cross-check.
 
-.. image:: img/image071.JPG
+.. image:: img/image071.jpg
 
 *Figure 2: Comparison of the intensity for uniformly distributed core-shell cylinders calculated from our 2D model and*
@@ -1538,9 +1538,9 @@
 *Solvent_sld* = 1e-6 |Ang^-2|, and *Background* = 0.0 |cm^-1|.
 
-.. image:: img/image061.JPG
+.. image:: img/image061.jpg
 
 *Figure. Definition of the angles for oriented core-shell cylinders.*
 
-.. image:: img/image062.JPG
+.. image:: img/image062.jpg
 
 *Figure. Examples of the angles for oriented pp against the detector plane.*
@@ -1591,9 +1591,9 @@
 All angle parameters are valid and given only for 2D calculation; ie, an oriented system.
 
-.. image:: img/image101.JPG
+.. image:: img/image101.jpg
 
 *Figure. Definition of angles for 2D*
 
-.. image:: img/image062.JPG
+.. image:: img/image062.jpg
 
 *Figure. Examples of the angles for oriented elliptical cylinders against the detector plane.*
@@ -1614,5 +1614,5 @@
 ==============  ========  =============
 
-.. image:: img/image102.JPG
+.. image:: img/image102.jpg
 
 *Figure. 1D plot using the default values (w/1000 data point).*
@@ -1625,5 +1625,5 @@
 and 76 degrees are taken for the angles of |theta|, |phi|, and |bigpsi| respectively).
 
-.. image:: img/image103.GIF
+.. image:: img/image103.gif
 
 *Figure. Comparison between 1D and averaged 2D.*
@@ -1632,5 +1632,5 @@
 the results of the averaging by varying the number of angular bins.
 
-.. image:: img/image104.GIF
+.. image:: img/image104.gif
 
 *Figure. The intensities averaged from 2D over different numbers of bins and angles.*
@@ -1660,5 +1660,5 @@
 *2.1.19.1. Definition*
 
-.. image:: img/image075.JPG
+.. image:: img/image075.jpg
 
 The chain of contour length, *L*, (the total length) can be described as a chain of some number of locally stiff
@@ -1683,5 +1683,5 @@
 ==============  ========  =============
 
-.. image:: img/image076.JPG
+.. image:: img/image076.jpg
 
 *Figure. 1D plot using the default values (w/1000 data point).*
@@ -1738,5 +1738,5 @@
 - The scattering function is negative for a range of parameter values and q-values that are experimentally accessible. A correction function has been added to give the proper behavior.
 
-.. image:: img/image077.JPG
+.. image:: img/image077.jpg
 
 The chain of contour length, *L*, (the total length) can be described as a chain of some number of locally stiff
@@ -1780,5 +1780,5 @@
 ==============  ========  =============
 
-.. image:: img/image078.JPG
+.. image:: img/image078.jpg
 
 *Figure. 1D plot using the default values (w/200 data points).*
@@ -1807,5 +1807,5 @@
 and SLDs.
 
-.. image:: img/image240.PNG
+.. image:: img/image240.png
 
 *(Graphic from DOI: 10.1039/C0NP00002G)*
@@ -1839,9 +1839,9 @@
 *Figure. 1D plot using the default values (w/200 data point).*
 
-.. image:: img/image061.JPG
+.. image:: img/image061.jpg
 
 *Figure. Definition of the angles for the oriented CoreShellBicelleModel.*
 
-.. image:: img/image062.JPG
+.. image:: img/image062.jpg
 
 *Figure. Examples of the angles for oriented pp against the detector plane.*
@@ -1869,5 +1869,5 @@
 The barbell geometry is defined as
 
-.. image:: img/image105.JPG
+.. image:: img/image105.jpg
 
 where *r* is the radius of the cylinder. All other parameters are as defined in the diagram.
@@ -1892,10 +1892,10 @@
 The volume of the barbell is
 
-.. image:: img/image108.JPG
+.. image:: img/image108.jpg
 
 
 and its radius-of-gyration is
 
-.. image:: img/image109.JPG
+.. image:: img/image109.jpg
 
 **The requirement that** *R* >= *r* **is not enforced in the model!** It is up to you to restrict this during analysis.
@@ -1916,5 +1916,5 @@
 ==============  ========  =============
 
-.. image:: img/image110.JPG
+.. image:: img/image110.jpg
 
 *Figure. 1D plot using the default values (w/256 data point).*
@@ -1923,13 +1923,13 @@
 |theta| = 45 deg and |phi| = 0 deg with default values for other parameters
 
-.. image:: img/image111.JPG
+.. image:: img/image111.jpg
 
 *Figure. 2D plot (w/(256X265) data points).*
 
-.. image:: img/image061.JPG
+.. image:: img/image061.jpg
 
 *Figure. Examples of the angles for oriented pp against the detector plane.*
 
-.. image:: img/image062.JPG
+.. image:: img/image062.jpg
 
 Figure. Definition of the angles for oriented 2D barbells.
@@ -1963,5 +1963,5 @@
 *2.1.23.1 Definition*
 
-.. image:: img/image079.GIF
+.. image:: img/image079.gif
 
 The scattering intensity *I(q)* is
@@ -2007,13 +2007,13 @@
 ==============  ========  =============
 
-.. image:: img/image085.JPG
+.. image:: img/image085.jpg
 
 *Figure. 1D plot using the default values (w/1000 data point).*
 
-.. image:: img/image086.JPG
+.. image:: img/image086.jpg
 
 *Figure. Examples of the angles for oriented stackeddisks against the detector plane.*
 
-.. image:: img/image062.JPG
+.. image:: img/image062.jpg
 
 *Figure. Examples of the angles for oriented pp against the detector plane.*
@@ -2038,5 +2038,5 @@
 This model provides the form factor, *P(q)*, for a 'pringle' or 'saddle-shaped' object (a hyperbolic paraboloid).
 
-.. image:: img/image241.PNG
+.. image:: img/image241.png
 
 *(Graphic from Matt Henderson, matt@matthen.com)*
@@ -2128,10 +2128,10 @@
 above.
 
-.. image:: img/image121.JPG
+.. image:: img/image121.jpg
 
 The *axis_theta* and *axis_phi* parameters are not used for the 1D output. Our implementation of the scattering
 kernel and the 1D scattering intensity use the c-library from NIST.
 
-.. image:: img/image122.JPG
+.. image:: img/image122.jpg
 
 *Figure. The angles for oriented ellipsoid.*
@@ -2143,5 +2143,5 @@
 software.
 
-.. image:: img/image123.JPG
+.. image:: img/image123.jpg
 
 *Figure 1: Comparison of the SasView scattering intensity for an ellipsoid with the output of the NIST SANS analysis*
@@ -2154,5 +2154,5 @@
 cross-check.
 
-.. image:: img/image124.JPG
+.. image:: img/image124.jpg
 
 *Figure 2: Comparison of the intensity for uniformly distributed ellipsoids calculated from our 2D model and the*
@@ -2186,5 +2186,5 @@
 all orientations for 1D.
 
-.. image:: img/image125.GIF
+.. image:: img/image125.gif
 
 The returned value is in units of |cm^-1|, on absolute scale.
@@ -2220,9 +2220,9 @@
 ==============  ========  =============
 
-.. image:: img/image127.JPG
+.. image:: img/image127.jpg
 
 *Figure. 1D plot using the default values (w/1000 data point).*
 
-.. image:: img/image122.JPG
+.. image:: img/image122.jpg
 
 *Figure. The angles for oriented CoreShellEllipsoid.*
@@ -2313,5 +2313,5 @@
 where the volume *V* = (4/3)\ |pi| (*Ra* *Rb* *Rc*), and the averaging < > is applied over all orientations for 1D.
 
-.. image:: img/image128.JPG
+.. image:: img/image128.jpg
 
 The returned value is in units of |cm^-1|, on absolute scale.
@@ -2349,5 +2349,5 @@
 ==============  ========  =============
 
-.. image:: img/image130.JPG
+.. image:: img/image130.jpg
 
 *Figure. 1D plot using the default values (w/1000 data point).*
@@ -2360,9 +2360,9 @@
 angles of |theta|, |phi|, and |psi| respectively).
 
-.. image:: img/image131.GIF
+.. image:: img/image131.gif
 
 *Figure. Comparison between 1D and averaged 2D.*
 
-.. image:: img/image132.JPG
+.. image:: img/image132.jpg
 
 *Figure. The angles for oriented ellipsoid.*
@@ -2399,5 +2399,5 @@
 The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 The returned value is in units of |cm^-1|, on absolute scale. In the parameters, *sld_bi* = SLD of the bilayer,
@@ -2414,5 +2414,5 @@
 ==============  ========  =============
 
-.. image:: img/image135.JPG
+.. image:: img/image135.jpg
 
 *Figure. 1D plot using the default values (w/1000 data point).*
@@ -2444,5 +2444,5 @@
 The form factor is
 
-.. image:: img/image137.JPG
+.. image:: img/image137.jpg
 
 where |delta|\ T = tail length (or *t_length*), |delta|\ H = head thickness (or *h_thickness*),
@@ -2451,5 +2451,5 @@
 The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 The returned value is in units of |cm^-1|, on absolute scale. In the parameters, *sld_tail* = SLD of the tail group,
@@ -2468,5 +2468,5 @@
 ==============  ========  =============
 
-.. image:: img/image138.JPG
+.. image:: img/image138.jpg
 
 *Figure. 1D plot using the default values (w/1000 data point).*
@@ -2519,5 +2519,5 @@
 The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 The returned value is in units of |cm^-1|, on absolute scale.
@@ -2535,5 +2535,5 @@
 ==============  ========  =============
 
-.. image:: img/image142.JPG
+.. image:: img/image142.jpg
 
 *Figure. 1D plot using the default values (w/6000 data point).*
@@ -2587,5 +2587,5 @@
 The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 The returned value is in units of |cm^-1|, on absolute scale. In the parameters, *sld_tail* = SLD of the tail group,
@@ -2607,5 +2607,5 @@
 ==============  ========  =============
 
-.. image:: img/image144.JPG
+.. image:: img/image144.jpg
 
 *Figure. 1D plot using the default values (w/6000 data point).*
@@ -2634,9 +2634,9 @@
 The scattering intensity *I(q)* is calculated as
 
-.. image:: img/image145.JPG
+.. image:: img/image145.jpg
 
 The form factor of the bilayer is approximated as the cross section of an infinite, planar bilayer of thickness *t*
 
-.. image:: img/image146.JPG
+.. image:: img/image146.jpg
 
 Here, the scale factor is used instead of the mass per area of the bilayer (*G*). The scale factor is the volume
@@ -2647,9 +2647,9 @@
 Non-integer numbers of stacks are calculated as a linear combination of the lower and higher values
 
-.. image:: img/image147.JPG
+.. image:: img/image147.jpg
 
 The 2D scattering intensity is the same as 1D, regardless of the orientation of the *q* vector which is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 The parameters of the model are *Nlayers* = no. of layers, and *pd_spacing* = polydispersity of spacing.
@@ -2668,5 +2668,5 @@
 ==============  ========  =============
 
-.. image:: img/image148.JPG
+.. image:: img/image148.jpg
 
 *Figure. 1D plot using the default values above (w/20000 data point).*
@@ -2695,5 +2695,5 @@
 The scattering intensity *I(q)* is calculated as
 
-.. image:: img/image149.JPG
+.. image:: img/image149.jpg
 
 where *scale* is the volume fraction of spheres, *Vp* is the volume of the primary particle, *V(lattice)* is a volume
@@ -2707,9 +2707,9 @@
 and nearest neighbor separation *D* is
 
-.. image:: img/image150.JPG
+.. image:: img/image150.jpg
 
 The distortion factor (one standard deviation) of the paracrystal is included in the calculation of *Z(q)*
 
-.. image:: img/image151.JPG
+.. image:: img/image151.jpg
 
 where *g* is a fractional distortion based on the nearest neighbor distance.
@@ -2717,14 +2717,14 @@
 The simple cubic lattice is
 
-.. image:: img/image152.JPG
+.. image:: img/image152.jpg
 
 For a crystal, diffraction peaks appear at reduced *q*\ -values given by
 
-.. image:: img/image153.JPG
+.. image:: img/image153.jpg
 
 where for a simple cubic lattice any *h*\ , *k*\ , *l* are allowed and none are forbidden. Thus the peak positions
 correspond to (just the first 5)
 
-.. image:: img/image154.JPG
+.. image:: img/image154.jpg
 
 **NB: The calculation of** *Z(q)* **is a double numerical integral that must be carried out with a high density of**
@@ -2748,5 +2748,5 @@
 default values.
 
-.. image:: img/image155.JPG
+.. image:: img/image155.jpg
 
 *Figure. 1D plot in the linear scale using the default values (w/200 data point).*
@@ -2756,7 +2756,7 @@
 computation.
 
-.. image:: img/image156.JPG
-
-.. image:: img/image157.JPG
+.. image:: img/image156.jpg
+
+.. image:: img/image157.jpg
 
 *Figure. 2D plot using the default values (w/200X200 pixels).*
@@ -2786,5 +2786,5 @@
 The scattering intensity *I(q)* is calculated as
 
-.. image:: img/image158.JPG
+.. image:: img/image158.jpg
 
 where *scale* is the volume fraction of spheres, *Vp* is the volume of the primary particle, *V(lattice)* is a volume
@@ -2798,9 +2798,9 @@
 *R* and nearest neighbor separation *D* is
 
-.. image:: img/image159.JPG
+.. image:: img/image159.jpg
 
 The distortion factor (one standard deviation) of the paracrystal is included in the calculation of *Z(q)*
 
-.. image:: img/image160.JPG
+.. image:: img/image160.jpg
 
 where *g* is a fractional distortion based on the nearest neighbor distance.
@@ -2808,14 +2808,14 @@
 The face-centered cubic lattice is
 
-.. image:: img/image161.JPG
+.. image:: img/image161.jpg
 
 For a crystal, diffraction peaks appear at reduced q-values given by
 
-.. image:: img/image162.JPG
+.. image:: img/image162.jpg
 
 where for a face-centered cubic lattice *h*\ , *k*\ , *l* all odd or all even are allowed and reflections where
 *h*\ , *k*\ , *l* are mixed odd/even are forbidden. Thus the peak positions correspond to (just the first 5)
 
-.. image:: img/image163.JPG
+.. image:: img/image163.jpg
 
 **NB: The calculation of** *Z(q)* **is a double numerical integral that must be carried out with a high density of**
@@ -2839,5 +2839,5 @@
 default values.
 
-.. image:: img/image164.JPG
+.. image:: img/image164.jpg
 
 *Figure. 1D plot in the linear scale using the default values (w/200 data point).*
@@ -2847,7 +2847,7 @@
 computation.
 
-.. image:: img/image165.GIF
-
-.. image:: img/image166.JPG
+.. image:: img/image165.gif
+
+.. image:: img/image166.jpg
 
 *Figure. 2D plot using the default values (w/200X200 pixels).*
@@ -2877,5 +2877,5 @@
 The scattering intensity *I(q)* is calculated as
 
-.. image:: img/image167.JPG
+.. image:: img/image167.jpg
 
 where *scale* is the volume fraction of spheres, *Vp* is the volume of the primary particle, *V(lattice)* is a volume
@@ -2889,9 +2889,9 @@
 *R* and nearest neighbor separation *D* is
 
-.. image:: img/image159.JPG
+.. image:: img/image159.jpg
 
 The distortion factor (one standard deviation) of the paracrystal is included in the calculation of *Z(q)*
 
-.. image:: img/image160.JPG
+.. image:: img/image160.jpg
 
 where *g* is a fractional distortion based on the nearest neighbor distance.
@@ -2899,14 +2899,14 @@
 The body-centered cubic lattice is
 
-.. image:: img/image168.JPG
+.. image:: img/image168.jpg
 
 For a crystal, diffraction peaks appear at reduced q-values given by
 
-.. image:: img/image162.JPG
+.. image:: img/image162.jpg
 
 where for a body-centered cubic lattice, only reflections where (\ *h* + *k* + *l*\ ) = even are allowed and
 reflections where (\ *h* + *k* + *l*\ ) = odd are forbidden. Thus the peak positions correspond to (just the first 5)
 
-.. image:: img/image169.JPG
+.. image:: img/image169.jpg
 
 **NB: The calculation of** *Z(q)* **is a double numerical integral that must be carried out with a high density of**
@@ -2930,5 +2930,5 @@
 default values.
 
-.. image:: img/image170.JPG
+.. image:: img/image170.jpg
 
 *Figure. 1D plot in the linear scale using the default values (w/200 data point).*
@@ -2938,7 +2938,7 @@
 computation.
 
-.. image:: img/image165.GIF
-
-.. image:: img/image171.JPG
+.. image:: img/image165.gif
+
+.. image:: img/image171.jpg
 
 *Figure. 2D plot using the default values (w/200X200 pixels).*
@@ -2967,5 +2967,5 @@
 For information about polarised and magnetic scattering, click here_.
 
-.. image:: img/image087.JPG
+.. image:: img/image087.jpg
 
 *2.1.37.1. Definition*
@@ -2991,9 +2991,9 @@
 parallel to the *x*-axis of the detector.
 
-.. image:: img/image090.JPG
+.. image:: img/image090.jpg
 
 *Figure. Definition of angles for 2D*.
 
-.. image:: img/image091.JPG
+.. image:: img/image091.jpg
 
 *Figure. Examples of the angles for oriented pp against the detector plane.*
@@ -3010,5 +3010,5 @@
 ==============  ========  =============
 
-.. image:: img/image092.JPG
+.. image:: img/image092.jpg
 
 *Figure. 1D plot using the default values (w/1000 data point).*
@@ -3021,5 +3021,5 @@
 angles of |theta|, |phi|, and |psi| respectively).
 
-.. image:: img/image093.GIF
+.. image:: img/image093.gif
 
 *Figure. Comparison between 1D and averaged 2D.*
@@ -3055,10 +3055,10 @@
 dimensions *A*, *B*, *C* such that *A* < *B* < *C*.
 
-.. image:: img/image087.JPG
+.. image:: img/image087.jpg
 
 There are rectangular "slabs" of thickness *tA* that add to the *A* dimension (on the *BC* faces). There are similar
 slabs on the *AC* (= *tB*) and *AB* (= *tC*) faces. The projection in the *AB* plane is then
 
-.. image:: img/image094.JPG
+.. image:: img/image094.jpg
 
 The volume of the solid is
@@ -3096,9 +3096,9 @@
 parallel to the *x*-axis of the detector.
 
-.. image:: img/image090.JPG
+.. image:: img/image090.jpg
 
 *Figure. Definition of angles for 2D*.
 
-.. image:: img/image091.JPG
+.. image:: img/image091.jpg
 
 *Figure. Examples of the angles for oriented cspp against the detector plane.*
@@ -3125,9 +3125,9 @@
 ==============  ========  =============
 
-.. image:: img/image096.JPG
+.. image:: img/image096.jpg
 
 *Figure. 1D plot using the default values (w/256 data points).*
 
-.. image:: img/image097.JPG
+.. image:: img/image097.jpg
 
 *Figure. 2D plot using the default values (w/(256X265) data points).*
@@ -3392,5 +3392,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -3402,5 +3402,5 @@
 ==============  ========  =============
 
-.. image:: img/image173.JPG
+.. image:: img/image173.jpg
 
 *Figure. 1D plot using the default values (w/200 data point).*
@@ -3429,5 +3429,5 @@
 The scattering intensity *I(q)* is calculated as
 
-.. image:: img/image174.JPG
+.. image:: img/image174.jpg
 
 Here the peak position is related to the d-spacing as *Q0* = 2|pi| / *d0*.
@@ -3435,5 +3435,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==================  ========  =============
@@ -3449,5 +3449,5 @@
 ==================  ========  =============
 
-.. image:: img/image175.JPG
+.. image:: img/image175.jpg
 
 *Figure. 1D plot using the default values (w/200 data point).*
@@ -3473,5 +3473,5 @@
 The scattering intensity *I(q)* is calculated as
 
-.. image:: img/image176.JPG
+.. image:: img/image176.jpg
 
 The first term describes Porod scattering from clusters (exponent = n) and the second term is a Lorentzian function
@@ -3484,5 +3484,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ====================  ========  =============
@@ -3497,5 +3497,5 @@
 ====================  ========  =============
 
-.. image:: img/image177.JPG
+.. image:: img/image177.jpg
 
 *Figure. 1D plot using the default values (w/500 data points).*
@@ -3524,5 +3524,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -3534,5 +3534,5 @@
 ==============  ========  =============
 
-.. image:: img/image179.JPG
+.. image:: img/image179.jpg
 
 *Â Figure. 1D plot using the default values (w/200 data point).*
@@ -3563,5 +3563,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -3573,5 +3573,5 @@
 ==============  ========  =============
 
-.. image:: img/image181.JPG
+.. image:: img/image181.jpg
 
 *Â Figure. 1D plot using the default values (w/200 data point).*
@@ -3606,5 +3606,5 @@
 ==============  ========  =============
 
-.. image:: img/image183.JPG
+.. image:: img/image183.jpg
 
 *Figure. 1D plot using the default values (w/200 data point).*
@@ -3629,5 +3629,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -3640,5 +3640,5 @@
 ==============  ========  =============
 
-.. image:: img/image185.JPG
+.. image:: img/image185.jpg
 
 *Figure. 1D plot using the default values (w/200 data point).*
@@ -3673,5 +3673,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -3687,5 +3687,5 @@
 ==============  ========  =============
 
-.. image:: img/image187.JPG
+.. image:: img/image187.jpg
 
 *Figure. 1D plot using the default values (w/200 data point).*
@@ -3705,5 +3705,5 @@
 *2.2.9.1. Definition*
 
-.. image:: img/mass_fractal_eq1.JPG
+.. image:: img/mass_fractal_eq1.jpg
 
 where *R* is the radius of the building block, *Dm* is the **mass** fractal dimension, |zeta| is the cut-off length,
@@ -3724,5 +3724,5 @@
 ==============  ========  =============
 
-.. image:: img/mass_fractal_fig1.JPG
+.. image:: img/mass_fractal_fig1.jpg
 
 *Figure. 1D plot using default values.*
@@ -3764,5 +3764,5 @@
 ==============  ========  =============
 
-.. image:: img/surface_fractal_fig1.JPG
+.. image:: img/surface_fractal_fig1.jpg
 
 *Figure. 1D plot using default values.*
@@ -3812,5 +3812,5 @@
 ==============  ========  =============
 
-.. image:: img/masssurface_fractal_fig1.JPG
+.. image:: img/masssurface_fractal_fig1.jpg
 
 *Figure. 1D plot using default values.*
@@ -3838,5 +3838,5 @@
 *2.2.12.1. Definition*
 
-.. image:: img/fractcore_eq1.GIF
+.. image:: img/fractcore_eq1.gif
 
 The form factor *P(q)* is that from CoreShellModel_ with *bkg* = 0
@@ -3854,5 +3854,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -3870,5 +3870,5 @@
 ==============  ========  =============
 
-.. image:: img/image188.JPG
+.. image:: img/image188.jpg
 
 *Figure. 1D plot using the default values (w/500 data points).*
@@ -3895,5 +3895,5 @@
 The scattering intensity *I(q)* is calculated as (eqn 5 from the reference)
 
-.. image:: img/image189.JPG
+.. image:: img/image189.jpg
 
 |bigzeta| is the length scale of the static correlations in the gel, which can be attributed to the "frozen-in"
@@ -3907,5 +3907,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ===================================  ========  =============
@@ -3919,5 +3919,5 @@
 ===================================  ========  =============
 
-.. image:: img/image190.JPG
+.. image:: img/image190.jpg
 
 *Figure. 1D plot using the default values (w/500 data points).*
@@ -3947,5 +3947,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -3988,5 +3988,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -4017,5 +4017,5 @@
 The following functional form is used
 
-.. image:: img/image193.JPG
+.. image:: img/image193.jpg
 
 This is based on the generalized Guinier law for such elongated objects (see the Glatter reference below). For 3D
@@ -4026,9 +4026,9 @@
 Enforcing the continuity of the Guinier and Porod functions and their derivatives yields
 
-.. image:: img/image194.JPG
+.. image:: img/image194.jpg
 
 and
 
-.. image:: img/image195.JPG
+.. image:: img/image195.jpg
 
 Note that
@@ -4054,5 +4054,5 @@
 ==============================  ========  =============
 
-.. image:: img/image196.JPG
+.. image:: img/image196.jpg
 
 *Figure. 1D plot using the default values (w/500 data points).*
@@ -4082,5 +4082,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -4110,5 +4110,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -4121,5 +4121,5 @@
 ==============  ========  =============
 
-.. image:: img/image199.JPG
+.. image:: img/image199.jpg
 
 *Figure. 1D plot using the default values (w/500 data points).*
@@ -4143,5 +4143,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -4154,5 +4154,5 @@
 ==============  ========  =============
 
-.. image:: img/image201.JPG
+.. image:: img/image201.jpg
 
 *Figure. 1D plot using the default values (w/500 data points).*
@@ -4190,5 +4190,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 This example dataset is produced using 200 data points, using 200 data points,
@@ -4204,5 +4204,5 @@
 ==============  ========  =============
 
-.. image:: img/image205.JPG
+.. image:: img/image205.jpg
 
 *Figure. 1D plot using the default values (w/200 data point).*
@@ -4230,5 +4230,5 @@
 The form factor  was originally presented in the following integral form (Benoit, 1957)
 
-.. image:: img/image206.JPG
+.. image:: img/image206.jpg
 
 where |nu| is the excluded volume parameter (which is related to the Porod exponent *m* as |nu| = 1 / *m*), *a* is the
@@ -4236,17 +4236,17 @@
 into an almost analytical form as follows (Hammouda, 1993)
 
-.. image:: img/image207.JPG
+.. image:: img/image207.jpg
 
 where |gamma|\ *(x,U)* is the incomplete gamma function
 
-.. image:: img/image208.JPG
+.. image:: img/image208.jpg
 
 and the variable *U* is given in terms of the scattering vector *Q* as
 
-.. image:: img/image209.JPG
+.. image:: img/image209.jpg
 
 The square of the radius-of-gyration is defined as
 
-.. image:: img/image210.JPG
+.. image:: img/image210.jpg
 
 Note that this model applies only in the mass fractal range (ie, 5/3 <= *m* <= 3) and **does not** apply to surface
@@ -4256,5 +4256,5 @@
 A low-*Q* expansion yields the Guinier form and a high-*Q* expansion yields the Porod form which is given by
 
-.. image:: img/image211.JPG
+.. image:: img/image211.jpg
 
 Here |biggamma|\ *(x)* = |gamma|\ *(x,inf)* is the gamma function.
@@ -4262,14 +4262,14 @@
 The asymptotic limit is dominated by the first term
 
-.. image:: img/image212.JPG
+.. image:: img/image212.jpg
 
 The special case when |nu| = 0.5 (or *m* = 1/|nu| = 2) corresponds to Gaussian chains for which the form factor is given
 by the familiar Debye_ function.
 
-.. image:: img/image213.JPG
+.. image:: img/image213.jpg
 
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 This example dataset is produced using 200 data points, *qmin* = 0.001 |Ang^-1|, *qmax* = 0.2 |Ang^-1| and the default
@@ -4285,5 +4285,5 @@
 ===================  ========  =============
 
-.. image:: img/image214.JPG
+.. image:: img/image214.jpg
 
 *Figure. 1D plot using the default values (w/500 data points).*
@@ -4369,5 +4369,5 @@
 =======================  ========  =============
 
-.. image:: img/image215.JPG
+.. image:: img/image215.jpg
 
 *Figure. 1D plot using the default values (w/500 data points).*
@@ -4399,5 +4399,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ===============================  ========  =============
@@ -4413,5 +4413,5 @@
 ===============================  ========  =============
 
-.. image:: img/image217.JPG
+.. image:: img/image217.jpg
 
 *Figure. 1D plot using the default values (w/500 data points).*
@@ -4435,5 +4435,5 @@
 The scattering intensity *I(q)* is calculated as
 
-.. image:: img/image218.JPG
+.. image:: img/image218.jpg
 
 where *qc* is the location of the crossover from one slope to the other. The scaling *coef_A* sets the overall
@@ -4445,5 +4445,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -4457,5 +4457,5 @@
 ==============  ========  =============
 
-.. image:: img/image219.JPG
+.. image:: img/image219.jpg
 
 *Figure. 1D plot using the default values (w/500 data points).*
@@ -4486,5 +4486,5 @@
 The empirical fit function isÂ 
 
-.. image:: img/image220.JPG
+.. image:: img/image220.jpg
 
 For each level, the four parameters *Gi*, *Rg,i*, *Bi* and *Pi* must be chosen.Â 
@@ -4497,5 +4497,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -4514,5 +4514,5 @@
 ==============  ========  =============
 
-.. image:: img/image221.JPG
+.. image:: img/image221.jpg
 
 *Figure. 1D plot using the default values (w/500 data points).*
@@ -4567,9 +4567,9 @@
 The scattered intensity *I(q)* is calculated as
 
-.. image:: img/image233.GIF
+.. image:: img/image233.gif
 
 where
 
-.. image:: img/image234.GIF
+.. image:: img/image234.gif
 
 Note that the first term reduces to the Ornstein-Zernicke equation when *D* = 2; ie, when the Flory exponent is 0.5
@@ -4587,5 +4587,5 @@
 ============================  ========  =============
 
-.. image:: img/image235.GIF
+.. image:: img/image235.gif
 
 *Figure. 1D plot using the default values (w/300 data points).*
@@ -4609,13 +4609,13 @@
 For a star with *f* arms:
 
-.. image:: img/star1.PNG
+.. image:: img/star1.png
 
 where
 
-.. image:: img/star2.PNG
+.. image:: img/star2.png
 
 and
 
-.. image:: img/star3.PNG
+.. image:: img/star3.png
 
 is the square of the ensemble average radius-of-gyration of an arm.
@@ -4646,10 +4646,10 @@
 Also see ReflectivityIIModel_.
 
-.. image:: img/image231.BMP
+.. image:: img/image231.bmp
 
 *Figure. Comparison (using the SLD profile below) with the NIST web calculation (circles)*
 http://www.ncnr.nist.gov/resources/reflcalc.html
 
-.. image:: img/image232.GIF
+.. image:: img/image232.gif
 
 *Figure. SLD profile used for the calculation (above).*
@@ -4675,13 +4675,13 @@
 1) Erf
 
-.. image:: img/image051.GIF
+.. image:: img/image051.gif
 
 2) Power-Law
 
-.. image:: img/image050.GIF
+.. image:: img/image050.gif
 
 3) Exp
 
-.. image:: img/image049.GIF
+.. image:: img/image049.gif
 
 The constant *A* in the expressions above (but the parameter *nu* in the model!) is an input.
@@ -4713,5 +4713,5 @@
 For a 2D plot, the wave transfer is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -4722,5 +4722,5 @@
 ==============  ========  =============
 
-.. image:: img/image224.JPG
+.. image:: img/image224.jpg
 
 *Figure. 1D plot using the default values (in linear scale).*
@@ -4754,5 +4754,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  =========  =============
@@ -4765,5 +4765,5 @@
 ==============  =========  =============
 
-.. image:: img/image226.JPG
+.. image:: img/image226.jpg
 
 *Figure. 1D plot using the default values (in linear scale).*
@@ -4806,5 +4806,5 @@
 ==============  ========  =============
 
-.. image:: img/image227.JPG
+.. image:: img/image227.jpg
 
 *Figure. 1D plot using the default values (in linear scale).*
@@ -4852,5 +4852,5 @@
 For 2D data: The 2D scattering intensity is calculated in the same way as 1D, where the *q* vector is defined as
 
-.. image:: img/image040.GIF
+.. image:: img/image040.gif
 
 ==============  ========  =============
@@ -4863,5 +4863,5 @@
 ==============  ========  =============
 
-.. image:: img/image230.JPG
+.. image:: img/image230.jpg
 
 *Figure. 1D plot using the default values (in linear scale).*
