function dy = glaze_97(z,y)

global phi_av phi_ad r_ab c_ab rho_ad rho_ab rho_av 

dy=zeros(5,1);

m = y(1)+y(2)+y(3)+m_s; % gesamtmasse in column
c_b = (y(1)*c_d+y(2)*c_v+y(3)*c_l+m_s*c_s)/m; % bulk specific heat


dy(1)=2*al*rho_ad*phi_ad*(a1*y(4)*c1)^0.5;
dy(2)=2*al*rho_av*phi_av*(a1*y(4)*c1)^0.5-omega*y(2)*d1;
dy(3)=omega*y(2)*d1;
dy(4)=e1*(rho_ab*(f1)-m^2);
dy(5)=2*al*rho_ab*c_ab*thetaa*(a1*y(4)*g1)^0.5+l*omega*y(2)*d1...
    -rho_ab*g*((h1)-(i1+j1));


end