This preview shows page 1. Sign up to view the full content.
Unformatted text preview: % L*x = s x = q*(u \ (l \ (p*s))); % dphi/ds = 1/n * (L*' \ ones) ss = (p'*(l' \ (u' \ (q'*ones(n,1)))))/n; % dphi/dwij = (x(j) - x(i))dphi/ds(i) Ws = spdiags(zeros(n,1),0,W); [i,j] = find(Ws); Ws = sparse(i,j,(x(j) - x(i)) .* ss(i)); % phi = dphi/ds dot s phi = sum(x)/n; end...
View Full Document
This note was uploaded on 03/12/2012 for the course CS 3220 taught by Professor Marschner during the Spring '09 term at Cornell.
- Spring '09