Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function Fx = F2(x,v0,x0)
- %funkciji glede na vpisno
- f1 = 1/(x(1)^2 + x(2)^2 +(x(3)-1)^2) + 1/(x(1)^2 + (x(2)-1)^2 + x(3)^2) + 1/((x(1)-1)^2 + x(2)^2 + x(3)^2) - 3; %primer 0
- f2 = x(1)^2 + cos(x(2))*x(3)^2 -1; %primer 1
- Fx = [f1;f2;v0*(x-x0)];
- endfunction
- function Jac = JF2(x,v0)
- %vrne Jacobijevo matriko za funkciji 0 in 1
- f1dx = -2*x(1)/(x(1)^2 + x(2)^2 +(x(3)-1)^2)^2;
- -2*x(1)/(x(1)^2 + (x(2)-1)^2 +x(3)^2)^2;
- -2*(x(1)-1)/((x(1)-1)^2 + x(2)^2 +x(3)^2)^2;
- f1dy = -2*x(2)/(x(1)^2 + x(2)^2 +(x(3)-1)^2)^2;
- -2*(x(2)-1)/(x(1)^2 + (x(2)-1)^2 +x(3)^2)^2 ;
- -2*x(2)/((x(1)-1)^2 + x(2)^2 +x(3)^2)^2;
- f1dz = -2*(x(3-1))/(x(1)^2 + x(2)^2 +(x(3)-1)^2)^2 ;
- -2*x(3)/(x(1)^2 + (x(2)-1)^2 +x(3)^2)^2 ;
- -2*x(3)/((x(1)-1)^2 + x(2)^2 +x(3)^2)^2;
- f2dx = 2*x(1);
- f2dy = -x(3)^2 * sin(x(2));
- f2dz = cos(x(2))*2*x(3);
- Jac = [f1dx f1dy f1dz;f2dx f2dy f2dz;v0];
- endfunction
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement