for i=1:p
%Campi nulli
B=[0 0 0]';
E=[0 0 0]';
% v0 = [0 0 v_H]'; %initial velocity column vector
% r0 = [0 0 0]'; %initial position of particle column vector
% y0 = [r0; v0]; %Concatena i due vettori colonna r0 e v0 e crea il vettore
%colonna delle condizioni iniziali
y0(i,:) = [Part_trasp(i).traiettoria(end,:),Part_trasp(i).velocita(end,:)];
%%%Risolvo il moto
f = @(t,y1) [y1(4:6); (q_over_m).*cross(y1(4:6),B)+(q_over_m).*E];
% the expression [y(4:6); (q/m)*cross(y(4:6),B)]
% combines two vectors of length 3 to make a vector of
% length 6(as long as y is a column vector).
options=odeset('RelTol',1e-7,'Events',@(t,y)Event_Stop_1(t,y,Deriva1)); % opzione per risoluzione equazione differenziale:
%set precisione calcolo---% Sintassi per utilizzare Event_Stop con Deriva1 come variabile in input
[t,y1] = ode23t(f,tspan,y0(i,:),options);
%traiettoria= y1; % aggiorno la traiettoria
end
Prima dell'end vorrei che aggiungesse alla traiettoria tutte le righe delle prime tre colonne di y1 e alla velocità tutte le righe delle altre 3 colonne di y1
gli ho detto:
Part_trasp(i).traiettoria(end+1,1:3)=y1(:,1:3);
e lui si mussia e dice:
??? Subscripted assignment dimension mismatch.
P.S.: mussia = storce il muso, si incazza...