Views
- State: visible
Ejemplo 2 Euler
Ejemplo simple de if
Size
1 kB
-
File type
text/x-objcsrc
Click here to get the file
File contents
% Este programa hace lo mismo que ej_euler, pero tiene un parámetro adicional: method.
% si method == 0, entonces el programa resuelve usando el método explícito. En caso contrario,
% utiliza el implícito.
function [t,y] = ej2_euler(k,tf,h,method) %devuelvo sólo una solución.
t = [0:h:tf]; %Defino t.
n = length(t); %cantidad de pasos temporales.
y = zeros(n,1); %Creo y.
y(1) = 1; %Dato inicial.
if method == 0 %Si me piden el método explícito
for i=1:n-1
y(i+1) = y(i)*(1+k*h);
end
else %En caso contrario
for i=1:n-1
y(i+1) = y(i)/(1-k*h);
end
end
