I think my program is complete What happen to this line? u(i)=sin(pi

[iwan89]

Homework Statement

I think my program is complete.. What happen to this line? u(i)=sin(pi*x(i)); ??

Homework Equations

n=10;
c1=0;
c2=0;
k= 0.0025;
L=1;
h=0.1;
alpha=1;
T=0.025;
n=L/h;m=T/k;
lambda =alpha*k/(h^2);
z=0:h:L;
disp('______________________________________________')
fprintf(' t x = ')
fprintf('%4.2f ',z)
fprintf('\n')
disp('____________________________________________ _')
fprintf('% 5.4f ',0)
% Compute the values of u at t=0
for i=1:n+1
u(i)=sin(pi*x(i));
fprintf('%10.6f ',u(i))
end
fprintf('\n')
% Compute the values of u t=jk, k=1,2,...,m
for j=1:m
t=j*k;
fprintf('% 5.4f ',t)
for i=1:n+1
if (i==1)
y(i)=c1;
elseif (i==n+1)
y(i)=c2;
else
y(i)=(1-2*lambda)*u(i)+lambda*(u(i+1)+u(i-1));
end;
fprintf('%10.6f ',y(i))
end;
fprintf('\n')
u=y;
end;

The Attempt at a Solution

There might be error somewhere. Please help me to locate it :(

 

[iwan89]

this is a MATLAB program