Advertisement
Guest User

Untitled

a guest
Jul 15th, 2018
85
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Scilab 0.53 KB | None | 0 0
  1. for n = 1 : npassos;
  2.    
  3.     if n ==1 then
  4.         //Passo UPWIND
  5.         for i = 2 : N + 1
  6.             Qnmenos1(i) = Q(i) - Courant*(Q(i) - Q(i - 1));
  7.             qntItera = qntItera + 1;
  8.         end
  9.     else    
  10.         //Calculo do LeapFrog
  11.         for i = 2 : N + 1
  12.             Qnmais1(i) = Qnmenos1(i) - Courant*(Q(i) - Q(i - 1));
  13.             qntItera = qntItera + 1;
  14.         end
  15.         Qnmenos1 = Qnmais1;
  16.         Q = Qnmais1;
  17.         Q(1) = Q(N+1);
  18.     end
  19.    
  20.     //Atualizando t e Q
  21.     t = t + dt;
  22.    
  23. end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement