home *** CD-ROM | disk | FTP | other *** search
/ Geek Gadgets 1 / ADE-1.bin / ade-dist / octave-1.1.1p1-src.tgz / tar.out / fsf / octave / test / octave.test / fsolve / fsolve-1.m next >
Text File  |  1996-09-28  |  689b  |  46 lines
  1. x_opt = [ 0.599054;
  2.           2.395931;
  3.           2.005014 ];
  4.  
  5. tol = 1.0e-5;
  6.  
  7. function retval = f (p) 
  8.  
  9.   x = p(1);
  10.   y = p(2);
  11.   z = p(3);
  12.  
  13.   retval = zeros (3, 1);
  14.  
  15.   retval(1) = sin(x) + y**2 + log(z) - 7;
  16.   retval(2) = 3*x + 2**y -z**3 + 1;
  17.   retval(3) = x + y + z - 5;
  18.  
  19. end
  20.  
  21. [x, info] = fsolve ("f", [ 0.5, 2.0, 2.5 ]);
  22.  
  23. val = f (x);
  24.  
  25. info_bad = (info != 1);
  26. solution_bad = sum (abs (x - x_opt) > tol);
  27. value_bad = sum (abs (val) > tol);
  28.  
  29. if (info_bad)
  30.   printf ("info bad\n");
  31. else
  32.   printf ("info good\n");
  33. endif
  34.  
  35. if (solution_bad)
  36.   printf ("solution bad\n");
  37. else
  38.   printf ("solution good\n");
  39. endif
  40.  
  41. if (value_bad)
  42.   printf ("value bad\n");
  43. else
  44.   printf ("value good\n");
  45. endif
  46.