Equazioni differenziali a coefficienti non costanti

di Anonimizzato9724 il
473 risposte

473 Risposte - Pagina 5

  • /
    Anonimizzato9724
    Utente anonimizzato
    Niente, non dice nulla solo "busy"... ma non capisco che cosa è occupato a fare...
  • /
    Anonimizzato9724
    Utente anonimizzato

    giug ha scritto:

    dà errore? che dice?

    Se la domanda rigurda è riferita a

    1keenan ha scritto:

    Per la cronaca: ho scritto il codice per risolvere l'equazione:

    f = @(z,y) [y(4:6); (q/m)* cross( cross( E, 1/y(4:6) ), 1/y(4:6) )+ (q/m)*cross( B, 1/y(4:6) ) ];
    [z,y] = ode23t(f,zspan,y0);

    fa cinque minuti di conti e poi fa una minchiata, quindi penso sia sbagliata l'equazione...
    allora non dà alcun errore, fa solo un plot che non c'entra niente con quello che mi aspetto, diciamo che fa una specie di "L" per questo secondo me è l'equazione che lo ricavata sbagliata.

    Nel frattempo avevo aggiornato il post per non farne un altro...
  • /
    Anonimizzato9528
    Utente anonimizzato
    Quant'è m_H? Altrimenti non riesco a far partire il codice...
  • /
    Anonimizzato9724
    Utente anonimizzato

    giug ha scritto:

    quant'è m_H? Altrimenti non riesco a far partire il codice...

    scusa, ho dimenticato di copiare e incollare le prime 5 righe:
    
clear all
close all

uma=1.66e-27; %[Kg], unità di massa atomica
q=1.602e-19; %[C], carica elementare
m_H = 1.00794*uma; %[kg], massa atomica di H
  • /
    Anonimizzato9528
    Utente anonimizzato
    Forse devi mettere && al posto di ||.
    Perché così gli stai dicendo
    se z>0.05 oppure z<0.13
    che se è sempre verificata...
  • /
    Anonimizzato9724
    Utente anonimizzato
    Così sembra che nemmeno ci entra nel ciclo...
  • /
    Anonimizzato9528
    Utente anonimizzato
    L'ultimo valore prima di entrare nel secondo ciclo è z= 2.3965e+004
    che non è compreso tra 0.05 e 0.13, quindi per forza non entra...
  • /
    Anonimizzato9724
    Utente anonimizzato
    Si, questo l'avevo capito... ma perchè non continua a iterare fino a che z=0.05 invece di uscire dal primo ciclo in cui c'è while z<=0.05?
  • /
    Anonimizzato9528
    Utente anonimizzato
    Perché entra con z=0 si fa i calcoli e gli viene subito quel numerone, quindi esce dal ciclo.
  • /
    Anonimizzato9724
    Utente anonimizzato
    Ma io perchè lo leggevo elevato alla -4?... vabbè sto fondendo.
    è da ieri che sono convinto che è un numero piccolo.
    C'è qualcosa che non funziona prima allora, non ha senso sto valore per quello che ci devo fare.
  • /
    Anonimizzato9724
    Utente anonimizzato
    Basta mettere:
    dt= 1*10^-10; %intervallo di tempo

    per fargli fare qualcosa di sensato
    Ora fa una retta fino a 0.05
    poi fa un semicerchio, che ci sta come risultato.
    Sono impazzito due giorni e mi bastava leggere con più attenzione.
    SONO UN CRETINO
  • /
    Anonimizzato9528
    Utente anonimizzato
    L'importante è arrivare (prima o poi) alla soluzione...
  • /
    Anonimizzato9724
    Utente anonimizzato
    GRANDE GIUG!!!
    GRAZIE GRAZIE GRAZIE!!!

    Ho fatto il terzo ciclo e sembra ok.
    Ora mi fumo una sigaretta e completo.

    Ti sono enormemente debitore.
Devi accedere o registrarti per scrivere nel forum
473 risposte