ご覧の通りの数値解法です。ソースはこれだけ、しかも初めの2行は文字列を書くだけ、
xset("font size",4);
xstring(1,4,"微分方程式 dy/dx=x-2y y(0)=2");
deff("ydot=f(x,y)", "ydot=x-2*y") //方程式の定義
x0=0 ; y0=2 ; //初期条件
x=0:0.1:5; //解の範囲
y=ode(y0 , x0 , x , f ); //解
plot2d(x,y,5,rect=[0,0,5,5]) //プロット
2つの方程式をまとめて解いてグラフを描くこともできます。
昨春、真中君（2人）の卒業後、castorも使わなくなってしまったが、Ver５．２用に書き換えています。