%%100_working code for conjugate gradient iterative algo for finding %%solution of Ax=b, where A is symmetric Positive definite %%written by vikram bhayyA on 18-6-26 ~13:19-13:59 %%general clc clear all A=input( 'enter a Symmetric Positive Definite matrix:' ); %A=[4,1,1;1,3,-1;1,-1,3]; %has to be symmetric Positve definite [m,n]=size(A); n b=input( 'enter the output vector of size n,1:' ); %b=[6,3 3]'; %x0=zeros(n,1); x0=[1 -1 0]'; %x0=[0 0 0]' ; % initial guess %initialize r0=b-A*x0; %residual if norm(r0)<0.1 display( 'solution of Ax=b is:' ) x0 ...
కామెంట్లు
కామెంట్ను పోస్ట్ చేయండి
దయచేసి మీ సలహాలను సూచనలను స్పష్టంగా పేర్కొనగలరు. plz see that ur comments are 'acceptable' in a value based society.