Lösung der Differentialgleichung:
restart: a(x):= ;# (* Ausdruck eingeben *) b(x):= ;# (* Ausdruck eingeben *) diffgl:=diff(y(x),x)+a(x)*y(x)=b(x): dsolve(diffgl,y(x));
Mit Randbedingung:
c:= ;# (* Randbedingung eingeben *) dsolve({diffgl,y(0)=c},y(x));
Kontrolle:
subs("[1],diffgl); simplify(");