Co=100;
k=0.05;
C=Co*exp(-k*time);
sfalma=1.*observed-C;
s=sfalma./C*100;
disp(s)
subplot(2,2,1)
plot(time,observed,'bo',time,C,'c-')
subplot(2,2,2)
plot(C,observed,'g*',observed,observed,'b')
subplot(2,2,3)
plot(time,sfalma,'r',time,zeros(size(time)),'k')
subplot(2,2,4)
plot(time,s,'m',observed,zeros(size(observed))+50,'k')