Page 71 - Chapra y Canale. Metodos Numericos para Ingenieros 5edición_Neat

2.6 OTROS LENGUAJES Y BIBLIOTECAS 47

Guarde este archivo como euler.m y después cree otro archivo-m para calcular la deri-
vada,
function dydt = dy(t, v, m, cd)
g = 9.8;
dydt = g – (cd / m) * v;

Guarde este archivo como dy.m y regrese al modo de comandos. Para llamar la función
y ver el resultado, teclee los siguientes comandos

>> m=68.1;
>> cd=12.5;
>> ti=0;
>> tf=2.;
>> vi=0;
>> dt=0.1;
>> euler(dt,ti,tf,vi,m,cd)

Una vez dado el último comando, se desplegará el resultado

ans =
16.5309

La combinación del ambiente de MATLAB con el lenguaje de programación para
los archivos-m nos abre un mundo de posibilidades para la solución de problemas en
ingeniería. En el siguiente capítulo veremos cómo se hace esto.

2.6 OTROS LENGUAJES Y BIBLIOTECAS

En la sección anterior mostramos cómo se escribe una función en Excel o MATLAB,
para el método de Euler, a partir de un algoritmo expresado en seudocódigo. Funciones
semejantes se escriben en los lenguajes de alto nivel como Fortran 90 y C++. Por ejem-
plo, una función en Fortran 90 para el método de Euler es

Function Euler(dt, ti, tf, yi, m, cd)
REAL dt, ti, tf, yi, m, cd
Real h, t, y, dydt
t = ti
y = yi
h = dt
Do
If (t + dt > tf) Then
h = tf – t
End If
dydt = dy(t, y, m, cd)
y = y + dydt * h
t = t + h
If (t >= tf) Exit
End Do

6/12/06 13:43:45
Chapra-02.indd 47 6/12/06 13:43:45
Chapra-02.indd 47