function F = parameterfun(x,hwein,cpWf,mL,CPr,lambda,volkoeff,g,rhoL,Tm,Tumg,L,eta,Pr,mW,hLein,K,a,V,Aaussen)
F = [mW*(hwein-cpWf*x(2)) + mL*(hLein-(0.00585*x(1)^3-0.497*x(1)^2+19.87*x(1)-207.61))...
    -CPr*lambda*(volkoeff*g*rhoL^2*(Tm-Tumg)*L^3*Pr/(eta^2))^(1/4)*Aaussen*(Tm-Tumg)/L ,...
    mL*((0.00585*x(1)^3-0.497*x(1)^2+19.87*x(1)-207.61)-hLein)-...
    K*a*V*(hwein-(0.00585*x(1)^3-0.497*x(1)^2+19.87*x(1)-207.61)-(cpWf*x(2)-hLein))/...
    (log((hwein-(0.00585*x(1)^3-0.497*x(1)^2+19.87*x(1)-207.61))/(cpWf*x(2)-hLein)))];
end