function y=F(U,gamma);
y(1,:)=U(2,:);
y(2,:)=U(2,:).^2./U(1,:)+(gamma-1).*(U(3,:)-0.5.*U(2,:).^2./U(1,:));
y(3,:)=(U(3,:)+(gamma-1).*(U(3,:)-0.5.*U(2,:).^2./U(1,:))).*U(2,:)./U(1,:);



