eq:'diff(y,x)-y=0;

ode2(eq,y,x);

sol1: ic1(%,x=1,y=4);

set_plot_option(['plot_format, 'xmaxima])$

plot2d(rhs(sol1),[x,-2,10],[y,-2,10]);

kill(all);

depends(u,x);

pde(u):=u*'diff(u,x)+'diff(u,x)^2=0$ pde(u);

    Solutions

u1(x,t):=c;

u2(x,t):=%e^(-x)*c;

u1(0,t)=sqrt(t);

u2(0,t)=sqrt(t);

u2(x,t):=%e^(-x)*sqrt(t);

u2(1,t)=0;

solve(%,t);

pde(%e^(-x)*sqrt(t));

ratsimp(%);

ratsimp(%), diff;