50/CH5/EX5.10/ex_5_10.sce


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15

// example 5.10;
// find the jacobian matrix;


// given two functions in x,y;
// and the point at which the jacobian has to be found out;

deff('[w]=f1(x,y)','w=x^2+y^2-x'); 

deff('[q]=f2(x,y)','q=x^2-y^2-y');

h=1;k=1;

J= jacobianmat (f1,f2,h,k);
disp(J);