charmed: Composite Hindered and Restricted Model for Diffusion

https://mybinder.org/badge_logo.svg

Contents

% This m-file has been automatically generated using qMRgenBatch(charmed)
% for publishing documentation.
% Command Line Interface (CLI) is well-suited for automatization
% purposes and Octave.
%
% Please execute this m-file section by section to get familiar with batch
% processing for charmed on CLI.
%
% Demo files are downloaded into charmed_data folder.
%
% Written by: Agah Karakuzu, 2017
% ==============================================================================

1. Print charmed information

qMRinfo('charmed');
 charmed: Composite Hindered and Restricted Model for Diffusion
Pulse Sequence Diagram


Assumptions:
Diffusion gradients are applied perpendicularly to the neuronal fibers.
Neuronal fibers model:
geometry                          cylinders
Orientation dispersion            NO
Permeability                      NO
Diffusion properties:
intra-axonal                      restricted in cylinder with Gaussian
Phase approximation
diffusion coefficient (Dr)       fixed by default. this assumption should have
little impact if the average
propagator is larger than
axonal diameter (sqrt(2*Dr*Delta)>8?m).
extra-axonal                      Gaussian
diffusion coefficient (Dh)       Constant by default. Time dependence (lc)
can be added

Inputs:
DiffusionData       4D DWI
(SigmaNoise)        map of the standard deviation of the noise per voxel. (OPTIONAL)
(Mask)              Binary mask to accelerate the fitting. (OPTIONAL)

Outputs:
fr                  Fraction of water in the restricted compartment.
Dh                  Apparent diffusion coefficient of the hindered compartment.
diameter_mean       Mean axonal diameter weighted by the axonal area --> biased toward the larger axons
fixed to 0 --> stick model (recommended if Gmax < 300mT/m).
fcsf                Fraction of water in the CSF compartment. (fixed to 0 by default)
lc                  Length of coherence. If > 0, this parameter models the time dependence
of the hindered diffusion coefficient Dh.
Els Fieremans et al. Neuroimage 2016.
Interpretation is not perfectly known.
Use option "Time-Dependent Models" to get different interpretations.
(fh)                Fraction of water in the hindered compartment, calculated as: 1 - fr - fcsf
(residue)           Fitting residuals

Protocol:
Various bvalues
diffusion gradient direction perpendicular to the fibers

DiffusionData       Array [NbVol x 7]
Gx                Diffusion Gradient x
Gy                Diffusion Gradient y
Gz                Diffusion Gradient z
Gnorm (T/m)         Diffusion gradient magnitude
Delta (s)         Diffusion separation
delta (s)         Diffusion duration
TE (s)            Echo time

Options:
Rician noise bias               Used if no SigmaNoise map is provided.
'Compute Sigma per voxel'     Sigma is estimated by computing the STD across repeated scans.
'fix sigma'                   Use scd_noise_std_estimation to measure noise level. Use 'value' to fix Sigma.
Display Type
'q-value'                     abscissa for plots: q = gamma.delta.G (?m-1)
'b-value'                     abscissa for plots: b = (2.pi.q)^2.(Delta-delta/3) (s/mm2)
S0 normalization
'Use b=0'                     Use b=0 images. In case of variable TE, your dataset requires a b=0 for each TE.
'Single T2 compartment'       In case of variable TE acquisition:
fit single T2 using data acquired at b<1000s/mm2 (assuming Gaussian diffusion))
Time-dependent models
'Burcaw 2015'                 XXX
'Ning MRM 2016'               XXX

Example of command line usage:
Model = charmed;  % Create class from model
Model.Prot.DiffusionData.Mat = txt2mat('Protocol.txt');  % Load protocol
data = struct;  % Create data structure
data.DiffusionData = load_nii_data('DiffusionData.nii.gz');  % Load data
data.Mask=load_nii_data('Mask.nii.gz');  % Load mask
FitResults = FitData(data,Model,1);  % Fit each voxel within mask
FitResultsSave_nii(FitResults,'DiffusionData.nii.gz');  % Save in local folder: FitResults/

For more examples: qMRusage(charmed)

Author: Tanguy Duval, 2016

References:
Please cite the following if you use this module:
Assaf, Y., Basser, P.J., 2005. Composite hindered and restricted model of diffusion (CHARMED) MR imaging of the human brain. Neuroimage 27, 48?58.
In addition to citing the package:
Karakuzu A., Boudreau M., Duval T.,Boshkovski T., Leppert I.R., Cabana J.F.,
Gagnon I., Beliveau P., Pike G.B., Cohen-Adad J., Stikov N. (2020), qMRLab:
Quantitative MRI analysis, under one umbrella doi: 10.21105/joss.02343

Reference page in Doc Center
doc charmed


2. Setting model parameters

2.a. Create charmed object

Model = charmed;

2.b. Set protocol and options

Protocol: MRI acquisition parameters that are accounted for by the respective model.

For example: TE, TR, FA FieldStrength. The assigned protocol values are subjected to a sanity check to ensure that they are in agreement with the data attributes.

Options: Fitting preferences that are left at user's discretion.

For example: linear fit, exponential fit, drop first echo.

2.b.1 Set protocol the CLI way

If you are using Octave, or would like to serialize your operations any without GUI involvement, you can assign protocol directly in CLI:

Gx = [0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071];
% Gx is a vector of [1791X1]
Gy = [0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.0000; 0.0000; 0.0000; 0.0000; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; -0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071; 0.7071];
% Gy is a vector of [1791X1]
Gz = [0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000; 0.0000];
% Gz is a vector of [1791X1]
Gnorm  = [0.0000; 0.0000; 0.0000; 0.0000; 0.8485; 0.8314; 0.8142; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6942; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6942; 0.7114; 0.7285; 0.7457; 0.7628; 0.7800; 0.7971; 0.8142; 0.8314; 0.8485; 0.8485; 0.8314; 0.8142; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6942; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6942; 0.7114; 0.7285; 0.7457; 0.7628; 0.0000; 0.0000; 0.0000; 0.0000; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.7800; 0.7971; 0.8143; 0.8314; 0.8485; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.0000; 0.0000; 0.0000; 0.0000; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.7800; 0.7971; 0.8143; 0.8314; 0.8485; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.0000; 0.0000; 0.0000; 0.0000; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.7800; 0.7971; 0.8143; 0.8314; 0.8485; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.0000; 0.0000; 0.0000; 0.0000; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.7800; 0.7971; 0.8143; 0.8314; 0.8485; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.0000; 0.0000; 0.0000; 0.0000; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.7800; 0.7971; 0.8143; 0.8314; 0.8485; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.0000; 0.0000; 0.0000; 0.0000; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.7800; 0.7971; 0.8143; 0.8314; 0.8485; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.0000; 0.0000; 0.0000; 0.0000; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.7800; 0.7971; 0.8143; 0.8314; 0.8485; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.0000; 0.0000; 0.0000; 0.0000; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628; 0.7800; 0.7971; 0.8143; 0.8314; 0.8485; 0.8485; 0.8314; 0.8143; 0.7971; 0.7800; 0.7628; 0.7457; 0.7285; 0.7114; 0.6943; 0.6771; 0.6600; 0.6428; 0.6257; 0.6085; 0.5914; 0.5743; 0.5571; 0.5400; 0.5228; 0.5057; 0.4885; 0.4714; 0.4543; 0.4371; 0.4200; 0.4028; 0.3857; 0.3685; 0.3514; 0.3343; 0.3171; 0.3000; 0.2828; 0.2657; 0.2486; 0.2314; 0.2143; 0.1971; 0.1800; 0.1628; 0.1457; 0.1286; 0.1114; 0.0943; 0.0771; 0.0600; 0.0429; 0.0257; 0.0086; 0.0086; 0.0257; 0.0429; 0.0600; 0.0771; 0.0943; 0.1114; 0.1286; 0.1457; 0.1628; 0.1800; 0.1971; 0.2143; 0.2314; 0.2486; 0.2657; 0.2828; 0.3000; 0.3171; 0.3343; 0.3514; 0.3685; 0.3857; 0.4028; 0.4200; 0.4371; 0.4543; 0.4714; 0.4885; 0.5057; 0.5228; 0.5400; 0.5571; 0.5743; 0.5914; 0.6085; 0.6257; 0.6428; 0.6600; 0.6771; 0.6943; 0.7114; 0.7285; 0.7457; 0.7628];
% Gnorm (T/m) is a vector of [1791X1]
Delta  = [0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0070; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0120; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0150; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0200; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0250; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0300; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0350; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400; 0.0400];
% Delta (s) is a vector of [1791X1]
delta  = [0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0030; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080; 0.0080];
% delta (s) is a vector of [1791X1]
TE  = [0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0362; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0462; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0473; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0523; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0573; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623; 0.0623];
% TE (s) is a vector of [1791X1]
Model.Prot.DiffusionData.Mat = [ Gx Gy Gz Gnorm  Delta  delta  TE ];

  See the generic notes section below for further information.

2.b.2 Set protocol and options the GUI way

The following command opens a panel to set protocol and options (if GUI is available to the user):

Model = Custom_OptionsGUI(Model);

If available, you need to close this panel for the remaining of the script to proceed.

  Using this panel, you can save qMRLab protocol files that can be used in both interfaces. See the generic notes section below for details.

3. Fit MRI data

3.a. Load input data

This section shows how you can load data into a(n) charmed object.

  • At the CLI level, qMRLab accepts structs containing (double) data in the fields named in accordance with a qMRLab model.

  See the generic notes section below for BIDS compatible wrappers and scalable
      qMRLab workflows.

%          |- charmed object needs 3 data input(s) to be assigned:
%          |-   DiffusionData
%          |-   SigmaNoise
%          |-   Mask

data = struct();
% DiffusionData.nii.gz contains [64    64     1  1791] data.
data.DiffusionData=double(load_nii_data('charmed_data/DiffusionData.nii.gz'));
% Mask.nii.gz contains [64  64] data.
data.Mask=double(load_nii_data('charmed_data/Mask.nii.gz'));

3.b. Execute fitting process

This section will fit the loaded data.

FitResults = FitData(data,Model,0);

Visit the generic notes section below for instructions to accelerate fitting by
      parallelization using ParFitData.

3.c. Display FitResults

You can display the current outputs by:

qMRshowOutput(FitResults,data,Model);

A representative fit curve will be plotted if available.

To render images in this page, we will load the fit results that had been saved before. You can skip the following code block;

% Load FitResults that comes with the example dataset.
FitResults_old = load('FitResults/FitResults.mat');
qMRshowOutput(FitResults_old,data,Model);

3.d. Save fit results

Outputs can be saved as *.nii.(gz) if NIfTI inputs are available:

% Generic function call to save nifti outputs
FitResultsSave_nii(FitResults, 'reference/nifti/file.nii.(gz)');

If not, FitResults.mat file can be saved. This file contains all the outputs as workspace variables:

% Generic function call to save FitResults.mat
FitResultsSave_mat(FitResults);

  FitResults.mat files can be loaded to qMRLab GUI for visualization and ROI
       analyses
.

The section below will be dynamically generated in accordance with the example data format (mat or nii). You can substitute FitResults_old with FitResults if you executed the fitting using example dataset for this model in section 3.b..

FitResultsSave_nii(FitResults_old, 'charmed_data/DiffusionData.nii.gz');

3.e. Re-use or share fit configuration files

qMRLab's fit configuration files (charmed_Demo.qmrlab.mat) store all the options and protocol in relation to the used model and the release version.

  *.qmrlab.mat files can be easily shared with collaborators to allow them fit their own
      data or run simulations using identical option and protocol configurations.

Model.saveObj('my_charmed_config.qmrlab.mat');

4. Simulations

4.a. Single Voxel Curve

Simulates single voxel curves
  1. Analytically generate synthetic MRI data
  2. Add rician noise
  3. Fit and plot the respective curve

      x = struct;
x.fr = 0.5;
x.Dh = 0.7;
x.diameter_mean = 6;
x.fcsf = 0;
x.lc = 0;
x.Dcsf = 3;
x.Dintra = 1.4;
% Set simulation options
Opt.SNR = 50;
% run simulation
figure('Name','Single Voxel Curve Simulation');
FitResult = Model.Sim_Single_Voxel_Curve(x,Opt);

4.b. Sensitivity Analysis

Simulates sensitivity to fitted parameters
  1. Iterate fitting parameters from lower (lb) to upper (ub) bound
  2. Run Sim_Single_Voxel_Curve for Nofruns times
  3. Compute the mean and std across runs

            fr            Dh            diameter_mean fcsf          lc            Dcsf          Dintra
      OptTable.st = [0.5           0.7           6             0             0             3             1.4]; % nominal values
OptTable.fx = [0             1             1             1             1             1             1]; %vary fr...
OptTable.lb = [0             0.3           3             0             0             1             0.3]; %...from 0
OptTable.ub = [1             3             10            1             8             4             3]; %...to 1
% Set simulation options
Opt.SNR = 50;
Opt.Nofrun = 5;
% run simulation
SimResults = Model.Sim_Sensitivity_Analysis(OptTable,Opt);
figure('Name','Sensitivity Analysis');
SimVaryPlot(SimResults, 'fr' ,'fr' );

5. Notes

5.a. Notes specific to charmed

Not provided.

5.b. Generic notes

5.b.1. Batch friendly option and protocol conventions

If you would like to load a desired set of options / protocols programatically, you can use *.qmrlab.mat files. To save a configuration from the protocol panel of charmed, first open the respective panel by running the following command in your MATLAB command window (MATLAB only):

Custom_OptionsGUI(charmed);

In this panel, you can arrange available options and protocols according to your needs, then click the save button to save my_charmed.qmrlab.mat file. This file can be later loaded into a charmed object in batch by:

Model = charmed;
Model = Model.loadObj('my_charmed.qmrlab.mat');

  Model.loadObj('my_charmed.qmrlab.mat') call won't update the fields in the Model object, unless the output is assigned to the object as shown above. This compromise on convenience is to retain Octave CLI compatibility.

If you don't have MATLAB, hence cannot access the GUI, two alternatives are available to populate options:

  1. Use qmrlab/mcrgui:latest Docker image to access GUI. The instructions are available here.
  2. Set options and protocols in CLI:
  • List available option fields using tab completion in Octave's command prompt (or window)
Model = charmed;
Model.option. % click the tab button on your keyboard and list the available fields.
  • Assign the desired field. For example, for a mono_t2 object:
Model = mono_t2;
Model.options.DropFirstEcho = true;
Model.options.OffsetTerm = false;

Some option fields may be mutually exclusive or interdependent. Such cases are handled by the GUI options panel; however, not exposed to the CLI. Therefore, manual CLI options assignments may be challenging for some involved methods such as qmt_spgr or qsm_sb. If above options are not working for you and you cannot infer how to set options solely in batch, please feel free to open an issue in qMRLab and request the protocol file you need.

Similarly, in CLI, you can inspect and assign the protocols:

Model = charmed;
Model.Prot. % click the tab button on your keyboard and list the available fields.

Each protocol field has two subfields of Format and Mat. The first one is a cell indicating the name of the protocol parameter (such as EchoTime (ms)) and the latter one contains the respective values (such as 30 x 1 double array containing EchoTimes).

The default Mat protocol values are set according to the example datasets served via OSF.

5.b.2 Parallelization

Beginning from release 2.5.0, you can accelerate fitting for the voxelwise models using parallelization.

Available in MATLAB only. Requires parallel processing toolbox.

In CLI, you can perform parallel fitting by:

parpool();
FitResults = ParFitData(data,Model);

If a parpool exists, the ParFitData will use it. If not, a new pool will be created using the local profile. By default, ParFitData saves outputs automatically every 5 minutes. You can disable this feature by:

FitResults = ParFitData(data, Model, 'AutosaveEnabled', false);

Alternatively, you can change the autosave interval (min 1 min) by:

FitResults = ParFitData(data,Model,'AutoSaveInterval',10);

If something went wrong during the fitting (e.g. your computer had to be restarted), you can recover the autosaved data by:

FitResults = ParFitData(data,Model,'RecoverDirectory','/ParFitTempResults_*/folder/from/the/previous/session');

GUI users will be prompted a question about whether they would like to use parallelization after clicking the Fit Data button, if the conditions are met. When called from GUI, ParFitData will be run with default options:

  • Save temporary results every 5 minutes or whenever a chunk has finished processing
  • Split data into chunks with a granularity factor of 3
  • Do not remove temporary fit results upon completion

For further information:

help ParFitData

The default parallelization options can be changed in the preferences.json file located at the root qMRLab directory.

6. Citations

qMRLab JOSS article

Karakuzu A., Boudreau M., Duval T.,Boshkovski T., Leppert I.R., Cabana J.F., Gagnon I., Beliveau P., Pike G.B., Cohen-Adad J., Stikov N. (2020), qMRLab: Quantitative MRI analysis, under one umbrella 10.21105/joss.02343

Reference article for charmed

Assaf, Y., & Basser, P.J. (2005). Composite hindered and restricted model of diffusion (CHARMED) MR imaging of the human brain. Neuroimage, 27(1), 48-58. 10.1016/j.neuroimage.2005.03.042


Quantitative MRI, under one umbrella.

| qMRPullseq | qMRFlow | Interactive Tutorials |

NeuroPoly Lab, Montreal, Canada