0001 function [out] = gradtest(v, model, expdata) 0002 0003 v = model.N\v; 0004 0005 Prob.user.expdata = expdata; 0006 Prob.user.model = model; 0007 0008 interval = -1:11; 0009 for i = 1:length(interval) 0010 i 0011 Prob.user.diff_interval = 10^(-interval(i)); 0012 grad = errorComputation2_grad(v, Prob); 0013 grads(i,:) = grad; 0014 end 0015 0016 for i = 1:size(grads,2) 0017 semilogy(interval, abs(grads(:,i))); 0018 hold on; 0019 end