//www.basicairdata.eu Example Uncertainty calculation
//Scilab 5.4.0
//(c)2013 JLJ  

vb=1;
ib=1;
ubv=0.1;
ubi=0.1;
ur2=((1/ib)^2*ubv^2+(-vb*ib^-2)^2*ubi^2)^0.5
printf("Total uncertainty symbolic manual derivation %f\n",ur2)

function y=F(x)
  y=[x(1)/x(2)];
endfunction
x=[1;1];
sd=[0.1;0.1]; 
[J]=derivative(F,x)
sd=(J^2*sd^2)^0.5
printf("Total uncertainty numeric derivation %f\n",sd) 
