File:Loworderpoly vs highorderpoly.png

X = linspace(0,1,100);

alpha = 1e-2;

A = [ 25; 
      -4;
      15;
      4;
      8;
      1;
      2;
      0.5;
      0.5;
      1;
      1;
      0.5; ];
  
lambda = [ 1; 2; 3; 4; 5; 8; 20; 21; 25; 28; 35; 100; ];

for i=1:length(X)
    % For each term in the temperature expression, 
    % at the given spatial location
    T(i) = 0;
    for j=1:length(lambda)
        B(j) = exp(-(lambda(j)*pi*alpha).^2);
        C(j) = sin( lambda(j)*pi*X(i) );
        T(i) = T(i) + A(j)*B(j)*C(j);
    end
end

U=30*X.^4-15*(X-1).^3+12*X.^2+X-5;

plot(X,T,'r',X,U,'b','LineWidth',2);