Upogib tanke opne na krožni zanki, osnutek
function [r,u] = upogib_opne(f,R,n)
% function [r,u] = upogib_opne(f,R,n)
%
% upogib_opne racuna obliko prereza opne, napete na krozno zanko
% r je delitev v radialni smeri
% u je vektor priblizkov za resitev
% f je desna stran enacbe u''+1/r u' = f(r)
% R je radij krozne zanke
% n je stevilo delilnih intervalov (indeksi: 0,1,2,...,n)
% uporabimo kompakten zapis matrike s 3 stolpci (resi3.m)
% diskretna mreza r-jev in velikost celic h
% desna stran
% diagonale
% resimo 3-diag sistem
% narisi u(r)
% narisi u(x,y) preko surf [neobvezni dodatek]
---------------------------------------------------------------------------
function x = resi3(a,b,c,f)
% function x = resi3(a,b,c,f)
%
% resi3: resevanje tridiagonalnega sistema, predstavljenega
% s tremi vektorji in desno stranjo
%
% x je resitev sistema,
% a,b,c so pod/glavna/nad diagonale dim. n-1, n, n-1,
% f je desna stran sistema