wenn ich mich da einmischen darf:
Du gehst von einer solchen PDGL aus:
[align=center]a(x,y)Zx+b(x,y)Zy+c(x,y)Z = d(x,y)
[/align]
Der Ansatz ist die charakteristische Gleichung:
[align=center]
y' = b(x,y) / a(x,y) oder x' = a(x,y) / b(x,y)
[/align]
Damit erhälst du eine gewöhnliche DGL, die es einfach zu lösen gilt.