clear all
clc
syms k1 k2 M
w0 = sqrt((k1+k2)/M);

dk1 = diff(w0,k1);
dk2 = diff(w0,k2);
dM = diff(w0, M);

ek1 = 0.0001;
ek2 = 0.0001;
eM = 0.0001;

err_sym = sqrt((dk1*ek1)^2 + (dk2*ek2)^2 + (dM*eM)^2);
err_theory_w0 = vpa( subs(err_sym,{k1, k2, M}, {2.7279, 2.6333, 0.2082+(11.8+17.54)*10^(-3)}) )
theory_w0 = vpa(subs(w0, {k1, k2, M}, {2.7279, 2.6333, 0.2082+(11.8+17.54)*10^(-3)}))