hold on
subplot(3,1,1)
Uds=[2.40 2.52 2.59 2.71 2.93 3.12];
k=0.34722;
Id=(k/2)*((Uds-2.1611).^2);
plot(Uds, Id)
xlim([1.5,3.5]);
grid on