...
...
solver.m
function [du, tol, its] = solver(A, b, du, min_tol, max_its);
du = -A\b; tol = 0.0 ; its = 0;
endfunction
gives:
error: sparse matrix type invalid as index value
error: run_lin.m: subscript indices must be either positive integers less than 2^31 or logicals
error: called from: main.c
if I put
du = -jac\res; tol = 0.0 ; its = 0;
It works ! why ? Thank you, Guido.