>> x=0:6.5:650;
y=0:27:2700;
[x,y]=meshgrid(x,y);
F=(669126134152931648175758856558983806728427764752864631.*35184372088832.^(1819./5000).*3695207012694701.^(3181./5000).*gamma(3181./5000))./642775217703596110216784836936465041008881197513117134120550400000000 - (6958236045425711425937142107971.*x.^(8181./5000).*exp(-(35184372088832.*x)./3695207012694701))./20769187434139310514121985316880384 - (669126134152931648175758856558983806728427764752864631.*35184372088832.^(1819./5000).*3695207012694701.^(3181./5000).*igamma(3181./5000, (35184372088832.*x)./3695207012694701))./642775217703596110216784836936465041008881197513117134120550400000000 - (210350875244555689461099923470287270269860976030451.*x.^(3181/5000).*exp(-(35184372088832.*x)./3695207012694701))./3653754093327257295509212081790707549139831357440000;
F2=(3926477034605225449330639167562343.*2.^(7551./10000).*5915813241695293.^(4243./1