二元一次方程 高斯消元法计算函数

2024-01-02 08:38

本文主要是介绍二元一次方程 高斯消元法计算函数,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!


   TRS = array [ 1..2 ]  of  Real ;

{二元一次方程 高斯消元法计算函数}
function  X1Y1_2 ( const  ia1 , ib1 , ic1 , ia2 , ib2 , ic2 : Real ): TRS ;
const  //为数组的行数与列数
   Rows  =  2 ;
   Cols  =  3 ;
var
   C :  array [ 1.. Rows ,  1.. Cols ]  of  Real ;  //系数及常数项
   R1 ,  R2 :  Real ;
   m , n , i :  Integer ;
begin
   C [ 1 , 1 ]:= ia1 ;
   C [ 1 , 2 ]:= ib1 ;
   C [ 1 , 3 ]:= ic1 ;
   C [ 2 , 1 ]:= ia2 ;
   C [ 2 , 2 ]:= ib2 ;
   C [ 2 , 3 ]:= ic2 ;

   //高斯消元
   for  n  :=  1  to  Rows  do
   begin
     R1  :=  C [ n , n ];
     for  m  :=  1  to  Cols  do  C [ n , m ]  :=  C [ n , m ]  /  R1 ;  //使每个方程对角线上的系数为1
     for  i  :=  1  to  Rows  do
       if  i  <>  n  then  //依次对其它方程的第n项进行消元
       begin
         R2  :=  C [ i , n ];
         for  m  :=  1  to  Cols  do  C [ i , m ]  :=  C [ i , m ]  -  C [ n , m ]* R2 ;
       end ;
   end ;
   Result [ 1 ]:= C [ 1 , 3 ];
   Result [ 2 ]:= C [ 2 , 3 ];
end ;

这篇关于二元一次方程 高斯消元法计算函数的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/561918

相关文章

hdu1171(母函数或多重背包)

题意:把物品分成两份,使得价值最接近 可以用背包,或者是母函数来解,母函数(1 + x^v+x^2v+.....+x^num*v)(1 + x^v+x^2v+.....+x^num*v)(1 + x^v+x^2v+.....+x^num*v) 其中指数为价值,每一项的数目为(该物品数+1)个 代码如下: #include<iostream>#include<algorithm>

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa

C++操作符重载实例(独立函数)

C++操作符重载实例,我们把坐标值CVector的加法进行重载,计算c3=c1+c2时,也就是计算x3=x1+x2,y3=y1+y2,今天我们以独立函数的方式重载操作符+(加号),以下是C++代码: c1802.cpp源代码: D:\YcjWork\CppTour>vim c1802.cpp #include <iostream>using namespace std;/*** 以独立函数

HDU 5833 高斯消元

n个数,任选>=1个数相乘,使得乘积是完全平方数。 其实就是开关,控制灯泡。 数 ----第i个质因子p的个数%2  = {1 , 0} == 开关----第i个灯泡 = {开 , 关} 最后使得所有灯泡都是灭着的方案数 = 2^自由变元个数   全部关着的情况     ==   一个数也不选   应省去 import java.io.BufferedReader;

函数式编程思想

我们经常会用到各种各样的编程思想,例如面向过程、面向对象。不过笔者在该博客简单介绍一下函数式编程思想. 如果对函数式编程思想进行概括,就是f(x) = na(x) , y=uf(x)…至于其他的编程思想,可能是y=a(x)+b(x)+c(x)…,也有可能是y=f(x)=f(x)/a + f(x)/b+f(x)/c… 面向过程的指令式编程 面向过程,简单理解就是y=a(x)+b(x)+c(x)

音视频入门基础:WAV专题(10)——FFmpeg源码中计算WAV音频文件每个packet的pts、dts的实现

一、引言 从文章《音视频入门基础:WAV专题(6)——通过FFprobe显示WAV音频文件每个数据包的信息》中我们可以知道,通过FFprobe命令可以打印WAV音频文件每个packet(也称为数据包或多媒体包)的信息,这些信息包含该packet的pts、dts: 打印出来的“pts”实际是AVPacket结构体中的成员变量pts,是以AVStream->time_base为单位的显

利用matlab bar函数绘制较为复杂的柱状图,并在图中进行适当标注

示例代码和结果如下:小疑问:如何自动选择合适的坐标位置对柱状图的数值大小进行标注?😂 clear; close all;x = 1:3;aa=[28.6321521955954 26.2453660695847 21.69102348512086.93747104431360 6.25442246899816 3.342835958564245.51365061796319 4.87