Matlab - Crear una matriz D a partir de varias jacobianas hechas con un bucle

Hola, mi problema es que no soy capaz de resolver lo siguiente:

Utiliza un bucle y la sentenza fjacobiana para obtener una aproximación de las derivadas parciales de primera
orden en los cuatro puntos (1.k, 2.0), k = 0, 1, 2, 3, 4, dentro de una matriz, D1, de forma que en cada fila k aparezca la jacobiana.

Gracias.

y qué función vas a optimizar ????, el jacobiano es una matriz , al menos que halles el determinante jacobiano , porque de lo contrario el bucle debes hacerlo en un array y no en una matriz

La función es sinh(x^2+x*y).
Meandan poner las jacobianas de esos puntoa dentro de una matriz llamada D1.