function [w1, w2, yA, MutlakHata]=quiz11(R1, R2, V1, V2, N) h = ... % h degeri nasil hesaplaniyor? p = @(r) ... % p(x)'i yazin. q = @(r) ... % q(x)'i yazin. r = @(r) ... % r(x)'i yazin. radi = R1:h:R2; % ilk denklemin çözümü u1(1) = ... % Baslangic kosullarini yazin. u2(1) = ... % Runge Kutta yönteminden hatirlayiniz. for i = 1:length(radi)-1 k11(i) = h*u2(i); k12(i) = h*(p(radi(i))*u2(i) + q(radi(i))*u1(i) + r(radi(i))); k21(i) = h*(u2(i) + (1/2)*k12(i)); k22(i) = h*(p(radi(i)+h/2)*(u2(i)+(1/2)*k12(i)) + q(radi(i)+h/2)*(u1(i)+(1/2)*k11(i)) + r(radi(i)+h/2)); k31(i) = h*(u2(i) + (1/2)*k22(i)); k32(i) = h*(p(radi(i)+h/2)*(u2(i)+(1/2)*k22(i)) + q(radi(i)+h/2)*(u1(i)+(1/2)*k21(i)) + r(radi(i)+h/2)); k41(i) = h*(u2(i) + k32(i)); k42(i) = h*(p(radi(i)+h)*(u2(i)+k32(i)) + q(radi(i)+h)*(u1(i)+k31(i)) + r(radi(i)+h)); u1(i+1) = u1(i) + (1/6)*(k11(i) + 2*k21(i) + 2*k31(i) + k41(i)); u2(i+1) = u2(i) + (1/6)*(k12(i) + 2*k22(i) + 2*k32(i) + k42(i)); end % ikinci denklemin cozumu v1(1) = ... % Baslangic kosullarini yazin. v2(1) = ... % Runge Kutta yonteminden hatirlayiniz. for i = 1:length(radi)-1 m11(i) = h*v2(i); m12(i) = h*(p(radi(i))*v2(i) + q(radi(i))*v1(i)); m21(i) = h*(v2(i) + (1/2)*m12(i)); m22(i) = h*(p(radi(i)+h/2)*(v2(i)+(1/2)*m12(i)) + q(radi(i)+h/2)*(v1(i)+(1/2)*m11(i))); m31(i) = h*(v2(i) + (1/2)*m22(i)); m32(i) = h*(p(radi(i)+h/2)*(v2(i)+(1/2)*m22(i)) + q(radi(i)+h/2)*(v1(i)+(1/2)*m21(i))); m41(i) = h*(v2(i) + m32(i)); m42(i) = h*(p(radi(i)+h)*(v2(i)+m32(i)) + q(radi(i)+h)*(v1(i)+m31(i))); v1(i+1) = v1(i) + (1/6)*(m11(i) + 2*m21(i) + 2*m31(i) + m41(i)); v2(i+1) = v2(i) + (1/6)*(m12(i) + 2*m22(i) + 2*m32(i) + m42(i)); end % w1 y çözümü, w2 y' çözümü olmak üzere for i = 1:length(radi) w1(i) = u1(i) + ((V2-u1(N+1))/v1(N+1))*v1(i); w2(i) = u2(i) + ((V2-u1(N+1))/v1(N+1))*v2(i); end yA = ... % Analitik cozumu yazin. MutlakHata = ... % Mutlak hatayi yazin.