
f=@(t,u) -5*t.*u.^2+5./t-1./(t.^2)
k=.1
[t,u]=fe(f,[1 25],1,k);
plot(t,u)
hold on
pause
k=.19
[t,u]=fe(f,[1 25],1,k);
plot(t,u)
pause
k=.21
[t,u]=fe(f,[1 25],1,k);
plot(t,u)
pause
k=.29
[t,u]=fe(f,[1 25],1,k);
plot(t,u)
pause

disp('Example 7.2')
figure
k=1e-3
f=@(t,u) -10*(u-cos(t)) - sin(t)
[t,u]=fe(f,[0 1],1,k);
plot(t,u)
hold on
pause
f=@(t,u) -2100*(u-cos(t)) - sin(t)
[t,u]=fe(f,[0 1],1,k);
plot(t,u)
pause

figure
k=.0004
[t,u]=fe(f,[0 1],1,k);
plot(t,u)
hold on
pause
k=.0008
[t,u]=fe(f,[0 1],1,k);
plot(t,u)
pause
k=.00095
[t,u]=fe(f,[0 1],1,k);
plot(t,u)
pause
k=.000976
[t,u]=fe(f,[0 1],1,k);
plot(t,u)
