Здравствуйте! Подскажите, пожалуйста, как можно решить задачу Гурса в Maple? Есть ли какие-нибудь примеры решения? Возможно ли решение с параметрми?
В качестве примера привожу код, который мне дал преподаватель. Что можно с этим сделать?
Код:
a:=1,-6,8,1,-3,0,0;
equ:=a[1]*diff(u(x,y),x,x)+a[2]*diff(u(x,y),x,y)+a[3]*diff(u(x,y),y,y)+a[4]*diff(u(x,y),x)+a[5]*diff(u(x,y),y)+a[6]*u(x,y)+a[7]=0;
eq:=lhs(equ);
A:=linalg[matrix](2,2,[coeff(eq,diff(u(x,y),x,x)),coeff(eq,diff(u(x,y),x,y))/2,coeff(eq,diff(u(x,y),x,y))/2,coeff(eq,diff(u(x,y),y,y))]);
Delta:=simplify(linalg[det](A));
A[1,1]*z^2-2*A[1,2]*z+A[2,2]=0;
res1:=solve(A[1,1]*z^2-2*A[1,2]*z+A[2,2],z);
res2:={seq(dsolve(diff(y(x),x)=res1[i],y(x)),i=1..2)};
res2:=subs(y(x)=y,res2);
{seq(solve(res2[i],_C1),i=1..nops(res2))};
if Delta >0 then itr:={xi=solve(res2[1],_C1)+solve(res2[2],_C1),eta=(solve(res2[1],_C1)-solve(res2[2],_C1))*I}; elif Delta < 0 then itr:={xi=solve(res2[1],_C1),eta=solve(res2[2],_C1)}; elif (Delta =0) and (Diff(xi,y)<>0) then itr:={xi=solve(res2[1],_C1),eta=x};else itr:={xi=solve(res2[1],_C1),eta=y};fi;
tr:=solve(itr,{x,y});PDEtools[dchange](tr,eq,itr,[eta,xi],simplify)=0;