si per gli ei ho capito che devo calcolarli dinamicamente, cioè in base al valore k che trovo, però io mi riferivo agli Ri. Ti faccio un esempio pratico:
Con il seed, m e a trovo Xo e da qui calcolo R0 e quindi K0. Immaginiamo che K0 = 5; tu ora qui mi consiglieresti di calcolarli solo ora gli altri 4 Ri necessari a calcolare l'erlangiano?
Questo intendevi?

Io ho fatto così...

Poco sopra spiegavo il mio approccio... lo so...non è che sia proprio chiarissimo quello che ho scritto :?


Però, in soldoni, potresti riuscire ad allocare lo spazio per un array solo per i K... e, se ne hai bisogno, degli Eri (la sommatoria degli Ei generati...) e, posso assicurarti, che funziona. (Provato pure a verificare con un foglio excel...)
Tutto sta a dimensionare correttamente l'array dei k (alle brutte, scrivi 3000 e non ci pensi più)...

si per gli ei ho capito che devo calcolarli dinamicamente, cioè in base al valore k che trovo, però io mi riferivo agli Ri. Ti faccio un esempio pratico:
Con il seed, m e a trovo Xo e da qui calcolo R0 e quindi K0. Immaginiamo che K0 = 5; tu ora qui mi consiglieresti di calcolarli solo ora gli altri 4 Ri necessari a calcolare l'erlangiano?
Questo intendevi?

scusa snoopy non riesco a cogliere l'accezione dinamica? praticamente che mi consigli di fare?

Di partire con k0 e generarti solamente le Expi di cui hai bisogno per calcolare T0, salvandoti in una lista i valori di ki ad essi relativi (che utilizzerai poi per T1, T2, ...)... e così via...

scusa snoopy non riesco a cogliere l'accezione dinamica? praticamente che mi consigli di fare?

ragazzi allora a nessuno è venuto in mente come evitare di allocare un array di 150*20 per gli Ri?

Si. Creando una lista, ad esempio. Cmq ti conviene generarti le exp dinamicamente e poi buttarle, salvandoti invece i vari k (risparmiando un sacco di memoria...).
Rinnovo l'appello: nessuno ha già riscritto le f1 e f2 per Mathematica (MUSA)...?? Non capisco dov'è che sbaglio...

ragazzi allora a nessuno è venuto in mente come evitare di allocare un array di 150*20 per gli Ri?

mi confermate la consistenza di questa formula:

ei = - 127/k * ln(Ri)

ho copiato male alla lavagna?

te la accendo

mi confermate la consistenza di questa formula:

ei = - 127/k * ln(Ri)

ho copiato male alla lavagna?

Abbiamo un problema con Mathematica: come avete specificato le funzioni f1 e f2??? Continua a darci errori sul determinante della matrice...

X pax22: a quali valori di n la media di Ti converge a 127??

Hai provato a dargli in pasto le funzioni già derivate?

Che modello vi è capitato, J&M o MUSA?

MUSA. Abbiamo provato con entrambe, ma niente...
Qualcuno sarebbe così gentile da postarmi la parte iniziale del listato (dove sono definite f1 e f2)???

Abbiamo un problema con Mathematica: come avete specificato le funzioni f1 e f2??? Continua a darci errori sul determinante della matrice...

X pax22: a quali valori di n la media di Ti converge a 127??

Hai provato a dargli in pasto le funzioni già derivate?

Che modello vi è capitato, J&M o MUSA?

A me non converge manco per sbaglio ...

Gli esponenziali che genero, sono sicuramente giusti ... tanto che, se non ci fosse la distribuzione K-erlangiana degli istanti di failure, ma fosse la solita esponenziale, sarei a cavallo... perché Phi0 viene 0.0006, come nei progetti passati....


Ora io fagianamente ho dichiarato tutti i valori del progetto come costanti...
[code type="markup"]#define A 1220703125UL
#define M 2147483648UL
...[/code] etc etc

Poi, per pura pigrizia, anziché fare il vettore da 3000 degli esponenziali... parto da k0, vedo quanto vale, genero i k0 esponenziali, tenendo traccia sia dei ki che ho calcolato nel frattempo, inserendoli in un vettore, sia dell'ultimo xi generato...
così, incremento il puntatore sul vettore dei k, vedo quanti esponenziali devo generare e riparto da quell'ultimo xi che avevo salvato (che mi indica l'ultimo esponenziale utilizzato per ricavare Ti)...

Casino, eh?

Boh, fatto sta... che anche facendoli a mano, i Ti sono giusti ... ma questa cacchio di media non converge :cry:

Abbiamo un problema con Mathematica: come avete specificato le funzioni f1 e f2??? Continua a darci errori sul determinante della matrice...

X pax22: a quali valori di n la media di Ti converge a 127??

per gli int...

Ma m=2^31

x pax: voi state usando solo interi? Noi lasciando i reali e approssimando per eccesso o difetto a seconda dei casi (se parte decimale > 0.5 per eccesso e viceversa) otteniamo una media pari a 115,...

Ci pare un valore sensato considerando tutte le varie approsimazioni. Poi ho provato ad aumentare il valore di numeri da generare(ad esempio 550) e mi sono avvicinato a 127 ottenendo 133. Magari con poche generazioni si rischia di approssimare molte più volte per difetto o per eccesso, mentre aumentando i numeri il tipo di approsimazione di bilancia e ci si avvicina di più a 127.

Intanto sto provando in mathematica, ho ottenuto dei valori ma mi sembrano un pò sballati....

si ma m = 2^32 quindi un intero non basta...

Potresti usare la classe BigInteger e vai più che tranquillo... 8)

si ma m = 2^32 quindi un intero non basta...

Per quanto riguarda il mio gruppo, lo stiamo facendo in C ....
Penso, cmq, che puoi usare i tipi primitivi...
A me int e float sono bastati ...

il progetto lo state facendo in java? per m usate la classe Long?

:shock:
Mi verrebbe un Fi(0) nell'ordine di 10^-5.... mi sa che non va bene ... devo aver sbagliato qualcosa... :? :cry:

anche a me Fi(0) ha quell'ordine di grandezza.