RMS error

n=500;
for j=1:n;
    rmsn(j)=sqrt((vibration_max_ansys(j)-Deformation_uy(j))^2);
end
% plot(t,Deformation_uy,t,vibration_max_ansys)

rms=sum(rmsn)/n;