clear;
clear all;
clear session;

format long;

global LAMBDA;

k0 = 1;
v0 = 3.0;

LAMBDA = 9.0;

fprintf('\n k=10:\n\n');
x10 = zeros(10,1);
x10(k0,1) = v0;
[t10,y10] = colloc3(x10);

fprintf('\n k=20:\n\n');
x20 = zeros(20,1);
x20(k0,1) = v0;
[t20,y20] = colloc3(x20);

fprintf('\n k=30:\n\n');
x30 = zeros(30,1);
x30(k0,1) = v0;
[t30,y30] = colloc3(x30);

fprintf('\n k=50:\n\n');
x50 = zeros(50,1);
x50(k0,1) = v0;
[t50,y50] = colloc3(x50);

plot(t50,y50);

fprintf(' y10-y50 = %e\n',norm(y10-y50,inf));
fprintf(' y20-y50 = %e\n',norm(y20-y50,inf));
fprintf(' y30-y50 = %e\n\n',norm(y30-y50,inf));