Question
Solution Preview
This material may consist of step-by-step explanations on how to solve a problem or examples of proper writing, including the use of citations, references, bibliographies, and formatting. This material is made available for the sole purpose of studying and learning - misuse is strictly forbidden.
function Question2()%clamped spline has endtype 2
%f(x)=x e^(-6x) f'(x)=(1-6x)e^(-6x)
%f''(x)=(-12+36x)e^(-6x)
f=@(x)(x.* exp(-6*x)); df=@(x)((1-6*x).*exp(-6*x));
endtype=2;
h=zeros(1,5);
enum=zeros(1,5);%vector of errors
alpha=zeros(1,5);%vector of orders of precision
X=0:0.001:1;
for k=1:5
n=2^k;
h(k)=1/n;
x=0:h(k):1;
Y2=[df(0),f(x),df(1)];
[pp2]=myspline(x,Y2,endtype);
[breaks2, coefs2]=unmkpp(pp2);
enum(k)=max(abs(ppval(pp2,X)-f(X)));
end
for k=2:5
alpha(k)=log(enum(k)/enum(k-1))/log(h(k)/h(k-1));
end
EnumClampedSpline=enum
alphaClampedSpline=alpha...
Return to homework library