Hej,
Jag vill integrera normalfördelningsfunktionen.
Så här långt har jag kommit:
clear all
syms x mu sigma
fx=exp(-(x-mu)^2/(2*sigma^2))/(sigma*sqrt(2*pi));
Fx=int(fx)
Fx blir då:
-(1125899906842624*2^(1/2)*pi^(1/2)*erf((2^(1/2)*(mu - x)*(1/sigma^2)^(1/2))/2))/(5644425081792261*sigma*(1/sigma^2)^(1/2))
Ur Fx vill jag lösa ut x, men det har jag inte fått att fungera. Mathematica fixar det men jag vill ha det i matlab.
Jag vet att det finns funktioner för detta (normcdf) men jag vill hitta x för att skapa en egen randn-funktion som bara bygger på rand.
Tack på förhand
Mvh
