facit12.mws

>    with(DEtools);

[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...
[DEnormal, DEplot, DEplot3d, DEplot_polygon, DFactor, DFactorLCLM, DFactorsols, Dchangevar, GCRD, LCLM, MeijerGsols, PDEchangecoords, RiemannPsols, Xchange, Xcommutator, Xgauge, abelsol, adjoint, auton...

Opg 1

>    ode1:=diff(y(x),x,x) + diff(y(x),x) - 6*y(x) = exp(2*x);

ode1 := diff(y(x),`$`(x,2))+diff(y(x),x)-6*y(x) = exp(2*x)

>    ans1:=dsolve(ode1);

ans1 := y(x) = exp(2*x)*_C2+exp(-3*x)*_C1+1/5*x*exp(2*x)

Opg 2

>    ode2:=diff(y(x),x,x) + 2*diff(y(x),x) + 2*y(x) = 2;

ode2 := diff(y(x),`$`(x,2))+2*diff(y(x),x)+2*y(x) = 2

>    ans2:=dsolve(ode2);

ans2 := y(x) = exp(-x)*sin(x)*_C2+exp(-x)*cos(x)*_C1+1

Opg 3

>    ode3:=diff(y(x),x,x) + y(x) = cos(x);

ode3 := diff(y(x),`$`(x,2))+y(x) = cos(x)

>    ans3:=dsolve(ode3);

ans3 := y(x) = sin(x)*_C2+cos(x)*_C1+1/2*cos(x)+1/2*sin(x)*x

Opg 4

>    ode4:=diff(y(x),x,x)+9*y(x)=80*cos(5*x);

ode4 := diff(y(x),`$`(x,2))+9*y(x) = 80*cos(5*x)

>    ans4:=dsolve({ode4,y(0)=0,D(y)(0)=0});

ans4 := y(x) = 5*cos(3*x)-5*cos(5*x)

Opg 5

>    ode5:=diff(y(x),x,x)-4*diff(y(x),x)+4*y(x)=0;

ode5 := diff(y(x),`$`(x,2))-4*diff(y(x),x)+4*y(x) = 0

>    ans5a:=dsolve(ode5);

ans5a := y(x) = _C1*exp(2*x)+_C2*exp(2*x)*x

>    ans5b:=dsolve({ode5,y(0)=0,D(y)(0)=2});

ans5b := y(x) = 2*exp(2*x)*x

Opg 6

>    ode6:=diff(y(x),x,x)-2*diff(y(x),x)+y(x)=x;

ode6 := diff(y(x),`$`(x,2))-2*diff(y(x),x)+y(x) = x

>    dsolve(ode6);

y(x) = exp(x)*_C2+exp(x)*x*_C1+2+x

Opg 7

>    ode7:=diff(y(x),x)=3*x^2*exp(-y(x));

ode7 := diff(y(x),x) = 3*x^2*exp(-y(x))

>    odeadvisor(ode7);

[_separable]

>    dsolve(ode7, [separable]);

y(x) = ln(x^3+3*_C1)

Opg 8

>    ode8:=(diff(y(x),x))=(x+1)/(y(x)+1);

ode8 := diff(y(x),x) = (x+1)/(y(x)+1)

>    odeadvisor(ode8);

[_separable]

>    ans8a:=dsolve({ode8, y(1)=-3});

ans8a := y(x) = -1-(1+x^2+2*x)^(1/2)

>    simplify(%);

y(x) = -1-csgn(x+1)*x-csgn(x+1)

>    ans8b:=dsolve(ode8, implicit);

ans8b := 1/2*x^2+x-1/2*y(x)^2-y(x)+_C1 = 0

>    Int(x*sin(x^2),x=0..1/2*sqrt(Pi)) =int(x*sin(x^2),x=0..1/2*sqrt(Pi));

Int(x*sin(x^2),x = 0 .. 1/2*Pi^(1/2)) = -1/4*2^(1/2)+1/2

>    Int(x*sin(x^2),x) = int(x*sin(x^2),x);

Int(x*sin(x^2),x) = -1/2*cos(x^2)

>