Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function [T,x]=generate_and_solve(N,k,h,Te)
- for i=1:N-1
- a(i)=h+2*k*N^2;
- end
- D1=diag(a);
- D=D1(1:N-1,1:N-1);
- for i=1:N-1
- b(i)=-k*N^2;
- end
- L1=diag(b,-1);
- L=L1(1:N-1,1:N-1);
- for i=1:N-1
- b(i)=-k*N^2;
- end
- U1=diag(b,1);
- U=U1(1:N-1,1:N-1);
- A=D+L+U;
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement