ODE45

[Barese1]

Salve ragazzi ho un altro piccolo GRANDE problema in matlab.
Dovrei risolvere il seguente sistema di equazioni con la funzione ODE45:
[r1
[m1 0 0 [x"1 [c1+c2 -c2 0 [x'1 [-k1 k2 -k2 0 x1 [F1(t)
0 m2 0 * x"2 + -c2 c2+c3 -c3 * x'2 + 0 -k2 k2+k3 -k3 * x2 = F2(t)
0 0 m3] x"3] 0 -c3 c3] x'3] 0 0 -k3 k3] x3] F3(t)]
essendo inoltre:
r'1(t)=x'1-beta*ass(x'1)*ass(r1)^(n-1)*r1-gamma*x'1*ass(r1)^n
come già detto in un altro post non sono un fenomeno con matlab (anzi sono proprio un novellino)
avevo pensato di creare le matrici mm cc kk e F, quest'ultima contenente come termini tre funzioni f(1), f(2) ed f(3)
che richiamavo da altri m.file, ma il resto??
come faccio a far capire al programma che il mio vettore incognito x ha 4 componenti mentre la sua derivata prima e
seconda ne hanno 3??
Se qualcuno sa come aiutarmi gli sarei immensamente grato

[Barese1]

scusate ragazzi non so perchè ma quando ho pubblicato il messaggio la spaziatura
è scomparsa rendendo incomprensibile il mio sistema di equazioni,
provo a riscriverlo diversamente:

[m1, 0, 0; 0, m2, 0;0, 0, m3]*[x"1, x"2, x"3] + [c1+c2, -c2, 0;-c2, c2+c3, -c3;0, -c3, c3]*[x'1, x'2, x'3]+[-k1, k2, -k2, 0;0, -k2, k2+k3, -k3;0, 0, -k3, k3]*[r1, x1, x2, x3]=[F1(t), F2(t), F3(t)]

essendo inoltre:
r'1(t)=x'1-beta*ass(x'1)*ass(r1)^(n-1)*r1-gamma*x'1*ass(r1)^n