function dxdt = nonlinear_restoring_k_mems_var_Uac(t,x)
% =====================================================================
% Differentialgleichung mit nichtlinearer Rueckstellkraft und
% Sinuserregung.
% Die Rueckstellkraft F_restoring(x) wird aus dem Look-up-Table 'Frest'
% durch lineare Interpolation an der Stelle x gewonnen.
% =====================================================================
global delta k0 m xF Frest Beta Udc fexc d0 zahl i j
w0 = sqrt(k0/m); % omega0
%Beta = 0.5*e0*pi*R^2;
Beta =  2.7914e-18;
% --- Differentialgleichung 2. Ordnung umgeformt in System 1. Ordnung
%dxdt = [ x(2) -2.*delta(j).*w0.*x(2) - interp1(xF,Frest,x(1))/m - ap*sin(2*pi*fexc(i)*t)]';
dxdt = [ x(2) -2*delta*x(2) - interp1(xF,Frest,x(1))/m + (1/m).*((Beta*(Udc + zahl(j)*sin(2*pi*fexc(i)*t)).^2)/((d0-x(1))^2))]';
end

