Esempio numerico pag.92

clear
clc
disp('Si consideri una particella che sia vincolata a muoversi su una retta')
disp('di un piano verticale xy, di equazione 3x+2y=0, sotto l''azione della')
disp('forza di gravità che ha la direzione dell''asse y positivo. Si vuole')
disp('trovare l''equazione del moto della particella')
disp('Condizioni iniziali')
dx=0
dy=0
m=1
disp(' La matrice dei coefficienti della equazione di vincolo è')
syms x y real
d=diff(sym('3*x+2*y'),x);
e=diff(sym('3*x+2*y'),y);
H=[d e];
A=sym2poly(H);
A
disp('-------------- premere un tasto per continuare -------------')
pause;
disp(' La matrice di massa del sistema è')
M=[m 0;0 m]
disp('-------------- premere un tasto per continuare -------------')
pause;
disp(' In virtù del sistema di riferimento adottato, il vettore forza è')
F=[0 9.81]
disp('-------------- premere un tasto per continuare -------------')
pause;
disp(' L''accelerazione del sistema in assenza di vincoli è')
a=[0 9.81]'
disp('-------------- premere un tasto per continuare -------------')
pause;
k=diff(d,x);
l=diff(e,y);
m=diff(d,y);
n=diff(e,x);
disp(' Il termine noto b è invece pari a')
c=-(k*(dx^2)+2*m*dx*dy+l*(dy^2));
b=sym2poly(c)
disp('-------------- premere un tasto per continuare -------------')
pause;
disp(' In questo caso il termine b è uno scalare')
disp(' Determiniamo ora la matrice dei vincoli B,sapendo che B=A*M^(-1/2)')
disp(' La matrice B è pari a')
D=M^(-1/2);
B=A*D
disp('-------------- premere un tasto per continuare -------------')
pause;
disp(' L''inversa generalizzata di Moore-Penrose della matrice B è')
MPB=pinv(B)
disp('-------------- premere un tasto per continuare -------------')
pause;
disp(' L''accelerazione effettiva del sistema è pari a')
ax=a+D*MPB*(b-A*a)

Interessante l'uso della matematica simbolica con Matlab di cui conoscevo l'esistenza ma che non mi è familiare.
Qualcuno ha voglia di confrontare i risultati ottenuti con il listato matlab proposto da campanile con quelli ottenuti, con gli setessi valori dei parametri ovviamente, con maple ?

non riesco a capire come mai il listato caricato mi dia il seguente errore:
Error using sym.svd
Too many input arguments.

Error in pinv (line 29)
[U,S,V] = svd(A,0);

Error in pinv (line 27)
X = pinv(A',varargin{:})';

Error in esercizio_nuovo (line 50)
MPB=pinv(B)