#include "mex.h"
#include <math.h>

static double addition(double *pulse[], int value)
{
    double summe; 
    int i, m, value_end;
    
    value_end=value+196;
    summe=0.0;
    
  /* Find the dimensions of the data */
    m = sizeof(pulse);
        
    for (i=value; i<=value_end; i++)
    {
        if (value_end<=m){
            summe=summe+((double)*pulse[i]);
            }
        else{
            summe=0.0;
            break;
            }
    }
    return summe;
}

static double difference(double *Temps[], int value)
{
    double diff;
    int m, value_end;
    
    value_end=value+196;
    diff=1.0;
    m = sizeof(Temps);
    
    if (value_end<m){
        diff=(*Temps[value_end])-(*Temps[value]);
        }
    else{
        diff=1.0;
    }
    return diff;
}


void mexFunction(int nlhs, mxArray *plhs[],
int nrhs, const mxArray *prhs[])
{
int i, n;
double data1, data2;
double *data3;

/* Find the dimensions of the data */
n=(int)*mxGetPr(prhs[2]);

/* Create an mxArray for the output data */
plhs[i] = mxCreateDoubleMatrix(n, 1, mxREAL);
/* Retrieve the input data */
/* Create a pointer to the output data */
data3 = mxGetPr(plhs[i]);
/* Put data in the output array */
for (i=0; i<n; i++){
        data1 = (addition(mxGetPr(prhs[0]),i));
        data2 = (difference(mxGetPr(prhs[1]),i));
        data3[i] = (data1/data2);
        }
}