用追赶法求解三对角方程组

如题：{2,1,0,0;1,3,0,0;0,1,1,1;0,0,2,1}{x1,x2,x3,x4}={1,2,2,0}

上面是矩阵相乘，会做的应该看得懂，在线等答案。。。

module tridiagonalmatriximplicit nonecontains !追赶法 function tri(a,b,c,f,n)integer n,ireal:: a(2:n),b(n),c(1:n-1),d(n-1),e(n-1),f(n)real,target::x(n)real,pointer::tri(:)tri=>x(1:n) !数组指针 !追赶法的核心 d(1)=f(1)/b(1)e(1)=c(1)/b(1)do i=2,n-1e(i)=c(i)/(b(i)-a(i)*e(i-1))d(i)=(f(i)-a(i)*d(i-1))/(b(i)-a(i)*e(i-1))enddox(n)=(f(n)-a(n)*d(n-1))/(b(n)-a(n)*e(n-1))do i=n-1,1,-1x(i)=d(i)-e(i)*x(i+1)end doreturnend functionend module!******************************************* !******************************************* program main !主函数 use tridiagonalmatriximplicit noneinteger,parameter::n=3 !变量x的个数 integer i,sreal:: a(2:n),b(n),c(1:n-1),f(n)a=(/2,3/)b=(/1,2,3/)c=(/2,3/) f=(/1,2,3/)write(*,*) 'written by chenke'write(*,*) 'The roots is:'write(*,*) tri(a,b,c,f,n)end

a=[-2,-2,-2,-2,-2,-2,-2]';       %由线性方程组得到的值，下面b,c,d均此道理 b=[2,5,5,5,5,5,5,5,]'; c=[-2,-2,-2,-2,-2,-2,-2]'; d=[220/27,0,0,0,0,0,0,0]'; n=length(b); %对A进行分解 u(1)=b(1); y(1)=d(1); for i=2:n         l(i-1)=a(i-1)/u(i-1);         %求l(i)         u(i)=b(i)-l(i-1)*c(i-1);   %求u(i),i=2,3,4...n end L=eye(n)+diag(l,-1); U=diag(u)+diag(c,1); for i=2:n y(i)=d(i)-y(i-1)*l(i-1);    %由Ly=d ,求解 y end %求解  Ux=y if (u(n)~=0)     x(n)=y(n)/u(n); end for i=n-1:-1:1     x(i)=(y(i)-c(i)*x(i+1))/u(i); end disp('L=');        %输出矩阵L disp(L); disp('U='); disp(U);          %输出矩阵U disp('x='); disp(x);

网上买礼物吧，我空间有详细的网购步骤和技巧总结。 我在网购经验比较多，网上的东西真的便宜很多，选购也很方便，我收集了很多热卖好评的商品和店铺，并做了统计排行，很多都是专家通过比较店铺信誉和销售记录以及网友的评价，做出的排行榜，当然也有很多是我通过购买和网友的交流统计出来的，都是热卖好评的，网购这么多年了，现在才知道，原来这样统计下，真的方便很多，现在分享给大家，当然主要是希望大家给我空间加加人气，还有采纳我的答案，让我赚赚分^_^，地址: taobibuy点cn（把“点”改成“.”访问） ,或者你可以百度搜"taobibuy"，第一个就是了，那里有我的超级经验分享，有我总结的详细购物步骤和购物心得，肯定对你购物有很大帮助！快去看看吧，登陆的人比较多，打不开，请多刷新几次. o(∩_∩)o希望对您有帮助，希望采纳我哦~

如以上回答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！

点此我要举报以上问答信息