کد روش نیوتون رافسون در متلب

[S_S_kavir]

سلام من سورس کد روش نیوتون رافسون را می خواستم که در محیط متلب پیاده سازی شده باشه. اگر کسی در اختیار داره ممنون میشم بفرسته. با تشکر

[saeed]

کد درخواستی شما در زیر قرار داده شده است:

  	 
  	
	clc;
	clear;
	
	syms x;
	f=input('enter Equation: ');
	f=char(f); % f(x)
	f_prime=char(diff(sym(f))); % f'(x) : Derivative of f(x)
	
	x0=input('Enter start point: ');
	disp(' ');
	
	maxiter=1000;
	iter=0;
	
	z=zeros(maxiter,1);
	z(1)=x0;
	
	emin=1e-10;
	e=inf;
	
	while iter<maxiter && e>emin
	iter=iter+1;
	fx=subs(f,z(iter));
	dfx=subs(f_prime,z(iter));
	z(iter+1)=z(iter)-fx/dfx;
	
	e=abs(z(iter+1)-z(iter));
	
	disp([&#39;Iteration: &#39; mat2str(iter) &#39;, x = &#39; mat2str(z(iter)) &#39;, f(x) = &#39; mat2str(fx)]);
	end
	disp(&#39; &#39;);
	disp([&#39;The answer is x = &#39; mat2str(z(iter))]);
	disp(&#39; &#39;);

[shabshekan_77]

سلام من چند سوال در مورد روش نيوتن و کرامر داشتم. چطور می تونم هم کد ها و هم توضيح کد ها رو پيدا کنم در این مورد یه روشی می خواستم که اونو بتونه به ساده ترین شکل ممکن توضیح بده مجکرم .