demo/eaton1.dem

   1 /* solution of simultaneous differential equations */
   2
   3 eq1:3*'diff(f(x),x,2)-2*'diff(g(x),x)=sin(x);
   4 eq2:a*'diff(g(x),x,2)+'diff(f(x),x)=a*cos(x);
   5 atvalue(g(x),x=0,1)$
   6 atvalue('diff(f(x),x),x=0,0)$
   7 atvalue('diff(g(x),x),x=0,1)$
   8 laplace_eq1:laplace(eq1,x,s);
   9 laplace_eq2:laplace(eq2,x,s);
  10 solve([laplace_eq1,laplace_eq2],
  11          ['laplace(f(x),x,s),'laplace(g(x),x,s)]);
  12 ilt(%,s,x);