#pragma rtGlobals=1		// Use modern global access method.

function deriv(waveA,waveB)

	wave waveA, waveB;
	duplicate /O waveA deriv_waveA
	duplicate /O waveA avg_waveB
	
	variable i;
	
	for(i=0;i<numpnts(waveA);i=i+1)	// Initialize variables;continue test
		deriv_waveA[i]=(waveA[i+1]-waveA[i])/(waveB[i+1]-waveB[i]);
		avg_waveB[i]=((waveB[i+1]+waveB[i])/2	);								// Condition;update loop variables
	endfor	// Execute body code until continue test is FALSE
	
										
	
end