d = daq.getDevices
s = daq.createSession('ni')
s1= daq.createSession('ni')
s.addAnalogInputChannel('Dev1','ai0', 'Voltage')
s.addAnalogInputChannel('Dev1','ai1', 'Voltage')
s.addAnalogInputChannel('Dev1','ai2', 'Voltage')
s.addAnalogInputChannel('Dev1','ai3', 'Voltage')
AO=s1.addAnalogOutputChannel('Dev1', 'ao0', 'Voltage')
s.Rate=1000;
s.DurationInSeconds =10
v.TerminalConfig = ' Differential';
lh = s.addlistener('DataAvailable', @plotData);
e= 18 - Vref
IE=sum(e);
u=C*er
for i=1:length(u)
t = timer('TimerFcn', ' putsample(AO,u)',...
'StartDelay',0.2);
start(t)
end
s.startBackground()
s.wait()
delete (lh)
s.release()