function [De,Dd,Du] = scaling(z1,z2)
%SCALING Summary of this function goes here
%   Detailed explanation goes here

De=diag([0.05 1 1 1 1 2 20 0.23]);
Dd=diag([0.001 3 3 1]);
Du=diag([min(z1,1-z1) min(z2,1-z2)]);
end

