Co=100;
k=0.05;
C=Co*exp(-k*time)
a=C-1.*observed;
sfalma=a./C*100;
disp(sfalma)