Question 665012
Obvious constraints are:
{{{x>=0}}} and {{{y>=0}}}
(no negative production)
Given constraints are:
{{{y-x<=200}}}
({{{y}}} cannot exceed {{{x}}} by more than {{{200}}}  ,
the difference has to be {{{200}}} or less),
and
{{{x+2y<=1600}}}.
 
Because of those constraints,
there is a feasibility region.
You only can work in that region of the x-y plane.
That region is bordered by the lines represented by
{{{x=0}}} , {{{y=0}}} , {{{y-x=200}}} , and {{{x+2y=1600}}}.
You can graph the lines.
{{{graph(300,300,-200,1800,-100,900,200+x,800-0.5x)}}}
{{{y-x=200}}} graphs as the red slanted line ,
and {{{x+2y=1600}}} graphs as the green line.
Your feasibility region is tue quadrilateral with 
parts of the x-axis, the y-axis and the green and red lines for sides.
You can find the intersection points for each pair of lines.
For example, solving {{{system(y-x=200,x+2y=1600)}}}
gives you the solution {{{x=400}}} with {{{y=600}}}
for point (400,600), where the red and green slanted lines intersect.
 
The vertices of your feasibility region are:
(0,0) , (0,200) , (400,600) and (1600,0).
The maximum for {{{P}}} will happen at 1 of those points.
(In some cases it could happen at 2 of vertices and the whole segment connecting them).
All you need to do is calculate {{{P}}} for each of those 4 points.
I will show you the calculation for 2 of them:
For point (400,600), with {{{x=400}}} and {{{y=600}}},
{{{P=14*400+22*600-900=17900}}}.
For point (1600,0), {{{P=14*1600+22*0-900=21500}}}.
The other points give you smaller values for {{{p}}},
so the solution is point (1600,0),
with {{{highlight(x=1600)}}} and {{{highlight(y=0)}}}.
 
THE REASON WHY IT WORKS THAT WAY:
The function to maximize, {{{P=14x+22y-200}}}
is a function of x and y,
which could be represented in 3 dimensions,
wit {{{z=14x+22y-200}}} being the third, dependent variable.
As with altitude as a function of 2-dimensional coordinates (latitude and longitude),
we can represent the function {{{P(x,y)=14x+22y-900}}} on paper by making a contour map.
The contour lines would be {{{14x+22y-200=K for various values of the constant {{{K}}}.
Luckily for us, the function is linear in {{{x}}} and {{{y}}},
so those contour lines will be straight lines, like the blue line below.
{{{graph(300,300,-200,1800,-100,900,200+x,800-0.5x,(10000-14x)/22)}}}
The blue line is the graph of
{{{14x+22y=10000}}} <---> {{{14x+22y-900=9100}}}
That is the line for {{{P=9100}}}.
As you change the constant, the line changes,
but all the other contour lines are parallel to that blue line.
As increase the value for the constant the lines moves away from the origin, until it moves out of your feasibility region.
You want the values (for {{{x}}}, {{{y}}} and {{{P}}} for the largest possible {{{P}}},
when you reach the end of the feasibility region.
In general, that will happen at one of the vertices, or at 2 of the vertices and the side that joins them.
In this case it happens at point (1600,0).
The maximum for P will be found at (1600,0)
where {{{P=14*1600-22*0-900=21500}}}