Ok... first of all, one should not post textbook exercises in this section. Anyway, i will show how the lambertW funtion appears from the second DE:
[tex](x^2-y)y'+x=0[/tex]
Using the substituition:
[tex]u=x^2-y[/tex]
we get:
[tex]y=x^2-u[/tex] and [tex]y'=2x-u'[/tex].
So, the DE becomes:
[tex]u(2x-u')+x=0[/tex]
Isolating u', we get:
[tex]u'=\frac{(2u+1)x}{u}[/tex]
Writing u' as du/dx and separating the variables:
[tex]\left[ \frac{u}{2u+1} \right] du = x dx[/tex]
Integrating both sides gives:
[tex]\int \frac{u du}{2u+1} = \frac{x^2}{2}+c[/tex]
The substituition v=2u+1 is useful, and gives:
[tex]\int \frac{1}{4} \left[ 1-\frac{1}{v} \right] dv = \frac{x^2}{2}+c[/tex]
Which, integrated, leads to:
[tex]v - ln(v)=2x^2+4c[/tex]
Ok... now we have to isolate v in this equation. And this is where the lambertW function appears. The http://en.wikipedia.org/wiki/Lambert_W_function" function W(x) is defined as the function that satisfies this equation:
[tex]x=W(x)e^{W(x)}[/tex]
Taking logarithms, we get:
[tex]ln(x) = ln[W(x)] + W(x)[/tex]
And this resembles very much the function v (which is actually v(x)) we got when solving the DE. I can't say that v(x)=W(x) because we have the - sign instead the + sign. Anyway, I am pretty sure the function v(x) is related to the lambertW function, as you mentioned that MATLAB gave it as a possible solution.
I didnt try to solve the first equation. But probably there is some substituition like the one i did that can make it a bit easier to solve.