x=0.01:0.05:0.98;
y=0.99-X;
[X,Y]=meshgrid(x,y);
Z=42*log(X/30797+Y/(3*15*100000/482)+(1-X-Y)/(100000*15));
figure;mesh(X,Y,Z);
figure;contour(X,Y,Z);
