clear all;
close all;
load ocean;
figure;
spy(A);

x = A\b;

% Plot the solution
figure;
x = [x; zeros(120,1)];
y = x(P);
y = reshape(y,121,60)';
contour(0:3:360,-89.5:3:89.5,y,20)
hold on;
load('topo.mat','topo','topomap1');
contour(0:359,-89:90,topo,[0 0],'b')

