Diciamo che io ora sta cosa la voglio usare come una funzione esterna per calcolare le componenti della velocità.
La funzione accetterebbe in input solo il modulo di v che lo calcolo dal main e mi restituisce v_x, v_y, v_z.
Però io ho 52000 e cocci puntini ma non così tante particelle, allora dò in input anche il numero di particelle.
Poi è inutile che lui usa i puntini più distanti dall'asse z perchè saranno particelle tutte sparate a popo, quindi voglio dirgli anche quali versori usare, e dargli quindi in input lo spread angolare. Questo lo posso fare passandogli un numero, che chiamo z e gli dico di usare i versori che hanno versore(:,3) compresa tra 1 e z che gli passo.
quindi l'intestazione della funzione dovrebbe essere:
function [v_x, v_y, v_z]=Event_Stop_1(n,z, v_ion)
e ancora non abbiamo finito.
a questo punto lui dovrebbe estrarre dei puntini a caso da:
versore_selezionato = find(versore(:,3)<=1 && versore(:,3)>=z)
n puntini a caso e usare le coordinate per calcolare le componenti della velocità.
Quindi devo capire come estrarre n riga a caso da versore_selezionato...