Advertisement
Guest User

Untitled

a guest
May 23rd, 2015
221
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 0.90 KB | None | 0 0
  1. function Fx = F2(x,v0,x0)
  2.     %funkciji glede na vpisno
  3.     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
  4.     f2 = x(1)^2 + cos(x(2))*x(3)^2 -1; %primer 1
  5.     Fx = [f1;f2;v0*(x-x0)];
  6. endfunction
  7.  
  8. function Jac = JF2(x,v0)
  9.     %vrne Jacobijevo matriko za funkciji 0 in 1
  10.     f1dx = -2*x(1)/(x(1)^2 + x(2)^2 +(x(3)-1)^2)^2;
  11.           -2*x(1)/(x(1)^2 + (x(2)-1)^2 +x(3)^2)^2;
  12.           -2*(x(1)-1)/((x(1)-1)^2 + x(2)^2 +x(3)^2)^2;
  13.     f1dy = -2*x(2)/(x(1)^2 + x(2)^2 +(x(3)-1)^2)^2;
  14.           -2*(x(2)-1)/(x(1)^2 + (x(2)-1)^2 +x(3)^2)^2 ;
  15.           -2*x(2)/((x(1)-1)^2 + x(2)^2 +x(3)^2)^2;
  16.     f1dz = -2*(x(3-1))/(x(1)^2 + x(2)^2 +(x(3)-1)^2)^2 ;
  17.           -2*x(3)/(x(1)^2 + (x(2)-1)^2 +x(3)^2)^2 ;
  18.           -2*x(3)/((x(1)-1)^2 + x(2)^2 +x(3)^2)^2;    
  19.     f2dx = 2*x(1);
  20.     f2dy = -x(3)^2 * sin(x(2));
  21.     f2dz = cos(x(2))*2*x(3);
  22.     Jac = [f1dx f1dy f1dz;f2dx f2dy f2dz;v0];
  23. endfunction
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement