From 0601e6f04e4cfc1b3a4ed250875b7b0154ede1a2 Mon Sep 17 00:00:00 2001 From: Ary Andrade Neto Date: Tue, 9 Apr 2019 10:39:31 -0300 Subject: [PATCH] minimo quadrados --- RESUL.txt | 9 +++----- lista2.f90 | 2 +- lista3 | Bin 0 -> 34776 bytes lista3.f90 | 28 ++++++++++++++++++++++++ methods.f90 | 59 +++++++++++++++++++++++++++++++++++++++++++++++++-- sistemaM.txt | 3 +++ 6 files changed, 92 insertions(+), 9 deletions(-) create mode 100755 lista3 create mode 100644 lista3.f90 create mode 100644 sistemaM.txt diff --git a/RESUL.txt b/RESUL.txt index bc4950e..ad94420 100644 --- a/RESUL.txt +++ b/RESUL.txt @@ -1,6 +1,3 @@ - 5.00000000 -4.00000000 1.00000000 0.00000000 - -4.00000000 6.00000000 -4.00000000 1.00000000 - 1.00000000 -4.00000000 6.00000000 -4.00000000 - 0.00000000 1.00000000 -4.00000000 5.00000000 - 5.79999876 10.1999979 10.7999983 7.19999886 - 0.00000000 + 1.00000000 2.00000000 3.00000000 4.00000000 + 1.00000000 2.50000000 3.50000000 4.30000019 + 0.100000381 1.09000015 diff --git a/lista2.f90 b/lista2.f90 index 84346e5..8164816 100644 --- a/lista2.f90 +++ b/lista2.f90 @@ -1,7 +1,7 @@ include 'methods.f90' program lista2 - integer n, i, j + integer n, m, i, j real, allocatable, dimension(:) :: b, x real, allocatable, dimension(:,:) :: a, d, c real lambda diff --git a/lista3 b/lista3 new file mode 100755 index 0000000000000000000000000000000000000000..bdfdd20136a8c3e7532cc7fe081e2a7c99efe8f0 GIT binary patch literal 34776 zcmeHw4SZD9weJ~{K-7>)C{_HbQ>!gnO;G_!g(pDn`F{tMgwgw z@ebo-9H&h`c%^;MHuu_(Xua1!#YXDK0FeorYonq@Ma5q2od|8Iw1};B-hZvV_nCbr zGlbsTzTf-(o(D2}@3q%nd+oi~UVH6z_Bm&RyKc!Ohr>`;k@4RQO0|^<$U>R1rCUZ2 zEHvgBC*k*O<1Aw`(!~rD`3psWTyY(fq*Yid;YA3k+`LHwkQ<&Q@l>cdq_9Nh%RK(W^otUMK4-)Gx;*|0-1D@R{bO)pKWm zrn%wb=BCzg=f$1#=3YE^c2!5P>QZhu(I?!imwLG->Q>cN#Pf9g5sdO@uPJ$9c=xh@ zy3O;gpHy@__4=n@uG~gt)gPiKd@V(xX?d`xqVjz3Vc_P2R{`$IL+AAZ^e-$xe_;Xm z|0+O-c*$4qX$A1FD?ooi0XlOF(76P3&dt-_hYH~DEdYNFH1hRlpa6cR0RB4#@H-2@ zKVN{(|0nY?zK)Q;J>>HV)clBwrjH1BtO>OH+5`S(-)-$pA&{|hJ3<6+ut9=t zsGz;6b*%<%54VPzS^~a6dwZ~5gUj;12C!%c2ipRz8pL9wsWlW>3wClju@H#Qd(JjK%mXnG~3q^Uae6Jc8qI+RdIczzujMtE+1Ekt@S>CNF#{)>zjie zw#ElSzF=FZDcIU!M=aQ6*(W+SX!zULc93|~lU&R-TKvt;!Fq^5eW%~Irm5B6+_at~ zf-s~+wD_A^jWz9ofN^zQ?V_uEmsVX`HOIJm$+GKKTvxNycR6Hm691#H2>&o?IVOHZ z-uPc2L6pHB#%Pgp29MjeJvKRX~LO8vIlZKA^!*)8I)B zK0|{KY4FoE_=pBSLxY=9;9O zA~M#^>0eQr#B!{S(=St+L}hFhr=O)XiT>DfPCrR$5}C1uoZdxg5|y!8oZdlcni69! zPH&?$iOg69r|+Y5Ii(Fw|1+gYbjA*U0Koa*qcn-k*a1#|lhPzAV#-!GGm9$o#`u0#~$;U4?;6Q@Qu}iU>OlCM6Y|nOF%QipD`l~ zeoiqj@l$5T7wkYPJrzBeHu1BkWGm7RVyBxJJF*)j(H=ryIkpBE6f)H}&A%-Exy zf5~Lb*!DIKHaBvRSVciRN88`%RU6w?^9xSf@Y29b+wqb$ z4@*3*D@6%ojq`WOBhE_qF{_%C`+XAjYWMrjAMA2I^lbFm@L5Nk6`o`GDf#T@&=Tj3 zJGn>%m2hC||ru)}T{oM9H7(C>b3^=cb|u@LLXeDmsE#H`q!=_aR0))J!mtinjj=lO+}HLCzVFU_^Tn zC0Rza4^fby=$0&WARA3)5&Bg$;T{1K-AJO?CItAOF%#PvA}Z!7wS%-hv2YHtkmYC= zCrd8QZH#l`#{B3#t1+jaP-7DA0Xwf?^|-utPz_mLH*&J%_G_snTE|yR%*_OO;a)T5g`V{Em~nTH)R`lQkJxqL z6E!U38qyHrLpqWf_x8l6nf1?tV$zK5GMx{2lR|?=pM~aqrZ)Cm%_A6iqc)bTeHWEv z@duC?u;V^cBZ$YlMs~1FR{v?kX=b7wTF|paDu?2&W@MM+;G2R}%DtZutG)Z3KiKJf zXm50{Ga4sBi+g(Fm!#aIoEsihx#6{O&wjIhubEiKL^0=zAC7xQ3FLfWs5UXbR3bQj z2^xLKpRV#nkzRnx9wzmWr9|4ESj5yXemH_k58l-i0cT?k7_pw7nwx7@*4(m^2LeMd ziXp(9@eW}KlIYRo2`pO{((%?T_mJ!(G#>@>Ei8*@1kr+0O<$6ts+hh>hmc03_2#Bs zqG?p7(%`jhYfq18i4dvlu{!ra{INRIAHj3eM_gh$4^f5!1Q{}1Xjgo9_Y=7;$`|L&gK)g>c zrrQd{yyUtUc)dCvq`fvaBF%-J2PO%DC3|Brs!pY6{E>RYy~SP~Pq}eqo|iK%51D8F zz+^U=@o;iL%!$Jr&JdOX_pzE5cbbu;rx>7YWJeV2MsyIFF-cMFV8%D&vCGOoH zo8C`idhlhi#b!p_KxWFlPt2Nqs-!cz1ll|9-iOl82X@ycrWcEZV;2h2xRbnKRlO(~ zb8kde%smq6+GseVeV~MNV!6YrA?Dp*?FLiAZgNI95{h>pDG}!b&($W%rh_v|IiA51 zsEFT*cJ6gn?R4Jh2EIx}=OF)bBvJCJn`DyLF|o#xFAMNz1e{i~RP+)2 z(xw6MT>?JD@DcpdW+N5tgmsPxVm*Mt72rrkzl|6#lL-dWV4rdXqm(FV%&qiRt7tdJ z@)SMr4i#leJ7{kQBYXjqNu{Gq)v9E&sH7;jlH06G5CNxoDtQH!JSM0wcpV0I`Wnz8 z0q~Bnp867+*@Fe(3im$e#z{0FJl2WjLQ;x*V(xBalkhCh$UL%5Aw|bBu=Me`jwP?AoEUK)x+ zQmQNld*etrV3A$u?$)tLd5?KXYfm7xpQcWlSZZ@oSj=_7F1r&JA6bUcucFpm^l|fm zv1pmO6R3&xm{mgdPsXxR?LqIVS#>-0t{3Z?oK^F!S@mLeR+a6@MG^dID8)Ol3P|4a zQiG8^UKq-js)0~s*KQ>gKLp!Sqq&Enu2>yq z7jdb`uFV-t(rWt1({_MHoA*y{2$5LK(9IQ|VYM|I=7qWs1ISh|}#%%;0#gO{WEw}|nS^;H)4L|S0VoF-cB-bDKDOTUHE!nCKj z5Ah!J_Cc}&zMz{r*TdFINoTYIz{PR*Ay0j>&iTNurHQi1GPx`<-I2{nB1a`X zj-7R|SB_joD=-+p7BxSjZAZbTj-N}a9AWB=4of2%Q(eRK0 zk%rn`ara=%os4wdWf^LN!cY@6#Jt1RP#)HjA`G=bYe^wPZP+r@6!GI5YD^i1+OTD) z&7RMOS|4~zMf>r4E=CTv)Jy>j`P*$Xbp*93kEK+2su}G1rrAA#EYN*?W#}$#GW`m9^ zO@^AF#DsnzW4m!rFgpDYd4Y7#C#gZb1pJv!Om_Ja^JrsZaEn$Dq8YJzH8lPOE zgwHE$gv30o|K;j23KFBbj<1c$l>-^Q5(yTD-7yA>qgPogurSzK8M6Dao|ShSl(}Nz zCs8=5n=6(J0ZaENyAbNCcahI)J?yJ+pIS2>vwU7474sg7bRB}lscijUqBoS${AS#H zi1!6x^3qsWx;mUz8a#Y$+?y8Kd&Er4uLLG*Z&>?Y$Uz0|1~{K9_x90JLDoN;&2Vdr z`vGr~O>huqvnb8ZId4DVZX;v!O)@q`mh%BrGux^LrWRLIye(VQGB&BU%HnO27FZeQ zkbtq%EIW$V;}&Hj>`7RhVT`(W6MSdg7@2N4N<)a-Jd^rRDkSv2J86x#vL`GW5R*~N zY|YoC-#SvpgoSRmP8ij?dlVkBOaP1}%@Xe>HI~9r1OFv801b(wLDH1|(NU~U+^Cp) zU!-fFBfG2k0jsfnk+fsnhGV_RQfi5OGs?bq!zGFJUhLQ{gMc^6u>Lu>4Ay|qZTIJ( z*@%=Kmk5~@H@Wi@u%EY0pRC-5BDl}3 zJkkQ^>>VIS*>cq+(U+~uLaY`^#4I5+RqcsFthHLOvt33pF?4EMtNXa3n4!|EHQvYN zzKGj9dEfd)o>#cD-^glm zycIvMh^(SVOW9W1Zg8;q*!yHTB@R;%k5V9~Tplall7Ag9XP(rK)MDdhqtsW2+sv4F zHY?tYMGsJ&sc0Gj`p+ON+(M9W8--qUeU5iFIwDGs%F^~KbIf~U3TjYiyQeO8jF{@N z6bAd0;q+D7J&(8c6R%Texi?W+X0?x)u+e3$WzkehNUJ&8iE$FDh@%#+s(+u@-I#gm zzvne78l^K5N?+-e8;No}N=2)JhhhCwgOn%45fCQiS zDkhUke}k!t8$zgLnyAE?Tgf)75=6jho=W}$l{_Y>FPMb#=^Kcal0LQT);GAHf&P;c zfZJ!@Vc0e)cXo>{4|-e?_ZHd%)4XFfjRlEMQQ${dP1wB+It~udULH0VFPMY! z=|<3!9c!(Y(pz!F`8n=Na-Q>uu{#~Zjg1yh`1r|i9`Z&@nb;>|hS@}8y(Jo5@V-dV zxJ4I@*_L|*V=XJHD1-5vl^#MFVjEmTAv;>w+&R*-8PY@8HtCsdNe_`w(!-{Bl1*`L zGxtDVhR7PXqezA<)7~nr_AyD zcQtiUm$g+@j&M<>+Nz!b26`~hzH&j%{Cm!l(r!yiKictZ1iB{X#F zB&tLVCbuRuwkoHqD(NPQtdh1iEIay56j3v87~tlZEM@1)Pvvy@YZG_)&-YOeZf+D^ zOC}m5~)#0M?m{)|N@lA@c5m8=xCFlXFtBvea9QczB%uM|lkr>XQT zPKw(|=OIqFj%eeN>$>3sOSJ7@W!v@_ITNggu@rP%bMjuSq?BNdOFi@uqT6pcMJT-tSMEox#C4uKu;F>gA>*@DyJJX7nYYu}P>LNl zD_?<0$AWN#r>OBQ@eWMqU>VwBtfCTWDfqqGfLkQ5}+ap23UG62<1pVs>uOr}iDIekVyI z+zY&%z3K5!4!t(QQh_b1-Zbjp@q5wZ2-kN6nN9LI55~0QLEM3WRAdh!Z-iX*SVDEC za#ay-NZ?B8I{=V-1ZQA3I=@3Vx_RRfntFTJzBJ~p^XtcG36uQh5i4i(RZYX=e!6U# zIv3*39?-i1r}3`v><>Y5<*B@LCP5=?6d8AWTo-oUDKnx~+A@oTH%pudm}3@{$2Pe5_M`A!GEdtqZxWw0dxcMQyahu8lG&w% z=-7r&XytfFP%ekqc)-dMYUwI5jDTqzNZOlTfEnJw znyEAv7U!w*oH&}TmP7b>l!hqlmsRmCdzSnD0;_pw$(5s-^t3WOq@y`|(uJ+h zoxeZ@@Oeu|GOKwTWxM@hSwkA-UTFf`tPjgdQ6?~NUnDQQSY;it$~COfgj<~a!g)-Q zQ*2H>sjx9+sgH)A8%*kK`Mr!c&=lO>&F8_u{UF}oBiEKN~G2M-_cDQkKJ;c=@W z%J-o=6C)%aV&#iZ(~nh%#}vXh$34SWc!t6k#y#-3>;)7520N!5PtN*`NiT|*@PIq1?khd#5=b5X`oEC$7!qQ{|T-L{iN+0!>D4$g98 zad6h|>AM}{CcSd(HKT*Gf;Bo7t60OWa`uy1D1tQ&MI4+Jtl3c|wehUs;4BWcw2Dz< zn<;*B%AF*iP_-v%c^NweLy(x~@x02ZI8k8T#;1jcZsQX-h4OXAh?Q}87GvV(A-GC9 z#`x3+yZzoEJxlKFu?L05mGTvH4@SBMwOiJw2+87`y*d7#U!zUhK^eNGOP&)xrKaQw z7Q1zlh7kDd3b|X<7{EOvPEKK3m&|e%Z+3)K53lWuv~+vf5rMU%Ey6zin%~*9POlDPO|yJyQ-%VC-a3YlmcG0zv3Vb8rZFl zJ0<0Dr<_Opq!Yq=#LvPXKw;Y>erQs6{3Kf7VxE4S>V_rZ>8D#r*zMo}JoD+=*i-6x z!RKq=#ck;wh+#A!jWY5 zf=szABFOac`SkB%_1u+3%c+g-UCJ6Oi@(d;eGdsXKi=_i5yOK4@6s&q$+Kn0lYyIP zqj#Jq1M9@yY7%5<2x)}E4`DxXR~xe)Mb+^z7WOk=n6ga9dLR{cYEm~XYh$~GblEF} znA#70@t#{fjGE_c!~93tp6v&jHt>x_w-kRRHFyDpy5S3H4UbqR>BodWXLGl}E;)@c zl^f?(DqJi;?$MO{eil$0ZFe6%`}Qr%_2!KF;CF*!pax@>ffe!KrFi0rR6@xhs|0LA zD}czS0-!-C%AG{}Zj6q2bXM+i=UZ+r^avF!-ei?j9%1%Hk>2AYsic$G2X^r1Jy}%r z!p;0t>D22g;2m5+*o~s;6`(qGK11<*y>L9CuT(ZlR^q5vqU7vylAP!eN{L6>>BvMX znnb!AT&JRgh{?Na1f&*k{sqDWw=5yEZog>O1@W?!vY<6 znSh098Vj;BUtuz-^v^lV)0e0I`?IB}3(gkxpOIVtBNNyE&4;rM#1NAR_+B5&{qMW6_wPP{)AJ_3VOx~cRjL?)Yswk z7m96>x{O?%+7aX05iJi#*+yt9Jt7Ye)2m8$S}F=BgZY7Zi_{Q`s&mGSiuG`dT1=BB z&KhGeRfmh^V!Bs8P#;T~@4krj5A3WQ;ib;3evQ_@-=olzSpQ1-L{eS#OT7Ml9BRLp zx0G`I!=WRre>ijmCQvCJtQXHqg%4Q|)_*4MIYbLP(<>_#sanTg6zC=Wb{_io$RenV zSOkR+4x)4UOykLvIfXLiQJF`y4EH0f_Q;fryyJ*I>x`0}C^x%qdI!K3+FT1}ycRRg zXfM#fLgLE~(Mz##8c$ncvEGBl`e^hhj?RbE^pqHu|GGr+KU(WJ47J?{HY>{R zAilnJ3ac-A_ecAJ6Dr{sB}k3E38Jy^0E)%J5Au^#>!1p-eBzArGK!7yeFIuL#aT>x zmVOJZm?<}T!Rbsu8}MY5JR%%Tl=NVBLpjiF76-lPO-E*;BqHFXe8fR7bSuGc6L8W- zLX%nrkIRxy+ykWYKPp zW5%6hR!ty&@XgtFu!ZQ z`yRaJ8>aBOd$5eB^T84hNZCUH@O7nG{lcwN*e*X+esbIim?>yZg<d^(dEKi6Uro50wI6*cj7$1XAuzOa>37%XXmi@4%5gP|8XI^!-Gge+$UR zSKgn(Iew~uyE<~9D00BDF&VC{UWo?=ls%xhgPnmo<>LX4C(WI|nuKE)6{bAHr!6I0 z(KLQS;3qgEV~ygUslh3Bt^dC;kRccjVj$azvG``1TWTQBqApF$*#aC4NCi>?_47X) zKaiIJhVJL$aS_OY+=}AqZ&NmA5vXfp$u0zdZJos4Q|b39>5c^xtUQze#JBTQq>9N@#KPYtTIm93wpm6m=R@l^k_|45 zJ=?+Buz+XvErbZYlA~AoWfIgs!dA5=%VYiAwLcX^vQwMh-!Q>Z%Wt9R-)(A7Mjwx# z+HHc{ztFm+E20?_bNkild27vkCS|+_W+O%y^EvqNEv%~XWtK4gZo~1K!RQ&RhKyb zsFPnccrk_7*Tnu(8#`Ke!CU+e$(?^GNivLqu|L(tjxLVFsok;zfJF${K86MPuKW%_hjokvKuX@cT&o?P&)9ARsEt=OCZ!3Z0M+3Q$5SL z&V9YNt}4_SY7ewE`|ATWt`st2$?1;9B5n_P4BV@FUWWKTs+Cri(iZD$6Jj%iBL@l^zrgtDv-?s5>>U%LT+%~>ThqNcPhI+b$-XEP|St* zkOy3sTuQ{(Gy$%0O?2tyyI>^P>bfl$Zf!$8RkD5f;`cFD_^~lGR7sE4+8g`=QA08_l?56 z)@tvts+;nkdm)plD%NOb`{m4H$s@_?OD|&~&5zefc*NsHz?=ViCbLuTXT5Lkr=AwF zb0ct{!GZ0y2(9*9ZSnYpNwV&4yF80`Mk&!Px-?3@_9Aovues`E-H5c5BQ79OZolY zH-kU$>Y`z0_@`N|Dl(r6-cn63lz z;~nIjrYp2kgyeOhOMN%xOp!Fm1f+``9(1Yir_jEHE|ves4>CbYledE|+N#jCP)2AQ zK-YlGmuwiVGEH|X=ptVgUHZ2;C{1QRT{D6%4rMXpm~61xf06ryF4}+6m6i-D_)M8r zd?y*p!1v+ArRx8fOe?%UXGHr@9fNRRhTEH#gs~6iwBifqGpGMm=EHpDv}%VW|2GS) z@4MwtJXLn=CuKNShD&64gACWm@Kzb#A;bSE!yn7=AsOzG;R`bSwG5BSaIzG*Q)T!` z8P1jA5*gkg!yL9JOsenmTz%D5pLNaju7;Z??7F<_vZ`4ZUwV0%lb8I%CG)Cg&Avc@ zC!|O+=CHo+r=!z%i`4hHZqy}8eD8d)z)I89@pkHFiza2i=nkRPn6hq069bRfE zxwFGhGL(GV;U`=3zz#pf(g$`p21Sme9X{1qm1n$)jB@Myi*|gcrC;pu3Zq-<9Xq_z zP~Q=>!>1XIt8|En7a1S1^p&1l1aFGcQ+D`A@m)lh4&r!`ajKxZJVj`7F$=N7>8@qdq`Q@=t`U|#{E_u%aNx-OC8$!ot@Cr?TD*C@<{Bq;kY5`&2FQRW0pmPfBTx!?) zPYHbXQzLpV!v#$LC0~BN`SJmfh@Tlb>&t3IUmEUwF95h92`^qGsIyNQ(XImYdjTi; zRO^K(g?liP&bkEx!ai(7=?p|ZK1U_~!#4|jSfuk)VGH^Lg6)Qs_KpG$3}gj z7W$lJYXD!a#1}PvYny|s@zv0VP_Vtj=MQ%p^}&|5=0GUWP<8p->bd!m=;M@4_Qu^`j(1vhJ%Q_U%Qa)hPhc2_I+FdtKP3=;Al zz@Jn#e2Z^dTH~p`%HWTJ;s$5PkAm5+rD`7)UAR#DuBQA7XMIPQ8?P6UAM&)JH=qaG zTkr`_-1M~3rw@k)T5Y&FpM3pdus*oDDW?#BSkzWZd@K|o?K`03J{4-uZV!aogZNr% zhtCMm?aWqxy+63tAMOAk+XcCwLbcU;+)vlmF2h$!8=6{u;f{d*!CCA3tJ()%ZB*%l zrgm1_g14dn_%1hiqwnSB{|K(H!H)ryLz^8%{gbh_np(omp{BNG-15FP0Ff7W%29*u zldE=C8Ui#Wa3>s{)g0D)N_`pDR)daS-xzETbbM)@&&a8lYQ$%MEdg)n2>KfRtqt@6 z;A|_hUzxSlkj1b%wn_iwt{qu?J(lE_CWpp|XnVjm_Uhxgb`lmxR(`{?!OW?ILVcG4820qpCF+3RKlBs=YYmUn^5OupIn)vyhOi}k$ZHQY`>BA8wKa!~Di*COM61>YkqiYo@k@&Ykb>OI zsz9S$78)BsM5RPgu|9|rDvlcbElu^{iF&CDf?!l(BDP??Q7h#lON+;g^t4K0%dM{uvc6l;w$k!cqK7u(2Ky>Qhx-`EgX3 zl!&T+RbI7!qb$Ej;-w{m3YDLf%F~g_-11)ojQT|RX_cQ?h031_6uxrDKL)5?-tr&O zpyFDU)8TZKmM#^38*#n7@-wMW@ymtqqtYt8JEy$b|Eh4(=Tu1%s(Kalo}BW^FRnu6 zA6EFe?Y~!+SL3Jloho#x3`roj{;fIXl^;xngQ{W~=GOmsPWi@F0#Sv%G&!_O^_v2b zUxj>63QgSyC?$sxu3UVTRw2D#(uqG+UisfDPO$tNC`+=Z$}9iN3|U@H3?-kcUe%A} zry}{W$_HIY(ftArjI;=>UsYa(ZzD%9Z~6V&L=jbwBB08t^kGB@rs`Mws|r~@7oY0V zuRkI|uvt0~aaCS!->qC#UJr9bjByvr`a7in(>il8{uDn-e^9pKoeRpJ0dfrDdQ~YN bkz`(Fl!HqA<}1H>gTSJFt$r!K73}{2J8g!o literal 0 HcmV?d00001 diff --git a/lista3.f90 b/lista3.f90 new file mode 100644 index 0000000..7cf2e37 --- /dev/null +++ b/lista3.f90 @@ -0,0 +1,28 @@ +include 'methods.f90' + +program lista3 + integer n, i, j + real, allocatable, dimension(:) :: b, x, y + real, allocatable, dimension(:,:) :: a, d, c + real lambda + + ! abertura do arquivo sistema.txt e leitura + open (1, file='sistemaM.txt', status='old', action='read') + read(1,*) n + allocate(x(n)) + read(1,*) ( x(j) , j= 1 , n) + + allocate(y(n)) + read(1,*) ( y(i) , i= 1 , n) + close(1) + + allocate(b(2)) + call minimosQuadrados(b, x, y, n) + + open(2, file='RESUL.txt', status='replace') + write(2,*) (x(j), j=1, n) + write(2, *) (y(i), i=1, n) + write(2, *) (b(i), i=1, 2) + close(2) + +end program lista3 \ No newline at end of file diff --git a/methods.f90 b/methods.f90 index f9b9695..1f7d02d 100644 --- a/methods.f90 +++ b/methods.f90 @@ -450,7 +450,7 @@ subroutine powerMethod(a, x, n, tol, lambda) subroutine jocobi(a, x, n, tol, lambda) - integer :: n, i, j, k, tol, vi, vj + integer :: n, i, j, k, vi, vj real :: a(n,n), tol, x(n,n), lambda(n,n), maior x = 0.0 @@ -464,13 +464,68 @@ subroutine jocobi(a, x, n, tol, lambda) maior = 0.0 do i=1, n do j=(i+1), n - if (abs(a(i,j)) >= maior) + if (abs(a(i,j)) >= maior) then maior = abs(a(i,j)) vi = i vj = j + end if end do end do ! FALTA 2.2, 3 e 4 do slide do professor ! PREGUICA DE PENSAR NO ALGORITMO DE 2.2 end do end subroutine + +subroutine minimosQuadrados(b, x, y, n) + integer n, i, j + real b(2), x(n), y(n), a(2,2), soma, c(2), aI(2,2), det + + b = 0.0 + soma = 0.0 + + do i=1, n + soma = soma + 1 + end do + a(1,1) = soma + + soma = 0.0 + do i=1, n + soma = soma + x(i) + end do + a(2,1) = soma + a(1,2) = soma + + soma = 0.0 + do i=1, n + soma = soma + (x(i))**2 + end do + a(2,2) = soma + + soma = 0.0 + do i=1, n + soma = soma + y(i) + end do + c(1) = soma + + soma = 0.0 + do i=1, n + soma = soma + x(i)*y(i) + end do + c(2) = soma + + det = a(1,1)*a(2,2) - (a(1,2)*a(2,1)) + if (det == 0) then + STOP + end if + aI(1,1) = a(2,2)/det + aI(1,2) = -a(1,2)/det + aI(2,1) = -a(2,1)/det + aI(2,2) = a(1,1)/det + + call multiplicaVetorMatriz(aI, c, b, 2) + write(*,*) a + write(*,*) c + write(*,*) b + + +end subroutine \ No newline at end of file diff --git a/sistemaM.txt b/sistemaM.txt new file mode 100644 index 0000000..cf0d4c4 --- /dev/null +++ b/sistemaM.txt @@ -0,0 +1,3 @@ +4 +1 2 3 4 +1 2.5 3.5 4.3 \ No newline at end of file