MAXIMAで微分方程式をときグラフを出力する方法です

Question

前提・実現したいこと

数理計算ソフトMAXIMAの使い方で質問です。
以下の微分方程式をとき、グラフをかく必要があります。
i' (t)+3i(t)=10 sin⁡(t)
初期値をｔ＝０、i=-1とするとエラーとなります。
記述方法を教え下さい。
よろしくお願いします。

発生している問題・エラーメッセージ

エラーメッセージ
#0: ic1(soln=i(t) = %e^-(3*t)*(%e^(3*t)*(3*sin(t)-cos(t))+%c),xc=t = 0,yc=i = -1) (ode2.mac line 323)
 -- an error. To debug this try: debugmode(true);

### 該当のソースコード

```ここに言語名を入力
eqn1:'diff(i(t),t)+3*i(t)=10*sin(t);
ode2(eqn1,i(t),t);
sol1:ic1(%,t=0,i=-1);

### 試したこと
plot2dまで行き着きません

### 補足情報（FW/ツールのバージョンなど）

微分方程式もMAXIMAも初めてなのであまりわかっていない中での質問で申し訳ありませんが、よろしくお願いします。

Answer 1

i=-1の個所おかしいかと思います。i(0)=-1として処理すればどうでしょうか。

eqn1:diff(i(t),t)+3*i(t)=10*sin(t);
ode2(eqn1,i(t),t);
/* i → i(0) */
sol1:ic1(%,t=0,i(0)=-1);
/* i(0)が消えないのでi(0)=-1を代入  */
sol1:subst(-1, i(0),sol1);
plot2d(rhs(sol1),[t,0,10])