function resid=fun(Put,Gmax,V1,Aut,gamma,P1,Pc,Vc,vc);

% Husk at ligningen har to løsninger, dvs at Put er en vektor med to elementer og 
% ./ og .^ må brukes når Put deles på eller opphøyes
vut     = Gmax*V1*(P1./Put).^(1/gamma)/Aut;
resid   = vut.^2-vc^2-(2*gamma/(gamma-1))*Pc*Vc*(1-(Put/Pc).^((gamma-1)/gamma));

