Message erreur maylab "attempt SCRIPT for as

excusez c MATLAB et le message 'attempt execute SCRIPT as a function"

je suis débutant, 59 ans.Comment réussir à faire tourner ce programme? il affiche "attempt execute SCRIPT for as a function"il bloque ligne 12.Comment aussi éviter l'affichage des résultats intermédiaires de la boucle ? puis je enlever les (n,l,m) à w1,w2,w3? Comment faire un graphique final 3d coloré des valeurs intermédiaires de min b et deltar en fonction de wi w2w3 ou 3 graphiques 2d en prenant w1,w2,w3 2 a 2 et impripmer les résultats ? merci d'avance
a=2;ww1=100;ww2=15/1000;ww3=4;
bb1=1;bb2=1;bb3=1;
tt1=1;tt2=1;tt3=1;
delta=10
Y1=a+bb1*sin(-3*ww1+tt1)+bb2*sin(-3*ww2+tt2)+bb3*sin(-3*ww3+tt3)
Y2=a+bb1*sin(-2*ww1+tt1)+bb2*sin(-2*ww2+tt2)+bb3*sin(-2*ww3+tt3)
Y3=a+bb1*sin(-1*ww1+tt1)+bb2*sin(-1*ww2+tt2)+bb3*sin(-1*ww3+tt3)
Y4=a+bb1*sin(0*ww1+tt1)+bb2*sin(0*ww2+tt2)+bb3*sin(0*ww3+tt3)
Y5=a+bb1*sin(1*ww1+tt1)+bb2*sin(1*ww2+tt2)+bb3*sin(1*ww3+tt3)
Y6=a+bb1*sin(2*ww1+tt1)+bb2*sin(2*ww2+tt2)+bb3*sin(2*ww3+tt3)
Y7=a+bb1*sin(3*ww1+tt1)+bb2*sin(3*ww2+tt2)+bb3*sin(3*ww3+tt3)
FOR n = 1:100:20000
FOR l = 1:100:20000
FOR m = 1:100:20000
w1(n)=n/100;cc=cos(w1(n));ss1=sin(w1(n));
w2(l)=l/100;cc2=cos(w2(l));ss2=sin(w2(l));
w3(m)=m/100;cc3=cos(w3(m));ss3=sin(w3(m));
D=[ss1*(-3+4*ss1*ss1),ss2*(-3+4*ss2*ss2),ss3*(-3+4*ss3*ss3),cc*(4*cc*cc-3)-1,cc2*(4*cc2*cc2-3)-1,cc3*(4*cc3*cc3-3)-1;-2*ss1*cc,-2*ss2*cc2,-2*ss3*cc3,2*cc*cc-2,2*cc2*cc2-2,2*cc3*cc3-2;-ss1,-ss2,-ss3,cc-1,cc2-1,cc3-1;ss1,ss2,ss3,cc-1,cc2-1,cc3-1;2*ss1*cc,2*ss2*cc2,2*ss3*cc3,2*cc*cc-2,2*cc2*cc2-2,2*cc3*cc3-2;ss1*(3-4*ss1*ss1),ss2*(3-4*ss2*ss2),ss3*(3-4*ss3*ss3),cc*(4*cc*cc-3)-1,cc2*(4*cc2*cc2-3)-1,cc3*(4*cc3*cc3-3)-1];
DD=det(D);
T=[y1-y4,ss2*(-3+4*ss2*ss2),ss3*(-3+4*ss3*ss3),cc*(4*cc*cc-3)-1,cc2*(4*cc2*cc2-3)-1,cc3*(4*cc3*cc3-3)-1;y2-y4,-2*ss2*cc2,-2*ss3*cc3,2*cc*cc-2,2*cc2*cc2-2,2*cc3*cc3-2;y3-y4,-ss2,-ss3,cc-1,cc2-1,cc3-1;y5-y4,ss2,ss3,cc-1,cc2-1,cc3-1;y6-y4,2*ss2*cc2,2*ss3*cc3,2*cc*cc-2,2*cc2*cc2-2,2*cc3*cc3-2;y7-y6,ss2*(3-4*ss2*ss2),ss3*(3-4*ss3*ss3),cc*(4*cc*cc-3)-1,cc2*(4*cc2*cc2-3)-1,cc3*(4*cc3*cc3-3)-1];
TT=det(T);
U=[ss1*(-3+4*ss1*ss1),y1-y4,ss3*(-3+4*ss3*ss3),cc*(4*cc*cc-3)-1,cc2*(4*cc2*cc2-3)-1,cc3*(4*cc3*cc3-3)-1;-2*ss1*cc,y2-y4,-2*ss3*cc3,2*cc*cc-2,2*cc2*cc2-2,2*cc3*cc3-2;-ss1,y3-y4,-ss3,cc-1,cc2-1,cc3-1;ss1,y5-y4,ss3,cc-1,cc2-1,cc3-1;2*ss1*cc,y6-y4,2*ss3*cc3,2*cc*cc-2,2*cc2*cc2-2,2*cc3*cc3-2;ss1*(3-4*ss1*ss1),y7-y4,ss3*(3-4*ss3*ss3),cc*(4*cc*cc-3)-1,cc2*(4*cc2*cc2-3)-1,cc3*(4*cc3*cc3-3)-1];
UU=det(U);
V=[ss1*(-3+4*ss1*ss1),ss2*(-3+4*ss2*ss2),y1-y4,cc*(4*cc*cc-3)-1,cc2*(4*cc2*cc2-3)-1,cc3*(4*cc3*cc3-3)-1;-2*ss1*cc,-2*ss2*cc2,y2-y4,2*cc*cc-2,2*cc2*cc2-2,2*cc3*cc3-2;-ss1,-ss2,y3-y4,cc-1,cc2-1,cc3-1;ss1,ss2,y5-y4,cc-1,cc2-1,cc3-1;2*ss1*cc,2*ss2*cc2,y6-y4,2*cc*cc-2,2*cc2*cc2-2,2*cc3*cc3-2;ss1*(3-4*ss1*ss1),ss2*(3-4*ss2*ss2),y7-y4,cc*(4*cc*cc-3)-1,cc2*(4*cc2*cc2-3)-1,cc3*(4*cc3*cc3-3)-1];
VV=det(V);
X=[ss1*(-3+4*ss1*ss1),ss2*(-3+4*ss2*ss2),ss3*(-3+4*ss3*ss3),y1-y4,cc2*(4*cc2*cc2-3)-1,cc3*(4*cc3*cc3-3)-1;-2*ss1*cc,-2*ss2*cc2,-2*ss3*cc3,y2-y4,2*cc2*cc2-2,2*cc3*cc3-2;-ss1,-ss2,-ss3,y3-y4,cc2-1,cc3-1;ss1,ss2,ss3,y5-y4,cc2-1,cc3-1;2*ss1*cc,2*ss2*cc2,2*ss3*cc3,y6-y4,2*cc2*cc2-2,2*cc3*cc3-2;ss1*(3-4*ss1*ss1),ss2*(3-4*ss2*ss2),ss3*(3-4*ss3*ss3),y7-y4,cc2*(4*cc2*cc2-3)-1,cc3*(4*cc3*cc3-3)-1];
XX=det(X);
S=[ss1*(-3+4*ss1*ss1),ss2*(-3+4*ss2*ss2),ss3*(-3+4*ss3*ss3),cc*(4*cc*cc-3)-1,y1-y4,cc3*(4*cc3*cc3-3)-1;-2*ss1*cc,-2*ss2*cc2,-2*ss3*cc3,2*cc*cc-2,y2-y4,2*cc3*cc3-2;-ss1,-ss2,-ss3,cc-1,y3-y4,cc3-1;ss1,ss2,ss3,cc-1,y5-y4,cc3-1;2*ss1*cc,2*ss2*cc2,2*ss3*cc3,2*cc*cc-2,y6-y4,2*cc3*cc3-2;ss1*(3-4*ss1*ss1),ss2*(3-4*ss2*ss2),ss3*(3-4*ss3*ss3),cc*(4*cc*cc-3)-1,y7-y4,cc3*(4*cc3*cc3-3)-1];
SS=det(S);
Z=[ss1*(-3+4*ss1*ss1),ss2*(-3+4*ss2*ss2),ss3*(-3+4*ss3*ss3),cc*(4*cc*cc-3)-1,cc2*(4*cc2*cc2-3)-1,y1-y4;-2*ss1*cc,-2*ss2*cc2,-2*ss3*cc3,2*cc*cc-2,2*cc2*cc2-2,y2-y4;-ss1,-ss2,-ss3,cc-1,cc2-1,y3-y4;ss1,ss2,ss3,cc-1,cc2-1,y5-y4;2*ss1*cc,2*ss2*cc2,2*ss3*cc3,2*cc*cc-2,2*cc2*cc2-2,y6-y4;ss1*(3-4*ss1*ss1),ss2*(3-4*ss2*ss2),ss3*(3-4*ss3*ss3),cc*(4*cc*cc-3)-1,cc2*(4*cc2*cc2-3)-1,y7-y4];
ZZ=det(Z);
B1r=SQRT((TT*TT)+(XX*XX))/abs(DD);C1r=ATAN(XX/TT);
B2r=SQRT((UU*UU)+(SS*SS))/abs(DD);C2r=ATAN(SS/UU);
B3r=SQRT((VV*VV)+(ZZ*ZZ))/abs(DD);C3r=ATAN(ZZ/VV);
Ar=Y4-B1r*SIN(C1r)-B2r*SIN(C2r)-B3r*SIN(C3r);
YF1=B1r*sin(-3*w1(n)+c1r)+B2r*sin(-3*w2(l)+c2r)+B3r*sin(-3*w3(m)+c3r);
YF2=B1r*sin(-2*w1(n)+c1r)+B2r*sin(-2*w2(l)+c2r)+B3r*sin(-2*w3(m)+c3r);
YF3=B1r*sin(-1*w1(n)+c1r)+B2r*sin(-1*w2(l)+c2r)+B3r*sin(-1*w3(m)+c3r);
YF4=B1r*sin(c1r)+B2r*sin(c2r)+B3r*sin(c3r);
YF5=B1r*sin(+1*w1(n)+c1r)+B2r*sin(+1*w2(l)+c2r)+B3r*sin(+1*w3(m)+c3r);
YF6=B1r*sin(+2*w1(n)+c1r)+B2r*sin(+2*w2(l)+c2r)+B3r*sin(+2*w3(m)+c3r);
YF7=B1r*sin(+3*w1(n)+c1r)+B2r*sin(+3*w2(l)+c2r)+B3r*sin(+3*w3(m)+c3r);
deltar =abs(Y1-YF1-Ar)+abs(Y2-YF2-Ar)+abs(Y3-YF3-Ar);
deltar=deltar+abs(Y4-YF4-Ar)+abs(Y5-YF5-Ar);
deltar=deltar+abs(Y6-YF6-Ar)+abs(Y7-YF7-Ar);
IF deltar+B1r+B2r+B3r<delta
delta=deltar+B1r+B2r+B3r;
minb=B1r+B2r+B3r;
Am=Ar;
B1m=B1r;w1m=w1;T1m=C1r;
B2m=B2r;w2m=w2;T2m=C2r;
B3m=B3r;w3m=w3;T3m=C3r;
END
END
END
END
YF1=B1m*sin(-3*w1m+T1m)+B2m*sin(-3*w2m+T2m)+B3m*sin(-3*w3m+T3m);
YF2=B1m*sin(-2*w1m+T1m)+B2m*sin(-2*w2m+T2m)+B3m*sin(-2*w3m+T3m);
YF3=B1m*sin(-1*w1m+T1m)+B2m*sin(-1*w2m+T2m)+B3m*sin(-1*w3m+T3m);
YF4=B1m*sin(T1m)+B2m*sin(T2m)+B3m*sin(T3m);
YF5=B1m*sin(+1*w1m+T1m)+B2m*sin(+1*w2m+T2m)+B3m*sin(+1*w3m+T3m);
YF6=B1m*sin(+2*w1m+T1m)+B2m*sin(+2*w2m+T2m)+B3m*sin(+2*w3m+T3m);
YF7=B1m*sin(+3*w1m+T1m)+B2m*sin(+3*w2m+T2m)+B3m*sin(+3*w3m+T3m);
delta =abs(Y1-YF1-Am)+abs(Y2-YF2-Am)+abs(Y3-YF3-Am);
delta=delta+abs(Y4-YF4-Am)+abs(Y5-YF5-Am);
delta=delta+abs(Y6-YF6-Am)+abs(Y7-YF7-Am);
Amm= Am;
B1mm=B1m;W1mm=W1m;T1mm=T1m;
B2mm=B2m;W2mm=W2m;T2mm=T2m;
B3mm=B3m;W3mm=W3m;T3mm=T3m;
deltafin=delta;minbfin=minnb;
comment réussir à faire tourner ce programme?